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.     

Structure and evolution of isolated giant gaseous protoplanets

The structure and evolution of isolated giant gaseous protoplanets in the mass range 0.3 to 4.5 Jovian masses is investigated. Under the assumptions of the calculations, the following properties are found: (1) The central region of protoplanets of mass less than about 1 Jovian mass is, at some evolutionary epoch, thermodynamically favorable to the liquification of major interstellar grain constituents. Grains in this region can grow and infall to form a planetary core in tens to hundreds of years. (2) All protoplanets studied are convective through-out most of their interior. This property is in contrast to Bodenheimer's fully radiative proto-Jupiter models. We attribute the difference to the use of improved opacities. The presence of convection has at least two important consequences. First, it can mix grains into the central regions during planetary core formation, possibly allowing a core of mass ～ 1 Earth mass to grow. Second, convection can transport angular momentum outward as the protoplanet quasi-statically contracts. (3) The thermal contraction time depends sensitively on the surface opacity (T < 200°K). This opacity is uncertain within a factor of 5. The contraction times imply that some protoplanets can remain stable against tidal disruption by the proto-Sun and solar nebula during core-forming stages.     

Accumulation processes in the primitive solar nebula

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

Temperature distribution in the solar nebula at successive stages of its evolution

We have constructed a model of the solar nebula that allows for the temperature and pressure distributions at various stages of its evolution to be calculated. The mass flux from the accretion envelope to the disk and from the disk to the Sun, the turbulent viscosity parameter α, the opacity of the disk material, and the initial angular momentum of the protosun are the input model parameters that are varied. We also take into account the changes in the luminosity and radius of the young Sun. The input model parameters are based mostly on data obtained from observations of young solar-type stars with disks. To correct the input parameters, we use the mass and chemical composition of Jupiter, as well as models of its internal structure and formation that allow constraints to be imposed on the temperature and surface density of the protoplanetary disk in Jupiter’s formation zone. Given the derived constraints on the input parameters, we have calculated models of the solar nebula at successive stages of its evolution: the formation inside the accretion envelope, the evolution around the young Sun going through the T Tauri stage, and the formation and compaction of a thin dust layer (subdisk) in the disk midplane. We have found the following evolutionary trend: an increase in the temperature of the disk at the stage of its formation, cooling at the T Tauri stage, and the subsequent internal heating of the dust subdisk by turbulence dissipation that causes a temperature rise in the formation zone of the terrestrial planets at the high subdisk density and the opacity in this zone. We have obtained the probable ranges of temperatures in the disk midplane, i.e., the temperatures of the protoplanetary material in the formation region of the terrestrial planets at the initial stage of their formation.     

The distributions and ages of refractory objects in the solar nebula

Refractory objects such as Calcium, Aluminum-rich Inclusions, Amoeboid Olivine Aggregates, and crystalline silicates, are found in primitive bodies throughout our Solar System. It is believed that these objects formed in the hot, inner solar nebula and were redistributed during the mass and angular momentum transport that took place during its early evolution. The ages of these objects thus offer possible clues about the timing and duration of this transport. Here we study how the dynamics of these refractory objects in the evolving solar nebula affected the age distribution of the grains that were available to be incorporated into planetesimals throughout the Solar System. It is found that while the high temperatures and conditions needed to form these refractory objects may have persisted for millions of years, it is those objects that formed in the first 10⁵ years that dominate (make up over 90%) those that survive throughout most of the nebula. This is due to two effects: (1) the largest numbers of refractory grains are formed at this time period, as the disk is rapidly drained of mass during subsequent evolution and (2) the initially rapid spreading of the disk due to angular momentum transport helps preserve this early generation of grains as opposed to later generations. This implies that most refractory objects found in meteorites and comets formed in the first 10⁵ years after the nebula formed. As these objects contained live ²⁶Al, this constrains the time when short-lived radionuclides were introduced to the Solar System to no later than 10⁵ years after the nebula formed. Further, this implies that the t=0 as defined by meteoritic materials represents at most, the instant when the solar nebula finished accreting significant amounts of materials from its parent molecular cloud.     

Protostars and the origin of the angular momentum of the solar system

Protostars in a group exert gravitational tidal torques on an aspherical nebula located in the group. The net torque transfers angular momentum from the orbital motions of the stars to rotation of the nebula. A relation can be derived between the parameters describing the protostars and the final angular momentum of the nebula. While the parameters concerned are uncertain, a conservative choice results in a value for the angular momentum equal to about 1/3 of that of the present solar system. This suggests that if the Sun formed in a group, tidal interactions with other protostars may account for a significant part of the angular momentum of the solar system.     

