Showing papers in "The Astrophysical Journal in 1994"

Toward a theory of interstellar turbulence. 2. Strong Alfvenic turbulence

Abstract: We continue to investigate the possibility that interstellar turbulence is caused by nonlinear interactions among shear Alfven waves. Here, as in Paper I, we restrict attention to the symmetric case where the oppositely directed waves carry equal energy fluxes. This precludes application to the solar wind in which the outward flux significantly exceeds the ingoing one. All our detailed calculations are carried out for an incompressible magnetized fluid. In incompressible magnetohydrodynamics (MHD), nonlinear interactions only occur between oppositely direct waves. Paper I contains a detailed derivation of the inertial range spectrum for the weak turbulence of shear Alfven waves. As energy cascades to high perpendicular wavenumbers, interactions become so strong that the assumption of weakness is no longer valid. Here, we present a theory for the strong turbulence of shear Alfven waves. It has the following main characteristics. (1) The inertial-range energy spectrum exhibits a critical balance beween linear wave periods and nonlinear turnover timescales. (2) The "eddies" are elongated in the direction of the field on small spatial scales; the parallel and perpendicular components of the wave vector, k_z and k_⊥, are related by k_z ≈ k^(2/3) _⊥L^(-1/3), where L is the outer scale of the turbulence. (3) The "one-dimensional" energy spectrum is proportional to k^(-5/3) _⊥-an anisotropic Kolmogorov energy spectrum. Shear Alfvenic turbulence mixes specific entropy as a passive contaminant. This gives rise to an electron density power spectrum whose form mimics the energy spectrum of the turbulence. Radio, wave scattering by these electron density fluctuations produces anisotropic scatter-broadened images. Damping by ion-neutral collisions restricts Alfvenic turbulence to highly ionized regions of the interstellar medium. We expect negligible generation of compressive MHD waves by shear Alfven waves belonging to the critically balanced cascade. Viscous and collisionless damping are also unimportant in the interstellar medium (ISM). Our calculations support the general picture of interstellar turbulence advanced by Higdon.

Advection-dominated Accretion: A Self-similar Solution

Abstract: We consider viscous rotating accretion flows in which most of the viscously dissipated energy is stored as entropy rather than being radiated. Such advection-dominated flows may occur when the optical depth is either very small or very large. We obtain a family of self-similar solutions where the temperature of the accreting gas is nearly virial and the flow is quasi-spherical. The gas rotates at much less than the Keplerian angular velocity; therefore, the central stars in such flows will cease to spin up long before they reach the break-up limit. Further, the Bernoulli parameter is positive, implying that advection-dominated flows are susceptible to producing outflows. Convection is likely in many of these flows and, if present, will tend to enhance the above effects. We suggest that advection-dominated accretion may provide an explanation for the slow spin rates of accreting stars and the widespread occurrence of outflows and jets in accreting systems.

Dust extinction of the stellar continua in starburst galaxies: The Ultraviolet and optical extinction law

Abstract: We analyze the International Ultraviolet Explorer (IUE) UV and the optical spectra of 39 starburst and blue compact galaxies in order to study the average properties of dust extinction in extended regions of galaxies. The optical spectra have been obtained using an aperture which matches that of IUE, so comparable regions within each galaxy are sampled. The data from the 39 galaxies are compared with five models for the geometrical distribution of dust, adopting as extinction laws both the Milky Way and the Large Magellanic Cloud laws. The commonly used uniform dust screen is included among the models. We find that none of the five models is in satisfactory agreement with the data. In order to understand the discrepancy between the data and the models, we have derived an extinction law directly from the data in the UV and optical wavelength range. The resulting curve is characterized by an overall slope which is more gray than the Milky Way extinction law's slope, and by the absence of the 2175 A dust feature. Remarkably, the difference in optical depth between the Balmer emission lines H(sub alpha) and H(sub beta) is about a factor of 2 larger than the difference in the optical depth between the continuum underlying the two Balmer lines. We interpret this discrepancy as a consequence of the fact that the hot ionizing stars are associated with dustier regions than the cold stellar population is. The absence of the 2175 A dust feature can be due either to the effects of the scattering and clumpiness of the dust or to a chemical composition different from that of the Milky Way dust grains. Disentangling the two interpretations is not easy because of the complexity of the spatial distribution of the emitting regions. The extinction law of the UV and optical spectral continua of extended regions can be applied to the spectra of medium- and high-redshift galaxies, where extended regions of a galaxy are, by necessity, sampled.

H II regions and the abundance properties of spiral galaxies

