Planetary growth with collisional fragmentation and gas drag

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct N-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about ∼0.1 Earth mass within 10⁷ years. This brings another difficulty in formation of gas giant planets, which requires cores with ∼10 Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos’ collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the Solar System.     

Orbital resonances and planetary formation sites

Orbital resonances may have played an important role in determining the locations where the planetesimal swarm eventually accreted into full-size planets. Several pairs of planets do indeed have commensurable orbital periods at present, but the case for control of planet formation by resonances is weakened by the fact that many pairs are not commensurable and that those which are do not necessarily exist at the strongest resonances. However, the mass loss and redistribution that occurred in the early solar system evolution can substantially alter the positions of planets and planetary embryos within the swarm. A cascaded resonance structure is hypothesized where planetesimal growth was accelerated at 2:1 interior and 1:2 exterior resonances with an early-formed Jupiter producing runaway growth of planetary embryos. These embryos produce their own resonances which, in turn, lead to additional embryos in a process that successively propagates inward and outward to generate a resonant configuration of embryos. In this manner, the early presence of Jupiter imposed a harmonic structure on the accumulating planetesimal swarm. For an accretion disk with surface density obeying a power law of index ?1.2 the positions of the planetary embryos can be moved into a reasonably good agreement with most of the present planetary positions that is as good as that given by the Titius-Bode law.     

From boulders to planetary systems

In this chapter, we will give a brief overview on our current theoretical understanding how planets form from the solid material in circumstellar disks in the core accretion-gas capture model. This chapter will not be as concise and complete as a review on this matter, yet will serve as an introductory text to generate interest in the subject. Students are referred to comprehensive text books and some important reviews.This chapter will discuss “dusty storms”, e.g. the dust transport in turbulent protoplanetary disks, followed by the latest model of planetesimal formation, e.g. gravoturbulent planetesimal formation, which deals with particle concentration in turbulence and N-body simulations thereof. We also briefly describe the core accretion-gas capture process and talk about nascent planets, e.g. the observability of planet–disk interaction concluding with the migration of young planets and the final arrangement of planetary systems.     

The full set of gas giant structures: I: On the origin of planetary masses and the planetary initial mass function

Context

Current planet search programs are detecting extrasolar planets at a rate of 60 planets per year. These planets show more diverse properties than was expected.

Aims

We try to get an overview of possible gas giant (proto-) planets for a full range of orbital periods and stellar masses. This allows the prediction of the full range of possible planetary properties which might be discovered in the near future.

Methods

We calculate the purely hydrostatic structure of the envelopes of proto-planets that are embedded in protoplanetary disks for all conceivable locations: combinations of different planetesimal accretion rates, host star masses, and orbital separations. At each location all hydrostatic equilibrium solutions to the planetary structure equations are determined by variation of core mass and pressure over many orders of magnitude. For each location we analyze the distribution of planetary masses.

Results

We get a wide spectrum of core-envelope structures. However, practically all calculated proto-planets are in the planetary mass range. Furthermore, the planet masses show a characteristic bimodal, sometimes trimodal, distribution. For the first time, we identify three physical processes that are responsible for the three characteristic planet masses: self-gravity in the Hill sphere, compact objects, and a region of very low adiabatic pressure gradient in the hydrogen equation of state. Using these processes, we can explain the dependence of the characteristic masses on the planet’s location: orbital period, host star mass, and planetesimal accretion rate (luminosity). The characteristic mass caused by the self-gravity effect at close proximity to the host star is typically one Neptune mass, thus producing the so-called hot Neptunes.

Conclusions

Our results suggest that hot Jupiters with orbital period less than 64 days (the exact location of the boundary depends on stellar type and accretion rate) have quite distinct properties which we expect to be reflected in a different mass distribution of these planets when compared to the “normal” planetary population. We use our theoretical survey to produce an upper mass limit for embedded planets: the maximum embedded equilibrium mass (MEEM). This naturally explains the lack of high mass planets between 3 and 64 days orbital period.     

Gravitational interaction between a planet and an optically thin disc

