Planetary migration to large radii

There is evidence for the existence of massive planets at orbital radii of several hundred au from their parent stars where the time-scale for planet formation by core accretion is longer than the disc lifetime. These planets could have formed close to their star and then migrated outwards. We consider how the transfer of angular momentum by viscous disc interactions from a massive inner planet could cause significant outward migration of a smaller outer planet. We find that it is in principle possible for planets to migrate to large radii. We note, however, a number of effects which may render the process somewhat problematic.     

The origin of Jovian and Saturnian satellites in accretion disks

The properties of gas-dust disks that surrounded Jupiter and Saturn during the final stage of their formation are analyzed. The sizes of the disks are determined by the total planetocentric angular momentum of the matter accreted by planets and correspond to the sizes of the orbits of their largest satellites. The mass of the solid component of disks is limited from below by the total mass of the Galilean satellites of Jupiter (no less than 4 × 10²⁶ g) and the mass of the largest Saturnian satellites (1.4 × 10²⁶ g), whereas the mass of the gaseous component is limited from above by the amount of hydrogen and helium that could have been later lost by the disks. Our analysis of the known mechanisms of dissipation of gas showed that its simultaneous content in the disks relative to the solid component was much lower than the corresponding gas-to-solid ratio in the Sun. A certain amount of solid compounds (including ice) could have been brought into the disks with planetesimals, which had undergone mutual collisions in the neighborhood of giant planets and served as germs of satellites. The bulk of solid matter appears to have been captured into disks when the latter were crossed by smaller and intermediate-sized planetesimals, which then became parts of the satellites.     

Computer simulation of the formation of asteroid belt structure

The system of two gravitational centers with variable separation between components one of which (the primary) loses its mass onto another (the secondary) is investigated under condition of total mass and angular momentum conservation. When the primary/secondary mass ratio becomes about that of Jupiter/Sun the small bodies ejected with the gaseous matter through the inner Lagrange point from the Roche lobe of the primary form a ring similar to the asteroid belt of the solar system. The formation of ring structure is calculated by numerical integration of Newtonian equations of N-body problem in orbital plane of the gravitational centers. The results are compared with the planar subsystem of the asteroid belt. The presence of the main gaps in the distribution of their mean motions at 2/1, 3/1, 5/2 and some other commensurabilities with the primary mean motion is found. More fine details of the belt structure are obtained, e.g. the gap asymmetry and a qualitative agreement with the eccentricity distribution. Within the scope of the same model the external part of the ring is investigated all the pairwise interactions being included. The clustering of bodies near 3/2 commensurability isolated from the main belt by the wide gap centered at 5/3 commensurability is obtained. It is supposed that the ring structure and the interplanetary spacing law for the terrestrial planets are due to the same mechanism.     

Distribution of Matter in the Solar System

A new formula for the distribution of matter in the solar system is derived by assuming that the planets were formed from trapped particles of a cosmic dust disk attached to the Sun. Contrary to Boltzmann's distribution which predicts thermal collapse of this cloud on the Sun, it is found that if the primeval particles move on circular orbits according to Kepler's law, then their velocities obey a 2-D global Maxwellian and their distribution in space is given by p ₀ (r)=(α r ²)\exp(-α r) (Km^-1); α = 888.73 × 10⁶ Km. The form ofp ₀ (r) agrees with the observed mass distribution of the planets and explains their present large angular momentum. PACS numbers: 96.35.Cp, 96.35.Fs This revised version was published online in July 2006 with corrections to the Cover Date.     

SOME MORPHOLOGIC SYSTEMATICS OF COMPLEX IMPACT STRUCTURES

