Protoplanetary core formation by rain-out of minerals

Models of giant gaseous protoplanets calculated by DeCampli and Cameron (1979) indicate that iron and probably other minerals in the interior of a planet would be in the liquid state during part of the protoplanet evolution. Liquid drops in a protoplanet would grow by coalescence much as cloud drops in the Earth's atmosphere grow to rain drops. We have modeled this process by using the stochastic collection equation (Slattery, 1978) for various initial conditions. In all of the cases considered, the growth time (to centimeter-sized droplets) is much shorter than the time, as estimated by detailed evolutionary calculations, that the drops are in the liquid state. Brownian collection is effective in quickly coalescing tiny liquid droplets to an average radius of about 0.005 cm with very few drops remaining with radii less than 0.001 cm. For radii larger than 0.005 cm gravitational collection is dominant. Since the particles are rapidly swept from interstellar grain sizes to much larger sizes, the opacity in the cloud layer is expected to drop sharply following melting of the grains.     

原行星盘的研究进展

原行星盘是环绕在年轻星天体(如T Tauri型星,HAe/Be星)周围的气体尘埃盘,是具有初始角动量的分子云核在塌缩形成恒星过程中的自然结果,是行星系统的起源地。原行星盘研究不仅是恒星形成理论的重要组成部分,而且是行星形成理论的基础。首先介绍了盘的形成与演化规律;然后介绍了年轻星天体的能谱分布,盘的模型和参数(质量吸积率、质量、尺度、温度、寿命);随后讨论了尘埃颗粒在盘中生长的观测证据以及行星在盘中形成的大致过程;最后对原行星盘研究的现状和未来做了总结与展望。     

Observational Diagnostics of Gas in Protoplanetary Disks

Protoplanetary disks are composed primarily of gas (99% of the mass). Nevertheless, relatively few observational constraints exist for the gas in disks. In this review, I discuss several observational diagnostics in the UV, optical, near-IR, mid-IR, and (sub)-mm wavelengths that have been employed to study the gas in the disks of young stellar objects. I concentrate in diagnostics that probe the inner 20 AU of the disk, the region where planets are expected to form. I discuss the potential and limitations of each gas tracer and present prospects for future research.     

Towards the Frequency and Evolutionary Timescale of Protoplanetary Disks

In this poster we will present preliminary results from our ongoing near-infrared survey of very young clusters, the majority of them still embedded in their parental molecular cloud. We will draw general conclusions on issues of the frequency and timescales for the initial stages of planet formation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Protoplanetary dust clouds in disks of Herbig Ae/Be stars

The very large brightness decrease of late-type Herbig Ae/Be stars is believed to be caused by obscuring dust clouds orbiting in the outer parts of their circumstellar disks. The distances of the dust clouds to the central stars have been estimated using the wavelength at maximum flux of the excess near-IR radiation, Wien's displacement law, and a formula derived by Rowan-Robinson (1980). The critical masses of these clouds were calculated employing Chandrasekhar's (1943) formula. The minimum size of the dust grains in the obscuring clouds was estimated using Aumannet al.'s (1984) formula they had applied to the star Lyr. However, it can be about ten times smaller if the dust grains are situated at the back of the cloud. The average size of these grains has been determined by assuming a size distribution similar to that in the asteroidal belt (Dohnanyi, 1969) and in the interstellar medium (Mathiset al., 1977). Their number density was determined by means of the extinction power of the dust cloud at theV pass-band. The results of our calculations show that above parameters are similar to those in our solar system. Therefore, we believe that most probably (a) the formation of planetesimals in the circumstellar disks of Herbig Ae/Be stars is on-going; and (b) the obscuring clouds will, in the long run, become planet-like objects.Paper presented at the Conference onPlanetary Systems: Formation, Evolution, and Detection held 7–10 December, 1992 at CalTech, Pasadena, California, U.S.A.     

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.     

The Earth's core formation due to the Rayleigh-Taylor instability

The recent theories of planetary formation lead to a gravitationally unstable structure of the proto-Earth in the accretion stage, which is composed of three layers: an innermost undifferentiated solid core, an intermediate metal-melt layer, and an outermost silicate-melt layer. Taking this configuration as an initial state, we investigate the Earth's core formation due to the Rayleigh-Taylor instability by using the quantitative results on the instability in a self-gravitating fluid sphere obtained from another paper (S. Ida, Y. Nakagawa, and K. Nakagawa, submitted). We find that the instability occurs through the translational mode on a time scale of about 10 hr if the thickness of the metal-melt layer ⪆1 km. This leads to the conclusion that the Earth's core began to form through the translation of the innermost undifferentiated solid core as soon as the outer layer was melted and differentiated in the late accretion stage. In addition, we examine the rotational effects of the instability; the translation occurs most often along the rotational axis. But this preference is weak, since the rotational energy is small compared to the gravitational one.     