We investigate the gravitational interaction between a planet and an optically thin protoplanetary disc, performing local three-dimensional hydrodynamical simulations. In the present study, we take account of radiative energy transfer in optically thin discs. Before the stage of planetary accretion, dust opacity is expected to decrease significantly because of grain growth and planetesimal formation. Thus, it would be reasonable to consider optically thin discs in the disc–planet interaction. Furthermore, we focus on small planets that can neither capture disc gas nor open a disc gap. The one-sided torque exerted on a planet by an optically thin disc is examined for various values of the disc optical thickness (<1). In optically thin discs, the temperature behind the density waves is lower than the unperturbed value because of radiative cooling. Heating due to shock dissipation is less effective than radiative cooling. Because of radiative cooling, the density distribution around the planet is not axisymmetric, which exerts an additional torque on the planet. The torque enhancement becomes maximum when the cooling time is comparable with the Keplerian period. The enhancement is significant for low-mass planets. For planets with  3 M_⊕  , the additional one-sided torque can be 40 per cent of the torque in the isothermal case. The radiative cooling is expected to change the differential torque and the migration speed of planets, too.     

“Horseshoe” Structures in the Debris Disks of Planet-Hosting Binary Stars

The formation of a planetary system from the protoplanetary disk leads to destruction of the latter; however, a debris disk can remain in the form of asteroids and cometary material. The motion of planets can cause the formation of coorbital structures from the debris disk matter. Previous calculations have shown that such a ring-like structure is more stable if there is a binary star in the center of the system, as opposed to a single star. To analyze the properties of the coorbital structure, we have calculated a grid of models of binary star systems with a circumbinary planet moving in a planetesimal disk. The calculations are performed considering circular orbits of the stars and the planet; the mass and position of the planet, as well as the mass ratio of the stars, are varied. The analysis of the models shows that the width of the coorbital ring and its stability significantly depend on the initial parameters of the problem. Additionally, the empirical dependences of the width of the coorbital structure on the parameters of the system have been obtained, and the parameters of the models with the most stable coorbital structures have been determined. The results of the present study can be used for the search of planets around binary stars with debris disks.     

Debris disks: seeing dust, thinking of planetesimals and planets

Debris disks are optically thin, almost gas-free dusty disks observed arounda significant fraction of main-sequence stars older than about 10 Myr. Since the circumstellar dust is short-lived, the very existence of these disks is considered as evi-dence that dust-producing planetesimals are still present in mature systems, in whichplanets have formed – or failed to form – a long time ago. It is inferred that theseplanetesimals orbit their host stars at asteroid to Kuiper-belt distances and continuallysupply ...     

内日冕的绿线强度分布研究

日冕是太阳大气活动的关键区域,是日地空间天气的源头.受观测限制,对日冕低层大气等离子体结构和磁场状态的研究非常欠缺,国际上对于可见光波段日冕低层大气的亮度分层研究很少.利用丽江日冕仪YOGIS(Yunnan Green-line Imaging System)的日冕绿线(Fe_ⅩⅣ5303?)观测资料,对内日冕区域(1.03R_☉-1.25R_☉,R_☉表示太阳半径)亮结构及其中冕环进行了有效的强度衰减分析.对亮结构的强度在太阳径向高度上进行了指数衰减拟合,比较这些拟合结果发现所得到的静态内冕环的衰减指数在一固定值附近.然后将比较明显的冕环提取出来,通过对不同高度的绿线强度进行指数拟合,得出的衰减指数与亮结构中也比较相近,这对进一步研究日冕中的各项物理参数演化提供了参考.     

The provenance and evolution of comets

