On Radiation Pressure in Static, Dusty HII Regions

Radiation pressure acting on gas and dust causes HII regions to have central densities that are lower than the density near the ionized boundary. HII regions in static equilibrium comprise a family of similarity solutions, parametrized by 3 parameters: beta, gamma, and the product (Q_0 n_rms); beta characterizes the stellar spectrum, gamma characterizes the dust/gas ratio, Q_0 is the ionizing output from the star (photons/s), and n_rms is the rms density within the ionized region. Adopting standard values for beta and gamma, varying (Q_0 n_rms) generates a one-parameter family of density profiles, ranging from nearly uniform density (small Q_0 n_rms), to hollow-sphere HII regions (large Q_0 n_rms). When (Q_0 n_rms) exceeds 10^{52} cm^{-3} s^{-1}, dusty HII regions have conspicuous central cavities, even if no stellar wind is present. For given beta, gamma and (Q_0 n_rms), a fourth quantity, which can be Q_0, determines the overall size and density of the HII region. Examples of density and emissivity profiles are given. We show how quantities of interest -- such as the peak-to-center emissivity ratio, the rms-to-mean density ratio, the edge-to-rms density ratio, and the fraction of the ionizing photons absorbed by the gas -- depend on the 3 parameters beta, gamma, and (Q_0 n_rms). For dusty HII regions, compression of the gas and dust into an ionized shell results in a substantial increase in the fraction of the >13.6 eV photons that actually ionize H (relative to a uniform density HII region with the same dust/gas ratio and density n=n_rms). We discuss the extent to which radial drift of dust grains in HII regions can alter the dust-to-gas ratio. The applicability of these solutions to real HII regions is discussed.Comment: New material and figures that were not in version 1. To appear in Ap

A Keck High Resolution Spectroscopic Study of the Orion Nebula Proplyds

