The Nobeyama Millimeter Array Survey for protoplanetary disks around protostar candidates and T Tauri stars in Taurus

The Nobeyama Millimeter Array Survey for protoplanetary disks has been made for 19 protostellar IRAS sources in Taurus; 13 of them were optically invisible protostars and 6 were young T Tauri stars. We observed 98-GHz continuum and CS(J = 2 – 1) line emissions simultaneously with spatial resolutions of 2 _. 8-8 _. 8 (360-1,200 AU). The continuum emission was detected from 5 out of 6 T Tauri stars and 2 out of 13 protostar candidates: the emission was not spatially resolved and was consistent with being originated from compact circumstellar disks. Extended CS emission was detected around 2 T Tauri stars and 11 protostar candidates. There is a remarkable tendency for the detectability of the 98-GHz continuum emission to be small for protostar candidates. This tendency is explained if the mass of protoplanetary disks around protostars is not as large as that around T Tauri stars; the disk mass may increase with the increase of central stellar mass by dynamical accretion in the course of evolution from protostars to T Tauri stars.Paper presented at the Conference onPlanetary Systems: Formation, Evolution, and Detection held 7–10 December, 1992 at CalTech, Pasadena, California, U.S.A.     

Wave propagation in protoplanetary disks: Formation of twin planets by ‘disk-brown dwarf’ collisions?

In this paper we derive a special linear non-vortical wave propagation solution in the shearing sheet, a model of a compressible two-dimensional fluid system with constant density, constant shear and constant Coriolis force, but without self-gravity. The linear analysis of the shearing sheet leads to a single differential equation for the azimuthal velocity perturbation. A detailed derivation of a special solution with a prescribed azimuthal wavenumber k is presented. More general wave solutions, eventually excited by large local ‘impacts’, can be derived by superimposing all k-modes. The special wave functions so obtained describe the formation of two independent spiral wave arms originating out of a ring-shaped structure. The motivation for this investigation lies in the fact that similar wave propagators can be excited by the transit of a solid or ‘clumpy’ object through a protoplanetary disk. We speculate that a disk-brown dwarf collision can produce in the disk a pair of two spiral density wave fragments triggering the rapid accretion of two giant planets by a gravitational shear instability simultaneously (Hypothesis of a mechanism for the production of giant planets in pairs).     

A submillimeter view of protoplanetary dust disks

We present some results from our submillimeter single-dish and aperture synthesis imaging surveys of protoplanetary disks using the JCMT, CSO, and Submillimeter Array (SMA) on Mauna Kea, Hawaii. Employing a simple disk model, we simultaneously fit the spectral energy distributions and spatially resolved submillimeter continuum emission from our SMA survey to constrain disk structure properties, including surface density profiles and sizes. The typical disk structure we infer is consistent with a fiducial accretion disk model with a viscosity parameter α≈0.01. Combined with a large, multiwavelength single-dish survey of similar disks, we show how these observations provide evidence for significant grain growth and rapid evolution in the outer regions of disks, perhaps due to an internal photoevaporation process. In addition, we discuss SMA observations of the disks in the Orion Trapezium (proplyds) in the context of disk evolution in a more extreme environment.     

Resonant character of the hall instability in protoplanetary disks

We consider nonaxisymmetric magnetosonic oscillations of a radially stratified, weakly ionized protoplanetary disk with a vertical magnetic field. The combined effect of the Hall electric field and the density and magnetic field inhomogeneities present in the disk has been previously predicted to lead to an instability of its small azimuthal perturbations. We revise the previous results and take into account the effect of inhomogeneous ionization of the protoplanetary material related to the inhomogeneity of the disk medium. We show that the instability criterion is governed by three parameters: the magnetic field and ionization fraction gradients and the plasma β. We have found that at high values of β typical of protoplanetary disks, the instability does not manifest itself if the gradients are directed oppositely. In the case of codirectional gradients, the interaction of magnetosonic fluctuations with inhomogeneities of a fixed size is resonant in character, giving rise to an instability in a narrow range of wave numbers.     

