Based on imaging and spectroscopic data, we develop a 3-D model for the Huygens Region of the Orion Nebula. Theta1OriC , the hottest star in the Trapezium, is surrounded by a wind-blown Central Bubble that opens SW into the Extended Orion Nebula. Outside of this feature lies a layer of ionized gas at about 0.4 pc from Theta1OriC. Both of these features are moving rapidly away from Theta1OriC with an expansion age for the Central Bubble of only 15,000 yrs.

(abridged) The ambiguous origin of [CII] 158um in the interstellar medium complicates its use for diagnostics concerning the star-formation rate and physical conditions in photodissociation regions (PDRs). We observed the giant HII region N11 in the Large Magellanic Cloud with SOFIA/GREAT in order to investigate the origin of [CII] to obtain the total H2 gas content, the fraction of CO-dark H2 gas, and the influence of environmental effects such as stellar feedback. We present an innovative spectral decomposition method that allows statistical trends to be derived. The [CII] line is resolved in velocity and compared to HI and CO, using a Bayesian approach to decompose the profiles. A simple model accounting for collisions in the neutral atomic and molecular gas was used in order to derive the H2 column density traced by C+. The profile of [CII] most closely resembles that of CO, but the integrated [CII] line width lies between that of CO and that of HI. Using various methods, we find that [CII] mostly originates from the neutral gas. We show that [CII] mostly traces the CO-dark H2 gas but there is evidence of a weak contribution from neutral atomic gas preferentially in the faintest components. Most of the molecular gas is CO-dark. The fraction of CO-dark H2 gas decreases with increasing CO column density, with a slope that seems to depend on the impinging radiation field from nearby massive stars. Finally we extend previous measurements of the photoelectric-effect heating efficiency, which we find is constant across regions probed with Herschel, with [CII] and [OI] being the main coolants in faint and diffuse, and bright and compact regions, respectively, and with PAH emission tracing the CO-dark H2 gas heating where [CII] and [OI] emit. Our study highlights the importance of velocity-resolved PDR diagnostics and higher spatial resolution for HI observations.

We have used existing optical emission and absorption lines, [C II] emission lines, and H I absorption lines to create a new model for a Central Column of material near the Trapezium region of the Orion Nebula. This was necessary because recent high spectral resolution spectra of optical emission lines and imaging spectra in the [C II] 158 micron line have shown that there are new velocity systems associated with the foreground Veil and the material lying between Theta 1 Ori C and the Main Ionization Front of the nebula. When a family of models generated with the spectral synthesis code Cloudy were compared with the surface brightness of the emission lines and strengths of the Veil absorption lines seen in the Trapezium stars, distances from Theta 1 Ori C, were derived, with the closest, highest ionization layer being 1.3 pc. The line of sight distance of this layer is comparable with the size of the inner Huygens Region in the plane of the sky. These layers are all blueshifted with respect to the Orion Nebula Cluster of stars, probably because of the pressure of a hot central bubble created by Theta 1 Ori C's stellar wind. We find velocity components that are ascribed to both sides of this bubble. Our analysis shows that the foreground [C II] 158 micron emission is part of a previously identified layer that forms a portion of a recently discovered expanding shell of material covering most of the larger Extended Orion Nebula.

The sensitive infrared telescopes, Spitzer and Herschel, have been used to target low-metallicity star-forming galaxies, allowing us to investigate the properties of their interstellar medium (ISM) in unprecedented detail. Interpretation of the observations in physical terms relies on careful modeling of those properties. We have employed a multiphase approach to model the ISM phases (HII region and photodissociation region) with the spectral synthesis code Cloudy. Our goal is to characterize the physical conditions (gas densities, radiation fields, etc.) in the ISM of the galaxies from the Herschel Dwarf Galaxy Survey. We are particularly interested in correlations between those physical conditions and metallicity or star-formation rate. Other key issues we have addressed are the contribution of different ISM phases to the total line emission, especially of the [CII]157um line, and the characterization of the porosity of the ISM. We find that the lower-metallicity galaxies of our sample tend to have higher ionization parameters and galaxies with higher specific star-formation rates have higher gas densities. The [CII] emission arises mainly from PDRs and the contribution from the ionized gas phases is small, typically less than 30% of the observed emission. We also find correlation - though with scatter - between metallicity and both the PDR covering factor and the fraction of [CII] from the ionized gas. Overall, the low metal abundances appear to be driving most of the changes in the ISM structure and conditions of these galaxies, and not the high specific star-formation rates. These results demonstrate in a quantitative way the increase of ISM porosity at low metallicity. Such porosity may be typical of galaxies in the young Universe.

Intermediate-mass black holes (IMBHs) with masses between $100 - 10^5M_{\odot}$ are crucial to our understanding of black hole seed formation and are the prime targets for LISA, yet black holes in this mass range have eluded detection by traditional optical spectroscopic surveys aimed at finding active galactic nuclei (AGNs). In this paper, we have modeled for the first time the dependence of the optical narrow emission line strengths on the black hole mass of accreting AGN over the range of $100-10^8M_{\odot}$. We show that as the black hole mass decreases, the hardening of the spectral energy distribution from the accretion disk changes the ionization structure of the nebula. The enhanced high energy emission from IMBHs results in a more extended partially ionized zone compared with models for higher mass black holes. This effect produces a net decrease in the predicted [OIII]/H$\beta$ and [NII]/H$\alpha$ emission line ratios. Based on this model, we demonstrate that the standard optical narrow emission line diagnostics used to identify massive black holes fail when black hole mass falls below $\approx10^4M_{\odot}$ for highly accreting IMBHs and for radiatively inefficient IMBHs with active star formation. Our models call into question the ability of common optical spectroscopic diagnostics to confirm AGN candidates in dwarf galaxies, and indicate that the low-mass black hole occupation fraction inferred from such diagnostics will be severely biased.

