Investigations of the dynamical evolution of protoplanetary nebulae

The spectral energy distributions (SED's) of all the known Class I (protostellar) sources in the Taurus molecular cloud were modelled. The Tereby, Shu, & Cassen (1984, TSC) density distribution for a rotating, infalling envelope was adopted. The radiative equilibrium temperature distribution from the spherical average of the TSC density distribution was calculated. The resulting spherically-symmetric temperature distribution then provides the source function to obtain the emergent spectrum at a given inclination angle i from the formal solution of the transfer equation, using the exact density (opacity) distribution. Results showed that the SED's of the protostar candidates in Taurus can be reproduced with TSC models having infall rates close to the values predicted by the theory of isothermal cloud collapse. Flat Spectrum T Tauri Stars were studied. The mid- to far-infrared fluxes of 'flat spectrum' T Tauri stars can be explained by radiative equilibrium emission from infalling dusty envelopes. This explanation indicated that models employing 'active' disks, in which the temperature distribution is a parameterized power law, should be invoked with caution. Infall also naturally explains the scattered light nebulae detected around many flat-spectrum sources. Results showed that the SED's of the protostar candidates in Taurus can be reproduced with TSC models having infall rates close to the values predicted by the theory of isothermal cloud collapse. Flattened infalling envelope models are also being studied

Infalling Envelopes and Pre-Main Sequence Disks

The goal of this project is to understand the observed infrared emission of young stellar objects, and explore the implications of this emission for the evolution of dusty envelopes and circumstellar disks. We are using sophisticated radiative transfer methods to compare models with observations, thereby making critical tests of the standard picture of low-mass star formation

A Preliminary Study of the Orion Nebula Cluster Structure and Dynamics

We use optical and near-infrared star counts to explore the structure and dynamics of the Orion Nebula Cluster (ONC). This very young (<1 Myr) cluster is not circularly symmetric in projection but is elongated north-south in a manner similar to the molecular gas distribution in the region, suggesting that the stellar system may still reflect the geometry of the protocluster cloud. Azimuthally averaged stellar source counts compare well with simple spherically symmetric, single-mass King cluster models. The model fits suggest that the inner Trapezium region should be regarded as the core of the ONC, not as a distinct entity as sometimes advocated. We estimate that the core radius of the cluster is 0.16-0.21 pc and that the central stellar density approaches 2 × 10^4 stars pc^(-3). Adopting the stellar velocity dispersion from published proper-motion studies, virial equilibrium would require a total mass within about 2 pc of the Trapezium of ~4500 M_☉, slightly more than twice the mass of the known stellar population and comparable to the estimated mass in molecular gas projected onto the same region of the sky. If ≳ 20% of the remaining molecular gas is converted into stars, thus adding to the binding mass, given that the present stellar population alone has a total energy close to zero, the ONC is likely to produce a gravitationally bound cluster. The ONC also exhibits mass segregation, with the most massive (Trapezium) stars clearly concentrated toward the center of the cluster and some evidence for the degree of central concentration to decrease with decreasing mass down to 1-2 M_☉, as would be expected for general mass segregation. Given the extreme youth of the stars compared with the estimated range of collisional relaxation times, the mass segregation is unlikely to be the result of cluster relaxation. Instead, we suggest that the mass segregation reflects a preference for higher mass stars to form in dense, central cluster regions

Advanced flight control system study

A fly by wire flight control system architecture designed for high reliability includes spare sensor and computer elements to permit safe dispatch with failed elements, thereby reducing unscheduled maintenance. A methodology capable of demonstrating that the architecture does achieve the predicted performance characteristics consists of a hierarchy of activities ranging from analytical calculations of system reliability and formal methods of software verification to iron bird testing followed by flight evaluation. Interfacing this architecture to the Lockheed S-3A aircraft for flight test is discussed. This testbed vehicle can be expanded to support flight experiments in advanced aerodynamics, electromechanical actuators, secondary power systems, flight management, new displays, and air traffic control concepts

Resolved Depletion Zones and Spatial Differentiation of N2H+ and N2D+

We present a study on the spatial distribution of N2D+ and N2H+ in thirteen protostellar systems. Eight of thirteen objects observed with the IRAM 30m telescope show relative offsets between the peak N2D+ (J=2-1) and N2H+ (J=1-0) emission. We highlight the case of L1157 using interferometric observations from the Submillimeter Array and Plateau de Bure Interferometer of the N2D+ (J=3-2) and N2H+ (J=1-0) transitions respectively. Depletion of N2D+ in L1157 is clearly observed inside a radius of ~2000 AU (7") and the N2H+ emission is resolved into two peaks at radii of ~1000 AU (3.5"), inside the depletion region of N2D+. Chemical models predict a depletion zone in N2H+ and N2D+ due to destruction of H2D+ at T ~ 20 K and the evaporation of CO off dust grains at the same temperature. However, the abundance offsets of 1000 AU between the two species are not reproduced by chemical models, including a model that follows the infall of the protostellar envelope. The average abundance ratios of N2D+ to N2H+ have been shown to decrease as protostars evolve by Emprechtinger et al., but this is the first time depletion zones of N2D+ have been spatially resolved. We suggest that the difference in depletion zone radii for N2H+ and N2D+ is caused by either the CO evaporation temperature being above 20 K or an H2 ortho-to-para ratio gradient in the inner envelope.Comment: Accepted to ApJ. 44 pages 13 Figure

Cooling, Gravity and Geometry: Flow-driven Massive Core Formation

We study numerically the formation of molecular clouds in large-scale colliding flows including self-gravity. The models emphasize the competition between the effects of gravity on global and local scales in an isolated cloud. Global gravity builds up large-scale filaments, while local gravity -- triggered by a combination of strong thermal and dynamical instabilities -- causes cores to form. The dynamical instabilities give rise to a local focusing of the colliding flows, facilitating the rapid formation of massive protostellar cores of a few 100 M$_\odot$. The forming clouds do not reach an equilibrium state, though the motions within the clouds appear comparable to ``virial''. The self-similar core mass distributions derived from models with and without self-gravity indicate that the core mass distribution is set very early on during the cloud formation process, predominantly by a combination of thermal and dynamical instabilities rather than by self-gravity.Comment: 13 pages, 12 figures, accepted by Ap