Models of Jupiter's growth incorporating thermal and hydrodynamic constraints

We model the growth of Jupiter via core nucleated accretion, applying constraints from hydrodynamical processes that result from the disk-planet interaction. We compute the planet's internal structure using a well tested planetary formation code that is based upon a Henyey-type stellar evolution code. The planet's interactions with the protoplanetary disk are calculated using 3-D hydrodynamic simulations. Previous models of Jupiter's growth have taken the radius of the planet to be approximately one Hill sphere radius, RH. However, 3-D hydrodynamic simulations show that only gas within ∼0.25RH remains bound to the planet, with the more distant gas eventually participating in the shear flow of the protoplanetary disk. Therefore in our new simulations, the planet's outer boundary is placed at the location where gas has the thermal energy to reach the portion of the flow not bound to the planet. We find that the smaller radius increases the time required for planetary growth by ∼5%. Thermal pressure limits the rate at which a planet less than a few dozen times as massive as Earth can accumulate gas from the protoplanetary disk, whereas hydrodynamics regulates the growth rate for more massive planets. Within a moderately viscous disk, the accretion rate peaks when the planet's mass is about equal to the mass of Saturn. In a less viscous disk hydrodynamical limits to accretion are smaller, and the accretion rate peaks at lower mass. Observations suggest that the typical lifetime of massive disks around young stellar objects is ∼3 Myr. To account for the dissipation of such disks, we perform some of our simulations of Jupiter's growth within a disk whose surface gas density decreases on this timescale. In all of the cases that we simulate, the planet's effective radiating temperature rises to well above 1000 K soon after hydrodynamic limits begin to control the rate of gas accretion and the planet's distended envelope begins to contract. According to our simulations, proto-Jupiter's distended and thermally-supported envelope was too small to capture the planet's current retinue of irregular satellites as advocated by Pollack et al. [Pollack, J.B., Burns, J.A., Tauber, M.E., 1979. Icarus 37, 587-611].     

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.     

Numerical experiments on the stability of preplanetary disks

Gravitational stability of gaseous protostellar disks is relevant to theories of planetary formation. Stable gas disks favor formation of planetesimals by the accumulation of solid material; unstable disks allow the possibility of direct condensation of gaseous protoplanets. We present the results of numerical experiments designed to test the stability of thin disks against large-scale, self-gravitational disruption. The disks are represented by a distribution of about 6 × 10⁴ point masses on a two-dimensional (r, φ) grid. The motions of the particles in the self-consistent gravity field are calculated, and the evolving density distributions are examined for instabilities. Two parameters that have major influences on stability are varied: the initial temperature of the disk (represented by an imposed velocity dispersion), and the mass of the protostar relative to that of the disk. It is found that a disk as massive as 1M_⊙, surrounding a 1M_⊙ protostar, can be stable against long-wavelength gravitational disruption if its temperature is about 300°K or greater. Stability of a cooler disk requires that it be less massive, but even at 100°K a stable disk can have an appreciable fraction ($\text{～}\text{1}\text{3}$) of a solar mass.     

On the Dynamics of Dust Grains in a Hierarchical Environment

Most main sequence stars are binaries or higher multiplicity Systems and it appears that at birth most stars have circumstellar disks. It is commonly accepted that planetary systems arise from the material of these disks; consequently, binary and multiple systems may have a main role in planet formation. In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimetre-to-metre sized primordial aggregates. We investigate the response of the particles, in a protoplanetary disk with radius R_D = 100 AU around a solar-like star, to the gravitational field of bound perturbing companions in a moderately wide (300–1600 AU) orbit. For this purpose, we have carried out a series of simulations of coplanar hierarchical configurations using a direct integration code that models gravitational and viscous forces. The massive protoplanetary disk is around one of the components of the binary. The evolution in time of the dust sub-disk depends mainly on the nature (prograde or retrograde) of the relative revolution of the stellar companion, and on the temperature and mass of the circumstellar disk. Our results show that for binary companions near the limit of tidal truncation of the disk, the perturbation leads to an enhanced accretion rate onto the primary, decreasing the lifetime of the particles in the protoplanetary disk with respect to the case of a single star. As a consequence of an enhanced accretion rate the mass of the disk decreases faster, which leads to a longer resultant lifetime for particles in the disk. On the other hand, binary companions may induce tidal arms in the dust phase of protoplanetary disks. Spiral perturbations with m = 1 may increase in a factor 10 or more the dust surface density in the neighbourhood of the arm, facilitating the growth of the particles. Moreover, in a massive disk (0.01M⊙) the survival time of particles is significantly shorter than in a less massive nebula (0.001M⊙) and the temperature of the disk severely influences the spiral-in time of particles. The rapid evolution of the dust component found in post T Tauri stars can be explained as a result of their binary nature. Binarity may also influence the evolution of circumpulsar disks. This revised version was published online in July 2006 with corrections to the Cover Date.     

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 local stability of accretion disk models with considering the role of various viscosity and cooling mechanisms