The process of comet formation through the hierarchical aggregation of originally submicron-sized interstellar grains to form micron-sized particles and then larger bodies in the protoplanetary disc, culminating in the formation of planetesimals in the disc extending from Jupiter to beyond Neptune, is briefly reviewed. The ‘planetesimal’ theory for the origin of comets implies the existence of distinct cometary reservoirs, with implications for the immediate provenance of observed comets (both long-period and short-period) and their evolution as a result of planetary perturbations and physical decay, for example splitting and sublimation. The principal mode of cometary decay and collisional interaction with the terrestrial planets is through the formation and evolution of streams of cometary debris and hitherto undiscovered ‘families’ of cometary asteroids. Recent dynamical results, in particular the sungrazing and sun-colliding end-state for short-period comet and asteroid orbits, are briefly discussed.     

Gap formation in a self-gravitating disk and the associated migration of the embedded giant planet

We present the results of our recent study on the interactions between a giant planet and a self-gravitating gas disk. We investigate how the disk's self-gravity affects the gap formation process and the migration of the giant planet. Two series of 1-D and 2-D hydrodynamic simulations are performed. We select several surface densities and focus on the gravitationally stable region. To obtain more reliable gravity torques exerted on the planet, a refined treatment of the disk's gravity is adopted in the vicinity of the planet. Our results indicate that the net effect of the disk's selfgravity on the gap formation process depends on the surface density of the disk. We notice that there are two critical values, ΣIand ΣII. When the surface density of the disk is lower than the first one, Σ0 ΣI, the effect of self-gravity suppresses the formation of a gap. When Σ0 ΣI, the self-gravity of the gas tends to benefit the gap formation process and enlarges the width/depth of the gap. According to our 1-D and2-D simulations, we estimate the first critical surface density to be ΣI≈ 0.8 MMSN.This effect increases until the surface density reaches the second critical value ΣII.When Σ0 ΣII, the gravitational turbulence in the disk becomes dominant and the gap formation process is suppressed again. Our 2-D simulations show that this critical surface density is around 3.5 MMSN. We also study the associated orbital evolution of a giant planet. Under the effect of the disk's self-gravity, the migration rate of the giant planet increases when the disk is dominated by gravitational turbulence. We show that the migration timescale correlates with the effective viscosity and can be up to 104yr.     

Planetary migration in a planetesimal disk: why did Neptune stop at 30 AU?

We study planetary migration in a gas-free disk of planetesimals. In the case of our Solar System we show that Neptune could have had either a damped migration, limited to a few AUs, or a forced migration up to the disk's edge, depending on the disk's mass density. We also study the possibility of runaway migration of isolated planets in very massive disk, which might be relevant for extra-solar systems. We investigate the problem of the mass depletion of the Kuiper belt in the light of planetary migration and conclude that the belt lost its pristine mass well before that Neptune reached its current position. Therefore, Neptune effectively hit the outer edge of the proto-planetary disk. We also investigate the dynamics of massive planetary embryos embedded in the planetesimal disk. We conclude that the elimination of Earth-mass or Mars-mass embryos originally placed outside the initial location of Neptune also requires the existence of a disk edge near 30AU.     

Late stages of planetary accretion

We derive first-order differential equations for the late stages of planetary accretion (planetesimal mass >10¹³ g). The effect of gravitational encounters, energy exchange, collisions, and gas drag has been included. Two simple models are discussed, namely, (i) when all planetesimals have the same mass and (ii) when there is one large planetesimal and numerous small planetesmals. Gravitational two-body encounters are modeled according to Chandrasekhar's classical theory from stellar dynamics. It is shown that the velocity increase due to mutual encounters can be modeled according to the simple theory of random flights. We find analytical equations for the average velocity decrease due to collisions. Gas drag, if present, is modeled in averaged form up to the first order in the eccentricities and inclinations of the planetesimals. Characteristic time scales for the formation of the terrestrial planets are found for the most favorable models to be of order 10⁸ year. The calculated mass of rock and ice of the giant planets is too low as compared to the observed one. This difficulty of our model could be overcome by assuming a several times larger surface density, an enlarged accretion cross section, and gas accretion during the final stages of accretion of the solid cores of the giant planets. Analytical and numerical results are presebted, the evolutionary tracks showing satisfactory agreement with observations for some models.     

Evolution of Planetary Chaotic Zones in Planetesimal Disks