Our objective is to investigate the physical properties of the ionised gas of the low-metallicity dwarf galaxy, IC 10, at various spatial scales: from individual HII regions to the entire galaxy scale and examine whether diagnostics for integrated measurements introduce bias in the results. We modeled the ionised gas combining the mid- and far-infrared fine-structure cooling lines observed with Spitzer/IRS and Herschel/PACS, with the photoionisation code Cloudy. The free parameters of the models are the age of the stellar cluster, the density and the ionisation parameter of the ionised gas as well as the depth of the cloud. The latter is used to investigate the leakage of the ionising photons from the analysed regions of IC 10. We investigate HII regions in the main star-forming body, on scales of ~25 pc, three in the main star-forming region in the center of the galaxy and two on the first arc. We then consider larger sizes on the scale of ~200 pc. We find that most clumps have nearly identical properties, density ~10$^{2.} $ - 10$^{2.6}$ cm$^{-3}$, ionisation parameter between 10$^{-2.2}$ and 10$^{-1.6}$ and age of the stellar cluster ~5.5 Myr. All of them are matter-bounded regions, allowing ionising photons to leak. The relatively uniform physical properties of the clumps suggest a common origin for their star formation activity, which could be related to the feedback from stellar winds or supernovae of a previous generation of stars. The properties derived for ~200 pc size "zones" have similar properties as the HII regions they encompass, but with the larger regions tending to be more radiation-bounded. Finally, we investigate the fraction of [CII] 157.7 {\mu}m, [SiII] 34.8 {\mu}m, and [FeII] 25.9 {\mu}m, emission arising from the ionised gas phase and we find that most of the emission originates from the neutral gas, not from the ionised gas.

Intermediate mass black holes (IMBHs), with masses between 100 to 10^5 M_\odot, represent the link between stellar mass black holes and the supermassive black holes that reside in galaxy centers. While IMBHs are crucial to our understanding of black hole seed formation, black holes of less than \approx 10^4 M_\odot eluded detection by traditional searches. Observations of the infrared coronal lines (CLs) offer us one of the most promising tools to discover IMBHs in galaxies. We have modeled the infrared emission line spectrum that is produced by gas photoionized by an AGN radiation field and explored for the first time the dependence of the infrared CL spectrum on black hole mass over the range of 10^2 M_\odot to 10^8 M_\odot. We show that infrared CLs are expected to be prominent in the spectra of accreting IMBHs and can potentially be a powerful probe of the black hole mass in AGNs. We identify key emission line ratios in the 1-30 \mu m range that are most sensitive to black hole mass. While variations in accretion rate and the physical parameters of the gas can also affect the CL spectrum, we demonstrate that the effect of black hole mass is likely to be the most dramatic over the mass range explored in our models. With the unprecedented sensitivity of JWST, a large number of CLs will be detectable for the first time, providing important insight into the existence and properties of IMBHs in the local universe, potentially revolutionizing our understanding of this class of object.

We conduct for the first time a theoretical investigation of the mid-infrared spectral energy distribution (SED) produced by dust heated by an active galactic nucleus (AGN) and an extreme starburst. These models employ an integrated modeling approach using photoionization and stellar population synthesis models in which both the line and emergent continuum is predicted from gas exposed to the ionizing radiation from a young starburst and an AGN. In this work, we focus on the infrared colors from the {\it Wide-field Infrared Survey Explorer}, predicting the dependence of the colors on the input radiation field, the ISM conditions, the obscuring column, and the metallicity. We find that an extreme starburst can mimic an AGN in two band mid-infrared color cuts employed in the literature. However, the three band color cuts employed in the literature require starbursts with extremely high ionization parameters or gas densities . We show that the extreme mid-IR colors seen in some blue compact dwarf galaxies are not due to metallicity but rather a combination of high ionization parameters and high column densities. Based on our theoretical calculations, we present a theoretical mid-infrared color cut that will exclude even the most extreme starburst that we have modeled in this work. The theoretical AGN demarcation region presented here can be used to identify elusive AGN candidates for future follow-up studies with the {\it James Webb Space Telescope (JWST)}. The full suite of simulated SEDs are available online.

The southwestern (SW) part of the Galactic H II region M17 contains an obscured ionization front that is most easily seen at infrared and radio wavelengths. It is nearly edge-on, thus offering an excellent opportunity to study the way in which the gas changes from fully ionized to molecular as radiation from the ionizing stars penetrates into the gas. M17 is also one of the very few H II regions for which the magnetic field strength can be measured in the photodissociation region ( PDR) that forms the interface between the ionized and molecular gas. Here we model an observed line of sight through the gas cloud, including the H+, H0 (PDR), and molecular layers, in a fully self-consistent single calculation. An interesting aspect of the M17 SW bar is that the PDR is very extended. We show that the strong magnetic field that is observed to be present inevitably leads to a very deep PDR, because the structure of the neutral and molecular gas is dominated by magnetic pressure, rather than by gas pressure, as previously had been supposed.We also show that a wide variety of observed facts can be explained if a hydrostatic geometry prevails, in which the gas pressure from an inner X-ray hot bubble and the outward momentum of the stellar radiation field compress the gas and its associated magnetic field in the PDR, as has already been shown to occur in the Orion Nebula. The magnetic field compression may also amplify the local cosmic-ray density. The pressure in the observed magnetic field balances the outward forces, suggesting that the observed geometry is a natural consequence of the formation of a star cluster within a molecular cloud.

We describe the 2017 release of the spectral synthesis code Cloudy. A major development since the previous release has been exporting the atomic data into external data files. This greatly simplifies updates and maintenance of the data. Many large datasets have been incorporated with the result that we can now predict well over an order of magnitude more emission lines when all databases are fully used. The use of such large datasets is not realistic for most calculations due to the time and memory needs, and we describe the limited subset of data we use by default. Despite the fact that we now predict significantly more lines than the previous Cloudy release, this version is faster because of optimization of memory access patterns and other tuning. The size and use of the databases can easily be adjusted in the command-line interface. We give examples of the accuracy limits using small models, and the performance requirements of large complete models. We summarize several advances in the H- and He-like iso-electronic sequences. We use our complete collisional-radiative models of the ionization of these one and two-electron ions to establish the highest density for which the coronal or interstellar medium (ISM) approximation works, and the lowest density where Saha or local thermodynamic equilibrium can be assumed. The coronal approximation fails at surprisingly low densities for collisional ionization equilibrium but is valid to higher densities for photoionized gas clouds. Many other improvements to the physics have been made and are described. These include the treatment of isotropic continuum sources such as the cosmic microwave background (CMB) in the reported output, and the ability to follow the evolution of cooling non-equilibrium clouds.