The standard thin accretion disk model can explain the soft X-ray spectra of Galactic black hole systems and AGN successfully. However, there are still a few observational documents for Radiation pressure theory in X-ray novae in black hole binary systems and AGN. The luminosity in accretion onto black holes is corresponds to L>0.01L _E . According to standard thin disk model, when the accretion rate is over a small fraction of the Eddington rate, L>0.01L _E , the inner region of the disk is radiation-pressure-dominated and thermally unstable. However, observations of the high/soft state of black hole X-ray binaries with luminosity within (0.01L _E <L<0.5L _E ) show that the disk is quite stable. Thus, this contradiction shows the objection of this model and maybe it is essential to change the standard viscosity law or one of the other basic assumptions in order to get a stable disk models. In this paper, we revisit and recalculate the thermal instability with a different models of viscosity and cooling functions and show that the choosing of an arbitrary cooling and viscosity functions can affect on the stability of a general disk model and hence maybe answer to a this problem in accretion disk theory. We choose an arbitrary functions of surface density Σ and half thickness of disk H for cooling and viscosity. Also, we discuss a general disk with thermal conduction, radial force and advection. Then, we solve the equations numerically. We obtain a fourth degree dispersions relation and discuss solutions and instability modes. This analysis shows the great sensitivity of stability of disk to the form of viscosity, so there are various effective factors to stabilize the disk. For example the exist of advection and thermal conduction can effect to stability of disks also.     

Cosmogony of the solar system

In Sections 1–6, we determine an approximate analytical model for the density and temperature distribution in the protoplanetary could. The rotation of the planets is discussed in Section 7 and we conclude that it cannot be determined from simple energy conservation laws.The velocity of the gas of the protoplanetary cloud is found to be smaller by about 5×10³ cm s^–1 in comparison to the Keplerian circular velocity. If the radius of the planetesimals is smaller than a certain limitr ₁, they move together with the gas. Their vertical and horizontal motion for this case is studied in Sections 8 and 9.As the planetesimals grow by accretion their radius becomes larger thanr ₁ and they move in Keplerian orbits. As long as their radius is betweenr ₁ and a certain limitr ₂ their gravitational interaction is negligible. In Section 10, we study the accretion for this case.In Section 11, we determine the change of the relative velocities due to close gravitational encounters. The principal equations governing the late stages of accretion are deduced in Section 12, In Section 13 there are obtained approximate analytical solutions.The effect of gas drag and of collisions is studied in Sections 14 and 15, respectively. Numerical results and conclusions concerning the last and principal stage of accretion are drawn in Section 16.     

Radial self-gravity of accretion disc around supermassive black holes

This paper studies the properties of self-gravity of accretion discs around supermassive black holes. With integration on the thin disc configuration, this paper has calculated the radial and vertical components of self-gravity of accretion discs. The discussion mainly concentrates on the radial component, and the results are briefly as follows: for accretion discs around supermassive black holes (M10 ⁸–10¹⁰ M). At the distance where (R/R _g)10⁵–10⁴, the radial component of self-gravity dominates over the central one where the dynamical structure of the accretion discs completely differs from that of Keplerian disc. A turbulence driven by radial self-gravity instability as a kind of energy source is proposed. This paper has two criteria for the comparison of magnitude between the self-gravity of accretion discs and the gravity of the central black hole, from which an analytic estimation for the outer radius of the accretion discs has been derived. The results of this paper may be used to explain the accretion discs of quasars and AGNs.This research was supported by the National Natural Science Foundation of China.     

