[Si II], [Fe II], [C II], and H2 emission from massive star-forming regions

We calculate the [Si II] 34.8 mu m, [Fe II] 26.0 mu m, and [C II] 158 mu m infrared fine-structure emission that may arise from ionized gas (H II regions) and/or associated neutral gas (photodissociation regions [PDRs]) in massive star-forming regions. Assuming thermal pressure balance between an H II region and a PDR, the relative amounts of fine-structure line emission from the H II region and PDR depend on the electron density and resulting thermal pressure in the H II region, gas-phase abundances of the emitting species, and the UV spectrum from the stellar population producing the H II region. For normal metallicity, we find that [C II] emission is always dominated by PDRs, while [Si II] and [Fe II] are dominated by PDRs for H II regions with electron density n(e) greater than or similar to 10 cm(-3). We also calculate the H(2) 0-0 S(0), 0-0 S(1), 0-0 S(2), and 0-0 S(3) pure rotational line emission arising from the PDR at the atomic-to-molecular interface. The overall intensity of the H(2) line emission directly traces warm molecular mass, while H(2) line ratios constrain the PDR temperature, gas density, and far-ultraviolet field strength. Models of the integrated emission of [Si II], [Fe II], [C II], and H(2) from Galactic and extragalactic star-forming regions are presented for use in interpreting observations with Spitzer, ISO, SOFIA, and the Herschel Space Observatory. We compare our results with observations of the Galactic source NGC 2023 (an individual H II/photodissociation region in Orion), the inner regions of the Milky Way, and the central regions of the nearby star-forming spiral galaxy NGC 7331. We also compare our results with recently published similar work by Abel and coworkers.