We present the results of spectroscopy of four bright proplyds in the Orion Nebula obtained at a velocity resolution of 6 km/s. After careful isolation of the proplyd spectra from the confusing nebular radiation, the emission line profiles are compared with those predicted by realistic dynamic/photoionization models of the objects. The spectral line widths show a clear correlation with ionization potential, which is consistent with the free expansion of a transonic, ionization-stratified, photoevaporating flow. Fitting models of such a flow simultaneously to our spectra and HST emission line imaging provides direct measurements of the proplyd size, ionized density and outflow velocity. These measurements confirm that the ionization front in the proplyds is approximately D-critical and provide the most accurate and robust estimate to date of the proplyd mass loss rate. Values of 0.7E-6 to 1.5E-6 Msun/year are found for our spectroscopic sample, although extrapolating our results to a larger sample of proplyds implies that 0.4E-6 Msun/year is more typical of the proplyds as a whole. In view of the reported limits on the masses of the circumstellar disks within the proplyds, the length of time that they can have been exposed to ionizing radiation should not greatly exceed 10,000 years - a factor of 30 less than the mean age of the proplyd stars. We review the various mechanisms that have been proposed to explain this situation, and conclude that none can plausibly work unless the disk masses are revised upwards by a substantial amount.Comment: 23 pages, 8 figures, uses emulateapj.sty, accepted for publication in The Astronomical Journal (scheduled November 1999

Estimating Electric Fields from Vector Magnetogram Sequences

Determining the electric field (E-field) distribution on the Sun's photosphere is essential for quantitative studies of how energy flows from the Sun's photosphere, through the corona, and into the heliosphere. This E-field also provides valuable input for data-driven models of the solar atmosphere and the Sun-Earth system. We show how Faraday's Law can be used with observed vector magnetogram time series to estimate the photospheric E-field, an ill-posed inversion problem. Our method uses a "poloidal-toroidal decomposition" (PTD) of the time derivative of the vector magnetic field. The PTD solutions are not unique; the gradient of a scalar potential can be added to the PTD E-field without affecting consistency with Faraday's Law. We present an iterative technique to determine a potential function consistent with ideal MHD evolution; but this E-field is also not a unique solution to Faraday's Law. Finally, we explore a variational approach that minimizes an energy functional to determine a unique E-field, similar to Longcope's "Minimum Energy Fit". The PTD technique, the iterative technique, and the variational technique are used to estimate E-fields from a pair of synthetic vector magnetograms taken from an MHD simulation; and these E-fields are compared with the simulation's known electric fields. These three techniques are then applied to a pair of vector magnetograms of solar active region NOAA AR8210, to demonstrate the methods with real data.Comment: 41 pages, 10 figure

Supernova 1987A: Rotation and a Binary Companion

In this paper we provide a possible link between the structure of the bipolar nebula surrounding SN1987A and the properties of its progenitor star. A Wind Blwon Bubble (WBB) scenario is emplyed, in which a fast, tenuous wind from a Blue Supergiant expands into a slow, dense wind, expelled during an earlier Red Supergiant phase. The bipolar shapre develops due to a pole-to-equator density contrast in the slow wind (ie, the slow wind forms a slow torus). We use the Wind Compressed Disk (WCD) model of Bjorkman & Cassinelli (1992) to determine the shape of the slow torus. In the WCD scenario, the shape of the torus is determined by the rotation of the progenitor star. We then use a self-similar semi-analytical method for wind blown bubble evolution to determine the shape of the resulting bipolar nebula. We find that the union of the wind-compressed-disk and bipolar-wind-blown- bubble models allows us to recover the salient properties of SN1987A's circumstellar nebula. In particular, the size, speed and density of SN1987A's inner ring are easily reproduced in our calculations. An exploration of parameter space shows the the red supergiant progenitor must be been rotating at > 0.3 of its breakup speed. We conclude that the progenitor was most likely spun up by a merger with a binary companion. Using a simple model for the binary merger we find that the companion is likely to have had a mass > 0.5 M_sun.Comment: 30 pages, 4 figure

Can Protostellar Jets Drive Supersonic Turbulence in Molecular Clouds?

Jets and outflows from young stellar objects are proposed candidates to drive supersonic turbulence in molecular clouds. Here, we present the results from multi-dimensional jet simulations where we investigate in detail the energy and momentum deposition from jets into their surrounding environment and quantify the character of the excited turbulence with velocity probability density functions. Our study include jet--clump interaction, transient jets, and magnetised jets. We find that collimated supersonic jets do not excite supersonic motions far from the vicinity of the jet. Supersonic fluctuations are damped quickly and do not spread into the parent cloud. Instead subsonic, non-compressional modes occupy most of the excited volume. This is a generic feature which can not be fully circumvented by overdense jets or magnetic fields. Nevertheless, jets are able to leave strong imprints in their cloud structure and can disrupt dense clumps. Our results question the ability of collimated jets to sustain supersonic turbulence in molecular clouds.Comment: 33 pages, 18 figures, accepted by ApJ, version with high resolution figures at: http://www.ita.uni-heidelberg.de/~banerjee/publications/jet_paper.pd

Three-dimensional chemically homogeneous and bi-abundance photoionization models of the "super-metal-rich" planetary nebula NGC 6153

Deep spectroscopy of the planetary nebula (PN) NGC\,6153 shows that its heavy element abundances derived from optical recombination lines (ORLs) are ten times higher than those derived from collisionally excited lines (CELs), and points to the existence of H-deficient inclusions embedded in the diffuse nebula. In this study, we have constructed chemically homogeneous and bi-abundance three-dimensional photoionization models, using the Monte Carlo photoionization code {\sc mocassin}. We attempt to reproduce the multi-waveband spectroscopic and imaging observations of NGC\,6153, and investigate the nature and origin of the postulated H-deficient inclusions, as well as their impacts on the empirical nebular analyses assuming a uniform chemical composition. Our results show that chemically homogeneous models yield small electron temperature fluctuations and fail to reproduce the strengths of ORLs from C, N, O and Ne ions. In contrast, bi-abundance models incorporating a small amount of metal-rich inclusions ($\sim 1.3$ per cent of the total nebular mass) are able to match all the observations within the measurement uncertainties. The metal-rich clumps, cooled down to a very low temperature ($\sim 800$~K) by ionic infrared fine-structure lines, dominate the emission of heavy element ORLs, but contribute almost nil to the emission of most CELs. We find that the abundances of C, N, O and Ne derived empirically from CELs, assuming a uniform chemical composition, are about 30 per cent lower than the corresponding average values of the whole nebula, including the contribution from the H-deficient inclusions. Ironically, in the presence of H-deficient inclusions, the traditional standard analysis of the optical helium recombination lines, assuming a chemically homogeneous nebula, overestimates the helium abundance by 40 per cent.Comment: 19 pages, 18 figures, accepted for publication in MNRA