Formation of solar nebula and mass distribution in the planetary system

The formation of the solar nebula and the distribution of mass in its planetary system is studied. The underlying idea is that the protosun, fragmented out from an interstellar cloud as a result of cluster formation, gathered the planetary material and, hence, spin angular momentum by gravitational accretion during its orbital motion around the centre of the Galaxy. The study gives the initial angular momentum of the solar nebula nearly equal to the present value of the solar system.     

Accretion disks around Jupiter and Saturn at the stage of regular satellite formation

Modern models of the formation of the regular satellites of giant planets, constructed with consideration for their structure and composition suggest that this process lasted for a considerable period of time (0.1–1 Myr) and developed in gas-dust circumplanetary disks at the final stage of giant planet formation. The parameters of protosatellite disks (e.g., the radial distribution of surface density and temperature) serve as important initial conditions for such models. Therefore, the development of protosatellite disk models that take into account currently known cosmochemical and physical restrictions remains a pressing problem. It is this problem that is solved in the paper. New models of the accretion disks of Jupiter and Saturn were constructed with consideration for the disk heating by viscous dissipation of turbulent motions, by accretion of material from the surrounding region of the solar nebula, and by radiation from the central planets. The influence of a set of input model parameters (the total rate of mass infall onto the disk, the turbulent viscosity and opacity of disk material, and the centrifugal radius of the disk) on thermal conditions in the accretion disks was studied. The dependence of opacity on temperature and the abundance and size of solid particles present in the disk was taken into account. Those constructed models that satisfy the existing constraints limit the probable values of input parameters (primarily rates of mass infall onto the disks of Jupiter and Saturn at the stage of regular satellite formation and, to a lesser extent, the disk opacities). Constraints on the location of the regions of formation of the major satellites of Jupiter and Saturn are suggested based on the constructed models and simple analytical estimates concerning the formation of satellites in the accretion disks. It is shown that Callisto and Titan could hardly be formed at significantly greater distances from their planets.     

Dust concentration and chondrule formation

Meteoritical and astrophysical models of planet formation make contradictory predictions for dust concentration factors in chondrule-forming regions of the solar nebula. Meteoritical and cosmochemical models strongly suggest that chondrules, a key component of the meteoritical record, formed in regions with solids-to-gas mass ratios orders above the solar nebula average. However, models of dust grain dynamics in protoplanetary disks struggle to surpass concentration factors of a few except during very short-lived stages in a dust grain's life. Worse, those models do not predict significant concentration factors for dust grains the size of chondrule precursors. We briefly develop the difficulty in concentrating dust particles in the context of nebular chondrule formation and show that the disagreement is sufficiently stark that cosmochemists should explore ideas that might revise the concentration factor requirements downward.     

Computer simulations of planetary accretion dynamics: Sensitivity to initial conditions

We have tested the implications and limitations of Program ACRETE, a scheme based merely on Newtonian physics and accretion with unit sticking efficiency, devised by Dole in 1970 to simulate the origin of the planets. The dependence of the results on a variety of radial and vertical density distribution laws, on the ratio of gas to dust in the solar nebula, on the total nebula mass, and on the orbital eccentricity, ?, of the accreting grains are explored. Only for a small subset of conceivable cases are planetary systems closely like our own generated. Many models have tendencies toward one of two preferred configurations: multiple-star systems, or planetary systems in which Jovian planets either have substantially smaller masses than in our system or are absent altogether. But for a wide range of cases recognizable planetary systems are generated, ranging from multiple-star systems with accompanying planets, to systems with Jovian planets at several hundred astronomical units, to single stars surrounded only by asteroids. Many systems exhibit planets like Pluto and objects of asteroidal mass, in addition to usual terrestrial and Jovian planets. No terrestrial planets were generated more massive than five Earth masses. The number of planets per system is for most cases of order 10, and, roughly, inversely proportional to ?. All systems generated obey a relation of the Titius-Bode variety for relative planetary spacing. The case with which planetary systems are generated using such elementary and incomplete physical assumptions supports the idea of abundant and morphologically diverse planetary systems throughout the Galaxy.     

The interaction between giant gaseous protoplanets and the primitive solar nebula