Similarities among impact structures on different planets and satellites suggest that the cratering process transcends variations in both target and impactor. In particular, impact may control the spacing of concentric rings, if not their actual emplacement. In at least four respects the scaled horizontal dimensions of complex meteorite-impact structures on Earth resemble those of multi-ring basins and large craters on the Moon, Mars, Mercury, and some outer satellites: (1) Base diameter of the (topographic) central peak is a constant 20% to 25% of the rim diameter in small complex craters; (2) it averages only half as much in large structures that also have concentric rings; (3) the inner ring of a two-ring crater lacking a central peak is half the diameter of the outer ring; (4) adjacent rings of complex craters that have more than two concentric rings are spaced at a constant interval of about (2.0 ± 0.2)^0.5 D, both inside and outside the main ring. Two minor differences in morphology suggest that uniquely terrestrial conditions may control some horizontal dimensions of meteorite craters: (1) the inner ring of a two-ringed structure that also has a central peak is 0.5X the diameter of the outer, not 0.4X as it is for peak-plus-ring basins on the planets; and (2) two-ring and multi-ring meteorite craters occupy the same size range, whereas on planets most two-ring basins are smaller than multi-ring basins.     

On the formation of Jupiter and Saturn

It is proposed that Jupiter and Saturn were initially formed as small rocks which grew into their present sizes as a result of accretion of matter from the gas-dust cloud surrounding the Sun. The energy released by the accretion of Jupiter and Saturn is computed. It is concluded that the ‘excess’ radiation from these planets is due to simple cooling and that the gravitational contraction from initially extended states most probably never occurred.     

Bulges versus discs: the evolution of angular momentum in cosmological simulations of galaxy formation

We investigate the evolution of angular momentum in simulations of galaxy formation in a cold dark matter universe. We analyse two model galaxies generated in the N -body/hydrodynamic simulations of Okamoto et al. Starting from identical initial conditions, but using different assumptions for the baryonic physics, one of the simulations produced a bulge-dominated galaxy and the other one a disc-dominated galaxy. The main difference is the treatment of star formation and feedback, both of which were designed to be more efficient in the disc-dominated object. We find that the specific angular momentum of the disc-dominated galaxy tracks the evolution of the angular momentum of the dark matter halo very closely: the angular momentum grows as predicted by linear theory until the epoch of maximum expansion and remains constant thereafter. By contrast, the evolution of the angular momentum of the bulge-dominated galaxy resembles that of the central, most bound halo material: it also grows at first according to linear theory, but 90 per cent of it is rapidly lost as pre-galactic fragments, into which gas had cooled efficiently, merge, transferring their orbital angular momentum to the outer halo by tidal effects. The disc-dominated galaxy avoids this fate because the strong feedback reheats the gas, which accumulates in an extended hot reservoir and only begins to cool once the merging activity has subsided. Our analysis lends strong support to the classical theory of disc formation whereby tidally torqued gas is accreted into the centre of the halo conserving its angular momentum.     

Embedded star clusters and the formation of the Oort cloud: III. Evolution of the inner cloud during the Galactic phase

In a previous publication [Brasser, R., Duncan, M.J., Levison, H.F., 2006. Icarus 184, 59-82], models of the inner Oort cloud were built which included the effect of an embedded star cluster on cometary orbits about the Sun. The main conclusions of that paper were that the formation efficiency is about 10% and the median distance of the cloud to the Sun only depends on the mean density of gas and stars the Sun encountered. Here we report on the results of simulations which followed the ensuing dynamical evolution of these comet clouds in the current Galactic environment once the Sun left the embedded star cluster. The goal is to determine whether or not the dynamical influence of passing Galactic field stars and the Galactic tidal field is sufficient to replenish the current outer cloud (semi-major axis a>20,000 AU) with enough material from the inner cloud (a<20,000 AU). Since visible new comets come directly from the outer cloud, a mass estimate only exists for the latter, with a lower limit of 1 M_⊕ [Francis, P.J., 2005. Astrophys. J. 635, 1348-1361]. Knowing the amount of expansion of the inner cloud may therefore yield an estimate of the mass of said (unseen) inner cloud. Our results indicate that typically only 10% of the comets from the inner cloud land in the outer cloud and are bound after 4.5 Gyr. If one assumes that in the extreme case all or the majority of the current population of the outer cloud has come from the inner cloud, then a typical value of the mass of the inner cloud is about 10 M_⊕. The results of [Brasser, R., Duncan, M.J., Levison, H.F., 2006. Icarus 184, 59-82] showed that ∼10% of comets from the Jupiter-Saturn region were implanted in the inner Oort cloud, which implies an uncomfortably large value of about 100 M_⊕ for the mass of solids in the primordial Jupiter-Saturn region. This extreme case might be remedied in two says: either the effect of Giant Molecular Cloud complexes on the inner Oort cloud must be much more severe than originally thought, or there was a two-stage formation process for the Oort cloud, in which the outer cloud was largely populated by comets scattered once the Sun had left its primordial birth cluster.     

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 Leonard Award Address Presented 1998 July 27, Dublin,Ireland: On the difficulties of making Earth-like planets

