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.     

Eccentricity growth of planetesimals in a self-gravitating protoplanetary disc

We investigate the orbital evolution of planetesimals in a self-gravitating circumstellar disc in the size regime (∼1–5000 km) where the planetesimals behave approximately as test particles in the disc's non-axisymmetric potential. We find that the particles respond to the stochastic, regenerative spiral features in the disc by executing large random excursions (up to a factor of 2 in radius in ∼1000 yr), although typical random orbital velocities are of the order of one tenth of the Keplerian speed. The limited time frame and small number of planetesimals modelled do not permit us to discern any net direction of planetesimal migration. Our main conclusion is that the high eccentricities (∼0.1) induced by interaction with spiral features in the disc is likely to be highly unfavourable to the collisional growth of planetesimals in this size range while the disc is in the self-gravitating regime. Thus if , as recently argued by Rice et al., the production of planetesimals gets under way when the disc is in the self-gravitating regime (either at smaller planetesimal size scales, where gas drag is important, or via gravitational fragmentation of the solid component), the planetesimals thus produced would not be able to grow collisionally until the disc ceases to be self-gravitating. It is unclear, however, given the large amplitude excursions undergone by planetesimals in the self-gravitating disc, whether they would be retained in the disc throughout this period, or whether they would instead be lost to the central star.     

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.     

Gravitational Lensing Effects in a Modified General Relativity Model

We investigate the effects of gravitational lensing in a modified generalrelativity flat cosmological model. We show that the predicted lensingprobability for low matter densities is larger than in the Einstein-deSitter model but considerably less than in constant- models. Wealso find that the optical depth at low density rises with the sourcered-shift more slowly than in constant- models and that the lensdistribution peaks at lower values of the lens red-shift.     

Kuiper Belt structure as a reflection of the migration process of a planet

One of the main particular features of the structure of the Kuiper Belt is that it contains clusters of objects of small orbital eccentricity and inclination (“cold population”). In order to solve the problem of the origin of the objects, we considered the process of the gravitational interaction of a comparatively small-mass planet with a planetesimal disk. We found that one particular property of the process is that the planet changes its direction of migration. The interaction with the planet results in the transportation of a considerable portion of planetesimals from the inner zone out to the Kuiper Belt. After such a transition of the objects, the planet returns to the inner regions of the planetesimal disk. Numerical simulations show that the reversible migration of a planet of a mass similar to that of the Earth can explain the main properties of the Kuiper Belt’s cold population orbit distribution.     

Gaseous drag and planetary formation by accretion

We have made numerical experiments of the collisional and gravitational interaction of a planetesimal swarm in the early Solar System. In particular we study the dynamical evolution of an initial population of kilometer-size planetesimals subject to collisions (accretion, rebound, cratering, and catastrophic fragmentation). This study is based on a Monte-Carlo statistical method and provides the mass and velocity distributions of the planetesimal swarm as a function of time as well as their distribution in heliocentric distance. Several experiments have been performed and three of them are presented here. They simulate the accretional growth of numerous planetesimals in the absence (or presence) of gaseous drag, with (or without) one larger embryo among them, and with (or without) a size gradient. The results show that (i) for a population of planetesimals submitted to a negative gradient in size as the heliocentric distance increases, the outer planetesimals spiral toward the Sun faster than inner ones, leading after some time to an accumulation of bodies inside the cloud which allows the formation of an embryo; (ii) the growth of one embryo among a population of planetesimals is accelerated by the presence of gas and is warranted as long as its feeding zone is fed by the inward flow of planetesimals due to gas drag. These results offer some complementary new insights in the understanding of the accretional formation of 4–5 terrestrial planets instead of the numerous Moon-size planets generally found in numerical experiments.     

10-µm images and spectra of T Tauri stars

T Tauri stars are young stars usually surrounded by dusty disks similar to the one from which we believe our own Solar System formed. Most T Tauri stars exhibit a broad emission or absorption band between 7.5 and 13.5µm which is attributed to silicate grains in the circumstellar environment. We imaged three spatially resolved T Tauri binaries through a set of broadband filters which include the spectral region occupied by the silicate band. Two of these objects (T Tauri and Haro 6–10) are infrared companion systems in which one component is optically much fainter but contributes strongly in the infrared. Both infrared companions exhibit a deep silicate absorption which is not present in their primaries, indicating that they suffer very strong local extinction which may be due to an edge-on circumstellar disk or to a dense shell. We also took low resolution spectra of the silicate feature of two unresolved T Tauris to look for narrow features in the silicate band which would indicate the presence of specific minerals such as olivine. We observed GK Tau, for which Cohen and Witteborn (1985) reported a narrow emission feature at 9.7µm, but do not find evidence for this feature, and conclude that it is either time-dependent or an artifact of absorption by telluric ozone.Paper presented at the Conference onPlanetary Systems: Formation, Evolution, and Detection held 7–10 December, 1992 at CalTech, Pasadena, California, U.S.A.     

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.     

