No Arabic abstract
The ionization parameter U is potentially useful for measuring radiation pressure feedback from massive star clusters, as it reflects the radiation-to-gas-pressure ratio and is readily derived from mid-infrared line ratios. We consider several effects which determine the apparent value of U in HII regions and galaxies. An upper limit is set by the compression of gas by radiation pressure. The pressure from stellar winds and the presence of neutral clumps both reduce U for a given radiation intensity. The most intensely irradiated regions are selectively dimmed by internal dust absorption of ionizing photons, inducing observational bias on galactic scales. We explore these effects analytically and numerically, and use them to interpret previous observational results. We find that radiation confinement sets the upper limit log_10 U = -1 seen in individual regions. Unresolved starbursts display a maximum value of ~ -2.3. While lower, this is also consistent with a large portion of their HII regions being radiation dominated, given the different technique used to interpret unresolved regions, and given the bias caused by dust absorption. We infer that many individual, strongly illuminated regions cannot be dominated by stellar winds, and that even when averaged on galactic scales, shocked wind pressures cannot be large compared to radiation pressure. Therefore, most HII regions cannot be adiabatic wind bubbles. Our models imply a metallicity dependence in the physical structure and dust attenuation of radiation-dominated regions, both of which should vary strongly across a critical metallicity of about one-twentieth solar.
We present our parameterizations of the log([NeIII]3869/[OII]3727) (Ne3O2) and log([OIII]5007/[OII]3727) ratios as comparable and effective diagnostics of ionization parameter in star-forming galaxies. Our calibrations are based on the most recent generations of the Starburst99/Mappings III photoionization models, which extend up to the extremely high values of ionization parameter found in high-redshift galaxies. While similar calibrations have been presented previously for O3O2, this is the first such calibration of Ne3O2. We illustrate the tight correlation between these two ratios for star-forming galaxies and discuss the underlying physics that dictates their very similar evolution. Based on this work, we propose the Ne3O2 ratio as a new and useful diagnostic of ionization parameter for star-forming galaxies. Given the Ne3O2 ratios relative insensitivity to reddening, this ratio is particularly valuable for use with galaxies that have uncertain amounts of extinction. The short wavelengths of the Ne3O2 ratio can also be applied out to very high redshifts, extending studies of galaxies ionization parameters out to z ~ 1.6 with optical spectroscopy and z ~ 5.2 with ground-based near-infrared spectra.
An ionization front (IF) surrounding an H II region is a sharp interface where a cold neutral gas makes transition to a warm ionized phase by absorbing UV photons from central stars. We investigate the instability of a plane-parallel D-type IF threaded by parallel magnetic fields, by neglecting the effects of recombination within the ionized gas. We find that weak D-type IFs always have the post-IF magnetosonic Mach number $mathcal{M}_{rm M2} leq 1$. For such fronts, magnetic fields increase the maximum propagation speed of the IFs, while reducing the expansion factor $alpha$ by a factor of $1+1/(2beta_1)$ compared to the unmagnetized case, with $beta_1$ denoting the plasma beta in the pre-IF region. IFs become unstable to distortional perturbations due to gas expansion across the fronts, exactly analogous to the Darrieus-Landau instability of ablation fronts in terrestrial flames. The growth rate of the IF instability is proportional linearly to the perturbation wavenumber as well as the upstream flow speed, and approximately to $alpha^{1/2}$. The IF instability is stabilized by gas compressibility and becomes completely quenched when the front is D-critical. The instability is also stabilized by magnetic pressure when the perturbations propagate in the direction perpendicular to the fields. When the perturbations propagate in the direction parallel to the fields, on the other hand, it is magnetic tension that reduces the growth rate, completely suppressing the instability when $mathcal{M}_{rm M2}^2 < 2/(beta_1 - 1)$. When the front experiences an acceleration, the IF instability cooperates with the Rayleigh-Taylor instability to make the front more unstable.
We employ ionization-parameter mapping (IPM) to infer the optical depth of HII regions in the northern half of M33. We construct [OIII]$lambda 5007$/[OII]$lambda 3727$ and [OIII]$lambda 5007$/[SII]$lambda 6724$ ratio maps from narrow-band images continuum-subtracted in this way, from which we classify the HII regions by optical depth to ionizing radiation, based on their ionization structure. This method works relatively well in the low metallicity regime, $12 + log(rm O/H) leq 8.4$, where [OIII]$lambdalambda4949,5007$ is strong. However, at higher metallicities, the method breaks down due to the strong dependence of the [OIII]$lambdalambda4959,5007$ emission lines on the nebular temperature. Thus, although O$^{++}$ may be present in metal-rich HII regions, these commonly used emission lines do not serve as a useful indicator of its presence, and hence, the O ionization state. In addition, IPM as a diagnostic of optical depth is limited by spatial resolution. We also report a region of highly excited [OIII] extending over an area $sim$ 1 kpc across and [OIII]$lambda5007$ luminosity of $4.9pm 1.5times10^{38}$ erg/s, which is several times higher than the ionizing budget of any potential sources in this portion of the galaxy. Finally, this work introduces a new method for continuum subtraction of narrow-band images based on the dispersion of pixels around the mode of the diffuse-light flux distribution. In addition to M33, we demonstrate the method on C III]$lambda$1909 imaging of Haro~11, ESO 338-IG004, and Mrk~71.
We compare the accuracy of various methods for determining the transfer of the diffuse Lyman continuum in HII regions, by comparing them with a high-resolution discrete-ordinate integration. We use these results to suggest how, in multidimensional dynamical simulations, the diffuse field may be treated with acceptable accuracy without requiring detailed transport solutions. The angular distribution of the diffuse field derived from the numerical integration provides insight into the likely effects of the diffuse field for various material distributions.
Interstellar bubbles around O stars are driven by a combination of the stars wind and ionizing radiation output. The wind contribution is uncertain because the boundary between the wind and interstellar medium is difficult to observe. Mid-infrared observations (e.g., of the H II region RCW 120) show arcs of dust emission around O stars, contained well within the H II region bubble. These arcs could indicate the edge of an asymmetric stellar wind bubble, distorted by density gradients and/or stellar motion. We present two-dimensional, radiation-hydrodynamics simulations investigating the evolution of wind bubbles and H II regions around massive stars moving through a dense (n=3000 cm^{-3}), uniform medium with velocities ranging from 4 to 16 km/s. The H II region morphology is strongly affected by stellar motion, as expected, but the wind bubble is also very aspherical from birth, even for the lowest space velocity considered. Wind bubbles do not fill their H II regions (we find filling factors of 10-20%), at least for a main sequence star with mass M~30 Msun. Furthermore, even for supersonic velocities the wind bow shock does not significantly trap the ionization front. X-ray emission from the wind bubble is soft, faint, and comes mainly from the turbulent mixing layer between the wind bubble and the H II region. The wind bubble radiates <1 per cent of its energy in X-rays; it loses most of its energy by turbulent mixing with cooler photoionized gas. Comparison of the simulations with the H II region RCW 120 shows that its dynamical age is <=0.4 Myr and that stellar motion <=4 km/s is allowed, implying that the ionizing source is unlikely to be a runaway star but more likely formed in situ. The regions youth, and apparent isolation from other O or B stars, makes it very interesting for studies of massive star formation and of initial mass functions.