Accretion-driven core collapse and the collisional formation of massive stars

We consider the conditions required for a cluster core to shrink, by adiabatic accretion of gas from the surrounding cluster, to densities such that stellar collisions are a likely outcome. We show that the maximum densities attained, and hence the viability of collisions, depend on the balance between core shrinkage (driven by accretion) and core puffing up (driven by relaxation effects). The expected number of collisions scales as     , where N _core is the number of stars in the cluster core and     is the free-fall velocity of the parent cluster (gas reservoir). Thus, whereas collisions are very unlikely in a relatively low-mass, low-internal-velocity system such as the Orion Nebula Cluster, they become considerably more important at the mass and velocity scales characteristic of globular clusters. Thus, stellar collisions in response to accretion-induced core shrinkage remain a viable prospect in more massive clusters, and may contribute to the production of intermediate-mass black holes in these systems.     

Collapse of primordial clouds and the formation of a central core

The collapse of rotating clouds is investigated using three-dimensional self-gravitating hydrodynamical simulations. We take into account the detailed non-equilibrium chemical reactions for primordial gas that consists of pure hydrogen. The parameters of the collapse and the condition of the fragmentation are compared with those of isothermal clouds. It is shown that the geometrical flatness of the central region of the disc is a good indicator for predicting whether the clouds fragment or not. If the flatness is greater than the critical value, ∼ 4π, a cloud fragments into filaments and blobs. On the other hand, if the flatness is smaller than the critical value, fragmentation is not expected before the central core formation even if the cooling is efficient and the total mass becomes much greater than the local Jeans mass at the centre. The critical mass is found to be 3 × 10⁶ M_⊙ for a typical initial condition. If the initial cloud mass is smaller than this critical value, fragmentation before the central core formation is not expected. For a typically estimated first collapsing cosmological baryonic object, M ≲ 10⁶ M_⊙, central core formation is expected and will have a significant effect on the later evolution of the whole system. This revised version was published online in July 2006 with corrections to the Cover Date.     

Unexpected formation modes of the first hard binary in core collapse

The conventional wisdom for the formation of the first hard binary in core collapse is that three-body interactions of single stars form many soft binaries, most of which are quickly destroyed, but eventually one of them survives. We report on direct N-body simulations to test these ideas, for the first time. We find that the assumptions are incorrect in the majority of the cases: (1) quite a few three-body interactions produce a hard binary from scratch; (2) in many cases there are more than three bodies directly and simultaneously involved in the production of the first binary. The main reason for the discrepancies is that the core of a star cluster, at the first deep collapse, contains typically only five or so stars. Therefore, the homogeneous background assumption, which still would be reasonable for, say, 25 stars, utterly breaks down. There have been some speculations in this direction, but we demonstrate this result here explicitly, for the first time.     

The formation of IIE iron meteorites investigated by the chondrule‐bearing Mont Dieu meteorite

A 435 kg piece of the Mont Dieu iron meteorite (MD) contains cm‐sized silicate inclusions. Based on the concentration of Ni, Ga, Ge, and Ir (8.59 ± 0.32 wt%, 25.4 ± 0.9 ppm, 61 ± 2 ppm, 7.1 ± 0.4 ppm, respectively) in the metal host, this piece can be classified as a IIE nonmagmatic iron. The silicate inclusions possess a chondritic mineralogy and relict chondrules occur throughout the inclusions. Major element analysis, oxygen isotopic analysis (Δ¹⁷O = 0.71 ± 0.02‰), and mean Fa and Fs molar contents (Fa_{15.7 ± 0.4} and Fs_{14.4 ± 0.5}) indicate that MD originated as an H chondrite. Because of strong similarities with Netschaëvo IIE, MD can be classified in the most primitive subgroup of the IIE sequence. ⁴⁰Ar/³⁹Ar ages of 4536 ± 59 Ma and 4494 ± 95 Ma obtained on pyroxene and plagioclase inclusions show that MD belongs to the old (~4.5 Ga) group of IIE iron meteorites and that it has not been perturbed by any subsequent heating event following its formation. The primitive character of MD sheds light on the nature of its formation process, its thermal history, and the evolution of its parent body.     

Structure and Instabilities of Irradiated Protoplanetary Disks: Effects of Dust Particle Growth

We have constructed self-consistent temperature and density profiles of irradiated active protoplanetary disks, using a two-dimensional radiative transfer calculation. By means of these profiles we have studied the stabilization of the convective instability by radiative heating and the magnetorotational instability (MRI) via ohmic dissipation, taking into account the effect of dust particle growth. Simple chemistry such as ionization by cosmic rays and recombination on dust grains are used to calculate the ionization degree of gas in the disks. Our results show that the dust growth stabilizes the convective instability due to the 2D effect of radiative transfer, while it enhances the MRI through the decrease in the recombination of ions on the dust grains. In addition, the influences of the dust settling toward the midplane of the disks on the instabilities are discussed.     