Properties of the cold interstellar medium of low-metallicity galaxies are not well-known due to the faintness and extremely small scale on which emission is expected. We present deep ALMA band 6 (230GHz) observations of the nearby, low-metallicity (12 + log(O/H) = 7.25) blue compact dwarf galaxy SBS0335-052 at an unprecedented resolution of 0.2 arcsec (52 pc). The 12CO J=2-1 line is not detected and we report a 3-sigma upper limit of LCO(2-1) = 3.6x10^4 K km/s pc^2. Assuming that molecular gas is converted into stars with a given depletion time, ranging from 0.02 to 2 Gyr, we find lower limits on the CO-to-H2 conversion factor alpha_CO in the range 10^2-10^4 Msun pc^-2 (K km/s)^-1. The continuum emission is detected and resolved over the two main super star clusters. Re-analysis of the IR-radio spectral energy distribution suggests that the mm-fluxes are not only free-free emission but are most likely also associated with a cold dust component coincident with the position of the brightest cluster. With standard dust properties, we estimate its mass to be as large as 10^5 Msun. Both line and continuum results suggest the presence of a large cold gas reservoir unseen in CO even with ALMA.

We investigate the physical conditions of the gas, atomic and molecular, in the filaments in the context of Photo-Dissociation Regions (PDRs) using the KOSMA-PDR mode of clumpy clouds. We also compare the [CII] vs. [NII] integrated intensity predictions in Abel et al. 2005 for HII regions and adjacent PDRs in the Galactic disk, and check for their applicability under the extreme physical conditions present in the GC. Our preliminary results show that observed integrated intensities are well reproduced by the PDR model. The gas is exposed to a relatively low Far-UV field between 10^2 - 10^3 Draine fields. The total volume hydrogen density is well constrained between 10^4 - 10^5 cm^-3. The hydrogen ionization rate due to cosmic-rays varies between 10^-15 and 4x10^-15 s^-1, with the highest value ~ 10^-14 s^-1 found towards G0.07+0.04. Our results show that the line-of-sight contribution to the total distance of the filaments to the Arches Cluster is not negligible. The spatial distribution of the [CII]/[NII] ratio shows that the integrated intensity ratios are fairly homogeneously distributed for values below 10 in energy units. Calculations including variation on the [C/N] abundance ratio show that tight constraints on this ratio are needed to reproduce the observations.

The [CII] 158 um fine structure line is one of the dominant cooling lines in the interstellar medium (ISM) and is an important tracer of star formation. Recent velocity-resolved studies with Herschel/HIFI and SOFIA/GREAT showed that the [CII] line can constrain the properties of the ISM phases in star-forming regions. The [CII] line as a tracer of star formation is particularly important in low-metallicity environments where CO emission is weak because of the presence of large amounts of CO-dark gas. The nearby irregular dwarf galaxy NGC 4214 offers an excellent opportunity to study an actively star-forming ISM at low metallicity. We analyzed the spectrally resolved [CII] line profiles in three distinct regions at different evolutionary stages of NGC 4214 with respect to ancillary HI and CO data in order to study the origin of the [CII] line. We used SOFIA/GREAT [CII] 158 um observations, HI data from THINGS, and CO(2-1) data from HERACLES to decompose the spectrally resolved [CII] line profiles into components associated with neutral atomic and molecular gas. We use this decomposition to infer gas masses traced by [CII] under different ISM conditions. Averaged over all regions, we associate about 46% of the [CII] emission with the HI emission. However, we can assign only around 9% of the total [CII] emission to the cold neutral medium (CNM). We found that about 79% of the total molecular hydrogen mass is not traced by CO emission. On average, the fraction of CO-dark gas dominates the molecular gas mass budget. The fraction seems to depend on the evolutionary stage of the regions: it is highest in the region covering a super star cluster in NGC 4214, while it is lower in a more compact, more metal-rich region.

We present newly processed archival Herschel images of molecular cloud MCLD 123.5+24.9 in the Polaris Flare. This cloud contains five starless cores. Using the spectral synthesis code Cloudy, we explore uncertainties in the derivation of column densities, hence, masses of molecular cores from Herschel data. We first consider several detailed grain models that predict far-IR grain opacities. Opacities predicted by the models differ by more than a factor of two, leading to uncertainties in derived column densities by the same factor. Then we consider uncertainties associated with the modified blackbody fitting process used by observers to estimate column densities. For high column density clouds (N(H) $\gg$ 10$^{22}$ cm$^{-2}$), this fitting technique can underestimate column densities by about a factor of three. Finally, we consider the virial stability of the five starless cores in MCLD 123.5+24.9. All of these cores appear to have strongly sub-virial masses, assuming, as we argue, that $^{13}$CO line data provide reliable estimates of velocity dispersions. Evidently, they are not self-gravitating, so it is no surprise that they are starless.

The far-infrared (FIR) lines are key tracers of the physical conditions of the interstellar medium (ISM) and are becoming workhorse diagnostics for galaxies throughout the universe. Our goal is to explain the differences and trends observed in the FIR line emission of dwarf galaxies compared to more metal-rich galaxies. We present Herschel PACS spectroscopic observations of the CII157um, OI63 and 145um, OIII88um, NII122 and 205um, and NIII57um fine-structure cooling lines in a sample of 48 low-metallicity star-forming galaxies of the guaranteed time key program Dwarf Galaxy Survey. We correlate PACS line ratios and line-to-LTIR ratios with LTIR, LTIR/LB, metallicity, and FIR color, and interpret the observed trends in terms of ISM conditions and phase filling factors with Cloudy radiative transfer models. We find that the FIR lines together account for up to 3 percent of LTIR and that star-forming regions dominate the overall emission in dwarf galaxies. Compared to metal-rich galaxies, the ratios of OIII/NII122 and NIII/NII122 are high, indicative of hard radiation fields. In the photodissociation region (PDR), the CII/OI63 ratio is slightly higher than in metal-rich galaxies, with a small increase with metallicity, and the OI145/OI63 ratio is generally lower than 0.1, demonstrating that optical depth effects should be small on the scales probed. The OIII/OI63 ratio can be used as an indicator of the ionized gas/PDR filling factor, and is found ~4 times higher in the dwarfs than in metal-rich galaxies. The high CII/LTIR, OI/LTIR, and OIII/LTIR ratios, which decrease with increasing LTIR and LTIR/LB, are interpreted as a combination of moderate FUV fields and low PDR covering factor. Harboring compact phases of low filling factor and a large volume filling factor of diffuse gas, the ISM of low-metallicity dwarf galaxies has a more porous structure than that in metal-rich galaxies.