The manner in which a giant gaseous protoplanet becomes embedded in the primitive solar nebula determines surface boundary conditions which must be used in studying the evolution of such objects. On the one hand, if the system resembles a contact binary system, then the envelope of the protoplanet should approach the entropy of the surrounding nebula. On the other hand angular. momentum transfer by resonance and tidal effects between the nebula and the protoplanet may cause the nebula to exhibit a zone of avoidance near the protoplanet, thus inhibiting exchange of material. This problem has been studied with a computer program developed by D. N. C. Lin which simulates disk hydrodynamics by particle motions with dissipation. These studies suggest that for expected values of the protoplanet/protosun mass ratios, significant inhibition of mass exchange is likely, so that it is a reasonable next step to undertake protoplanet evolution studies with the imposition of minimum protoplanet surface temperatures.     

Physics of the primitive solar accretion disk

The theory of viscous accretion disks developed by Lynden-Bell and Pringle has been applied to the evolution of the primitive solar nebula. The additional physical input needed to determine the structure of the disk is described. A series of calculations was carried out using a steady flow approximation to explore the effects on the disk properties of variations in such parameters as the angular momentum and accretion rate of the infalling material from a collapsing interstellar cloud fragment. The more detailed evolutionary calculations involved five cases with various combinations of parameters. It was concluded that the late stages of evolution of the disks would be dominated by the effects of mass loss from the expansion of a hot disk corona into space, and the effects of this were included in the evolutionary calculations. A new theory of comet formation is formulated upon these results. The most important result is the conclusion, which appears to be inescapable, that the primitive solar accretion disk was repeatedly unstable against axisymmetric perturbations, in which rings would form and collapse upon themselves, with the subsequent formation of giant gaseous protoplanets.     

Calculations of the effects of angular momentum on the early evolution of Jupiter

The planet Jupiter is assumed to have formed as a subcondensation in the solar nebula. The initial phase of its evolution is one of hydrostatic contraction with radiative energy transport. Calculations of evolutionary sequences through this phase are presented, including the effects of angular momentum. The calculations are carried out in two space dimensions under the assumptions of axial symmetry, constancy of angular velocity on cylindrical surfaces about the rotation axis, a pressure-density relation given by the polytrope of index 3, conservation of angular momentum, and a homogeneous composition. The results show that under certain physically reasonable initial distributions of density and angular momentum the formation of a central planet and a rotating circumplanetary envelope is possible, while under assumptions a point of instability is reached that probably results in the breakup of the condensation by fission into two or more parts. The models are discussed with reference to the present angular momenta of Jupiter and its regular satellites.     

Evolution of grains in a turbulent solar nebula

In accretion disk models of the solar nebula, turbulence is driven by convective instability. This mechanism requires high opacity, which must be provided by solid grains. Evolution of the grain size distribution in a turbulent disk is computed numerically, using realistic collisional outcomes and strengths of grain aggregates, rather than an arbitrary “sticking efficiency.” The presence of turbulence greatly increases the rate of grain collisions; the coagulation rate is initially much greater than in a nonturbulent disk. Aggregates quickly reach sizes ～0.1–1 cm, but erosion and breakup in collisions prevent growth of larger bodies for plausible aggregate impact strengths. These aggregates are too small to settle to the plane of the disk, and planetesimal formation is impossible as long as the turbulence persists. However, the opacity of the disk is reduced by aggregate formation; some combinations of opacity law and surface density produce an optically thin disk, cutting off turbulent convection. The disk may experience alternating periods of turbulence and quiescence, as grains are depleted by coagulation and replenished by infall from the presolar cloud. Planetesimals can form only during the quiescent intervals; it is argued that such episodes were rare during the lifetime of the accretion disk.     

太阳系起源和演化理论的研究

通过角动量守恒计算,证明了原始星云角动量不足,单纯靠星云自转惯性离心力无法抗衡中心部位星云的吸引力,无法在星云赤道处形成星云盘.原始星云角动量不足,同时星云收缩时径向方向速度不等,内快外慢,结果中心部位星云形成太阳,外部赤道部位星云物质因赶不上内部星云物质收缩而掉队形成星云盘.再由星云盘分裂、掉队形成星云环;星云环形成行星、卫星.对太阳系一些主要特征,作了分析和说明.     

A large-eddy simulation of turbulent compressible convection: differential rotation in the solar convection zone

We present the results of two simulations of the convection zone, obtained by solving the full hydrodynamic equations in a section of a spherical shell. The first simulation has cylindrical rotation contours (parallel to the rotation axis) and a strong meridional circulation, which traverses the entire depth. The second simulation has isorotation contours about mid-way between cylinders and cones, and a weak meridional circulation, concentrated in the uppermost part of the shell.
We show that the solar differential rotation is directly related to a latitudinal entropy gradient, which pervades into the deep layers of the convection zone. We also offer an explanation of the angular velocity shear found at low latitudes near the top. A non-zero correlation between radial and zonal velocity fluctuations produces a significant Reynolds stress in that region. This constitutes a net transport of angular momentum inwards, which causes a slight modification of the overall structure of the differential rotation near the top. In essence, the thermodynamics controls the dynamics through the Taylor–Proudman momentum balance . The Reynolds stresses only become significant in the surface layers, where they generate a weak meridional circulation and an angular velocity 'bump'.     