Abstract: We investigate the relationships between the characteristic oxygen abundance, the radial abundance gradient, and the macroscopic properties of spiral galaxies by examining the properties of individual H II regions within those galaxies Our observations of the line flux ratio (O II) lambda lambda 3726, 3729 + (O III) lambda lambda 4959, 5007)/H beta for 159 H II regions in 14 spiral galaxies are combined with published data to provide a sample of 39 disk galaxies for which (O II) + (O III)/H beta has been measured for at least five H II regions We find that the characteristic gas-phase abundances and luminosities of spiral galaxies are strongly correlated This relationship maps almost directly onto the luminosity-metallicity relationship of irregular galaxies and is also quite similar to that found for elliptical and dwarf spheroidal galaxies Within our sample of spirals, a strong correlation between characteristic abundance and Hubble type also exists The correlation between luminosity and Hubble type complicates the issue, but we discuss several interpretations of the correlations The relationship between circular velocity and characteristic abundance is also discussed We find that the slopes of the radial abundance gradients, when expressed in units of dex/isophotal radius, do not significantly correlate with either luminosity or Hubble type However, the hypothesis that both early and very late type spirals have shallower gradients than intermediate spirals is consistent with the data We find suggestive evidence that the presence of a bar induces a flatter gradient and also briefly discuss whether abundance gradients are exponential, as is usually assumed We investigate the properties of individual H II regions in a subset of 42 regions for which we have spectra that cover almost the entire spectral range from 3500 to 9800 A We use those data to estimate the densitites and ionizing spectra within the H II regions We confirm that the ionizing spectrum hardens with increasing radius and decreasing abundance We find no correlation between the ionization parameter and either radius or abundance, but this may be due to significant scatter introduced by the simple conversion of line ratios to ionization parameter

Dynamics of binary-disk interaction. 1: Resonances and disk gap sizes

Abstract: We investigate the gravitational interaction of a generally eccentric binary star system with circumbinary and circumstellar gaseous disks. The disks are assumed to be coplanar with the binary, geometrically thin, and primarily governed by gas pressure and (turbulent) viscosity but not self-gravity. Both ordinary and eccentric Lindblad resonances are primarily responsible for truncating the disks in binaries with arbitrary eccentricity and nonextreme mass ratio. Starting from a smooth disk configuration, after the gravitational field of the binary truncates the disk on the dynamical timescale, a quasi-equilibrium is achieved, in which the resonant and viscous torques balance each other and any changes in the structure of the disk (e.g., due to global viscous evolution) occur slowly, preserving the average size of the gap. We analytically compute the approximate sizes of disks (or disk gaps) as a function of binary mass ratio and eccentricity in this quasi-equilibrium. Comparing the gap sizes with results of direct simulations using the smoothed particle hydrodynamics (SPH), we obtain a good agreement. As a by-product of the computations, we verify that standard SPH codes can adequately represent the dynamics of disks with moderate viscosity, Reynolds number R approximately 10(exp 3). For typical viscous disk parameters, and with a denoting the binary semimajor axis, the inner edge location of a circumbinary disk varies from 1.8a to 2.6a with binary eccentricity increasing from 0 to 0.25. For eccentricities 0 less than e less than 0.75, the minimum separation between a component star and the circumbinary disk inner edge is greater than a. Our calculations are relevant, among others, to protobinary stars and the recently discovered T Tau pre-main-sequence binaries. We briefly examine the case of a pre-main-sequence spectroscopic binary GW Ori and conclude that circumbinary disk truncation to the size required by one proposed spectroscopic model cannot be due to Linblad resonances, even if the disk is nonviscous.

Magnetocentrifugally driven flows from young stars and disks. 1: A generalized model

Abstract: We propose a generalized model for stellar spin-down, disk accretion, and truncation, and the origin of winds, jets, and bipolar outflows from young stellar objects. We consider the steady state dynamics of accretion of matter from a viscous and imperfectly conducting disk onto a young star with a strong magnetic field. For an aligned stellar magnetosphere, shielding currents in the surface layers of the disk prevent stellar field lines from penetrating the disk everywhere except for a range of radii about pi = R(sub x), where the Keplerian angular speed of rotation Omega(sub x) equals the angular speed of the star Omega(sub *). For the low disk accretion rates and high magnetic fields associated with typical T Tauri stars, R(sub x) exceeds the radius of the star R(sub *) by a factor of a few, and the inner disk is effectively truncated at a radius R(sub t) somewhat smaller than R(sub x). Where the closed field lines between R(sub t) and R(sub x) bow sufficiently inward, the accreting gas attaches itself to the field and is funneled dynamically down the effective potential (gravitational plus centrifugal) onto the star. Contrary to common belief, the accompanying magnetic torques associated with this accreting gas may transfer angular momentum mostly to the disk rather than to the star. Thus, the star can spin slowly as long as R(sub x) remains significantly greater than R(sub *). Exterior to R(sub x) field lines threading the disk bow outward, which makes the gas off the mid-plane rotate at super-Keplerian velocities. This combination drives a magnetocentrifugal wind with a mass-loss rate M(sub w) equal to a definite fraction f of the disk accretion rate M(sub D). For high disk accretion rates, R(sub x) is forced down to the stellar surface, the star is spun to breakup, and the wind is generated in a manner identical to that proposed by Shu, Lizano, Ruden, & Najita in a previous communication to this journal. In two companion papers (II and III), we develop a detailed but idealized theory of the magnetocentrifugal acceleration process.

