No Arabic abstract
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.
In this paper we focus on understanding the physical processes that lead to stable or unstable ionization fronts (I-fronts) observed in simulations of moving black holes (BHs). The front instability may trigger bursts of gas accretion, rendering the BH significantly more luminous than at steady-state. We perform a series of idealized three dimensional radiation hydrodynamics simulations resolving the I-fronts around BHs of mass $M_mathrm{BH}$ and velocity $v_infty$ accreting from a medium of density $n_mathrm{H}$. The I-front, with radius $R_mathrm{I}$, transitions from D-type to R-type as the BH velocity becomes larger than a critical value $v_mathrm{R}sim 40,mathrm{km/s}$. The D-type front is preceded by a bow-shock of thickness $Delta R_mathrm{I}$ that decreases as $v_infty$ approaches $v_mathrm{R}$. We find that both D-type and R-type fronts can be unstable given the following two conditions: i) for D-type fronts the shell thickness must be $Delta R_mathrm{I}/R_mathrm{I}<0.05$ (i.e., $v_infty gtrsim 20,mathrm{km/s}$.), while no similar restriction holds for R-type fronts; ii) the temperature jump across the I-front must be $T_mathrm{II}/T_mathrm{I}>3$. This second condition is satisfied if $T_mathrm{I}<5000,mathrm{K}$ or if $n_mathrm{H},M_mathrm{BH} gtrsim 10^6,M_odot,mathrm{cm^{-3}}$. Due to X-ray pre-heating typically $T_mathrm{I} sim 10^4,mathrm{K}$, unless the D-type shell is optically thick to X-rays, which also happens when $n_mathrm{H},M_mathrm{BH}$ is greater than a metallicity-dependent critical value. We thus conclude that I-fronts around BHs are unstable only for relatively massive BHs moving trough very dense molecular clouds. We briefly discuss the observational consequences of the X-ray luminosity bursts likely associated with this instability.
We present measurements of the singly ionized helium to hydrogen ratio ($n_{He^+}/n_{H^+}$) toward diffuse gas surrounding three Ultra-Compact HII (UCHII ) regions: G10.15-0.34, G23.46-0.20 & G29.96-0.02. We observe radio recombination lines (RRLs) of hydrogen and helium near 5 GHz using the GBT to measure the $n_{He^+}/n_{H^+}$ ratio. The measurements are motivated by the low helium ionization observed in the warm ionized medium (WIM) and in the inner Galaxy diffuse ionized regions (DIR). Our data indicate that the helium is not uniformly ionized in the three observed sources. Helium lines are not detected toward a few observed positions in sources G10.15-0.34 & G23.46-0.20 and the upper limits of the $n_{He^+}/n_{H^+}$ ratio obtained are 0.03 and 0.05 respectively. The selected sources harbor stars of type O6 or hotter as indicated by helium line detection toward the bright radio continuum emission from the sources with mean $n_{He^+}/n_{H^+}$ value 0.06$pm$0.02. Our data thus show that helium in diffuse gas located a few pc away from the young massive stars embedded in the observed regions is not fully ionized.We investigate the origin of the non-uniform helium ionization and rule out the possibilities : (a) that the helium is doubly ionized in the observed regions and (b) that the low $n_{He^+}/n_{H^+}$ values are due to additional hydrogen ionizing radiation produced by accreting low-mass stars (Smith 2014). We find that selective absorption of ionizing photons by dust can result in low helium ionization but needs further investigation to develop a self-consistent model for dust in HII regions.
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.
Sgr B1 is a luminous H II region in the Galactic Center immediately next to the massive star-forming giant molecular cloud Sgr B2 and apparently connected to it from their similar radial velocities. In 2018 we showed from SOFIA FIFI-LS observations of the [O III] 52 and 88 micron lines that there is no central exciting star cluster and that the ionizing stars must be widely spread throughout the region. Here we present SOFIA FIFI-LS observations of the [O I] 146 and [C II] 158 micron lines formed in the surrounding photodissociation regions (PDRs). We find that these lines correlate neither with each other nor with the [O III] lines although together they correlate better with the 70 micron Herschel PACS images from Hi-GAL. We infer from this that Sgr B1 consists of a number of smaller H II regions plus their associated PDRs, some seen face-on and the others seen more or less edge-on. We used the PDR Toolbox to estimate densities and the far-ultraviolet intensities exciting the PDRs. Using models computed with Cloudy, we demonstrate possible appearances of edge-on PDRs and show that the density difference between the PDR densities and the electron densities estimated from the [O III] line ratios is incompatible with pressure equilibrium unless there is a substantial pressure contribution from either turbulence or magnetic field or both. We likewise conclude that the hot stars exciting Sgr B1 are widely spaced throughout the region at substantial distances from the gas with no evidence of current massive star formation.
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.