Identification of a collapsing protostar

The globular molecular cloud B335 contains a single, deeply embedded, far-infrared source. Our recent observations of H₂CO and CS lines toward this source provide direct kinematic evidence for collapse. Both the intensity and detailed shape of the line profiles match those expected from inside-out collapse inside a radius of 0.036 pc. The collapse began about 1.5 × 10⁵ years ago, similar to the onset of the outflow. The mass accretion rate is about 10 times the outflow rate, and about 0.4M should have now accumulated in the star and disk. Because B335 rotates only very slowly, any disk would still be very small (about 3 AU). The accretion luminosity should be adequate to power the observed luminosity. Consequently, we believe that B335 is indeed a collapsing protostar.Paper presented at the Conference onPlanetary Systems: Formation, Evolution, and Detection held 7–10 December, 1992 at CalTech, Pasadena, California, U.S.A.     

Molecular Line Emission from Accretion Disks Around YSOs

In this work, we model the expected molecular emission from protoplanetary disks, modifying different physical parameters, such as dust grain size, mass accretion rate, viscosity, and disk radius, to obtain observational signatures in these sources. Having in mind possible future observations, we study correlations between physical parameters and observational characteristics. Our aim is to determine the kind of observations that will allow us to extract information about the physical parameters of disks. We also present prospects for molecular line observations of protoplanetary disks, using millimeter and submillimeter interferometers (e.g., SMA or ALMA), based on our results.     

Radial migration of planetesimals

Planetesimals orbiting a protostar in a circumstellar disk are affected by gravitational interaction among themselves and by gas drag force due to disk gas. Within the Kyoto model of planetesimal accretion, the migration rate is interpreted as the inverse of the planetary formation time scale. Here, we study time scales of gravitational interaction and gas drag force and their influence on planetesimal migration in detail. Evaluating observations of 86 T Tauri stars (Beckwithet al., 1990), we find the mean radial temperature profile of circumstellar disks. The disk mass is taken to be 0.01M in accordance with minimum mass models and observed T Tauri disks. The time scale of gravitational interaction between planetesimals is studied analogously to Chandrasekhar's stellar dynamics. Hence, Chandrasekhar's coefficient , defined as the fraction between the mean separation of planetesimals and the impact parameter, plays an important role in determining the migration rate. We find ln to lie between 5 and 10 within the protosolar disk. Our result is that, at the stage of disk evolution considered here, gas drag force affects the radial migration of planetesimals by a few orders of magnitude more than gravitational interaction.Paper presented at the Conference on Planetary Systems: Formation, Evolution, and Detection held 7–10 December, 1992 at CalTech, Pasadena, California, U.S.A.     

The role of toroidal magnetic field in the stability of accretion disks with alpha viscosity and other viscosities

The standard thin accretion disk model predicts that the inner regions of alpha model disks, where radiation pressure is dominant, are thermally and viscously unstable. However, observations show that the bright X-ray binaries and AGN accretion disks, corresponding to radiation-pressure thin disks, are stable. In this paper, we reconsider the linear and local instability of accretion disks in the presence of a toroidal magnetic field. In the basic equations, we consider physical quantities such as advection, thermal conduction, arbitrary viscosity, and an arbitrary cooling function also. A fifth order diffusion equation is obtained and is solved numerically. The solutions are compared to non-magnetic cases. The results show that the toroidal magnetic field can make the thermal instability in radiation pressure-dominated slim disks disappear if ? _m ≥0.3. However, it causes a more thermal instability in radiation pressure alpha disks without advection. Also, we consider the thermal instability in accretion disks with other values of the viscosity and obtain a general criterion for thermal instability in the long-wavelength limit and in the presence of a toroidal magnetic field.