Unsteady outflow models for cosmological gamma-ray bursts

Abstract: The 'event' that triggers a gamma-ray burst cannot last for more than a few seconds. This is, however, long compared with the dynamical timescale of a compact stellar-mass object (approximately 10 (exp-3) s). Energy is assumed to be released as an outflow with high mean Lorentz factor Gamma. But a compact stellar-mass collapse or merger is, realistically, likely to generate a mass (or energy) flux that is unsteady on some timescales in the range 10(exp -3) - 10 s. If Gamma fluctuates by a factor of approximately 2 around its mean value, relative motions within the outflowing material will themselves (in the comoving frame) be relativistic, and can give rise to internal shocks. For Gamma approximately 10(exp 2), the resultant dissipation occurs outside the 'photosphere' and can convert a substantial fraction of the overall outflow energy into nonthermal radiation. This suggests a mechanism for cosmological bursts that demands less extreme assumptions (in respect of Gamma-values, freedom from baryonic contamination, etc.) than earlier proposals.

Comptonization of Diffuse Ambient Radiation by a Relativistic Jet: The Source of Gamma Rays from Blazars?

Abstract: Recent Energy Gamma Ray Experiment Telescope (EGRET) observations of blazars have revealed strong, variable gamma-ray fluxes with no signatures of gamma-ray absorption by pair production. This radiation probably originates from the inner parts of relativistic jets which are aimed nearly toward us. On sub-parsec scales, the jet will be pervaded by radiation from the broad-line region, as well as by photons from the central continuum source (some of which will be scattered by thermal plasma). In a frame moving with the relativistic outflow, the energy of this ambient radiation would be enhanced. This radiation would be Comptonized by both cold and relativistic electrons in the jet, yielding (in the observer's frame) a collimated beam of X-rays and gamma rays. On the assumption that this process dominates self-Comptonization of synchrotron radiation, we develop a self-consistent model for variable gamma-ray emission, involving a single population of relativistic electrons accelerated by a disturbance in the jet. The spectral break between the X-ray and gamma-ray band, observed in 3C 279 and deduced for other blazars, results from inefficient radiative cooling of lower energy electrons. The existence of such a break strongly favors a model involving Comptonization of an external radiation field over a synchrotron self-Compton model. We derive constraints on such model parameters as the location and speed of the source, its dimensions and internal physical parameters, the maximum photon energies produced in the source, and the density and distribution of ambient radiation. Finally, we discuss how observations might discriminate between our model and alternative ones invoking Comptonization of ambient radiation.

Composition and radiative properties of grains in molecular clouds and accretion disks

Abstract: We define a model of the compositon and abundances of grains and gases in molecular cloud cores and accretion disks around young stars by employing a wide range of astronomical data and theory, the composition of primitive bodies in the solar system, and solar elemental abundances. In the coldest portions of these objects, we propose that the major grain species include olivine (Fe, Mg, 2SiO4), orthopyroxene (Fe, Mg, SiO3), volatile and refractory organics, water ice, troilite (FeS), and metallic iron. This compositional model differs from almost all previous models of the interstellar medium (ISM) by having organics as the major condensed C species, rather than graphite; by including troilite as a major grain species; and by specifying the mineralogical composition of the condensed silicates. Using a combination of laboratory measurements of optical constants and asymptotic theory, we derive values of the real and imaginary indices of refraction of these grain species over a wavelength range that runs from the vacuum ultraviolet (UV) to the radio domain. The above information on grain properties is used to estimate the Rosseland mean opacity of the grains and their monochromatic opacity.

Low-Temperature Rosseland Opacities

Abstract: A new, comprehensive set of low-temperature opacity data has been assembled. From this basic data set, Rosseland and Planck mean opacities have been computed for temperatures between 12,500 and 700 K. In addition to the usual continuous absorbers, atomic line absorption (with more than 8 million lines), molecular line absorption (with nearly 60 million lines), and grain absorption and scattering (by silicates, iron, carbon, and SiC) have been accounted for. The absorption due to lines is computed monochromatically and included in the mean with the opacity sampling technique. Grains are assumed to form in chemical equilibrium with the gas and to form into a continuous distribution of ellipsoids. Agreement of these opacities with other recent tabulations of opacities for temperatures above 5000 K is excellent. It is shown that opacities which neglect molecules become unreliable for temperatures below 5000 K. Triatomic molecules become important absorbers at 3200 K. Similarly, grains must be included in the computation for temperatures below 1700 K.

Determining structure in molecular clouds

Abstract: We descibe an automatic, objective routine for analyzing the clumpy structure in a spectral line position-position-velocity data cube. The algorithm works by first contouring the data at a multiple of the rms noise of the observations, then searches for peaks of emission which locate the clumps, and then follows them down to lower intensities. No a proiri clump profile is assumed. By creating simulated data, we test the performance of the algorithm and show that a contour map most accurately depicts internal structure at a contouring interval equal to twice the rms noise of the map. Blending of clump emission leads to small errors in mass and size determinations and in severe cases can result in a number of clumps being misidentified as a single unit, flattening the measured clump mass spectrum. The algorithm is applied to two real data sets as an example of its use. The Rosette molecular cloud is a 'typical' star-forming cloud, but in the Maddalena molecular cloud high-mass star formation is completely absent. Comparison of the two clump lists generated by the algorithm show that on a one-to-one basis the clumps in the star-forming cloud have higher peak temperatures, higher average densities, and are more gravitationally bound than in the non-star-forming cloud. Collective properties of the clumps, such as temperature-size-line-width-mass relations appear very similar, however. Contrary to the initial results reported in a previous paper (Williams & Blitz 1993), we find that the current, more thoroughly tes ted analysis finds no significant difference in the clump mass spectrum of the two clouds.

Generalized comptonization models and application to the recent high-energy observations

Abstract: The theory of spectral formation in thermal X-ray sources, where the effects of Comptonization and Klein-Nishina corrections are important, is presented. Analytical expressions are obtained for the produced spectrum as a function of such input parameters as the plasma temperature, the optical depth of the plasma cloud and the injected soft photon spectrum. The analytical theory developed here takes into account the dependence of the scattering opacity on the photon energy. It is shown that the plasma temperature as well as the asymptotic rate of photon escape from the plasma cloud determine the shape of the upscattered hard tail in the emergent spectra, even in the case of very small optical depths. The escape distributions of photons are given for any optical depth of the plasma cloud and their asymptotic dependence for very small and large optical depths are examined. It is shown that this new generalized approach can fit spectra for a large variety of hard X-ray sources and determine the plasma temperature in the region of main energy release in Cyg X-1 and the Seyfert galaxy NGC 4151.

Using fu orionis outbursts to constrain self-regulated protostellar disk models

Abstract: 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.

R(sub nu)-dependent optical and near-ultraviolet extinction

Abstract: We have derived extinctions A(lambda)/A(V) at the wavelengths of the uvby filters for 22 stars, with a range of values of R(sub nu), from the sample of Cardelli, Clayton, & Mathis (1989, hereafter CCM). We have fit these extinctions, and also UBVRIJHKL, IUE and ANS extinction measurements, with linear relations A(lambda)/A/(V) = a+b/R(sub nu) and fit a and b as a function of x(=1/lambda) with polynomials to obtain an R(sub nu)-dependent mean extinction law (A(x)/A(V) = a(x) + b(x)/R(sub nu))in the optical and near-ultraviolet (1.1/micrometer less than or equal to 3.3/micrometer). This law is virtually identical to the CCM extinction law for large values of R(sub nu)(R(sub nu) approximately 5) but is slightly lower in the near-ultraviolet for smaller R(sub nu) (R(sub nu) approximately 3). The extinction law presented here agrees much better with a high-resolution extinction curve for the diffuse interstellar medium (R(sub nu) approximately 3.1), presented by Bastiaansen (1992), than CCM. The deviations of individual extinction curves from the mean are dominated by observational errors. The wavelength resolution of this work is not high enough to show evidence for or against the existence of very broad structure in optical extinction curves.

Power spectrum analysis of three-dimensional redshift surveys

Abstract: We develop a general method for power spectrum analysis of three dimensional redshift surveys. We present rigorous analytical estimates for the statistical uncertainty in the power and we are able to derive a rigorous optimal weighting scheme under the reasonable (and largely empirically verified) assumption that the long wavelength Fourier components are Gaussian distributed. We apply the formalism to the updated 1-in-6 QDOT IRAS redshift survey, and compare our results to data from other probes: APM angular correlations; the CfA and the Berkeley 1.2Jy IRAS redshift surveys. Our results bear out and further quantify the impression from e.g.\ counts-in-cells analysis that there is extra power on large scales as compared to the standard CDM model with $\Omega h\simeq 0.5$. We apply likelihood analysis using the CDM spectrum with $\Omega h$ as a free parameter as a phenomenological family of models; we find the best fitting parameters in redshift space and transform the results to real space. Finally, we calculate the distribution of the estimated long wavelength power. This agrees remarkably well with the exponential distribution expected for Gaussian fluctuations, even out to powers of ten times the mean. Our results thus reveal no trace of periodicity or other non-Gaussian behavior.

Measurement of the Cosmic Microwave Background spectrum by the COBE FIRAS instrument

Abstract: The cosmic microwave background radiation (CMBR) has a blackbody spectrum within 3.4 x 10(exp -8) ergs/sq cm/s/sr cm over the frequency range from 2 to 20/cm (5-0.5 mm). These measurements, derived from the Far-Infrared Absolute Spectrophotomer (FIRAS) instrument on the Cosmic Background Explorer (COBE) satellite, imply stringent limits on energy release in the early universe after t approximately 1 year and redshift z approximately 3 x 10(exp 6). The deviations are less than 0.30% of the peak brightness, with an rms value of 0.01%, and the dimensionless cosmological distortion parameters are limited to the absolute value of y is less than 2.5 x 10(exp -5) and the absolute value of mu is less than 3.3 x 10(exp -4) (95% confidence level). The temperature of the CMBR is 2.726 +/- 0.010 K (95% confidence level systematic).

Observational and theoretical constraints on singular dark matter halos

Abstract: The distribution of dark matter around galactic or cluster halos has usually been assumed to be approximately isothermal with a non-zero core radius, which is expected to be of the order of the size of the visible matter distribution. Recently, the possibility has been raised that dark matter halos might be singular in the sense that the dark matter density $\rho$ could increase monotonically with radius $r$ down to a very small distance from the center of galaxies or clusters. Such central cusps in the dark matter density could lead to a high flux of gamma rays from WIMP dark matter annihilation. Here we analyze two possibilities that have been discussed in the literature, $\rho \propto r^{-n}$ with $n \approx 1\ {\rm or}\ 2$, and point out that such density profiles are excluded by gravitional lensing analyses on cluster scales and by the rotation curves of gas-rich, halo-dominated dwarf spirals on small scales. We also point out that if spiral galaxies form by gas infall inside dark matter halos, as they are expected to do in any hierarchical clustering model, such profiles almost always lead to falling rotation curves after infall, contrary to observations.

The photoelectric heating mechanism for very small graphitic grains and polycyclic aromatic hydrocarbons

Abstract: We have theoretically modeled the gas heating associated with the photoelectric ejection of electrons from a size distribution of interstellar carbon grains which extends into the molecular domain. We have considered a wide range of physical conditions for the interstellar gas (1 less than G(sub 0) less than 10(exp 5), with G(sub 0) being the intensity of the incident far-UV field in units of the Habing interstellar radiation field; 2.5 x 10( exp -3) less than n(sub e) less than 75/cu cm, with n(sub e) being the electron density; 10 less than T less than 10,000 K, with T being the gas temperature). The results show that about half of the heating is due to grains less than 1500 C atoms (less than 15 A). The other half originates in somewhat larger grains (1500-4.5 x 10(exp 5) C atoms; 15 less than 100 A). While grains larger than this do absorb about half of the available far-UV photons, they do not contribute appreciably to the gas heating. This strong dependence of gas heating on size results from the decrease in yield and from the increased grain charge (hence larger Coulomb losses) with increasing grain size. We have determined the net photoelectric heating rate and evaluated a simple analytical expression for the heating efficiency, dependent only on G(sub 0), T, and n(sub e). This expression is accurate to 3% over the whole parameter range and is valid up to gas temperatures of 10(exp 4) K, at which point the dominant gas-dust heat exchange mechanism becomes the recombination of electrons with grains rather than photoelectric ejection. The calculated heating efficiency for neutral grains is in good agreement with that derived from observations of the diffuse interstellar clouds. Our results also agree well with the Far Infrared Absolute Spectrometer (FIRAS) observations on the Cosmic Background Explorer Satellite. Finally, our photoelectric heating efficiency is compared to previous studies.

Old stellar populations. 5: Absorption feature indices for the complete LICK/IDS sample of stars

Abstract: Twenty-one optical absorption features, 11 of which have been previously defined, are automatically measured in a sample of 460 stars Following Gorgas et al, the indices are summarized in fitting functions that give index strengths as functions of stellar temperature, gravity, and (Fe/H) This project was carried out with the purpose of predicting index strengths in the integrated light of stellar populations of different ages and metallicities, but the data should be valuable for stellar studies in the Galaxy as well Several of the new indices appear to be promising indicators of metallicity for old stellar populations A complete list of index data and atmospheric parameters is available in computer-readable form

On the hydrodynamic interaction of shock waves with interstellar clouds. 1: Nonradiative shocks in small clouds

Abstract: The interstellar medium (ISM) is inhomogeneous, with clouds of various temperatures and densities embedded in a tenuous intercloud medium. Shocks propagating through the ISM can ablate or destroy the clouds, at the same time significantly altering the properties of the intercloud medium. This paper presents a comprehensive numerical study of the simplest case of the interaction between a shock wave and a spherical cloud, in which the shock far from the cloud is steady and planar, and in which radiative losses, thermal conduction, magnetic fields, and gravitational forces are all neglected. As a result, the problem is completely specified by two numbers: the Mach number of the shock, M, and the ratio of the density of the cloud to that of the intercloud medium, Chi. For strong shocks we show that the dependence on M scales out, so the primary independent parameter is Chi. Variations from this simple case are also considered: the potential effect of radiative losses is assessed by calculations in which the ratio of specific heats in the cloud is 1.1 instead of 5/3; the effect of the initial shape of the cloud is studied by using a cylindrical cloud instead of a spherical one; and the role of the initial shock is determined by considering the case of a cloud embedded in a wind. Local adaptive mesh refinement techniques with a second-order, two-fluid, two-dimensional Godunov hydrodynamic scheme are used to address these problems, allowing heretofore unobtainable numerical resolution. Convergence studies to be described in a subsequent paper demonstrate that about 100 zones per cloud radius are needed for accurate results; previous calculations have generally used about a third of this number. The results of the calculations are analyzed in terms of global quantities which provide an overall description of te shocked cloud: the size and shape of the cloud, the mean density, the mean pressure, the mean velocity, the velocity dispersion, and the total circulation.

The r-process and neutrino-heated supernova ejecta

Abstract: As a neutron star is formed by the collapse of the iron core of a massive star, its Kelvin-Helmholtz evolution is characterized by the release of gravitational binding energy as neutrinos. The interaction of these neutrinos with heated material above the neutron star generates a hot bubble in an atmosphere that is nearly in hydrostatic equilibrium and heated, after approximately 10 s, to an entropy of S/N(sub AS)k greater than or approximately = 400. The neutron-to-proton ratio for material moving outward through this bubble is set by the balance between neutrino and antineutrino capture on nucleons. Because the electron antineutrino spectrum at this time is hotter than the electron neutrino spectrum, the bubble is neutron-rich (0.38 less than or approximately = Y(sub e) less than or approximately = 0.47). Previous work using a schematic model has shown that these conditions are well suited to the production of heavy elements by the r-process. In this paper we have advanced the numerical modeling of a 20 solar mass 'delayed' supernova explosion to the point that we can follow the detailed evolution of material moving through the bubble at the late times appropiate to r-process nucleosynthesis. The supernova model predicts a final kinetic energy for the ejecta of 1.5 x 10(exp 51) ergs and leaves behind a remnant with a baryon mass of 1.50 solar mass (and a gravitational mass of 1.445 solar mass). We follow the thermodynamic and compositional evolution of 40 trajectories in rho(t), T(t), Y(sub e)(t) for a logarithmic grid of mass elements for the last approximately = 0.03 solar mass to be ejected by the proto-neutron star down to the last less than 10(exp -6) solar mass of material expelled at up to approximately = 18 s after core collapse. We find that an excellent fit to the solar r-process abundance distribution is obtained with no adjustable parameters in the nucleosynthesis calculations. Moreover, the abundances are produced in the quantities required to account for the present Galactic abundances. However, at earlier times, this one-dimensional model ejects too much material with entropies S/N(sub A)k approximately 50 and Y(sub e) approximately 0.46. This leads to an acceptable over production of N = 50 nuclei, particularly Sr-88, Y-89, and Zr-90, relative to their solar abundances. We speculate on various means to avoid the early overproduction and/or ejection of N = 50 isotonic nuclei while still producing and ejecting the correct amount of r-process material.

Proton and Electron Mean Free Paths: The Palmer Consensus Revisited

Abstract: We present experimental and theoretical evidence suggesting that the mean free path of cosmic-ray electrons and protons may be fundamentally different at low to intermediate (less than 50 MV) rigidities. The experimental evidence is from Helios observations of solar energetic particles, which show that the mean free path of 1.4 MV electrons is often similar to that of 187 MV protons, even though proton mean free paths continue to decrease comparatively rapidly with decreasing rigidty down to the lowest channels (about 100 MV) observed. The theoretical evidence is from computations of particle scattering in dynamical magnetic turbulence, which predict that electrons will have a larger mean free path than protons of the same rigidity. In the light of these new results, 'consensus' ideas about cosmic-ray mean free paths may require drastic revision.

Photoevaporation of Disks around Massive Stars and Application to Ultracompact H II Regions

Abstract: Young massive stars produce sufficient Lyman continuum photon luminosity Phi(sub i) to significantly affect the structure and evolution of the accretion disks surrounding them. A nearly static, ionized, isothermal 10(exp 4) K atmosphere forms above the neutral disk for disk radii r less than r(sub g) = 10(exp 15) M(sub 1) cm, where M(sub *) = 10 solar mass M(sub 1) is the stellar mass. For r approximately greater than r(sub g) the diffuse field created by hydrogen recombinations to the ground state in the photoionized gas above the disk produces a steady evaporation at the surface of the disk, and this H II gas flows freely out to the ISM (the 'disk wind'). The detailed structure depends on the mass-loss rate dot-M(sub w) of the fast, approximately greater than 1000 km/sec, stellar wind from the massive star. A critical mass-loss rate dot-M(sub cr) is defined such that the ram pressure of the stellar wind equals the thermal pressure of the H II atmosphere at r(sub g). In the weak stellar wind solution, dot-M(sub w) less than dot-M(sub cr), the diffuse photons from the atmosphere above r(sub g) produce a photoevaporative mass-loss rate from the disk at r approximately greater than r(sub g) of order 1 x 10(exp -5)(Phi(sub 49))(exp 1/2)(M(sub 1))(exp 1/2) solar mass/year, where Phi(sub i) = 10(exp 49) Phi(sub 49)/sec. The resulting slow (10 to 50 km/sec) ionized outflow, which persists for approximately greater than 10(exp 5) year for disk masses M(sub d) approximately 0.3 M(sub *), may explain the observational characteri stics of unresolved, ultracompact H II regions. In the strong stellar wind solution, dot-M(sub w) greater than dot-M(sub cr), the ram pressure of the stellar wind blows down the atmosphere for r less than r(sub g) and allows the stellar photons to penetrate to greater radii and smaller heights. A slow, ionized outflow produced mainly by diffuse photons is again created for r less than r(sub g); however, it is now dominated by the flow at r(sub w)(greater than r(sub g)), the radius at which the stellar wind ram pressure equals the thermal pressure in the evaporating flow. The mass-loss rate from the disk is of order 6 x 10(exp -5)dot-M(sub w-6) v(sub w8)(Phi (sub 49))(exp -1/2) solar mass/year, where dot-M(sub w-6) = M(sub w)/10(exp -6) solar mass/year and v(sub w8) = v(sub w)/1000 km/sec is the stellar wind velocity. The resulting outflow, which also persists for approximately greater than 10(exp 5) year may explain many of the more extended (r approximately greater than 10(exp 16) cm) ultracompact H II regions. Both the weak-wind and the strong-wind models depend entirely on stellar parameters Phi(sub i), M(sub *), dot-M(sub w)) and are independent of disk parameters as long as an extended r much greater than (r(sub g)), neutral disk exists. We compare both weak-wind and strong-wind model results to the observed radio free-free spectra and luminosities of ultracompact H II regions and to the interesting source MWC 349.