Theoretical,observational, and isotopic estimates of the lifetime of the solar nebula

Abstract— There are a variety of isotopic data for meteorites which suggest that the protostellar nebula existed and was involved in making planetary materials for some 10⁷ yr or more. Many cosmochemists, however, advocate alternative interpretations of such data in order to comply with a perceived constraint, from theoretical considerations, that the nebula existed only for a much shorter time, usually stated as ≤ 10⁶ yr. In this paper, we review evidence relevant to solar nebula duration which is available through three different disciplines: theoretical modelling of star formation, isotopic data from meteorites, and astronomical observations of T Tauri stars. Theoretical models based on observations of present star-forming regions indicate that stars like the Sun form by dynamical gravitational collapse of dense cores of cold molecular clouds in the interstellar medium. The collapse to a star and disk occurs rapidly, on a time scale of the order 10⁵ yr. Disks evolve by dissipating energy while redistributing angular momentum, but it is difficult to predict the rate of evolution, particularly for low mass (compared to the star) disks which nonetheless still contain enough material to account for the observed planetary system. There is no compelling evidence, from available theories of disk structure and evolution, that the solar nebula must have evolved rapidly and could not have persisted for more than 1 Ma. In considering chronologically relevant isotopic data for meteorites, we focus on three methodologies: absolute ages by U-Pb/Pb-Pb, and relative ages by short-lived radionuclides (especially ²⁶Al) and by evolution of ⁸⁷Sr/⁸⁶Sr. Two kinds of meteoritic materials-refractory inclusions such as CAIs and differentiated meteorites (eucrites and angrites)—appear to have experienced potentially dateable nebular events. In both cases, the most straightforward interpretations of the available data indicate nebular events spanning several Ma. We also consider alternative interpretations, particularly the hypothesis of radically heterogeneous distribution of ²⁶Al, which would avoid these chronological interpretations. The principal impetus for such alternative interpretations seems to be precisely the obviation of the chronological interpretation (i.e., the presumption rather than the inference of a short (≤1 Ma) lifetime of the nebula). Astronomical observations of T Tauri stars indicate that the presence of dusty disks is a common if not universal feature, that the disks are massive enough to accomodate a planetary system such as ours, and that at least some persist for 10⁷ yr or more. The results are consistent with the time scales inferred from the meteorite isotopic data. They cannot be considered conclusive with regard to solar nebula time scales, however, in part because it is difficult to relate disk observations to processes that affect meteorites, and in part because the ages assigned for these stars could be wrong by a factor of several in either direction. We conclude that the balance of available evidence favors the view that the nebula existed and was active for at least several Ma. However, because the evidence is not definitive, it is important that the issue be perceived to be an open question, whose answer should be sought rather than presumed.     

Simulation of turbulent convection in a slowly rotating red giant star

The first 3-D non-linear hydrodynamical simulation of the inner convective envelope of a rotating low mass red giant star is presented. This simulation, computed with the ASH code, aims at understanding the redistribution of angular momentum and heat in extended convection zones. The convection patterns achieved in the simulation consist of few broad and warm upflows surrounded by a network of cool downflows. This asymmetry between up and downflows leads to a strong downward kinetic energy flux, that must be compensated by an overluminous enthalpy flux in order to carry outward the total luminosity of the star. The influence of rotation on turbulent convection results in the establishment of largescale mean flows: a strong radial differential rotation and a single cell poleward meridional circulation per hemisphere. A detailed analysis of angular momentum redistribution reveals that the meridional circulation transports angular momentum outward in the radial direction and poleward in the latitudinal direction, with the Reynolds stresses acting in the opposite direction. This simulation indicates that the classical hypothesis of mixing length theory and solid-body rotation in the envelope of red giants assumed in 1-D stellar evolution models are unlikely to be realized and thus should be reconsidered. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The first ten million years in the solar nebula