Astronomy Letters - Extensive numerical experiments on the long-term dynamics of planetesimal disks with planets in systems of single stars have been carried out. The planetary chaotic zone...     

Effects of Planetary Migration on Natural Satellites of the Outer Planets

Numerous studies in the past few years have analyzed possible effects of planetary migration on the small bodies of the Solar System (mainly asteroids and KBOs), with the double aim of explaining certain dynamical structures in these systems, as well as placing limits on the magnitude of the radial migration of the planets. Here we undertake a similar aim, only this time concentrating on the dynamical stability of planetary satellites in a migration scenario. However, different from previous works, the strongest perturbations on satellite systems are not due to the secular variation of the semimajor axes of the planets, but from the planetesimals themselves. These perturbations result from close approaches between the planetesimals and satellites.We present results of several numerical simulations of the dynamical evolution of real and fictitious satellite systems around the outer planets, under the effects of multiple passages of a population of planetesimals representing the large-body component of a residual rocky disk. Assuming that this component dominated the total mass of the disk, our results show that the present systems of satellites of Uranus and Neptune do not seem to be compatible with a planetary migration larger than even one quarter that suggested by previous studies, unless these bodies were originated during the late stage of evaporation of the planetesimal disk. For larger variations of the semimajor axes of the planets, most of the satellites would either be ejected from the system or suffer mutual collisions due to excitation in their eccentricities. For the systems of Jupiter and Saturn, these perturbations are not so severe, and even large migrations do not introduce large instabilities.Nevertheless, even a small number of 1000-km planetesimals in the region may introduce significant excitation in the eccentricities and inclinations of satellites. Adequate values of this component may help explain the present dynamical distribution of distant satellites, including the highly peculiar orbit of Nereid.     

On the diversity of giant planets—Simulating the evolution of solids in protoplanetary disks

We describe a model designed to track simultaneously the evolution of gas and solids in protoplanetary disks from an early stage, when all solids are in the dust form, to the stage when most solids are in the form of a planetesimal swarm. The model is computationally efficient and allows for a global, comprehensive approach to the evolution of solid particles due to gas–solid coupling, coagulation, sedimentation, and evaporation/condensation. We have used it to calculate the co-evolution of gas and solids starting from a comprehensive domain of initial conditions. Then based on the core accretion-gas capture scenario, we have estimated the planet-bearing capability of the environment defined by the final planetesimal swarm and the still evolving gaseous component of the disk. We describe how the disk's capability of formation of giant planets depends on the initial mass and size of a protoplanetary disk, its thermal structure, mass of the central star and properties of the material forming solid grains.     

Astronomical evidence for the rapid growth of millimeter‐sized particles in protoplanetary disks

Abstract– I summarize recent surveys of protoplanetary disks at millimeter wavelengths and show that the distribution of luminosity, equivalent to the mass in small dust grains, declines rapidly. This contrasts with statistics on the lifetime of disks from infrared observations and the high occurrence of planets from radial velocity and transit surveys. I suggest that these disparate results can be reconciled if most of the dust in a disk is locked up in millimeter and larger‐sized particles within about 2 Myr. This statistical result on disk evolution agrees with detailed modeling of a small number of individual disks and with cosmochemical measurements of rapid planetesimal formation.     

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.     

The origin of the atmospheres of the earth and the planets

From a systematic analysis of the whole history of the protoplanetary cloud and the observational facts of the Earth's atmosphere, we propose a new theory of the origin of the atmospheres of the Earth and the planets. For the Earth-like planets, there were extended primordial atmospheres originating from the protoplanetary cloud by the accretion of the planetary embryoes. These primordial atmospheres existed on a time scale of 10³–10⁷ years and were composed of chemically reducing gases. The presence of such a reducing atmosphere may be of great significance to the theories of cosmogony and the origin of life.The contents are 1. The escape of the nebula and the planetary atmospheres. 2. The thermal dissipation of the atmospheres and their blowing-off by the solar wind. 3. The accretion of gases by the planetary enbryoes. 4. The primordial atmospheres.     

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.     

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_⊕.