The Hubble Space Telescope Extragalactic Distance Scale Key Project. 1: The discovery of Cepheids and a new distance to M81

Abstract: We report on the discovery of 30 new Cepheids in the nearby galaxy M81 based on observations using the Hubble Space Telescope (HST). The periods of these Cepheids lie in the range of 10-55 days, based on 18 independent epochs using the HST wide-band F555W filter. The HST F555W and F785LP data have been transformed to the Cousins standard V and I magnitude system using a ground-based calibration. Apparent period-luminosity relations at V and I were constructed, from which apparent distance moduli were measured with respect to assumed values of mu(sub 0) = 18.50 mag and E(B - V) = 0.10 mag for the Large Magellanic Cloud. The difference in the apparent V and I moduli yields a measure of the difference in the total mean extinction between the M81 and the LMC Cepheid samples. A low total mean extinction to the M81 sample of E(B - V) = 0.03 +/- 0.05 mag is obtained. The true distance modulus to M81 is determined to be 27.80 +/- 0.20 mag, corresponding to a distance of 3.63 +/- 0.34 Mpc. These data illustrate that with an optimal (power-law) sampling strategy, the HST provides a powerful tool for the discovery of extragalactic Cepheids and their application to the distance scale. M81 is the first calibrating galaxy in the target sample of the HST Key Project on the Extragalactic Distance Scale, the ultimate aim of which is to provide a value of the Hubble constant to 10% accuracy.