We present a framework for the interpretation of the far-infrared spectra of galaxies in which we have expanded the model parameters compared with previous work by varying the ionization parameter $U$, column density $N$(H), and gas density at the cloud face $n$(H$^{+}$) for a central starburst or AGN. We compare these models carried out with the $Cloudy$ spectral synthesis code to trends in line-to-total far-infrared luminosity ratios, far-infrared fine-structure line ratios, IRAS colors, and OH and H$_{2}$O column densities with trends found in the well-studied sample of ten nearby galaxies from the IRAS Bright Galaxy Sample with infrared luminosities greater than 10$^{10}$ L$_{\odot}$ and IRAS 60 micron fluxes equal to or greater than that of the nearby ULIRG Arp 220. We find that the spectral sequence extending from normal starburst-type emission line spectra to ULIRG-type absorption-dominated spectra with significant absorption from excited levels, can be best explained by simultaneously increasing the hydrogen column density, from as low as 10$^{21}$ cm$^{-2}$ to as high as 10$^{24.8}$ cm$^{-2}$ or greater, and the ionization parameter, from as low as 10$^{-4}$ to as high as 1. The starburst models best reproduce most of the sequence, while AGN models are somewhat better able to produce the high OH and H$_{2}$O column densities in Arp 220. Our results suggest that the molecular ISM in ULIRG-like, absorption-dominated systems is located close to and at least partially obscures the source of power throughout much of the far-infrared, which must be taken into account in order to properly interpret diagnostics of both their sources of power and of feedback.

We aim to quantify the molecular gas reservoir in a subset of 6 low-metallicity galaxies from the Herschel Dwarf Galaxy Survey with newly acquired CO data, and link this reservoir to the observed star formation activity. We present CO(1-0), CO(2-1), and CO(3-2) observations obtained at the ATNF Mopra 22-m, APEX, and IRAM 30-m telescopes, as well as [CII] 157um and [OI] 63um observations obtained with the Herschel/PACS spectrometer in the 6 galaxies: Haro11, Mrk1089, Mrk930, NGC4861, NGC625, and UM311. We derive molecular gas mass from several methods including the use of the CO-to-H2 conversion factor Xco (both Galactic and metallicity-scaled values) and of dust measurements. The molecular and atomic gas reservoirs are compared to the star formation activity. We also constrain the physical conditions of the molecular clouds using the non-LTE code RADEX and the spectral synthesis code Cloudy. We detect CO in 5 of the 6 galaxies, including first detections in Haro11 (Z~0.4 Zsun), Mrk930 (0.2 Zsun), and UM311 (0.5 Zsun), but CO remains undetected in NGC4861 (0.2 Zsun). The CO luminosities are low while [CII] is bright in these galaxies, resulting in [CII]/CO(1-0)>10000. Our dwarf galaxies are in relatively good agreement with the Schmidt-Kennicutt relation for total gas. They show short molecular depletion time scales, even when considering metallicity-scaled Xco factors. Those galaxies are dominated by their HI gas, except Haro11 which has high star formation efficiency and is dominated by ionized and molecular gas. We determine the mass of each ISM phase in Haro11 using Cloudy and estimate an equivalent Xco factor which is 10 times higher than the Galactic value. Overall, our results confirm the emerging picture that CO suffers from significant selective photodissociation in low-metallicity dwarf galaxies.

This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime. This paper summarizes advances made since the last major review in 1998. Much of the recent development has emphasized dusty molecular environments, improvements to the ionization / chemistry solvers, and how atomic and molecular data are used. We present two types of simulations to demonstrate the capability of the code. We consider a molecular cloud irradiated by an X-ray source such as an Active Nucleus and show how treating EUV recombination lines and the full SED affects the observed spectrum. A second example illustrates the very wide range of particle and radiation density that can be considered.

(Abridged) The Magellanic Clouds provide a nearby laboratory for metal-poor dwarf galaxies. The low dust abundance enhances the penetration of UV photons into the interstellar medium (ISM), resulting in a relatively larger filling factor of the ionized gas. Furthermore, there is likely a hidden molecular gas reservoir probed by the [CII]157um line. We present Herschel/PACS maps in several tracers, [CII], [OI]63um,145um, [NII]122um, [NIII]57um, and [OIII]88um in the HII region N11B in the Large Magellanic Cloud. Halpha and [OIII]5007A images were used as complementary data to investigate the effect of dust extinction. Observations were interpreted with photoionization models to infer the gas conditions and estimate the ionized gas contribution to the [CII] emission. Photodissociation regions (PDRs) are probed through polycyclic aromatic hydrocarbons (PAHs). We first study the distribution and properties of the ionized gas. We then constrain the origin of [CII]157um by comparing to tracers of the low-excitation ionized gas and of PDRs. [OIII] is dominated by extended emission from the high-excitation diffuse ionized gas; it is the brightest far-infrared line, ~4 times brighter than [CII]. The extent of the [OIII] emission suggests that the medium is rather fragmented, allowing far-UV photons to permeate into the ISM to scales of >30pc. Furthermore, by comparing [CII] with [NII], we find that 95% of [CII] arises in PDRs, except toward the stellar cluster for which as much as 15% could arise in the ionized gas. We find a remarkable correlation between [CII]+[OI] and PAH emission, with [CII] dominating the cooling in diffuse PDRs and [OI] dominating in the densest PDRs. The combination of [CII] and [OI] provides a proxy for the total gas cooling in PDRs. Our results suggest that PAH emission describes better the PDR gas heating as compared to the total infrared emission.