Abstract— Here I discuss the series of events that led to the formation and evolution of our planet to examine why the Earth is unique in the solar system. A multitude of factors are involved: These begin with the initial size and angular momentum of the fragment that separated from a molecular cloud; such random factors are crucial in determining whether a planetary system or a double star develops from the resulting nebula. Another requirement is that there must be an adequate concentration of heavy elements to provide the 2% “rock” and “ice” components of the original nebula. An essential step in forming rocky planets in the inner nebula is the loss of gas and depletion of volatile elements, due to early solar activity that is linked to the mass of the central star. The lifetime of the gaseous nebula controls the formation of gas giants. In our system, fine timing was needed to form the gas giant, Jupiter, before the gas in the nebula was depleted. Although Uranus and Neptune eventually formed cores large enough to capture gas, they missed out and ended as ice giants. The early formation of Jupiter is responsible for the existence of the asteroid belt (and our supply of meteorites) and the small size of Mars, whereas the gas giant now acts as a gravitational shield for the terrestrial planets. The Earth and the other inner planets accreted long after the giant planets, from volatile-depleted planetesimals that were probably already differentiated into metallic cores and silicate mantles in a gas-free, inner nebula. The accumulation of the Earth from such planetesimals was essentially a stochastic process, accounting for the differences among the four rocky inner planets—including the startling contrast between those two apparent twins, Earth and Venus. Impact history and accretion of a few more or less planetesimals were apparently crucial. The origin of the Moon by a single massive impact with a body larger than Mars accounts for the obliquity (and its stability) and spin of the Earth, in addition to explaining the angular momentum, orbital characteristics, and unique composition of the Moon. Plate tectonics (unique among the terrestrial planets) led to the development of the continental crust on the Earth, an essential platform for the evolution of Homo sapiens. Random major impacts have punctuated the geological record, accentuating the directionless course of evolution. Thus a massive asteroidal impact terminated the Cretaceous Period, resulted in the extinction of at least 70% of species living at that time, and led to the rise of mammals. This sequence of events that resulted in the formation and evolution of our planet were thus unique within our system. The individual nature of the eight planets is repeated among the 60-odd satellites—no two appear identical. This survey of our solar system raises the question whether the random sequence of events that led to the formation of the Earth are likely to be repeated in detail elsewhere. Preliminary evidence from the “new planets” is not reassuring. The discovery of other planetary systems has removed the previous belief that they would consist of a central star surrounded by an inner zone of rocky planets and an outer zone of giant planets beyond a few astronomical units (AU). Jupiter-sized bodies in close orbits around other stars probably formed in a similar manner to our giant planets at several astronomical units from their parent star and, subsequently, migrated inwards becoming stranded in close but stable orbits as “hot Jupiters”, when the nebula gas was depleted. Such events would prevent the formation of terrestrial-type planets in such systems.     

Resonant structures in the solar system

The resonance theory is discussed with respect to the Solar System with a view to show that every triad of successive planets in the Solar System follows Laplace's resonance relation. With rings now known to exist around three of the four major planets, scientists have begun to speculate about the possible existence of ring structure and one or two small planets going around the Sun itself. It is also believed that the ring systems may exist around the planets Neptune and Mars. In this paper an attempt is made to provide a basis to these beliefs using Laplace's resonance relation. The triads of successive innermost objects (rings and/or satellites) in the satellite — systems of Jupiter, Saturn and Uranus are also shown to follow Laplace's resonance relation.     

Numerical models of the primitive solar nebula