The soft x-ray properties of a complete sample of optically selected quasars. 1: First results

Abstract: We present the results of ROSAT position sensitive proportional counter (PSPC) observations of 10 quasars. These objects are part of our ROSAT program to observe a complete sample of optically selected quasars. This sample includes all 23 quasars from the bright quasar survey with a redshift z less than or = 0.400 and a Galactic H I column density N(sup Gal sub H I) less than 1.9 x 10(exp 20)/sq cm. These selection criteria, combined with the high sensitivity and improved energy resolution of the PSPC, allow us to determine the soft (approximately 0.2-2 keV) X-ray spectra of quasars with about an order of magnitude higher precision compared with earlier soft X-ray observations. The following main results are obtained: Strong correlations are suggested between the soft X-ray spectral slope alpha(sub x) and the following emission line parameters: H beta Full Width at Half Maximum (FWHM), L(sub O III), and the Fe II/H beta flux ratio. These correlations imply the following: (1) The quasar's environment is likely to be optically thin down to approximately 0.2 keV. (2) In most objects alpha(sub x) varies by less than approximately 10% on timescales shorter than a few years. (3) alpha(sub x) might be a useful absolute luminosity indicator in quasars. (4) The Galactic He I and H I column densities are well correlated. Most spectra are well characterized by a simple power law, with no evidence for either significant absorption excess or emission excess at low energies, to within approximately 30%. We find mean value of alpha(sub x) = -1.50 +/- 0.40, which is consistent with other ROSAT observations of quasars. However, this average is significantly steeper than suggested by earlier soft X-ray observations of the Einstein IPC. The 0.3 keV flux in our sample can be predicted to better than a factor of 2 once the 1.69 micrometer(s) flux is given. This implies that the X-ray variability power spectra of quasars flattens out between f approximately 10(exp -5) and f approximately 10(exp -8) Hz. A steep alpha(sub x) is mostly associated with a weak hard X-ray component, relative to the near-IR and optical emission, rather than a strong soft excess, and the scatter in the normalized 0.3 keV flux is significantly smaller than the scatter in the normalized 2 keV flux. This argues against either thin or thick accretion disks as the origin of the soft X-ray emission. Further possible implications of the results found here are briefly discussed.

