Evidence for Outflow-Infall Interactions Terminating Accretion in Protostellar Disks
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.
We report the detection of emission from methanol in a compact source coincident with the position of the L1157 infrared source, which we attribute to molecules in the disk surrounding this young, class 0 protostellar object.
One-dimensional, convective, vertical structure models and one dimensional time-dependent, radial diffusion models are combined to create a self-consistent picture in which FU Orionis outbursts occur in young stellar objects (YSOs) as the result of a large-scale, self-regulated, thermal ionization instability in the surrounding protostellar accretion disk. Although active accretion disks have long been postulated to be ubiqitous among low-mass YSOs, few constraints have until now been imposed on physical conditions in these disks. By fitting the results of time-dependent disk models to observed timescales of FU Orionis events, we estimate the magnitude of the effective viscous stress in the inner disk (r approximately less than 1 AU) to be, in accordance with an ad hoc 'alpha' prescription, the product of the local sound speed, pressure scale height, and an efficiency factor alpha of 10(exp -4) where hydrogen is neutral and 10(exp 3) where hydrogen is ionized. We hypothesize that all YSOs receive infall onto their outer disks which is steady (or slowly declining with time) and that FU Orionis outbursts are self-regulated, disk outbursts which occur only in systems which transport matter inward at a rate sufficiently high to cause hydrogen to be ionized in the inner disk. We estimate a critical mass flux of dm(sub crit)/dt = 5 x 10(exp 7) solar mass/yr independent of the magnitude of alpha for systems with one solar mass, three solar radius central objects. Infall accretion rates in the range of dm(sub in)/dt = 1-10) x 10(exp -6) solar mass/yr produce observed FU Orionis timescales consistent with estimates of spherical molecular cloud core collapse rates. Modeled ionization fronts are typically initiated near the inner edge of the disk and propogate out to a distance of several tens of stellar radii. Beyond this region, the disk transports mass steadily inward at the supplied constant infall rate. Mass flowing through the innermost disk annulus is equal to dm(sub in)/dt only in a time-averaged sense and is regulated by the ionization of hydrogen in the inner disk such that long intervals (approximately 1000 yr) of low-mass flux: (1-30) x 10(exp -8) solar mass/yr are punctuated by short intervals (approximately 100 yr) of high-mass flux: (1-30) x 10(exp -5) solar mass/yr. Timescales and mass fluxes derived for quiescent and outburst stages are consistent with estimates from observations of T Tauri and FU Orionis systems, respectively.
Recently we have shown (sup 1,2) from (sup 12) CO and C (sup 18) O OVRO-MMA observations evidence of infall-outflow interactions in the Young Stellar Object IRS1 in B5. We present new data obtained with the OVRO-MMA in (sup 13) CO(2-1) and (1-0) to trace the full extent of this infall-outflow interaction region.
Hydrodynamical calculations of the collapse of an axisymmetric, rotating one solar mass protostellar cloud, including the effects of radiative transfer and radiative acceleration but without magnetic fields, are presented. The results include calculations of infrared protostellar spectra as a function of time and viewing angle. A numerical algorithm involving explicit nested grids is used to resolve the region of initial disk formation and at the same time to include the outer regions in the calculation. The central part of the protostar is modeled approximately. Initial conditions are systematically varied to investigate their influence on the evolution and final configuration of central star plus circumstellar disk. The initial state for the standard case is a centrally condensed molecular cloud core of one solar mass with a mean density of 8 x 10 exp -18 g/cu cm and a specific angular momentum at the outer edge of 7 x 10 exp 20 sq cm/s. The collapse is followed for 8 x 10 exp 4 yr, at which point 0.45 solar mass is contained in a rapidly rotating central object and most of the remainder in a surrounding equilibrium disk. The stability of this final structure is qualitatively analyzed.
We consider the origin and intital growth of the disks that form around protostars during the collapse of rotating molecular cloud cores. These disks are assumed to be inviscid and pressure free, and to have masses small compared to those of their central stars. We find that there exist three distinct components-an outer disk, in which shocked gas moves with comparable azimuthal and radical velocities; and inner disk, where material follows nearly circular orbits, but spirals slowly toward the star because of the drag exerted by adjacent onfalling matter, and a turbulent ring adjoining the first two regions. Early in the evolution, i.e., soon after infalling matter begins to miss the star, only the outer disk is present, and the total mass acceration rate onto the protostar is undiminished. Once the outer disk boundary grows to more than 2.9 times the stellar radius, first the ring, and then the inner disk appear. Thereafter, the radii of all three components expand as t(exp 3). The mass of the ring increase with time and is always 13% of the total mass that has fallen from the cloud. Concurrently with the buildup of the inner disk and ring, the accretion rate onto the star falls off. However, the protostellar mass continue to rise, asymptotically as t(exp 1/4). We calculated the radiated flux from the inner and outer disk components due to the release of gravitational potential energy. The flux from the inner disk is dominant and rises steeply toward the stellar surface. We also determine the surface temperature of the inner disk as a function of radius. The total disk luminosity decreases slowly with time, while the contributions from the ring and inner disk both fall as t(exp -2).
The objective is to obtain quantitative information on the turbulent transport of mass, angular momentum, and energy under the conditions that characterize the solar nebula, by direct numerical calculations. These calculations were made possible by research conducted on supercomputers (Cray XMP and Cray 2) by the Ames Computational Fluid Dynamics Branch. Techniques were developed that permitted the accurate representation of turbulent flows over the full range of important eddy sizes. So far, these techniques were applied (and verified) primarily in mundane laboratory situations, but they have a strong potential for astrophysical applications. A sequence of numerical experiments were conducted to evaluate the Reynold's stress tensor, turbulent heat transfer rate, turbulent dissipation rate, and turbulent kinetic energy spectrum, as functions of position, for conditions relevant to the solar nebula. Emphasis is placed on the variation of these properties with appropriate nondimensional quantities, so that relations can be derived that will be useful for disk modeling under a variety of hypotheses and initial conditions.
The formation and evolution of an accretion disk formed during the collapse of a rotating cloud core are considered. The effect of the usual 'turbulent' alpha-viscosity as well as the effective viscosity due to self-gravitation in the disk are taken into account. For observed values of the cloud's initial rotation rates, and for reasonable estimates of the efficiency of the various processes, it is found that the disks so formed are large (radii several hundred to a few thousand AU), relatively massive (comparable in mass to the central star), and long-lived (tens of infall time scales).
Models of the solar nebula reveal that it might have been gravitationally unstable, both early and later in its evolution. Such instabilities produce density waves and associated gravitational torques, which are potent agents of angular momentum transport. In previous work, we conducted a series of numerical simulations designed to quantify the effects of gravitational instabilities in a generalizable way (Tomley, Cassen, & Steiman-Cameron 1991). Here we present a second series of simulations in which we examine disks of greater size, increased star/disk mass ratio, and flatter surface density distribution than those in our initial study. The purpose is to represent disks at a later stage of evolution than those already studied, to test the quantitative relations derived in our earlier work and to explore the effects of mass ratio on the results. The new results indicate that the tendencies for unstable, uncooled disks to heat to stability and for dynamical evolution rates to be proportional to cooling rates are general characteristics of the behavior of gravitationally unstable disks. Nevertheless, there are quantitative, and (for strong cooling) even qualitative differences that are revealed in the new simulations, particularly with regard to the cooling rates at which clumping tends to occur.
Recent observations of binary, pre-main-sequence, solar-type stars provide evidence that such systems may coexist with circumstellar disks. The binary disk systems, besides being of general interest for the study of star formation, potentially provide useful tests of companion-disk interaction theories prominent in current hypotheses of planet formation. In this paper, we apply an asymptotic analysis of the nonlinear, resonant interaction of a stellar companion with a disk to understand the dependence of such interactions on the properties of the system: the binary mass ratio, the physical properties of the disk, and the effective dissipation (treated herein as viscosity). The method is based on a WKBJ approximation and exploits the conditions that the disk is thin and much less massive than the primary, but does not require that the companion-induced disturbance be small. Both isothermal and adiabatic responses are treated. Only circular orbit resonances are considered in this paper. It is demonstrated that the temperature of the disk as well as the relative mass of the companion affects the degree of nonlinearity, and that nonlinearity promotes high wave compression ratios, long wavelengths, and increased propagation distances. Nevertheless, the total torque exerted between the companion and the disk is well represented by linear theory. The amplitudes of density disturbances are reduced by viscosity and nonisothermality. Because resonant interactions are generally strong and capable of driving rapid evolution, one might expect observations of systems undergoing strong, resonant-driven evolution to be rare. In this connection, it is pointed out that the m = 1 resonance is distinguished by being anomalously weaker than the others and is therefore of observational interest. It is speculated that, in conditions of intrinsically small dissipation, the propagation of resonant-driven density waves is limited by the tendency of their wavelength to diminish with distance, and that the propagation distance (and therefore the region of the disk to which angular momentum is redistributed) is set by the distance at which the wavelength becomes comparable to the disk thickness.
We have examined the physical and chemical processes which accompany the propagation of a J-type shock through molecular astrophysical gas of density 10(exp 7.5) - 10(exp 12) H nuclei per cu cm. Our study is primarily relevant to the high-density accretion shocks which are associated with the supersonic infall of material during the collapse of a molecular cloud core to form a protostar. We have carried out a general parameter study to determine the characteristic infrared emission spectrum of dense molecular shocks and the extent of grain destruction within such shocks as a function of relevant shock parameters. Accretion shocks with mass accretion rates approximately greater than 10(exp -5) solar masses per y are plausible sources of the CO v = 2-0 vibrational band emission observed in many protostars, at least for sources of moderate CO v = 2-0 luminosity (approximately less than few x 10(exp 23) W). Rovibrational emissions from H2O and OH are also predicted to be luminous in dense molecular accretion shocks. Vaporization is the dominant mechanism of grain destruction within dense molecular shocks, leading to the efficient removal of those grains which are heated to their vaporization temperature by the combined effects of radiative heating, drag heating, and thermal heating by the hot postshock gas. The criteria for grain vaporization are obtained for grains composed of several different materials: metallic iron, silicate, troilite, refractory organic materials, volatile organic materials, and water ice.
Hydrodynamical calculations of the evolution of a collapsing, rotating axisymmetric 10 solar masses molecular clump, including the effects of radiative acceleration but without magnetic fields, are represented. The initial cloud is assumed to be uniformly rotating, centrally condensed sphere with rho is proportional to r(exp -2). Several cases are considered, in which both the overall clump size and the total amount of angular momentum are varied. The calculations show how a warm, quasi-hydrostatic disk surrounding a central unresolved core of only a few solar masses forms and grows in size and mass. The disk is encased in two distinct accretion shock fronts, both of which are several scale heights above the equatorial plane. At the end of the calculation of our standard case, the central unresolved region is found to have a mass of 2.7 solar masses and a ratio of rotational to gravitational energy of approximately 0.45, sufficiently large to be unstable to nonaxisymmetric perturbations. In addition, the inner portions of the disk containing most of the mass are unstable according to the local Toomre criterion, implying that also in this region nonaxisymmetric perturbations will lead to rapid evolution. Under the assumption that gravitational torques would transport angular momentum out of this region, a central core of less than or approximately 8 solar masses with a stable disk of greater than or approximately = 2 solar masses should result. Frequency-dependent radiative transfer calculations of the standard case at selected ages show how the continuum spectrum of the structure depends on the disk's orientation and age and how the observed isophotal contours vary with wavelength. Because of the strong dependence on viewing angle, continuum spectra alone should not be used to estimate the evolutionary stage of development of these objects. Comparable results were obtained for the other cases considered.
The results of the time-dependent disk models developed in Bell & Lin are compared with observed properties of FU Orionis variables. Specific models are fit to the light curves of Fu Ori, V1515 Cyg, and V1057 Cyg. The slow risetime of V1515 Cyg can be matched by a self-regulated outburst model. The rapid risetimes of FU Ori and V1057 Cyg can be fitted with the application of modest perturbations to the disk surface density. Model disks display spectral features characteristic of observed objects. The color evolution of V1057 Cyg is naturally explained if mass flux drops in the inner disk (r less than 1/4 AU) while remaining steady in the outer disk. The decrease in optical line width (rotational velocity) observed during the decay of V1057 Cyg may be accounted for by an outward-propagating ionization front. We predict that before final decay to the quiescent phase, short-wavelength line widths (lambda less than 1.5 microns) will again increase. It is suggested that FU Orionis outbursts primarily occur to systems during the embedded phase with ages less than several times 10(exp 5) yr.
We examine the infrared emission of the Herbig Ae/Be stars and show that some possess characteristics indicative of partially crystalline grains similar to those seen in Beta Pictoris and some solar system comets.
We model the thermal balance, the chemistry, and the radiative transfer in dusty disks orbiting young, low mass stars. These models are motivated by observations of infrared and ultraviolet transitions of H2 from protoplanetary disks, as well as millimeter and submillimeter observations of other molecules such as CO, and infrared continuum observations of the dust. The dust grains are heated primarily by the stellar radiation and the infrared radiation field produced by the dust itself. The gas is heated by collisions with warmer dust grains, X-rays from the region close to the stellar surface, UV pumping of hydrogen molecules, and the grain photoelectric heating mechanism initiated by UV photons from the central star. We treat cases where the gas to dust ratio is high, because the dust has settled to the midplane and coagulated into relatively large objects. We discuss situations in which the infrared emission from H2 can be detected, and how the comparison of the observations with our models can deduce physical parameters such as the mass and the density and temperature distribution of the gas.
This grant was intended for, and has been used for, the support of the specialized CalFUSE data reduction required for these non-standard observations. The goal of the program was to search for FUV continuum radiation from T-Tauri stars, which would have affected the chemistry in the proto-planetary disk. Because of the low flux expected, the standard background subtraction method would not be appropriate. Rather on-chip background determination was expected to be required.