Numerical models have been constructed to represent probable conditions in the primitive solar nebula. A two solar mass fragment of a collapsing interstellar gas cloud has been represented by a uniformly rotating sphere. Two cases have been considered: one in which the internal density of the sphere is uniform and the other in which the density falls linearly from a central value to zero at the surface (the uniform and linear models). These assumptions served to define the distribution of angular momentum per unit mass with mass fraction. The spheres were flattened into disks, and models of the disks were found in which there was a force balance in the radial and vertical directions, subject to certain approximations, and with everywhere the assigned values of angular momentum per unit mass. The radial pressure gradient of the gas was included in the force balance. The energy transport in the vertical direction involved convection and radiative equilibrium; the principal contributors to opacity at lower temperatures were metallic iron grains and ice. The models contained two convection zones, an inner one due to the dissociation of hydrogen molecules, and an outer one in which there was a high opacity due to metallic iron grains. The characteristic semithickness of the disks ranged from about 0.1 astronomical units near the center to about one astronomical unit near the exterior. Characteristic angular momentum transport times and radiation lifetimes for these models of the initial solar nebula were estimated. Both types of characteristic lifetime were as short as a few years near the inner part of the models, and became about 10⁴ years or longer at distances greater than ten astronomical units.     

On the formation of protostellar disks

An analysis is presented of the hydrodynamic aspects of the growth of protostellar disks from the accretion (or collapse) of a rotating gas cloud. The size, mass, and radiative properties of protostellar disks are determined by the distribution of mass and angular momentum in the clouds from which they are formed, as well as from the dissipative processes within the disks themselves. The angular momentum of the infalling cloud is redistributed by the action of turbulent viscosity on a shear layer near the surface of the disk (downstream of the accretion shock) and on the radial shear across cylindrical surfaces parallel to the rotation axis. The fraction of gas that is fed into a central core (protostar) during accretion depends on the ratio of the rate of viscous diffusion of angular momentum to the accretion rate; rapid viscous diffusion (or a low accretion rate) promotes a large core-to-disk mass ratio. The continuum radiation spectrum of a highly viscous disk is similar to that of a steady-state accretion disk without mass addition. It is possible to construct models of the primitive solar nebula as an accretion disk, formed by the collapse of a slowly rotating protostellar cloud, and containing the minimum mass required to account for the planets. Other models with more massive disks are also possible.     

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.     

Reassessing the formation of the inner Oort cloud in an embedded star cluster

We re-examine the formation of the inner Oort comet cloud while the Sun was in its birth cluster with the aid of numerical simulations. This work is a continuation of an earlier study (Brasser, R., Duncan, M.J., Levison, H.F. [2006]. Icarus 184, 59–82) with several substantial modifications. First, the system consisting of stars, planets and comets is treated self-consistently in our N-body simulations, rather than approximating the stellar encounters with the outer Solar System as hyperbolic fly-bys. Second, we have included the expulsion of the cluster gas, a feature that was absent previously. Third, we have used several models for the initial conditions and density profile of the cluster – either a Hernquist or Plummer potential – and chose other parameters based on the latest observations of embedded clusters from the literature. These other parameters result in the stars being on radial orbits and the cluster collapses. Similar to previous studies, in our simulations the inner Oort cloud is formed from comets being scattered by Jupiter and Saturn and having their pericentres decoupled from the planets by perturbations from the cluster gas and other stars. We find that all inner Oort clouds formed in these clusters have an inner edge ranging from 100 AU to a few hundred AU, and an outer edge at over 100,000 AU, with little variation in these values for all clusters. All inner Oort clouds formed are consistent with the existence of (90377) Sedna, an inner Oort cloud dwarf planetoid, at the inner edge of the cloud: Sedna tends to be at the innermost 2% for Plummer models, while it is 5% for Hernquist models. We emphasise that the existence of Sedna is a generic outcome. We define a ‘concentration radius’ for the inner Oort cloud and find that its value increases with increasing number of stars in the cluster, ranging from 600 AU to 1500 AU for Hernquist clusters and from 1500 AU to 4000 AU for Plummer clusters. The increasing trend implies that small star clusters form more compact inner Oort clouds than large clusters. We are unable to constrain the number of stars that resided in the cluster since most clusters yield inner Oort clouds that could be compatible with the current structure of the outer Solar System. The typical formation efficiency of the inner Oort cloud is 1.5%, significantly lower than previous estimates. We attribute this to the more violent dynamics that the Sun experiences as it rushes through the centre of the cluster during the latter’s initial phase of violent relaxation.     