Simulation of the Ion Modification Processes in the Silicate Dust of Protoplanetary Disks

The modification of crystalline olivine under irradiation with 1.43-MeV deuterons was studied experimentally. For the first time, the partial destruction of olivine in the surface (100–200 nm) layer was found at a current density of 5–10 A/cm². This effect is caused by ionization and desorption of Mg and Fe atoms at deuteron fluences (3–8 × 10¹⁶cm^–2) which are substantially below the amorphization threshold. It was suggested that the negative charge of the anion group SiO₄ ^4–is neutralized by the interaction with the deuterium ions and Frenkel pairs. The mechanism of sputtering of the crystalline olivine surfaces under irradiation with H, D, He, and Ar ions was studied using of the TRIM 98 program.     

Jeans Instability of a Protoplanetary Gas Cloud with Radiation in Nonextensive Tsallis Kinetics

Solar System Research - In the framework of Tsallis statistics, we study the effect of medium nonextensivity on the Jeans gravitational instability criterion for a self-gravitating protoplanetary...     

Crystallization of magmatic iron meteorites: The role of mixing in the molten core

Abstract— The IIIAB group is the largest of the magmatic iron meteorite groups and consequently is commonly used to test models of asteroid core crystallization. Simple fractional crystallization calculations appear to reproduce the general shape of the elemental trends observed in the IIIAB group when these trends are plotted vs. Ni, as is traditionally done. However, when the elemental trends are examined vs. another element (such as Ge vs. Ir), simple fractional crystallization fails to match a significant portion of the trend, specifically meteorites formed during the final stages of crystallization. Our simple mixing model, which attempts to account for the possibility of inhomogeneities in the molten metallic core, is able to reproduce the entire IIIAB trend observed. This model is a variant of simple fractional crystallization and involves mixing between a zone of liquid involved in the crystallization process and a second zone too far from the crystallizing solid to be actively involved in crystallization. This model does not suggest one unique solution for the method by which an asteroidal core crystallizes; rather it demonstrates that including the effects of mixing in the molten core can account for the observed IIIAB elemental trends, particularly the late-stage crystallizing members, which other models have difficulty explaining.     

On the opaque core appearing in the early phase of star formation

The spherical collapse of a protostar with one solar mass is calculated from a gravitationally unstable initial stage, by solving the equations of hydrodynamics and pherical radiative transfer. The temperature dependence of the dust opacity is taken into account in contrast to the earlier calculations with temperature independent opacities. It is shown that the opaque core is equivalent to the adiabatic core in the purely adiabatic case. The blanketing effect of dust grains strengthens the shock of the opaque core and may result in raising the central entropy of protostars.Paper presented at the Conference on Protostars and Planets, held at the Planetary Science Institute, University of Arizona, Tucson, Arizona, between January 3 and 7, 1978.     

Gravitational collapse of a rotating iron stellar core: A model with Fermi-Dirac neutron degeneracy

The effect of Fermi-Dirac neutron degeneracy on the collapse process is studied. Two computations of a spherically symmetric quasi-one-dimensional model are considered: the case with a full allowance for neutrino deposition and the case with total transparency of the star to neutrino emission. The hindering character of neutron degeneracy with respect to collapse has been revealed quantitatively. A weak sensitivity of both thermodynamics and neutrino characteristics of a rotating collapsar to this effect has been established. However, the hydrostatically equilibrium neutron star obtained is composed of a strongly degenerate neutron gas by more than half. This undoubtedly confirms that the performed computations of a quasi-one-dimensional model with neutron degeneracy are topical.     

Gravitational collapse of a rotating iron stellar core and physical properties of the accompanying neutrino emission

We have ascertained an important role of rotation effects in a collapsing stellar core using a quasi-one-dimensional hydrodynamic model with a rigorous allowance for the neutrino energy losses including the neutrino opacity stage. However, the neutrino scattering processes are not considered in the neutrino emission kinetics as secondary compared to the absorption processes. The quasi-one-dimensional approximation (with averaging of the expression for the centrifugal force over the polar angle) allows numerical calculations to be performed relatively easily up to the formation of a hydrostatically equilibrium neutron star after a very long stage of collapsar cooling by neutrino emission (about 2 s). We present detailed results of our numerical solution, including the neutrino spectra, with electron neutrinos making a dominant contribution to them and the contribution from electron antineutrinos being smaller by an order of magnitude. In the model under consideration, we solve the equation of matter neutronization kinetics by taking into account the main process of nuclear reactions on free nucleons, although the contribution from iron and helium nuclei is included in the equation of state.