SEARCH · Engineering Papers
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.
In recent years the nearby cold, dark clouds have been shown to possess a rich chemistry, with interesting differences with respect to warmer massive-star-forming regions and also among the cold clouds themselves. Thirty-nine molecular species are now known in these regions. Recent molecular detections and upper limits in dark clouds are discussed, with particular emphasis on the tricarbon species C3O, C3H, and C3H2.
The collapse and fragmentation of rotating protostellar clouds is explored, starting from nonuniform density and nonuniform rotation initial conditions. Whether binary fragmentation occurs during the first dynamic collapse phase depends strongly on the initial density profile. Exponential clouds are only somewhat more resistant to fragmentation than uniform-density clouds, but power-law clouds do not undergo fragmentation for likely values of a relevant parameter. Because binary fragments start from profiles intermediate between uniform density and exponential clouds, minimum protostellar mass for population I stars should be increased to approximately 0.02 solar mass. The axisymmetric Terey et al. (1984) model should be stable with respect to nonaxisymmetric perturbations. Considering the observed binary frequency, collapse from power-law initial conditions appears to be less common than collapse from more uniform initial conditions.
The chemical evolution of oxygen and sulfur species in shocked dense clouds is studied. Reaction rate constants for several important neutral reactions are examined, and revised values are suggested. The one-fluid magnetohydrodynamic shock structure and postshock chemical evolution are calculated for shocks of velocity v(s) = 10 km/s through clouds of initial number density n(0) = 100,000/cu cm and of molecule/atom ratios H2/H = 10, 1000, and 100,000 with most sulfur contained initially in molecules SO2 and SO. Abundances of SO2, SO, CS, and OCS remain near their preshock values, except in clouds containing substantial amounts of atomic hydrogen, where significant destruction of sulfur-oxygen species occurs. Abundances of shock-enhanced molecules HS and H2O are sensitive to the molecule/atom ratio. Nonthermal oxygen-hydrogen chemistry has a minor effect on oxygen-sulfur molecules in the case H2/H = 10.
The 1(10)-1(01) transition of C3HD was detected at 19.418 GHz at twelve positions in cold, dark clouds, and the D hyperfine components were resolved in two sources (L1498 and TMC-1C) well enough to derive values for the D quadrupole coupling constants. Simultaneous observations of C3H2 in each source yield relative integrated line intensities in the range 0.10-0.18, from which relative C3HD/C3H2 abundances in the range 0.05-0.15 were derived. These are among the highest deuteration ratios yet observed. Within the limits of the observational and modeling uncertainties, it is possible to explain the derived C3HD/C3H2 ratios by ion-molecule chemistry if e is equal to about 3 x 10 to the -7th.
New data for C2 toward X Per, HD 206267, HD 207198, and Gamma Cep, for CH and CN toward X Per, and for CO toward HD 207198 have been obtained. The column densities of CH, C2, CN, and CO toward the stars in the Cepheus OB2 association are similar to reddened directions in Perseus and in Ophiuchus, indicating a similarity in physical conditions for the foreground clouds. The available data for other directions have been analyzed and the resulting data applied to study the transition from a photochemical regime to a chemical regime. The data for N(CN), N(C2), and N(CO) have been plotted against N(CH) to elucidate the chemistry of carbon-bearing molecules more clearly. The observed trends for CN and C2 suggest a change in slope at N(CH) of roughly 5 x 10 to the 13th/sq cm. Below this value, photodestruction is predicted to dominate and the slope is determined by the photochemistry. For directions with more N(CH), a linear correlation consistent with destruction by chemical reactions is expected.
The effects of the presence of a substantial component of large molecules on the chemistry of diffuse molecular clouds are explored, and detailed models of the zeta Persei and zeta Ophiuchi clouds are constructed. The major consequence is a reduction in the abundances of singly charged atomic species. The long-standing discrepancy between cloud densities inferred from rotational and fine-structure level populations and from the ionization balance can be resolved by postulating a fractional abundance of large molecules of 1 x 10 to the -7th for zeta Persei and 6 x 10 to the -7th for zeta Ophiuchi. If the large molecules are polycyclic aromatic hydrocarbons (PAH) containing about 50 carbon atoms, they contain 1 percent of the carbon in zeta Persei and 7 percent in zeta Ophiuchi. Other consequences of the possible presence of PAH molecules are discussed.
This paper investigates the effect of variations in the Rosseland mean opacity of dust grains on numerical models of three-dimensional protostellar collapse and fragmentation. In particular, it is found that increasing the dust grain opacity by factors of three to four has little effect upon the gross characteristics of protostellar fragmentation. Consequently, theoretical quantities such as the estimated minimum protostellar mass for Population I star formation are insensitive to the precise value of the opacity.
This paper explores fluorescence and reflection models of the companion to SN 1987A obseved by speckle interferometry, recalling a 1901 precedent. The apparent small angular size of the companion is a severe constraint. A fluorescence model cannot reach the observed brightness unless the ultraviolet burst from the supernova contained as many as 2 x 10 to the 58th ionizing photons. This is about 25 times stronger than generous current models. Even then, the expected line ratios and widths do not fit the observations. The absence of narrow H-alpha and H-beta lines in the supernova spectrum, the ratio of fluxes of the companion in H-alpha and forbidden N II line filters, the invisibility of the companion at 4861 (H-beta), and its detection at 5330 fail to agree with theory. A dust-reflection model is more promising, and the color can be reddened by the evaporation of small grains, but the model still falls more than about 1 mag short in brightness. Furthermore, a dust reflection should have increased in relative brightness in May-June 1987, rather than disappearing as the mystery spot did. If all the observations are correct, neither model is likely to work.
The results are presented of a variety of spherically symmetric one-dimensional (1D) calculations intended to determine the robustness of the dynamical hiccup phenomenon in protostellar cores. The 1D models show that the phenomenon is relatively insensitive to changes in the equations of state, numerical resolution, initial density and temperature, and the radiative transfer approximation. In 1D, the hiccup results in an explosive destruction of the entire inner protostellar core. Inner core formation is studied with a sequence of three-dimensional models which show that rapid inner core rotation stabilizes the hiccup instability. Instead, the inner core becomes quite flat and undergoes a cycle of binary fragmentation, binary decay into a single object surrounded by a bar, breakup of the bar into a binary, etc. When lesser amounts of rotation are involved, the inner core does hiccup somewhat, but mass is ejected in only a few directions, leading to several broad streams of ejecta.
Evidence suggesting the existence of old stars in young galactic clusters, has motivated an investigation of a model in which such stars are captured from the field during the gas cloud collapse phase of star formation. A hyperbolic field star cannot be gravitationally captured by a static cloud since the star's total energy is conserved. However, if the cloud's potential is an explicit function of time, the energy of a penetrating star can be decreased. As a first step in modeling the dynamical capture of field stars, capture by freely falling, homogeneous spherical cloud is considered. This simplified model is solved, and a prediction is obtained for number of captured field stars in terms of the cloud parameters and dispersion velocities of the field stars.
Explore the source record for details and available documents.
The 1(10)-1(11) transition of ortho-H2D(+) at 372 GHz has been sought in several dark clouds. The transition was not detected; the best upper limits obtained are about 0.3 K (3 sigma). We derive upper limits for the ortho-H2D(+) column density and briefly discuss their meaning in comparison with a simple chemical model we have developed (Pagani et al., 1992).
High-resolution observations of gamma(sup 2) Vel with the Goddard High-Resolution Spectrograph (GHRS) echelle on the Hubble Space Telescope reveal the presence of seven narrow absorption components, with LSR velocities between -23 and +9 km s(exp -1). Three of these show column density ratios N(S(++))/N(S(+)) and N(P(++))/N(P+)) of about 1 or more, and can be identified as H II regions, while the other four are H I regions, consistent with the O I profile and with the overall H(sup 0) column density of 5.9 x 10(exp 19) cm(exp -2), given the usual assumptions that S is undepleted while O has a depletion D(O) = -0.3 dex. The depletions of Fe, Si, and Mn, which could be measure accurately for two of the four H I regions (components 6 and 7), differ somewhat from the values of D(sub ws) found for slowly moving warm clouds in HD 93521; in particular, for the component at 4.0 km s(exp -1) (No. 6), abosolute of D exceeds absolute of D(sub ws) by 0.1-0.4 dex, while for that at 9.3 km s(exp -1) (No. 7), absolute of D equals absolute of D(sub ws) on the average. The observed ratio of Fe + Mg atoms to Si atoms in the grains of component 6 is 2.04 +/-0.10, consistent with an olivine grain composition; the Fe/Mg ratio is 1.5 +/- 0.2. The electron density in component 6, determined from the C II(sup *) feature, is 0.075 +/- 0.013 cm (exp -3), about two-thirds of that found for clouds of this velocity in HD 93521. In the two conspicuous H II regions, components 3 and 4, n(sub e), determined from the Si II(sup *) feature, is about 1 cm(exp -3). From the column density of S(+) + S(++) in these two components, the total H II path length is about 40 pc. With the radius of a wind-blown bubble around gamma(sup 2) Vel set equal to 60 pc, the effective Stromgren radius is about 100 pc, requiring that T approx. equal to 50,000 K for the Wolf-Rayet component of the gamma(sup 2) Vel binary. Since zeta Pup is a comparable source of ionizing radiation, this temperature is an upper limit. The profiles of the strongest H2 absorption features, from Copernicus archives, indicate that the absorbing molecules have a mean velocity identical with that of the strongest H II component (No. 4). We have no explanation for the possible presence of these H2 molecules in a region of ionized H. Alternatively, the H2 profiles can be explained by molecules in the two adjacent (in velocity) H I regions, components 2 and 5, provided their H I gas has densities and temperatures typical of normal cold clouds. The GHRS data show absorption by highly ionized atoms Si(3+) and C(3+), N(4+) in broad features, in addition to the narrow-line absorption by Si(3+) and C(3+) observed in the dominant H II components, Nos. 3 and 4. The broad C(3+) and N(4+) features have widths corresponding to T in the range (4-8) x 10(exp 5) K, consistent with the broad O(5+) line shown in Copernicus data. Despite some observational uncertainties, the ratios of column densities in the broad C(3+), N(4+), and O(5+) features agree to +/- 0.1 dex with theoretical values for warm gas, heating and evaporating by thermal conduction from an adjacent hot region. Outward evaporation from an isolated cloud in a hot ambient gas cannot be distinguished, on the basis of these data, from inward evaporation of a warm shell, compressed by an expanding, hot stellar-wind bubble. For several halo stars, the C IV/O VI ratio has a quite different average value, perhaps consistent with cooling of infalling hot gas instead of conductive heating and evaporation.
The spectrum of HD 149881, a halo star 1300 pc above the Galactic plane, has been analyzed, based on high-resolution GHRS data from Hubble Space Telescope (HST) together with 21 cm emission observations and profiles of Na I and Ca II features. The results reveal 11 absorption components, of which nine are produced in H I regions, with log N(H(sup 0) in each region mostly approximately 19.5; one component clearly originates in an H II region. Detailed column densities show that (D(Zn)), the absolute value of the depletion D, for Zn, like that for S, does not much exceed 0.1 dex. With Zn used as a standard of reference, D(Fe) varies from -1.1 to -0.6 dex; (D(Cr)) and (D(Mn)) average less than (D(Fe)) by about 0.1 and 0.2 dex, respectively, with dispersions of about 0.1 dex in this difference. As in Papers I and II, the measured D(Si) is consistent with a 2 to 1 ratio of Fe to Si atoms in grains. The kinetic temperature is these H I regions varies widely among components. Maximum temperatures, found from the widths of 21 cm emission components, range from about 100 to 10,000 K, while actual temperatures, found by comparing the H I widths with the b-values for the ultraviolet features, are about equal or less than 1000 K for six of the H I components. In contrast, toward the halo star HD 93521, temperatures of about 6000 K were found for most of the components. With use of these HD 149881 temperatures, values of n(sub e) were obtained from the ionization equilibrium of Ca(+) and also of Na(0). These values range over an order of magnitude, from about 4.9 x 10(exp -2) to 1.3 x 10(exp -3)/cu cm, systematically smaller (but also less certain) than in HD 93521. For the low values of n(sub e) found in the coldest components, the electrons can be accounted for by ionization of C atoms, if n(sub H) approximately equals 10/cu cm. For the warmest component, an ionization probability of beta(H(sup 0)) approximately equals 3 x 10(exp -15)/s is required, with lower values for other H I components.
We have used the ROSAT high-resolution imager and optical emission line data to study the structure of a bright optical knot on the southeastern rim of the Cygnus Loop supernova remnant. This knot has been identified as an encounter between the blast wave and a small isolated cloud. The knot appears in projection just behind the blast wave, which is traced by Balmer line filaments which bound the X-ray emission. The knot is a prominent X-ray feature, consisting of a numerb of filaments which are correlated with the optical line emission. These data permit a detailed view of the blast wave interaction. By combining the optical and X-ray data it is possible to trace the blast wave as a continuous surface from its southern edge around the western side of the cloud that then continues on to the north. The southeastern knot is an indentation on the surface of the blast wave. Thus the southeastern knot is not a small cloud that has been overrun by the blast wave, but the tip of a larger cloud. Bright X-ray emission is associated with bright radiative filaments. The location of this emission, upstream of the radiative shocks, implies that the enhanced X-rays come from a reverse shock. The presence of a reverse shock is further evidence that the southeastern knot represents an early stage of a blast wave encountering a large cloud. A consistent picture of the Cygnus Loop as an explosion within a preexisting cavity is emerging. These observations agree with that picture.
A heuristic criterion, based on linear perturbation analysis, is applied to the initial growth of density perturbations in isothermal or adiabatic gas clouds, with initially uniform density and uniform rotation. The heuristic criterion is shown to be consistent with the available results from numerical calculations of cloud collapse. The criterion predicts that perturbations varying as cos (m(phi)) will be most likely to grow when )pi is small, unless the cloud is nearly pressureless.
We modeled the long-term and small-scale evolution of molecular clouds using direct 2D and 3D magnetohydrodynamic (MHD) simulations. This work followed up on previous research by our group under auspices of the ATP in which we studied the energetics of turbulent, magnetized clouds and their internal structure on intermediate scales. Our new work focused on both global and smallscale aspects of the evolution of turbulent, magnetized clouds, and in particular studied the response of turbulent proto-cloud material to passage through the Galactic spiral potential, and the dynamical collapse of turbulent, magnetized (supercritical) clouds into fragments to initiate the formation of a stellar cluster. Technical advances under this program include developing an adaptive-mesh MHD code as a successor to ZEUS (ATHENA) in order to follow cloud fragmentation, developing a shearing-sheet MHD code which includes self-gravity and externally-imposed gravity to follow the evolution of clouds in the Galactic potential, and developing radiative transfer models to evaluate the internal ionization of clumpy clouds exposed to external photoionizing UV and CR radiation. Gammie's work at UIUC focused on the radiative transfer aspects of this program.
The recent confirmation by Ziurys and Apponi of the detection of HOC+ toward Sgr B2 (OH), and their identification of the ion in Orion-KL and several other sources show that HOC+ is far more abundant than predicted by previous ion-molecule models. In these models, the reaction HOC(+) + H2 --> HCO(+) + H2 is assumed to rapidly destroy HOC+, based on the results of a prior calculation. We have recalculated the rate of this reaction as a function of temperature using a new ab initio potential surface and a phase space approach to the dynamics which includes tunneling. The newly calculated rate is small (< or = 1 x 10(-10) cm3 s-1) at temperatures under 100 K.