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.     

Topical Problems in the Theory of Planet Formation: Formation of Planetesimals

Izvestiya, Physics of the Solid Earth - We consider the history of the development of the theory of the formation of the Earth and planets over 70 years with special emphasis placed on the...     

Gas Dissipation from the Protosatellite Disk of Jupiter

We consider the dissipation of the gaseous component from the gas–dust accretion disk of Jupiter in which the Galilean satellites were formed. The thermal dissipation of hydrogen and helium is shown to be ineffective. It could ensure the loss of gas only for a low-mass disk and only if the rarefied outer layers of the disk are heated to 10⁴ K. Such a high disk temperature is not reached through Jupiter's radiation in existing models of its formation, but it could be provided by UV radiation of the early Sun after the dissipation of the protoplanetary disk. The viscous dissipation (with a viscosity parameter 10^–3 in the -disk model) related to disk accretion onto Jupiter could disperse a low-mass disk in 10⁷ years. A magnetocentrifugal mechanism, which produced a disk wind during accretion capable of carrying away 0.1 of the accreted gas mass, was probably also involved in the dispersal of the Jovian disk. Differential dispersion, with the loss of only hydrogen and helium and the retention of water vapor and heavier gases in the disk, is possible only in a low-mass disk model. We conclude that the water contained in the Galilean satellites was brought in mainly by solid planetesimals captured into the disk during mutual inelastic collisions in Jupiter's sphere of influence.     

The specific features of the evolution of the viscous protoplanetary circumsolar disk

The problem of angular-momentum and mass transport in the disk is discussed and the disk viscosity is estimated. The evolution of the gas-dust protoplanetary disk at the stage of its formation inside the protostellar (protosolar) accretion envelope is considered. The conditions for the radial growth of the disk are estimated. For the subsequent period, when the central star (young Sun) is in the T Tauri phase, the temporal variations of the radius, mass, and the surface density of the disk, as well as the total mass flux from the disk onto the star (Sun), i.e., the mass accretion rate, are evaluated. The constraints on the initial value of the angular momentum of the protoplanetary circumsolar disk (that is, on the angular momentum of the protosolar cloud) are discussed with due regard for cosmochemical data.Translated from Astronomicheskii Vestnik, Vol. 38, No. 6, 2004, pp. 559–576.Original Russian Text Copyright © 2004 by Makalkin.     

Formation of Planetesimals in the Trans-Neptunian Region of the Protoplanetary Disk

We consider the formation of cometlike and larger bodies in the trans-Neptunian region of the protoplanetary gas–dust disk. Once the particles have reached 1–10 cm in size through mutual collisions, they compact and concentrate toward the midplane of the disk to form a dust subdisk there. We show that after the subdisk has reached a critical density, its inner, equatorial layer that, in contrast to the two subsurface layers, contains no shear turbulence can be gravitationally unstable. The layer breaks up into 10¹²-cm clumps whose small fragments (10⁹ cm) can rapidly contract to form bodies 10 km in size. We consider the sunward drift of dust particles at a velocity that decreases with decreasing radial distance as the mechanism of radial contraction and compaction of the layer that contributes to its gravitational instability and the formation of larger (100 km) planetesimals. Given all of the above processes, it takes 10⁶ yr for planetesimals to form, which is an order of magnitude shorter than the lifetime of the gas–dust protoplanetary disk. We discuss peculiarities of the structure of planetesimals.     

Gravitational Instability in the Dust Layer of a Protoplanetary Disk with Interaction between the Layer and the Surrounding Gas

We consider gravitational instability of the dust layer in the midplane of a protoplanetary disk with turbulence and shear stresses between the gas in the disk and that in the dust layer. We solve a linearized system of hydrodynamic equations for perturbations of dust (monodisperse) and gas phases in the incompressible gas approximation. We take into account the gas drag of solid particles (dust aggregates), turbulent diffusion and the velocity dispersion of particles, and the perturbation of the azimuthal velocity of gas in the layer upon the transfer of angular momentum from solid particles to it and from this gas to the surrounding gas in the disk. We obtain and solve the dispersion equation for the layer with the ratio of surface densities of the dust phase and gas being well above unity. The following parameters of gravitational instability in the dust layer are calculated: the critical surface density of solid matter and the Stokes number of particles corresponding to the onset of instability, the wavelength range in which instability occurs, and the rate of its growth as a function of the perturbation wavelength in the circumsolar disk at radial distances of 1 and 10 AU. We show that at 10 AU, the maximum instability growth rate increases due to the transfer of angular momentum of gas in the layer to gas outside it, a new maximum emerges at a longer wavelength, a long-wavelength instability “tail” forms, and the critical surface density initiating instability decreases relative to that determined without the transfer of angular momentum to gas outside the layer. None of these effects are observed at 1 AU, since instability in this region probably develops faster than the transfer of angular momentum to the surrounding gаs of a protoplanetary disk occurs.     

Gravitational instability in the dust layer of a protoplanetary disk: Interaction of solid particles with turbulent gas in the layer