(abridged) Our goal is to describe the multi-phase ISM of the IR bright low-metallicity galaxy Haro 11, dissecting the photoionised and photodissociated gas components. We present observations of the mid- and far-IR fine-structure cooling lines obtained with the Spitzer/IRS and Herschel/PACS spectrometers. We use the spectral synthesis code Cloudy to methodically model the ionised and neutral gas from which these lines originate. We find that the mid- and far-IR lines account for ~1% of the total IR luminosity L_TIR. Haro 11 is undergoing a phase of intense star formation, as traced by the brightest line [OIII] 88um, with L_[OIII]/L_TIR ~0.3%, and high ratios of [NeIII]/[NeII] and [SIV]/[SIII]. Due to their different origins, the observed lines require a multi-phase modeling comprising: a compact HII region, dense fragmented photodissociation regions (PDRs), a diffuse extended low-ionisation/neutral gas which has a volume filling factor of at least 90%, and porous warm dust in proximity to the stellar source. For a more realistic picture of the ISM of Haro 11 we would need to model the clumpy source and gas structures. We combine these 4 model components to explain the emission of 17 spectral lines, investigate the global energy balance of the galaxy through its spectral energy distribution, and establish a phase mass inventory. While the ionic emission lines of Haro 11 essentially originate from the dense HII region component, a diffuse low-ionisation gas is needed to explain the [NeII], [NII], and [CII] line intensities. The [OIII] 88um line intensity is not fully reproduced by our model, hinting towards the possible presence of yet another low-density high-ionisation medium. The [OI] emission is consistent with a dense PDR of low covering factor, and we find no evidence for an X-ray dominated component. The PDR component accounts for only 10% of the [CII] emission.

We have conducted a high-resolution spectroscopic study using Spitzer of 18 bulgeless (Sd/Sdm) galaxies that show no definitive signatures of nuclear activity in their optical spectra. This is the first systematic mid-IR search for weak or hidden AGNs in a statistically significant sample of bulgeless disk galaxies. Based on the detection of the high-ionization [NeV] line, we report the discovery of an AGN in one out of the 18 galaxies in the sample. This galaxy, NGC 4178, is a nearby edge-on Sd galaxy, which likely hosts a prominent nuclear star cluster (NSC). The bolometric luminosity of the AGN inferred from the [NeV] luminosity is ~ 8e41 ergs/s. This is almost two orders of magnitude greater than the luminosity of the AGN in NGC 4395, the best studied AGN in a bulgeless disk galaxy. Assuming that the AGN in NGC 4178 is radiating below the Eddington limit, the lower mass limit for the black hole is ~ 6e3M_sun. The fact that none of the other galaxies in the sample shows any evidence for an AGN demonstrates that while the AGN detection rate based on mid-IR diagnostics is high (30-40%) in optically quiescent galaxies with pseudobulges, it drops drastically in Sd/Sdm galaxies. Our observations therefore confirm that AGNs in completely bulgeless disk galaxies are not hidden in the optical but truly are rare. Of the three Sd galaxies with AGNs known so far, all have prominent NSCs, suggesting that in the absence of a well-defined bulge, the galaxy must possess a NSC in order to host an AGN. While the presence of a NSC appears to be a requirement for hosting an AGN in bulgeless galaxies, neither the properties of the NSC nor those of the host galaxy appear exceptional in late-type AGN hosts. The recipe for forming and growing a central black hole in a bulgeless galaxy therefore remains unknown.

The Orion Bar is one of the nearest and best-studied photodissociation or photon-dominated regions (PDRs). Observations reveal the presence of H2 lines from vibrationally or rotationally excited upper levels that suggest warm gas temperatures (400 to 700 K). However, standard models of PDRs are unable to reproduce such warm rotational temperatures. In this paper we attempt to explain these observations with new comprehensive models which extend from the H+ region through the Bar and include the magnetic field in the equation of state. We adopt the model parameters from our previous paper which successfully reproduced a wide variety of spectral observations across the Bar. In this model the local cosmic-ray density is enhanced above the galactic background, as is the magnetic field, and which increases the cosmic-ray heating elevating the temperature in the molecular region. The pressure is further enhanced above the gas pressure in the H+ region by the momentum transferred from the absorbed starlight. Here we investigate whether the observed H2 lines can be reproduced with standard assumptions concerning the grain photoelectric emission. We also explore the effects due to the inclusion of recently computed H2 + H2, H2 + H and H2 + He collisional rate coefficients.

The observed faintness of infrared fine-structure line emission along with the warm far-infrared (FIR) colors of ultraluminous infrared galaxies (ULIRGs) is a long-standing problem. In this work, we calculate the line and continuum properties of a cloud exposed to an Active Galactic Nucleus (AGN) and starburst spectral energy distribution (SED). We use an integrated modeling approach, predicting the spectrum of ionized, atomic, and molecular environments in pressure equilibrium. We find that the effects of high ratios of impinging ionizing radiation density to particle density (i.e. high ionization parameters, or U) can reproduce many ULIRG observational characteristics. Physically, as U increases, the fraction of UV photons absorbed by dust increases, corresponding to fewer photons available to photoionize and heat the gas, producing what is known as a "dust-bounded" nebula. We show that high U effects can explain the "[C II] deficit", the ~1 dex drop in the [C II] 158 micron /FIR ratio seen in ULIRGs when compared to starburst or normal galaxies. Additionally, by increasing U through increasing the ionizing photon flux, warmer dust and thus higher IRAS F(60)/F(100) ratios result. High U effects also predict an increase in [O I]63 micron /[C II] 158 micron and a gradual decline in [O III] 88 micron /FIR, similar to the magnitude of the trends observed, and yield a reasonable fit to [Ne V]14 micron /FIR ratio AGN observations.

The three dimensional structure of the brightest part of the Orion Nebula is assessed in the light of published and new data. We find that the widely accepted model of a concave blister of ionized material needs to be altered in the southwest direction from the Trapezium, where we find that the Orion-S feature is a separate cloud of very optically thick molecules within the body of ionized gas, which is probably the location of the multiple embedded sources that produce the outflows that define the Orion-S star formation region. Evidence for this cloud comes from the presence of H2CO lines in absorption in the radio continuum and discrepancies in the extinction derived from radio-optical and optical only emission. We present an equilibrium Cloudy model of the Orion-S cloud, which successfully reproduces many observed properties of this feature. We also report the discovery of an open-sided shell of [O III] surrounding the Trapezium stars, revealed through emission line ratio images and the onset of radiation shadows beyond some proplyds. We show that the observed properties of the shell are consistent with it being a stationary structure, produced by shock interactions between the ambient nebular gas and the high-velocity wind from theta^1 Ori C. We examine the implications of the recently published evidence for a large blueshifted velocity of theta^1 Ori C with respect to the Orion Molecular Cloud, which could mean that this star has only recently begun to photoionize the Orion Nebula. We show that current observations of the Nebula do not rule out such a possibility, so long as the ionization front has propagated into a pre-existing low-density region. In addition, a young age for the Nebula would help explain the presence of nearby proplyds with a short mass-loss timescale to photoablation.