Dust to planetesimals: Settling and coagulation in the solar nebula

The behavior of solid particles in a low-mass solar nebula during settling to the central plane and the formation of planetesimals is examined. Gravitational instability in a dust layer and collisional accretion are considered as possible mechanisms of planetesimal formation. Non-Keplerian rotation of the nebula results in shear between the gas and a dust layer. This shear produces turbulence within the layer which inhibits gravitational instability, unless the mean particle size exceeds a critical value, ～1 cm at 1 AU. The size requirement is less stringent at larger heliocentric distances, suggesting a possible difference in planetesimal formation mechanisms between the inner and outer nebula. Coagulation of grains during settling is expected in the solar nebula environment. Van der Waals forces appear adequate to produce centimeter-sized aggregates. Growth is primarily due to sweepup of small particles by larger ones due to size-dependent settling velocities. A numerical model for computing simultaneous coagulation and settling is described. Relative velocities are determined by gas drag and the non-Keplerian rotation of the nebula. The settling is very nonhomologous. Most of the solid matter reaches the central plane as centimeter-sized aggregates in a few times 10³ revolutions, but some remains suspended in the form of fine dust. Drag-induced relative velocities result in collisions. The growth of bodies in the central plane is initially rapid. After sizes reach ～10³ cm, relative velocities decrease and the growth rate declines. Gas drag rapidly damps the out-of-plane motions of these intermediate-sized bodies. They settle into a thin layer which is subject to gravitational instability. Kilometer-sized planetesimals are formed by this composite process.     

Planetesimal formation by turbulent concentration

The formation of 1-1000 km diameter planetesimals from dust grains in a protoplanetary disk is a key step in planet formation. Conventional models for planetesimal formation involve pairwise sticking of dust grains, or the sedimentation of dust grains to a thin layer at the disk midplane followed by gravitational instability. Each of these mechanisms is likely to be frustrated if the disk is turbulent. Particles with stopping times comparable to the turnover time of the smallest eddies in a turbulent disk can become concentrated into dense clumps that may be the precursors of planetesimals. Such particles are roughly millimeter-sized for a typical protoplanetary disk. To survive to become planetesimals, clumps need to form in regions of low vorticity to avoid rotational breakup. In addition, clumps must have sufficient self gravity to avoid break up due to the ram pressure of the surrounding gas. Given these constraints, the rate of planetesimal formation can be estimated using a cascade model for the distribution of particle concentration and vorticity within eddies of various sizes in a turbulent disk. We estimate planetesimal formation rates and planetesimal diameters as a function of distance from a star for a range of protoplanetary disk parameters. For material with a solar composition, the dust-to-gas ratio is too low to allow efficient planetesimal formation, and most solid material will remain in small particles. Enhancement of the dust-to-gas ratio by 1-2 orders of magnitude, either vertically or radially, allows most solid material to be converted into planetesimals within the typical lifetime of a disk. Such dust-to-gas ratios may occur near the disk midplane as a result of vertical settling of short-lived clumps prior to clump breakup. Planetesimal formation rates are sensitive to the assumed size and rotational speed of the largest eddies in the disk, and formation rates increase substantially if the largest eddies rotate more slowly than the disk itself. Planetesimal formation becomes more efficient with increasing distance from the star unless the disk surface density profile has a slope of −1.5 or steeper as a function of distance. Planetesimal formation rates typically increase by an order-of-magnitude or more moving outward across the snow line for a solid surface density increase of a factor of 2. In all cases considered, the modal planetesimal size increases with roughly the square root of distance from the star. Typical modal diameters are 100 km and 400 km in the regions corresponding to the asteroid belt and Kuiper belt in the Solar System, respectively.     

The application of a non-static spherically-symmetric metric to clusters of galaxies

An approximate metric is found which represents a sphere of matter embedded in a background of dust. The use of this metric in conjunction with the Friedmann equations gives values of for the three possible values ofk as +6×10^–36 (k=+1), +3×10^–35 (k=0), +10^–36 (k=–1). These values depend on data regarding clusters of galaxies, and are probably accurate to within an order of magnitude given the correctness of the assumptions on which their derivation rests.     

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.     

The circumstellar gas of sigma orionis E