Inside the Supernova: A Powerful Convective Engine

Abstract: We present an extensive study of the inception of supernova explosions by following the evolution of the cores of two massive stars (15 and 25 Solar mass) in multidimension. Our calculations begin at the onset of core collapse and stop several hundred milliseconds after the bounce, at which time successful explosions of the appropriate magnitude have been obtained. Similar to the classical delayed explosion mechanism of Wilson, the explosion is powered by the heating of the envelope due to neutrinos emitted by the protoneutron star as it radiates the gravitational energy liberated by the collapse. However, as was shown by Herant, Benz, & Colgate, this heating generates strong convection outside the neutrinosphere, which we demonstrate to be critical to the explosion. By breaking a purely stratified hydrostatic equilibrium, convection moves the nascent supernova away from a delicate radiative equilibrium between neutrino emission and absorption, Thus, unlike what has been observed in one-dimensional calculations, explosions are rendered quite insensitive to the details of the physical input parameters such as neutrino cross sections or nuclear equation of state parameters. As a confirmation, our comparative one-dimensional calculations with identical microphysics, but in which convection cannot occur, lead to dramatic failures. Guided by our numerical results, we have developed a paradigm for the supernova explosion mechanism. We view a supernova as an open cycle thermodynamic engine in which a reservoir of low-entropy matter (the envelope) is thermally coupled and physically connected to a hot bath (the protoneutron star) by a neutrino flux, and by hydrodynamic instabilities. This paradigm does not invoke new or modified physics over previous treatments, but relies on compellingly straightforward thermodynamic arguments. It provides a robust and self-regulated explosion mechanism to power supernovae that is effective under a wide range of physical parameters.