Our understanding of physical processes in Photodissociation regions or Photon Dominated Regions (PDRs) largely depends on the ability of spectral synthesis codes to reproduce the observed infrared emission-line spectrum. In this paper, we explore the sensitivity of a single PDR model to microphysical details. Our calculations use the Cloudy spectral synthesis code, recently modified to include a wealth of PDR physical processes. We show how the chemical/thermal structure of a PDR, along with the calculated spectrum, changes when the treatment of physical processes such as grain physics and atomic/molecular rates are varied. We find a significant variation in the intensities of PDR emission lines, depending on different treatments of the grain physics. We also show how different combinations of the cosmic-ray ionization rate, inclusion of grain-atom/ion charge transfer, and the grain size distribution can lead to very similar results for the chemical structure. Additionally, our results show the utility of Cloudy for the spectral modeling of molecular environments.

(Abridged) We carried out a comprehensive far-ultraviolet (UV) survey of ^12CO and H_2 column densities along diffuse molecular Galactic sight lines in order to explore in detail the relationship between CO and H_2. We measured new CO abundances from HST spectra, new H_2 abundances from FUSE data, and new CH, CH^+, and CN abundances from the McDonald and European Southern Observatories. A plot of log N(CO) versus log N(H_2) shows that two power-law relationships are needed for a good fit of the entire sample, with a break located at log N(CO, cm^-2) = 14.1 and log N(H_2) = 20.4, corresponding to a change in production route for CO in higher-density gas. Similar logarithmic plots among all five diatomic molecules allow us to probe their relationships, revealing additional examples of dual slopes in the cases of CO versus CH (break at log N = 14.1, 13.0), CH^+ versus H_2 (13.1, 20.3), and CH^+ versus CO (13.2, 14.1). These breaks are all in excellent agreement with each other, confirming the break in the CO versus H_2 relationship, as well as the one-to-one correspondence between CH and H_2 abundances. Our new sight lines were selected according to detectable amounts of CO in their spectra and they provide information on both lower-density (< 100 cm^-3) and higher-density diffuse clouds. The CO versus H_2 correlation and its intrinsic width are shown to be empirically related to the changing total gas density among the sight lines of the sample. We employ both analytical and numerical chemical schemes in order to derive details of the molecular environments. In the low-density gas, where equilibrium-chemistry studies have failed to reproduce the abundance of CH^+, our numerical analysis shows that nonequilibrium chemistry must be employed for correctly predicting the abundances of both CH^+ and CO.

Archival Spitzer observations of 41 starburst galaxies that span a wide range in metallicity reveal for the first time a correlation between the [FeII]/[NeII] 26.0/12.8 micron ratio and the electron gas density as traced by the 18.7/33.4 micron [SIII] ratio, with the [FeII]/[NeII] ratio decreasing with increasing gas density. We also find a strong correlation between the gas density and the PAH peak to continuum strength. Using shock and photoionization models, we see the driver of the observed [FeII]/[NeII] ratios is metallicity. The majority of [FeII] emission in low metallicity galaxies may be shock-derived, whilst at high metallicity, the [FeII] emission may be instead dominated by contributions from HII and in particular from dense PDR regions. However, the observed [FeII]/[NeII] ratios may instead be following a metallicity-abundance relationship, with iron being less depleted onto grains in low metallicity galaxies - a result that would have profound implications for the use of iron emission lines as unambiguous tracers of shocks.

Archival Spitzer observations of 41 starburst galaxies that span a wide range in metallicity reveal for the first time a correlation between the [FeII]/[NeII] 26.0/12.8 micron ratio and the electron gas density as traced by the 18.7/33.4 micron [SIII] ratio, with the [FeII]/[NeII] ratio decreasing with increasing gas density. The correlations of the [FeII]/[NeII] ratio, the PAH peak to continuum strength & metallicity found in an earlier paper were confirmed for a larger sample. We also find a strong correlation between the gas density and the PAH peak to continuum strength. Using shock and photoionization models, we see the driver of the observed [FeII]/[NeII] ratios is metallicity. The majority of [FeII] emission in low metallicity galaxies may be shock-derived, whilst at high metallicity, the [FeII] emission may be instead dominated by contributions from HII and in particular from dense PDR regions. However, the observed [FeII]/[NeII] ratios may instead be following a metallicity-abundance relationship, with iron being less depleted onto grains in low metallicity galaxies - a result that would have profound implications for the use of iron emission lines as unambiguous tracers of shocks.

Recent laboratory measurements for the S^+2 + H2 reaction find a total rate coefficient significantly larger than previously used in theoretical models of X-ray dominated regions (XDRs). While the branching ratio of the products is unknown, one energetically possible route leads to the SH+ molecule, a known XDR diagnostic. In this work, we study the effects of S^+2 on the formation of SH+ and the destruction of S^+2 in XDRs. We find the predicted SH+ column density for molecular gas surrounding an Active Galactic Nucleus (AGN) increases by as much as 2 dex. As long as the branching ratio for S^+2 + H2 -> SH+ + H+ exceeds a few percent, doubly ionized chemistry will be the dominant pathway to SH+, which then initiates the formation of other sulfur-bearing molecules. We also find that the high rate of S^+2 + H2 efficiently destroys S^+2 once H2 forms, while the S^+2 abundance remains high in the atomic hydrogen region. We discuss the possible consequences of S^+2 in the atomic hydrogen region on mid-infrared diagnostics. The enhanced SH+ abundance has important implications in the study of XDRs, while our conclusions for S^+2 could potentially impact the interpretation of Spitzer and SOFIA observations.

