We present the first fully calibrated H$_2$, 1-0 S(1) image of the entire 30 Doradus nebula. The observations were conducted using the NOAO Extremely Wide-Field Infrared Imager on the CTIO 4-meter Blanco Telescope. Together with a NEWFIRM Br$gamma$ i
mage of 30 Doradus, our data reveal the morphologies of the warm molecular gas and ionized gas in 30 Doradus. The brightest H$_2$-emitting area, which extends from the northeast to the southwest of R136, is a photodissociation region viewed face-on, while many clumps and pillar features located at the outer shells of 30 Doradus are photodissociation regions viewed edge-on. Based on the morphologies of H$_2$, Br$gamma$, $^{12}$CO, and 8$mu$m emission, the H$_2$ to Br$gamma$ line ratio and Cloudy models, we find that the H$_2$ emission is formed inside the photodissociation regions of 30 Doradus, 2 - 3 pc to the ionization front of the HII region, in a relatively low-density environment $<$ 10$^4$ cm$^{-3}$. Comparisons with Br$gamma$, 8$mu$m, and CO emission indicate that H$_2$ emission is due to fluorescence, and provide no evidence for shock excited emission of this line.
The emission line ratios [OIII]5007/H-beta and [NII]6584/H-alpha have been adopted as an empirical way to distinguish between the fundamentally different mechanisms of ionization in emission-line galaxies. However, detailed interpretation of these di
agnostics requires calculations of the internal structure of the emitting HII regions, and these calculations depend on the assumptions one makes about the relative importance of radiation pressure and stellar winds. In this paper we construct a grid of quasi-static HII region models to explore how choices about these parameters alter HII regions emission line ratios. We find that, when radiation pressure is included in our models, HII regions reach a saturation point beyond which further increases in the luminosity of the driving stars does not produce any further increase in effective ionization parameter, and thus does not yield any further alteration in an HII regions line ratio. We also show that, if stellar winds are assumed to be strong, the maximum possible ionization parameter is quite low. As a result of this effect, it is inconsistent to simultaneously assume that HII regions are wind-blown bubbles and that they have high ionization parameters; some popular HII region models suffer from this inconsistency. Our work in this paper provides a foundation for a companion paper in which we embed the model grids we compute here within a population synthesis code that enables us to compute the integrated line emission from galactic populations of HII regions.
Optical and infrared emission lines from HII regions are an important diagnostic used to study galaxies, but interpretation of these lines requires significant modeling of both the internal structure and dynamical evolution of the emitting regions. M
ost of the models in common use today assume that HII region dynamics are dominated by the expansion of stellar wind bubbles, and have neglected the contribution of radiation pressure to the dynamics, and in some cases also to the internal structure. However, recent observations of nearby galaxies suggest that neither assumption is justified, motivating us to revisit the question of how HII region line emission depends on the physics of winds and radiation pressure. In a companion paper we construct models of single HII regions including and excluding radiation pressure and winds, and in this paper we describe a population synthesis code that uses these models to simulate galactic collections of HII regions with varying physical parameters. We show that the choice of physical parameters has significant effects on galactic emission line ratios, and that in some cases the line ratios can exceed previously claimed theoretical limits. Our results suggest that the recently-reported offset in line ratio values between high-redshift star-forming galaxies and those in the local universe may be partially explained by the presence of large numbers of radiation pressured-dominated HII regions within them.
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 effect
s 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 have mapped the proto-binary source IRAS 16293-2422 in CO 2-1, 13CO 2-1, and CO 3-2 with the Submillimeter Array (SMA). The maps with resolution of 1.5-5 reveal a single small scale (~3000 AU) bipolar molecular outflow along the east-west directio
n. We found that the blueshifted emission of this small scale outflow mainly extends to the east and the redshifted emission to the west from the position of IRAS 16293A. A comparison with the morphology of the large scale outflows previously observed by single-dish telescopes at millimeter wavelengths suggests that the small scale outflow may be the inner part of the large scale (~15000 AU) E-W outflow. On the other hand, there is no clear counterpart of the large scale NE-SW outflow in our SMA maps. Comparing analytical models to the data suggests that the morphology and kinematics of the small scale outflow can be explained by a wide-angle wind with an inclination angle of ~30-40 degrees with respect to the plane of the sky. The high resolution CO maps show that there are two compact, bright spots in the blueshifted velocity range. An LVG analysis shows that the one located 1 to the east of source A is extremely dense, n(H_2)~10^7 cm^-3, and warm, T_kin >55 K. The other one located 1 southeast of source B has a higher temperature of T_kin >65 K but slightly lower density of n(H_2)~10^6 cm^-3. It is likely that these bright spots are associated with the hot core-like emission observed toward IRAS 16293. Since both two bright spots are blueshifted from the systemic velocity and are offset from the protostellar positions, they are likely formed by shocks.