Gravitational instability of the dust layer formed after the aggregates of dust particles settle toward the midplane of a protoplanetary disk under turbulence is considered. A linearized system of hydrodynamic equations for perturbations of dust (monodisperse) and gas phases in the incompressible gas approximation is solved. Turbulent diffusion and the velocity dispersion of solid particles and the perturbation of gas azimuthal velocity in the layer upon the transfer of angular momentum from the dust phase due to gas drag are taken into account. Such an interaction of the particles and the gas establishes upper and lower bounds on the perturbation wavelength that renders the instability possible. The dispersion equation for the layer in the case when the ratio of surface densities of the dust phase and the gas in the layer is well above unity is obtained and solved. An approximate gravitational instability criterion, which takes the size-dependent stopping time of a particle (aggregate) in the gas into account, is derived. The following parameters of the layer instability are calculated: the wavelength range of its subsistence and the dependence of the perturbation growth rate on the perturbation wavelength in the circumsolar disk at a radial distance of 1 and 10 AU. It is demonstrated that at a distance of 1 AU, the gas–dust disk should be enriched with solids by a factor of 5–10 relative to the initial abundance as well as the particle aggregates should grow to the sizes higher than about 0.3 m in order for the instability to emerge in the layer in the available turbulence models. Such high disk enrichment and aggregate growth is not needed at a distance of 10 AU. The conditions under which this gravitational instability in the layer may be examined with no allowance made for the transfer of angular momentum from the gas in the layer to the gas in a protoplanetary disk outside the layer are discussed.     

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.     

Models of the protosatellite disk of Saturn: Conditions for Titan’s formation

Models of the protosatellite accretion disk of Saturn are developed that satisfy cosmochemical constraints on the volatile abundances in the atmospheres of Saturn and Titan with due regard for the data obtained with the Cassini orbiter and the Huygens probe, which landed on Titan in January 2005. All basic sources of heating of the disk and protosatellite bodies are taken into account in the models, namely, dissipation of turbulence in the disk, accretion of gaseous and solid material onto the disk from the feeding zone of Saturn in the solar nebula, and heating by the radiation of young Saturn and thermal radiation of the surrounding region of the solar nebula. Two-dimensional (axisymmetric) temperature, pressure, and density distributions are calculated for the protosatellite disk. The distributions satisfy the cosmochemical constraints on the disk temperature, according to which the temperature at the stage of the satellite formation ranged from 60–65 K to 90–100 K at pressures from 10^?7 to ?10^?4 bar in the zone of Titan’s formation (according to estimates, r = 20–35R _Sat). Variations of the basic input parameters (the accretion rate onto the protosatellite disk of Saturn from the feeding zone of the planet ?; the parameter α characterizing turbulent viscosity of the disk; and the mass concentration ratio in the solid/gas system) satisfying the aforementioned temperature constraint are found. The spectrum of models satisfying the cosmochemical constraints covers a considerable range of consistent parameters. A model with a rather small flux of ? = 10^?8 M _Sat/ yr and a tenfold depletion of Saturn’s disk in gas due to gas scattering from the solar nebula is at one side of this range. A model with a much higher flux of ? = 10^?6 M _Sat/yr and a hundredfold decrease in opacity of the disk matter owing to decreased concentration of dust particles and/or their agglomeration into large aggregates and sweeping up by planetesimals is at the other side of the range.     

Hydrated Silicates on Edgeworth-Kuiper Objects – Probable Ways of Formation

Visible-range absorption bands at 600–750 nm were recently detected on two Edgeworth-Kuiper Belt (EKB) objects (Boehnhardt et al., 2002). Most probably the spectral features may be attributed to hydrated silicates originated in the bodies. We consider possibilities for silicate dressing and silicate aqueous alteration within them. According to present models of the protoplanetary disk, the temperatures and pressures at the EKB distances (30–50 AU) at the time of formation of the EKB objects (10⁶ to 10⁸ yr) were very low (15–30 K and 10^-9–10^-10 bar). At these thermodynamic conditions all volatiles excluding hydrogen, helium and neon were in the solid state. An initial mass fraction of silicates (silicates/(ices + dust)) in EKB parent bodies may be estimated as 0.15–0.30. Decay of the short-lived ²⁶Al in the bodies at the early stage of their evolution and their mutual collisions (at velocities ≥1.5 km s^-1) at the subsequent stage were probably two main sources of their heating, sufficient for melting of water ice. Because of the former process, large EKB bodies (R ≥ 100 km) could contain a large amount of liquid water in their interiors for the period of a few 10⁶ yr. Freezing of the internal ocean might have begun at ≈ 5 × 10⁶ yr after formation of the solar nebula (and CAIs). As a result, aqueous alteration of silicates in the bodies could occur. A probable mechanism of silicate dressing was sedimentation of silicates with refractory organics, resulting in accumulation of large silicate-rich cores. Crushing and removing icy covers under collisions and exposing EKB bodies' interiors with increased silicate content could facilitate detection of phyllosilicate spectral features.