We report the discovery using Spitzers high resolution spectrograph of 7 Active Galactic Nuclei (AGN) in a sample of 32 late-type galaxies that show no definitive signatures of AGN in their optical spectra. Our observations suggest that the AGN detection rate in late-type galaxies is possibly 4 times larger than what optical spectroscopic observations alone suggest. We demonstrate using photoionization models with an input AGN and an extreme EUV-bright starburst ionizing radiation field that the observed mid-infrared line ratios cannot be replicated unless an AGN contribution, in some cases as little as 10% of the total galaxy luminosity, is included. These models show that when the fraction of the total luminosity due to the AGN is low, optical diagnostics are insensitive to the presence of the AGN. In this regime of parameter space, the mid-infrared diagnostics offer a powerful tool for uncovering AGN missed by optical spectroscopy. The AGN bolometric luminosities in our sample range from ~3 X 10^41 - ~2 X 10^43 ergs s^-1, which, based on the Eddington limit, corresponds to a lower mass limit for the black hole that ranges from ~3 X 10^3Mdot to as high as ~1.5 X 10^5Mdot. These lower mass limits however do not put a strain on the well-known relationship between the black hole mass and the host galaxy's stellar velocity dispersion established in predominantly early-type galaxies. Our findings add to the growing evidence that black holes do form and grow in low-bulge environments and that they are significantly more common than optical studies indicate.

Detecting Active Galactic Nuclei (AGN) in galaxies dominated by powerful nuclear star formation and extinction effects poses a unique challenge. Due to the longer wavelength emission and the ionization potential of Ne^3+, infrared [Ne V] emission lines are thought to be excellent AGN diagnostics. However, stellar evolution models predict Wolf-Rayet stars in young stellar clusters emit significant numbers of photons capable of creating Ne^4+. Recent observations of [Ne V] emission in optically classified starburst galaxies require us to investigate whether [Ne V] can arise from star-formation activity and not an AGN. In this work, we calculate the optical and IR spectrum of gas exposed to a young starburst and AGN SED. We find: 1) a range of parameters where [Ne V] emission can be explained solely by star-formation and 2) a range of relative AGN to starburst luminosities that reproduce the [Ne V] observations, yet leaves the optical spectrum looking like a starburst. We also find infrared emission-line diagnostics are much more sensitive to the AGN than optical diagnostics, particularly for weak AGN. We apply our model to the optically classified, yet [Ne V] emitting, starburst galaxy NGC 3621. We find, when taking the IR and optical spectrum into account, about 30 - 50% of the galaxy's total luminosity is due to an AGN. Our calculations show that [Ne V] emission is almost always the result of AGN activity. The models presented in this work can be used to determine the AGN contribution to a galaxy's power output.

The spatial distribution of the cosmic-ray flux is important in understanding the Interstellar Medium (ISM) of the Galaxy. This distribution can be analyzed by studying different molecular species along different sight lines whose abundances are sensitive to the cosmic-ray ionization rate. Recently several groups have reported an enhanced cosmic-ray ionization rate in diffuse clouds compared to the standard value, zeta(average)=2.5e-17 s^-1, measured toward dense molecular clouds. In an earlier work we reported an enhancement in the cosmic ray rate of 20 towards HD185418. McCall et al. have reported enhancements of 48 towards zeta Persei based on the observed abundance of H$_{3}^+ while Le Petit et al. found a cosmic ray enhancement of ~10 to be consistent with their models for this same sight line. Here we revisit zeta Persei and perform a detailed calculation using a self-consistent treatment of the hydrogen chemistry, grain physics, energy and ionization balance, and excitation physics. We show that the value of zeta deduced from the H$_{3}^+$ column density in the diffuse region of the sightline depends strongly on the properties of the grains because they remove free electrons and change the hydrogen chemistry. The observations are largely consistent with a cosmic ray enhancement of 40, with several diagnostics indicating higher values. This underscores the importance of a full treatment of grain physics in studies of interstellar chemistry.

This work focuses on [O I] and [C II] emission towards NGC 6334 A, an embedded H+ region/PDR only observable at infrared or longer wavelengths. A geometry where nearly all the emission escapes out the side of the cloud facing the stars, such as Orion, is not applicable to this region. Instead, we find the geometry to be one where the H+ region and associated PDR is embedded in the molecular cloud. Constant density PDR calculations are presented which predict line intensities as a function of AV (or N(H)), hydrogen density (nH), and incident UV radiation field (G0). We find that a single component model with AV ~650 mag, nH = 5x10^5 cm-3, and G0 = 7x10^4 reproduces the observed [O I] and [C II] intensities, and that the low [O I] 63 to 146 micron ratio is due to line optical depth effects in the [O I] lines, produced by a large column density of atomic/molecular gas. We find that the effects of a density-law would increase our derived AV, while the effects of an asymmetric geometry would decrease AV, with the two effects largely canceling. We conclude that optically selected H+ regions adjacent to PDRs, such as Orion, likely have a different viewing angle or geometry than similar regions detected through IR observations. Overall, the theoretical calculations presented in this work have utility for any PDR embedded in a molecular cloud.

We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.

The [C II] 158 micron line is an important emission line diagnostic in Photodissociation Regions (PDRs), but this emission line can also emerge from ionized gas. This work calculates the contribution of [C II] emission from ionized gas over a wide range of parameter space by considering the simplified case of an H+ region and PDR in pressure equilibrium. Additionally, these calculations also predict the strong correlation observed between [N II] 205 micron emission and [C II] discussed by previous authors. Overall, the results of these calculations have wide-ranging applications to the interpretation of [C II] emission in astrophysical environments.

The emission line regions of starburst galaxies and active nuclei reveal a wealth of spectroscopic information. A unified picture of the relationship between ionized, atomic, and molecular gas makes it possible to better understand these observations. We performed a series of calculations designed to determine the equation of state, the relationship between density, temperature, and pressure, through emission-line diagnostic ratios that form in the H+ region and PDR. We consider a wide range of physical conditions in the H+ region. We connect the H+ region to the PDR by considering two constant pressure cases, one with no magnetic field and one where the magnetic field overwhelms the thermal pressure. We show that diagnostic ratios can yield the equation of state for single H+ regions adjacent to single PDRs, with the results being more ambiguous when considering observations of entire galaxies. As a test, we apply our calculations to the Orion H+/PDR region behind the Trapezium. We find the ratio of thermal to magnetic pressure in the PDR to be \~1.2. If magnetic and turbulent energy are in equipartition, our results mean the magnetic field is not the cause of the unexplained broadening in M 42, but may significantly affect line broadening in the PDR. Since Orion is often used to understand physical processes in extragalactic environments, our calculations suggest magnetic pressure should be considered in modeling such regions.