Rings Around Planets, Atmospheric Superrotation, and Their Great Spots

A common explanation is offered of the facts that the (outer) planets Jupiter, Saturn, Uranus, and Neptune (i) have ring systems inside their magnetospheres, (ii) show alternating atmospheric super- and sub-rotation in latitude belts, and (iii) have great (coloured) whirling spots in their atmospheres, at latitudes of maximal shear flow. The common reason, so we argue, is the action of magnetic torques between the various ring systems in non-synchronous rotation which drive electric currents, help ionize the orbiting gas, and redistribute angular momentum. Very similar reasoning has been used earlier – though incorrect in detail – to explain the complicated system of torsional oscillations in the solar convection zone.     

The circumnuclear material in the Galactic Centre: a clue to the accretion process

On the basis of 'sticky particle' calculations, it is argued that the gas features observed within 10 pc of the Galactic Centre — the circumnuclear disc (CND) and the ionized gas filaments — as well as the newly formed stars in the inner 1 pc can be understood in terms of tidal capture and disruption of gas clouds on low angular momentum orbits in a potential containing a point mass. The calculations demonstrate that a dissipative component forms a 'dispersion ring', an asymmetric elliptical torus precessing counter to the direction of rotation, and that this shape can be maintained for many orbital periods. For a range of plausible initial conditions, such a structure can explain the morphology and kinematics of the CND and of the most conspicuous ionized filament. While forming the dispersion ring, a small cloud with low specific angular momentum is drawn into a long filament which repeatedly collides with itself at high velocity. The compression in strong shocks is likely to lead to star formation even in the near tidal field of the point mass. This process may have general relevance to accretion on to massive black holes in normal and active galactic nuclei.     

Collapse and fragmentation of rotating magnetized clouds – II. Binary formation and fragmentation of first cores

Subsequent to Paper I, the evolution and fragmentation of a rotating magnetized cloud are studied with use of three-dimensional magnetohydrodynamic nested grid simulations. After the isothermal runaway collapse, an adiabatic gas forms a protostellar first core at the centre of the cloud. When the isothermal gas is stable for fragmentation in a contracting disc, the adiabatic core often breaks into several fragments. Conditions for fragmentation and binary formation are studied. All the cores which show fragmentation are geometrically thin, as the diameter-to-thickness ratio is larger than 3. Two patterns of fragmentation are found. (1) When a thin disc is supported by centrifugal force, the disc fragments into a ring configuration (ring fragmentation). This is realized in a rapidly rotating adiabatic core as  Ω > 0.2τ⁻¹_ff  , where Ω and  τ_ff  represent the angular rotation speed and the free-fall time of the core, respectively. (2) On the other hand, the disc is deformed to an elongated bar in the isothermal stage for a strongly magnetized or rapidly rotating cloud. The bar breaks into 2–4 fragments (bar fragmentation). Even if a disc is thin, the disc dominated by the magnetic force or thermal pressure is stable and forms a single compact body. In either ring or bar fragmentation mode, the fragments contract and a pair of outflows is ejected from the vicinities of the compact cores. The orbital angular momentum is larger than the spin angular momentum in the ring fragmentation. On the other hand, fragments often quickly merge in the bar fragmentation, since the orbital angular momentum is smaller than the spin angular momentum in this case. Comparison with observations is also shown.     

Conservation laws and mass distribution in the planet formation process

Within the framework of the nebular theory of the origin of the solar system, conservation laws are applied to the condensation of a ring shaped cloud of orbiting particles. The final configuration is assumed to be a point-like planet in a circular orbit around the Sun. On this ground, it is possible to relate the masses of the planets with the interplanetary distances. This relation is confirmed satisfactorily by the observed masses and orbital radii of several planets and satellites of the solar system.     

Shock fragmentation model for gravitational collapse

A cloud of gas collapsing under gravity will fragment.We present a new theory for this process,in which layers of shocked gas fragment due to their gravitational instability.Our model explains why angular momentum does not inhibit the collapse process.The theory predicts that the fragmentation process produces objects which are significantly smaller than most stars,implying that accretion onto the fragments plays an essential role in determining the initial masses of stars.This prediction is also consistent...