Abstract— The ²⁶Al/²⁷Al ratio in a large number of calcium-aluminum inclusions (CAIs) is a rather uniform 5 × 10^?5, whereas in chondrules the ratio is either undetectable or has a much lower value; the simplest interpretation of this is that there was an interval of a few million years between the times that these two meteoritic constituents formed stable solids. The present investigation was undertaken as an exploration of the physics of the processes in the solar nebula during and after the accumulation of the Sun. Understanding the time scales of events in this nebular model, to see if this would cast light on this apparent CAI to chondrule time interval, was the major motivation for the exploration. There were four stages in the history of the solar nebula; in stage 1, a fragment of an interstellar molecular cloud collapsed to form the Sun and solar nebula; in stage 2, the nebula was in approximate steady state balance between infall from the cloud and accretion onto the Sun and was in its FU Orionis accumulation stage; in stage 3, the Sun had been mainly accumulated and there was a slow residual mass flow into the Sun while it was in its classical T Tauri stage; and in stage 4, the nebula had finished accreting material onto the Sun (now a weak-lined T Tauri star) and was in a static condition with no significant dissipation or motions, other than removal at the inner edge due to the T Tauri solar wind and photoevaporation beyond 9 astronomical units (AU). It is found that the energy source keeping the nebula warm during stages 3 and 4 is recombination of ionized H in the ionized bipolar jets and the T Tauri coronal expansion solar wind. The parameters of the heating model were adjusted to locate the ice sublimation line at 5.2 AU. In this work, a nebular model is used with a surface density of 4.25 × 10³ gm/cm² at 1 AU and a variation with radial distance as the inverse first power. Under normal conditions in the nebula, there is a negative pressure gradient that provides partial radial support for the gas, which thus circles the Sun more slowly than large solid objects do. Large objects undergo a slow inward spiral due to the gas drag; very small objects move essentially with the gas but have a slow inward drift; and intermediate objects (e.g., 1 m) have a fairly large inward drift velocity that traverses the full radial extent of the nebula in considerably less than the CAI to chondrule time interval. Such objects are thus lost unless they can grow rapidly to larger sizes. Near the inner edge (bow) of the nebula during stage 4, the pressure gradient becomes positive, creating a narrow zone of zero gas drag toward which solids drift from both directions, facilitating planetesimal formation in the inner solar nebula. Recent theoretical and experimental results on sticking probabilities of solids show that icy surfaces have the best sticking properties, but icy interstellar grains can only stick together when subjected to impact velocities of less than 2000 cm/sec. However, if the solid objects are very underdense, then a collision leads to interpenetration and many points at which the small constituent grains can adhere to one another, and thus coagulation becomes possible for such underdense objects. Simulations were made of such coagulation in the outer solar nebula, and it was found that the central plane of the nebula quickly becomes filled with meter-sized and larger bodies that rapidly accumulated near the top of the nebula and rapidly descended; in a few thousand years this quickly leads to gravitational instabilities that can form planetesimals. These processes led to the rapid formation of Jupiter in the nebula (and the slightly less rapid formation of the other giant planets). The early formation of Jupiter opens an annular gap in the nebula, and thus a second region is created in the nebula with zero gas drag. It is concluded that CAIs were formed at the end of stage 2 of the nebula history and moved out into the nebula for long-term storage, and that most chondrules were formed by magnetic reconnection flares in the bow region of the nebula during stage 4, several million years later. Carbonaceous meteorites should be formed on the far side of the Jovian gap, with the chondrules being heated by flares on the early Jupiter irradiating materials in the nearby zone of zero gas drag, and they should have essentially the same ²⁶Al ages as the CAIs (this will be very hard to confirm owing to scarcity of Al mineral phases in these chondrules).     

The effect of multiple particle sizes on cooling rates of chondrules produced in large‐scale shocks in the solar nebula

Chondrules represent one of the best probes of the physical conditions and processes acting in the early solar nebula. Proposed chondrule formation models are assessed based on their ability to match the meteoritic evidence, especially experimental constraints on their thermal histories. The model most consistent with chondrule thermal histories is passage through shock waves in the solar nebula. Existing models of heating by shocks generally yield a good first‐order approximation to inferred chondrule cooling rates. However, they predict prolonged heating in the preshock region, which would cause volatile loss and isotopic fractionation, which are not observed. These models have typically included particles of a single (large) size, i.e., chondrule precursors, or at most, large particles accompanied by micron‐sized grains. The size distribution of solids present during chondrule formation controls the opacity of the affected region, and significantly affects the thermal histories of chondrules. Micron‐sized grains evaporate too quickly to prevent excessive heating of chondrule precursors. However, isolated grains in chondrule‐forming regions would rapidly coagulate into fractal aggregates. Preshock heating by infrared radiation from the shock front would cause these aggregates to melt and collapse into intermediate‐sized (tens of microns) particles. We show that inclusion of such particles yields chondrule cooling rates consistent with petrologic and isotopic constraints.