In order to explain the variable H emission and the eclipse-like light variation of Ori E, we investigated the circumstellar gas trapped by the stellar magnetic field and corotating with the star. By considering the potential along the magnetic field line, we found that the gas concentrates to a potential minimum. The circumstellar gas forms either two condensations or a disk, depending on the inclination of the magnetic dipole to the stellar rotation axis. The geometrical thickness of the circumstellar disk, of about 0.2 stellar radii, and the distance from the center of the star to the inner edge of the disk, of about 3 stellar radii, were obtained. The H emission line profile at its maximum phase and the amplitude of light variation were calculated by assuming the isothermal gas in LTE with the maximum gas density which the magnetic field can hold. The model gives good agreement with observation in the low obliquity case, and also explains the phase correlation among the H emission maximum, the light minimum, and the magnetic extreme. The model, however, failed to explain the large IR excess in theM band.     

Orbital Migration Of Giant Planets: Using Numerical Integration To Investigate Consequences For Other Bodies

A number of extrasolar planets have been detected in close orbits around nearby stars. It is probable that these planets did not form in these orbits but migrated from their formation locations beyond the ice line. Orbital migration mechanisms involving angular momentum transfer through tidal interactions between the planets and circumstellar gas-dust disks or by gravitational interaction with a residual planetesimal disk together with several means of halting inward migration have been identified. These offer plausible schemes to explain the orbits of observed extrasolar giant planets and giant planets within the Solar System. Recent advances in numerical integration methods and in the power of computer workstations have allowed these techniques to be applied to modelling directly the mechanisms and consequences of orbital migration in the Solar System. There is now potential for these techniques also to be applied to modelling the consequences of the orbital migration of planets in the observed exoplanetary systems. In particular the detailed investigation of the stability of terrestrial planets in the habitable zone of these systems and the formation of terrestrial planets after the dissipation of the gas disk is now possible. The stability of terrestrial planets in the habitable zone of selected exoplanetary systems has been established and the possibility of the accretion of terrestrial planets in these systems is being investigated by the author in collaboration with Barrie W. Jones (Open University), and with John Chambers (NASA-Ames) and Mark Bailey of Armagh Observatory, using numerical integration. The direct simulation of orbital migration by planetesimal scattering must probably await faster hardware and/or more efficient algorithms. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

Dust in the disk winds from young stars as a source of the circumstellar extinction

We consider the problem of dust grain survival in the disk winds from T Tauri and Herbig Ae stars. For our analysis, we have chosen a disk wind model in which the gas component of the wind is heated through ambipolar diffusion to a temperature of ～10⁴ K. We show that the heating of dust grains through their collisions with gas atoms is inefficient compared to their heating by stellar radiation and, hence, the grains survive even in the hot wind component. As a result, the disk wind can be opaque to the ultraviolet and optical stellar radiation and is capable of absorbing an appreciable fraction of it. Calculations show that the fraction of the wind-absorbed radiation for T Tauri stars can be from 20 to 40% of the total stellar luminosity at an accretion rate ? _a = 10^?8-10^?6 M _⊙yr^?1. This means that the disk winds from T Tauri stars can play the same role as the puffed-up inner rim in current accretion disk models. In Herbig Ae stars, the inner layers of the disk wind (r ≤ 0.5 AU) are dust-free, since the dust in this region sublimates under the effect of stellar radiation. Therefore, the fraction of the radiation absorbed by the disk wind in this case is considerably smaller and can be comparable to the effect from the puffed-up inner rim only at an accretion rate of the order of or higher than 10^?6 M _⊙yr^?1. Since the disk wind is structurally inhomogeneous, its optical depth toward the observer can be variable, which should be reflected in the photometric activity of young stars. For the same reason, moving shadows from gas and dust streams with a spiral-like shape can be observed in high-angular-resolution circumstellar disk images.     

Wind dynamics and circumstellar extinction variations in the T Tauri star RY Tau

The wind interaction with the dusty environment of the classical T Tauri star RY Tau has been investigated. During two seasons from 2013 to 2015, we carried out a spectroscopicmonitoring of this star with simultaneous BV R photometry. A correlation between the stellar brightness and the radial velocity of the wind determined from the Hα and Na D line profiles has been found. The irregular stellar brightness variations are shown to be caused by extinction in a dusty disk wind at a distance of about 0.2 AU from the star. We hypothesize that the circumstellar extinction variations result from a cyclic rearrangement of the magnetosphere and coronal mass ejections, which affect the dusty disk wind near the inner boundary of the circumstellar disk.     

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.     

Anisotropic viscous universe with variableG and Λ

An anisotropic model with variableG and and bulk viscosity is considered. The model exhibits an inflationary behavior during which the coefficient of bulk viscosity varies lineraly with the energy density. This allows the anisotropy energy to decrease exponentially with time. Other results overlap with our earlier work with a different ansatz for . The gravitational constant was found to increase during the radiation and matter epochs.