Orion's Veil is an absorbing screen that lies along the line of sight to the Orion H II region. It consists of two or more layers of gas that must lie within a few parsecs of the Trapezium cluster. Our previous work considered the Veil as a whole and found that the magnetic field dominates the energetics of the gas in at least one component. Here we use high-resolution STIS UV spectra that resolve the two velocity components in absorption and determine the conditions in each. We derive a volume hydrogen density, 21 cm spin temperature, turbulent velocity, and kinetic temperature, for each. We combine these estimates with magnetic field measurements to find that magnetic energy significantly dominates turbulent and thermal energies in one component, while the other component is close to equipartition between turbulent and magnetic energies. We observe molecular hydrogen absorption for highly excited v, J levels that are photoexcited by the stellar continuum, and detect blueshifted S III and P III. These ions must arise from ionized gas between the mostly neutral portions of the Veil and the Trapezium and shields the Veil from ionizing radiation. We find that this layer of ionized gas is also responsible for He I absorption in the Veil, which resolves a 40-year-old debate on the origin of He I absorption towards the Trapezium. Finally, we determine that the ionized and mostly atomic layers of the Veil will collide in less than 85,000 years.

We have developed a complete model of the hydrogen molecule as part of the spectral simulation code Cloudy. Our goal is to apply this to spectra of high-redshift star-forming regions where H2 absorption is seen, but where few other details are known, to understand its implication for star formation. The microphysics of H2 is intricate, and it is important to validate these numerical simulations in better-understood environments. This paper studies a well-defined line-of-sight through the Galactic interstellar medium (ISM) as a test of the microphysics and methods we use. We present a self-consistent calculation of the observed absorption-line spectrum to derive the physical conditions in the ISM towards HD185418, a line-of-sight with many observables. We deduce density, temperature, local radiation field, cosmic ray ionization rate, chemical composition and compare these conclusions with conditions deduced from analytical calculations. We find a higher density, similar abundances, and require a cosmic ray flux enhanced over the Galactic background value, consistent with enhancements predicted by MHD simulations.

We have performed a series of calculations designed to reproduce infrared diagnostics used to determine physical conditions in star forming regions. We self-consistently calculate the thermal and chemical structure of an H II region and photodissociation region (PDR) that are in pressure equilibrium. This differs from previous work, which used separate calculations for each gas phase. Our calculations span a wide range of stellar temperatures, gas densities, and ionization parameters. We describe improvements made to the spectral synthesis code Cloudy that made these calculations possible. These include the addition of a molecular network with ~1000 reactions involving 68 molecular species and improved treatment of the grain physics. Data from the Spitzer First Look Survey, along with other archives, are used to derive important physical characteristics of the H II region and PDR. These include stellar temperatures, electron densities, ionization parameters, UV radiation flux, and PDR density. Finally, we calculate the contribution of the H II region to PDR emission line diagnostics, which allows for a more accurate determination of physical conditions in the PDR.

Conditions within the Homunculus nebula around Eta Car are determined by many of the same physical processes that occur in molecular clouds in the interstellar medium. But there is one major exception -- we know when the ejection occurred and something about its composition and initial state. The gas was warm, ionized, and dust-free when it was located within the star's atmosphere and it is currently cold, molecular, and dusty. It undertook this transformation in a bit over 150 years. It offers a laboratory for the study of physical processes in a well-constrained environment. We derive a photoionization model of the Homunculus nebula that reproduces many of its observed properties. We conclude by outlining how observations of the Homunculus could address basic problems in the physics of the interstellar medium.

Much of the baryonic matter in the Universe is in the form of H2 which includes most of the gas in galactic and extragalactic interstellar clouds. Molecular hydrogen plays a significant role in establishing the thermal balance in many astrophysical environments and can be important as a spectral diagnostic of the gas. Modeling and interpretation of observations of such environments requires a quantitatively complete and accurate treatment of H2. Using this micro-physical model of H2, illustrative calculations of prototypical astrophysical environments are presented. This work forms the foundation for future investigations of these and other environments where H2 is an important constituent.

Orion's veil consists of several layers of largely neutral gas lying between us and the main ionizing stars of the Orion nebula. It is visible in 21cm H I absorption and in optical and UV absorption lines of H I and other species. Toward the Trapezium, the veil has two remarkable properties, high magnetic field (~100 microGauss) and a surprising lack of molecular hydrogen given its total hydrogen column density. Here we compute photoionization models of the veil to establish its gas density and its distance from the Trapezium. We use a greatly improved model of the hydrogen molecule that determines level populations in ~1e5 rotational/vibrational levels and provides improved estimates of molecular hydrogen destruction via the Lyman-Werner bands. Our best fit photoionization models place the veil 1-3 pc in front of the star at a density of 1e3-1e4 cubic centimeters. Magnetic energy dominates the energy of non-thermal motions in at least one of the 21cm H I velocity components. Therefore, the veil is the first interstellar environment where magnetic dominance appears to exist. We find that the low ratio of molecular to atomic hydrogen (< 1e-4) is a consequence of high UV flux incident upon the veil due to its proximity to the Trapezium stars and the absence of small grains in the region.

It has long been known that infrared fine structure lines of abundant ions, like the [O III] 88 micron line, can become optically thick in H II regions under certain high luminosity conditions. This could mitigate their potential as diagnostic tools, especially if the source is too dusty for optical spectroscopy to otherwise determine the system's parameters. We examined a series of photoionization calculations which were designed to push the nebulae into the limit where many IR lines should be quite optically thick. We find that radiative transfer effects do not significantly change the observed emission line spectrum. This is due to a combination of grain absorption of the hydrogen ionizing continuum and the fact that the correction for stimulated emission in these lines is large. Given these results, and the likelihood that real objects have non-thermal line broadening, it seems unlikely that line optical depth presents a problem in using these lines as diagnostics of the physical conditions or chemical composition.