We present observations of velocity-resolved [CII] 158 $\mu$m emission from both a dense and a more diffuse photodissociation region (PDR) in the Perseus giant molecular cloud (GMC) using the Heterodyne Instrument for the Far-Infrared onboard the Herschel Space Telescope. We detect [CII] emission from 80\% of the total positions, with a 95\% detection rate from the dense boundary region. The integrated intensity of the [CII] emission remains relatively constant across each boundary, despite the observed range in optical extinction between 1 mag and 10 mag. This flat profile indicates a constant heating and cooling rate within both regions observed. The integrated intensity of [CII] emission is reasonably well correlated with the neutral hydrogen (HI) column density, as well as total gas column density. This, in addition to the 80$'$ (7 pc) extent of the [CII] emission from cloud center, suggests that the HI envelope plays a dominant role in explaining the [CII] emission emanating from Perseus. We compare the [CII] and $^{12}$CO integrated intensities with predictions from a 1-D, two-sided slab PDR model and show that a simple core $+$ envelope, equilibrium model without an additional "CO-dark" H$_2$ component can reproduce observations well. Additional observations are needed to disentangle how much of the [CII] emission is associated with the "CO-dark" H$_2$ gas, as well as constrain spatial variations of the dust-to-gas ratio across Perseus.

Strong winds and ultraviolet (UV) radiation from O-type stars disrupt and ionize their molecular core birthplaces, sweeping up material into parsec-size shells. Owing to dissociation by starlight, the thinnest shells are expected to host low molecular abundances and therefore little star formation. Here, we expand previous maps taken with the IRAM 30m telescope and present square-degree 12CO and 13CO (J=2-1) maps of the wind-driven "Veil bubble'' that surrounds the Trapezium cluster and its natal Orion molecular core (OMC). Although widespread and extended CO emission is largely absent from the Veil, we show that several CO "globules'' exist and are embedded in the [CII]158um-bright shell that confines the bubble. This includes the first detection of quiescent CO at negative LSR velocities in Orion. Given the harsh UV irradiation conditions in this translucent material, the detection of CO globules is surprising. These globules are small (R=7,100 AU), not massive (M=0.3M_Sun), and are moderately dense: n_ H=4x10^4 cm^-3 (median values). They are confined by the external pressure of the shell, P_ext/k~10^7 cm^-3 K, and are likely magnetically supported. They are either transient objects formed by instabilities or have detached from pre-existing molecular structures, sculpted by the passing shock associated with the expanding shell and by UV radiation from the Trapezium. Some represent the first stages in the formation of small pillars, others of isolated small globules. Although their masses do not suggest they will form stars, one globule matches the position of a known YSO. The lack of extended CO in the "Veil shell'' demonstrates that feedback from massive stars expels, agitates, and reprocesses most of the disrupted molecular cloud gas, thereby limiting the star-formation rate in the region. The presence of globules is a result of this feedback.

We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4-3) and [CI](2-1) lines; we hypothesize this is due to the similar upper level temperature of these two lines. We modeled the [CI] and CO line emission using large velocity gradient models combined with an empirical template. According to this modeling, the [CI](1-0) line is clearly dominated by the low-excitation component. We determine [CI] to molecular mass conversion factors for both the [CI](1-0) and [CI](2-1) lines, with mean values of $\alpha_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $\alpha_{\mathrm{[CI](2-1)}} = 34 $ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ with logarithmic root-mean-square spreads of 0.20 and 0.32 dex, respectively. The similar spread of $\alpha_{\mathrm{[CI](1-0)}}$ to $\alpha_{\mathrm{CO}}$ (derived using the CO(2-1) line) suggests that [CI](1-0) may be just as good a tracer of cold molecular gas as CO(2-1) in galaxies of this type. On the other hand, the wider spread of $\alpha_{\mathrm{[CI](2-1)}}$ and the tight relation found between [CI](2-1) and CO(4-3) suggest that much of the [CI](2-1) emission may originate in warmer molecular gas.

Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies.

The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that originates from the ionized and neutral phases of the ISM. Sub-kiloparsec measurements of the [CII] 158 micron and [NII] 205 micron line in nearby galaxies have recently become available as part of the Key Insights in Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) programs. With the information from these two far-infrared lines along with the multi-wavelength suite of KINGFISH data, a calibration of the [CII] emission line as a star formation rate indicator and a better understanding of the [CII] deficit are pursued. [CII] emission is also compared to PAH emission in these regions to compare photoelectric heating from PAH molecules to cooling by [CII] in the neutral and ionized phases of the ISM. We find that the [CII] emission originating in the neutral phase of the ISM does not exhibit a deficit with respect to the infrared luminosity and is therefore preferred over the [CII] emission originating in the ionized phase of the ISM as a star formation rate indicator for the normal star-forming galaxies included in this sample.

Supersonic turbulence results in strong density fluctuations in the interstellar medium (ISM), which have a profound effect on the chemical structure. Particularly useful probes of the diffuse ISM are the ArH$^+$, OH$^+$, H$_2$O$^+$ molecular ions, which are highly sensitive to fluctuations in the density and the H$_2$ abundance. We use isothermal magnetohydrodynamic (MHD) simulations of various sonic Mach numbers, $\mathcal{M}_s$, and density decorrelation scales, $y_{\rm dec}$, to model the turbulent density field. We post-process the simulations with chemical models and obtain the probability density functions (PDFs) for the H$_2$, ArH$^+$, OH$^+$ and H$_2$O$^+$ abundances. We find that the PDF dispersions increases with increasing $\mathcal{M}_s$ and $y_{\rm dec}$, as the magnitude of the density fluctuations increases, and as they become more coherent. Turbulence also affects the median abundances: when $\mathcal{M}_s$ and $y_{\rm dec}$ are high, low-density regions with low H$_2$ abundance become prevalent, resulting in an enhancement of ArH$^+$ compared to OH$^+$ and H$_2$O$^+$. We compare our models with Herschel observations. The large scatter in the observed abundances, as well as the high observed ArH$^+$/OH$^+$ and ArH$^+$/H$_2$O$^+$ ratios are naturally reproduced by our supersonic $(\mathcal{M}_s=4.5)$, large decorrelation scale $(y_{\rm dec}=0.8)$ model, supporting a scenario of a large-scale turbulence driving. The abundances also depend on the UV intensity, CR ionization rate, and the cloud column density, and the observed scatter may be influenced by fluctuations in these parameters.

We report the results from a new, highly sensitive ($\Delta T_{mb} \sim 3 $mK) survey for thermal OH emission at 1665 and 1667 MHz over a dense, 9 x 9-pixel grid covering a $1\deg$ x $1\deg$ patch of sky in the direction of $l = 105\deg, b = +2.50\deg$ towards the Perseus spiral arm of our Galaxy. We compare our Green Bank Telescope (GBT) 1667 MHz OH results with archival CO J=1-0 observations from the Five College Radio Astronomy Observatory (FCRAO) Outer Galaxy Survey within the velocity range of the Perseus Arm at these galactic coordinates. Out of the 81 statistically-independent pointings in our survey area, 86% show detectable OH emission at 1667 MHz, and 19% of them show detectable CO emission. We explore the possible physical conditions of the observed features using a set of diffuse molecular cloud models. In the context of these models, both OH and CO disappear at current sensitivity limits below an A$_{\rm v}$ of 0.2, but the CO emission does not appear until the volume density exceeds 100-200 cm$^{-3}$. These results demonstrate that a combination of low column density A$_{\rm v}$ and low volume density $n_{H}$ can explain the lack of CO emission along sight lines exhibiting OH emission. The 18-cm OH main lines, with their low critical density of $n^{*}$ $ \sim 1 $ cm$^{-3}$, are collisionally excited over a large fraction of the quiescent galactic environment and, for observations of sufficient sensitivity, provide an optically-thin radio tracer for diffuse H$_2$.

We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $\mu$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a velocity resolution of 0.17 km s$^{-1}$. The observations cover a 0.25$^\circ$ by 0.28$^\circ$ area with central position {\it l} = 297.34$^\circ$, {\it b} = -0.60$^\circ$. The kinematics show that bright [C\,{\sc ii}] structures are spatially and spectrally correlated with the surfaces of CO clouds, tracing the photodissociation region and ionization front of each molecular cloud. Along 7 lines of sight that traverse Tr 14 into the dark ridge to the southwest, we find that the [C\,{\sc ii}] luminosity from the HII region is 3.7 times that from the PDR. In same los we find in the PDRs an average ratio of 1:4.1:5.6 for the mass in atomic gas:dark-CO gas: molecular gas traced by CO. Comparing multiple gas tracers including HI 21cm, [C\,{\sc ii}], CO, and radio recombination lines, we find that the HII regions of the Carina Nebula Complex are well-described as HII regions with one-side freely expanding towards us, consistent with the champagne model of ionized gas evolution. The dispersal of the GMC in this region is dominated by EUV photoevaporation; the dispersal timescale is 20-30 Myr.

Chemistry plays an important role in the interstellar medium (ISM), regulating heating and cooling of the gas, and determining abundances of molecular species that trace gas properties in observations. Although solving the time-dependent equations is necessary for accurate abundances and temperature in the dynamic ISM, a full chemical network is too computationally expensive to incorporate in numerical simulations. In this paper, we propose a new simplified chemical network for hydrogen and carbon chemistry in the atomic and molecular ISM. We compare results from our chemical network in detail with results from a full photo-dissociation region (PDR) code, and also with the Nelson & Langer (1999) (NL99) network previously adopted in the simulation literature. We show that our chemical network gives similar results to the PDR code in the equilibrium abundances of all species over a wide range of densities, temperature, and metallicities, whereas the NL99 network shows significant disagreement. Applying our network in 1D models, we find that the $\mathrm{CO}$-dominated regime delimits the coldest gas and that the corresponding temperature tracks the cosmic ray ionization rate in molecular clouds. We provide a simple fit for the locus of $\mathrm{CO}$ dominated regions as a function of gas density and column. We also compare with observations of diffuse and translucent clouds. We find that the $\mathrm{CO}$, $\mathrm{CHx}$ and $\mathrm{OHx}$ abundances are consistent with equilibrium predictions for densities $n=100-1000~\mathrm{cm^{-3}}$, but the predicted equilibrium $\mathrm{C}$ abundance is higher than observations, signaling the potential importance of non-equilibrium/dynamical effects.

Massive stars inject mechanical and radiative energy into the surrounding environment, which stirs it up, heats the gas, produces cloud and intercloud phases in the interstellar medium, and disrupts molecular clouds (the birth sites of new stars). Stellar winds, supernova explosions and ionization by ultraviolet photons control the lifetimes of molecular clouds. Theoretical studies predict that momentum injection by radiation should dominate that by stellar winds, but this has been difficult to assess observationally. Velocity-resolved large-scale images in the fine-structure line of ionized carbon ([C II]) provide an observational diagnostic for the radiative energy input and the dynamics of the interstellar medium around massive stars. Here we report observations of a one-square-degree region (about 7 parsecs in diameter) of Orion molecular core -- the region nearest to Earth that exhibits massive-star formation -- at a resolution of 16 arcseconds (0.03 parsecs) in the [C II] line at 1.9 terahertz (158 micrometres). The results reveal that the stellar wind originating from the massive star ${\theta}^{1}$ Orionis C has swept up the surrounding material to create a bubble roughly four parsecs in diameter with a 2,600-solar-mass shell, which is expanding at 13 kilometres per second. This finding demonstrates that the mechanical energy from the stellar wind is converted very efficiently into kinetic energy of the shell and causes more disruption of the Orion molecular core 1 than do photo-ionization and evaporation or future supernova explosions.

Cosmic rays are predominantly accelerated in shocks associated with star formation such as supernova remnants and stellar wind bubbles, so the cosmic-ray flux and thus cosmic-ray ionization rate, $\zeta_{\rm H}$, should correlate with the star-formation rate in a galaxy. Submillimeter bright galaxies (SMGs) are some of the most prolific star forming galaxies in the Universe, and gravitationally lensed SMGs provide bright continuum sources suitable for absorption line studies. Abundances of OH$^+$ and H$_2$O$^+$ are useful for inferring $\zeta_{\rm H}$ when combined with chemical models, and have been used for this purpose within the Milky Way. At redshifts $z\gtrsim2$ transitions out of the ground rotational states of OH$^+$ and H$_2$O$^+$ are observable with ALMA, and we present observations of both molecules in absorption toward the lensed SMGs SMM J2135$-$0102 and SDP 17b. These detections enable an exploration of $\zeta_{\rm H}$ in galaxies with extreme star formation and high supernova rates, both of which should significantly enhance cosmic-ray production. The observed OH$^+$ and H$_2$O$^+$ absorption is thought to arise in massive, extended halos of cool, diffuse gas that surround these galaxies. Using a chemical model designed to focus on the reaction network important to both species, we infer cosmic-ray ionization rates of $\zeta_{\rm H}\sim10^{-16}$-$10^{-14}$ s$^{-1}$ in these extended gaseous halos. As our estimates come from gas that is far away from the sites of cosmic-ray acceleration, they imply that cosmic-ray ionization rates in the compact regions where star formation occurs in these galaxies are orders of magnitude higher.

Quantitative understanding of the interstellar medium requires knowledge of its physical conditions. Low frequency carbon radio recombination lines (CRRLs) trace cold interstellar gas, and can be used to determine its physical conditions (e.g., electron temperature and density). In this work we present spatially resolved observations of the low frequency ($\leq390$ MHz) CRRLs centered around C$268\alpha$, C$357\alpha$, C$494\alpha$ and C$539\alpha$ towards Cassiopeia A on scales of $\leq1.2$ pc. We compare the spatial distribution of CRRLs with other ISM tracers. This comparison reveals a spatial offset between the peak of the CRRLs and other tracers, which is very characteristic for photodissociation regions and that we take as evidence for CRRLs being preferentially detected from the surfaces of molecular clouds. Using the CRRLs we constrain the gas electron temperature and density. These constraints on the gas conditions suggest variations of less than a factor of two in pressure over $\sim1$ pc scales, and an average hydrogen density of $200$-$470$ cm$^{-3}$. From the electron temperature and density maps we also constrain the ionized carbon emission measure, column density and path length. Based on these, the hydrogen column density is larger than $10^{22}$ cm$^{-2}$, with a peak of $\sim4\times10^{22}$ cm$^{-2}$ towards the South of Cassiopeia A. Towards the southern peak the line of sight length is $\sim40$ pc over a $\sim2$ pc wide structure, which implies that the gas is a thin surface layer on a large (molecular) cloud that is only partially intersected by Cassiopeia A. These observations highlight the utility of CRRLs as tracers of low density extended HI and CO-dark gas halo's around molecular clouds.

The Small Magellanic Cloud (SMC) provides the only laboratory to study the structure of molecular gas at high resolution and low metallicity. We present results from the Herschel Spectroscopic Survey of the SMC (HS$^{3}$), which mapped the key far-IR cooling lines [CII], [OI], [NII], and [OIII] in five star-forming regions, and new ALMA 7m-array maps of $^{12}$CO and $^{13}$CO $(2-1)$ with coverage overlapping four of the five HS$^{3}$ regions. We detect [CII] and [OI] throughout all of the regions mapped. The data allow us to compare the structure of the molecular clouds and surrounding photodissociation regions using $^{13}$CO, CO, [CII], and [OI] emission at $<10$" ($<3$ pc) scales. We estimate Av using far-IR thermal continuum emission from dust and find the CO/[CII] ratios reach the Milky Way value at high A$_{V}$ in the centers of the clouds and fall to $\sim{1/5-1/10}\times$ the Milky Way value in the outskirts, indicating the presence of translucent molecular gas not traced by bright CO emission. We estimate the amount of molecular gas traced by bright [CII] emission at low A$_{V}$ and bright CO emission at high A$_{V}$. We find that most of the molecular gas is at low A$_{V}$ and traced by bright [CII] emission, but that faint CO emission appears to extend to where we estimate the H$_{2}$-to-HI transition occurs. By converting our H$_{2}$ gas estimates to a CO-to-H$_{2}$ conversion factor ($X_{CO}$), we show that $X_{CO}$ is primarily a function of A$_{V}$, consistent with simulations and models of low metallicity molecular clouds.

Observations towards L1630 in the Orion B molecular cloud, comprising the iconic Horsehead Nebula, allow us to study the interplay between stellar radiation and a molecular cloud under relatively benign conditions, that is, intermediate densities and an intermediate UV radiation field. Contrary to the well-studied Orion Molecular Cloud 1 (OMC1), which hosts much harsher conditions, L1630 has little star formation. We aim to relate the [CII] fine-structure line emission to the physical conditions predominant in L1630 and compare it to studies of OMC1. The [CII] $158\,\mu\mathrm{m}$ emission from an area of $12' \times 17'$ in L1630 was observed using the upGREAT instrument onboard SOFIA. Of the [CII] emission from the mapped area 95%, $13\,L_{\odot}$, originates from the molecular cloud; the adjacent HII region contributes only 5%, that is, $1\,L_{\odot}$. From comparison with other data (CO (1-0)-line emission, far-infrared (FIR) continuum studies, emission from polycyclic aromatic hydrocarbons (PAHs)), we infer a gas density of the molecular cloud of $n_{\mathrm{H}}\sim 3\cdot 10^3\,\mathrm{cm^{-3}}$, with surface layers, including the Horsehead Nebula, having a density of up to $n_{\mathrm{H}}\sim 4\cdot 10^4\,\mathrm{cm^{-3}}$. The temperature of the surface gas is $T\sim 100\,\mathrm{K}$. The average [CII] cooling efficiency within the molecular cloud is $1.3\cdot 10^{-2}$. The fraction of the mass of the molecular cloud within the studied area that is traced by [CII] is only $8\%$. Our PDR models are able to reproduce the FIR-[CII] correlations and also the CO (1-0)-[CII] correlations. Finally, we compare our results on the heating efficiency of the gas with theoretical studies of photoelectric heating by PAHs, clusters of PAHs, and very small grains, and find the heating efficiency to be lower than theoretically predicted, a continuation of the trend set by other observations.

The [CII] 158um fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, [CII] balances the heating, including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of [CII] emission remains unclear, because C+ can be found in multiple phases of the interstellar medium. Here we measure the fractions of [CII] emission originating in the ionized and neutral gas phases of a sample of nearby galaxies. We use the [NII] 205um fine-structure line to trace the ionized medium, thereby eliminating the strong density dependence that exists in the ratio of [CII]/[NII] 122um. Using the FIR [CII] and [NII] emission detected by the KINGFISH and Beyond the Peak Herschel programs, we show that 60-80% of [CII] emission originates from neutral gas. We find that the fraction of [CII] originating in the neutral medium has a weak dependence on dust temperature and the surface density of star formation, and a stronger dependence on the gas-phase metallicity. In metal-rich environments, the relatively cooler ionized gas makes substantially larger contributions to total [CII] emission than at low abundance, contrary to prior expectations. Approximate calibrations of this metallicity trend are provided.

We have obtained estimates for the cosmic-ray ionization rate (CRIR) in the Galactic disk, using a detailed model for the physics and chemistry of diffuse interstellar gas clouds to interpret previously-published measurements of the abundance of four molecular ions: ArH$^+$, OH$^+$, H$_2$O$^+$ and H$_3^+$. For diffuse $atomic$ clouds at Galactocentric distances in the range $R_g \sim 4 - 9$ kpc, observations of ArH$^+$, OH$^+$, and H$_2$O$^+$ imply a mean primary CRIR of $(2.2 \pm 0.3) \exp [(R_0-R_g)/4.7\,\rm{kpc}] \times 10^{-16} \rm \, s^{-1}$ per hydrogen atom, where $R_0=8.5$ kpc. Within diffuse $molecular$ clouds observed toward stars in the solar neighborhood, measurements of H$_3^+$ and H$_2$ imply a primary CRIR of $(2.3 \pm 0.6) \times 10^{-16}\,\,\rm s^{-1}$ per H atom, corresponding to a total ionization rate per H$_2$ molecule of $(5.3 \pm 1.1) \times 10^{-16}\,\,\rm s^{-1},$ in good accord with previous estimates. These estimates are also in good agreement with a rederivation, presented here, of the CRIR implied by recent observations of carbon and hydrogen radio recombination lines along the sight-line to Cas A. Here, our best-fit estimate for the primary CRIR is $2.9 \times 10^{-16}\,\,\rm s^{-1}$ per H atom. Our results show marginal evidence that the CRIR in diffuse molecular clouds decreases with cloud extinction, $A_{\rm V}({\rm tot})$, with a best-fit dependence $\propto A_{\rm V}({\rm tot})^{-1}$ for $A_{\rm V}({\rm tot}) \ge 0.5$.

We present 5-20 micron spectral maps of the reflection nebula NGC2023 obtained with the Infrared Spectrograph SL and SH modes on board the Spitzer Space Telescope which reveal emission from polycyclic aromatic hydrocarbons (PAHs), C60, and H2 superposed on a dust continuum. We show that several PAH emission bands correlate with each other and exhibit distinct spatial distributions revealing a spatial sequence with distance from the illuminating star. We explore the distinct morphology of the 6.2, 7.7 and 8.6 micron PAH bands and find that at least two spatially distinct components contribute to the 7--9 micron PAH emission in NGC2023. We report that the PAH features behave independently of the underlying plateaus. We present spectra of compact oval PAHs ranging in size from C_66 to C_210, determined computationally using density functional theory, and investigate trends in the band positions and relative intensities as a function of PAH size, charge and geometry. Based on the NASA Ames PAH database, we discuss the 7--9 micron components in terms of band assignments and relative intensities. We assign the plateau emission to very small grains with possible contributions from PAH clusters and identify components in the 7--9 micron emission that likely originates in these structures. Based on the assignments and the observed spatial sequence, we discuss the photochemical evolution of the interstellar PAH family as they are more and more exposed to the radiation field of the central star in the evaporative flows associated with the PDRs in NGC2023.

Dynamic and thermal processes regulate the structure of the multi-phase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $\mu$m, HI, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of $\sim$1 kiloparsec. We find a distribution of thermal pressures in the $P_{\rm th}/k\sim10^3-10^5$ K cm$^{-3}$ range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of $P_{\rm th}/k\approx3600$ K cm$^{-3}$. These results are consistent with thermal pressure measurements using other observational methods. We find that $P_{\rm th}$ increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral medium could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about $\sim30$% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks.

We present [CII] 158um measurements from over 15,000 resolved regions within 54 nearby galaxies of the KINGFISH program to investigate the so-called [CII] "line cooling deficit" long known to occur in galaxies with different luminosities. The [CII]/TIR ratio ranges from above 1% to below 0.1% in the sample, with a mean value of 0.48+-0.21%. We find that the surface density of 24um emission dominates this trend, with [CII]/TIR dropping as nuInu{24um} increases. Deviations from this overall decline are correlated with changes in the gas phase metal abundance, with higher metallicity associated with deeper deficits at a fixed surface brightness. We supplement the local sample with resolved [CII] measurements from nearby luminous infrared galaxies and high redshift sources from z=1.8-6.4, and find that star formation rate density drives a continuous trend of deepening [CII] deficit across six orders of magnitude in SFRD. The tightness of this correlation suggests that an approximate star formation rate density can be estimated directly from global measurements of [CII]/TIR, and a relation is provided to do so. Several low-luminosity AGN hosts in the sample show additional and significant central suppression of [CII]/TIR, but these deficit enhancements occur not in those AGN with the highest X-ray luminosities, but instead those with the highest central starlight intensities. Taken together, these results demonstrate that the [CII] cooling line deficit in galaxies likely arises from local physical phenomena in interstellar gas.

We present a general parameter study, in which the abundance of interstellar argonium (ArH$^+$) is predicted using a model for the physics and chemistry of diffuse interstellar gas clouds. Results have been obtained as a function of UV radiation field, cosmic-ray ionization rate, and cloud extinction. No single set of cloud parameters provides an acceptable fit to the typical ArH$^+$, OH$^+$ and $\rm H_2O^+$ abundances observed in diffuse clouds within the Galactic disk. Instead, the observed abundances suggest that ArH$^+$ resides primarily in a separate population of small clouds of total visual extinction of at most 0.02 mag per cloud, within which the column-averaged molecular fraction is in the range $10^{-5} - 10^{-2}$, while OH$^+$ and $\rm H_2O^+$ reside primarily in somewhat larger clouds with a column-averaged molecular fraction $\sim 0.2$. This analysis confirms our previous suggestion that the argonium molecular ion is a unique tracer of almost purely atomic gas.

The [NII] 122 and 205 \mu m transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in the electron density within individual galaxies can be as a high as a factor of ~50, frequently with strong radial gradients. We find that $n_e$ increases as a function of infrared color, dust-weighted mean starlight intensity, and star formation rate surface density ($\Sigma_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between [NII] surface brightness and $\Sigma_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code.

We present spectral data cubes of the [CI] 809GHz, 12CO 115GHz, 13CO 110GHz and HI 1.4GHz line emission from an 1 square degree region along the l = 328{\deg} (G328) sightline in the Galactic Plane. Emission arises principally from gas in three spiral arm crossings along the sight line. The distribution of the emission in the CO and [CI] lines is found to be similar, with the [CI] slightly more extended, and both are enveloped in extensive HI. Spectral line ratios per voxel in the data cubes are found to be similar across the entire extent of the Galaxy. However, towards the edges of the molecular clouds the [CI]/13CO and 12CO/13CO line ratios rise by ~50%, and the [CI]/HI ratio falls by ~10$%. We attribute this to these sightlines passing predominantly through the surfaces of photodissociation regions (PDRs), where the carbon is found mainly as C or C+, while the H2 is mostly molecular, and the proportion of atomic gas also increases. We undertake modelling of the PDR emission from low density molecular clouds excited by average interstellar radiation fields and cosmic-ray ionization to quantify this comparison, finding that depletion of sulfur and reduced PAH abundance is needed to match line fluxes and ratios. Roughly one-third of the molecular gas along the sightline is found to be associated with this surface region, where the carbon is largely not to be found in CO. ~10% of the atomic hydrogen along the sightline is cold gas within PDRs.

We combine Herschel observations of a total of 12 sources to construct the most uniform survey of HF and H2O in our Galactic disk. Both molecules are detected in absorption along all sight lines. The high spectral resolution of the Heterodyne Instrument for the Far-Infrared (HIFI) allows us to compare the HF and H2O distributions in 47 diffuse cloud components sampling the disk. We find that the HF and H2O velocity distributions follow each other almost perfectly and establish that HF and H2O probe the same gas-phase volume. Our observations corroborate theoretical predictions that HF is a sensitive tracer of H2 in diffuse clouds, down to molecular fractions of only a few percent. Using HF to trace H2 in our sample, we find that the N(H2O)-to-N(HF) ratio shows a narrow distribution with a median value of 1.51. Our results further suggest that H2O might be used as a tracer of H2 -within a factor 2.5- in the diffuse interstellar medium. We show that the measured factor of ~2.5 variation around the median is driven by true local variations in the H2O abundance relative to H2 throughout the disk. The latter variability allows us to test our theoretical understanding of the chemistry of oxygen-bearing molecules in the diffuse gas. We show that both gas-phase and grain-surface chemistry are required to reproduce our H2O observations. This survey thus confirms that grain surface reactions can play a significant role in the chemistry occurring in the diffuse interstellar medium n_H < 1000 cm^-3.

We report additional detections of the chloronium molecular ion, H$_2$Cl$^+$, toward four bright submillimeter continuum sources: G29.96, W49N, W51, and W3(OH). With the use of the HIFI instrument on the Herschel Space Observatory, we observed the $2_{12}-1_{01}$ transition of ortho-H$_2^{35}$Cl$^+$ at 781.627 GHz in absorption toward all four sources. Much of the detected absorption arises in diffuse foreground clouds that are unassociated with the background continuum sources and in which our best estimates of the $N({\rm H_2Cl^+})/N({\rm H})$ ratio lie in the range $(0.9 - 4.8) \times 10^{-9}$. These chloronium abundances relative to atomic hydrogen can exceed the predictions of current astrochemical models by up to a factor of 5. Toward W49N, we have also detected the $2_{12}-1_{01}$ transition of ortho-H$_2^{37}$Cl$^+$ at 780.053 GHz and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ at 485.418 GHz. These observations imply $\rm H_2^{35}Cl^+/H_2^{37}Cl^+$ column density ratios that are consistent with the solar system $^{35}$Cl/$^{37}$Cl isotopic ratio of 3.1, and chloronium ortho-to-para ratios consistent with 3, the ratio of spin statistical weights.

Photodissociation regions (PDRs) contain a large fraction of all of the interstellar matter in galaxies. Classical examples include the boundaries between ionized regions and molecular clouds in regions of massive star formation, marking the point where all of the photons energetic enough to ionize hydrogen have been absorbed. In this paper we determine the physical properties of the PDRs associated with the star forming regions IRAS 23133+6050 and S 106 and present them in the context of other Galactic PDRs associated with massive star forming regions. We employ Herschel PACS and SPIRE spectroscopic observations to construct a full 55-650 {\mu}m spectrum of each object from which we measure the PDR cooling lines, other fine- structure lines, CO lines and the total far-infrared flux. These measurements are then compared to standard PDR models. Subsequently detailed numerical PDR models are compared to these predictions, yielding additional insights into the dominant thermal processes in the PDRs and their structures. We find that the PDRs of each object are very similar, and can be characterized by a two-phase PDR model with a very dense, highly UV irradiated phase (n $\sim$ 10^6 cm^(-3), G$_0$ $\sim$ 10^5) interspersed within a lower density, weaker radiation field phase (n $\sim$ 10^4 cm^(-3), G$_0$ $\sim$ 10^4). We employed two different numerical models to investigate the data, firstly we used RADEX models to fit the peak of the $^{12}$CO ladder, which in conjunction with the properties derived yielded a temperature of around 300 K. Subsequent numerical modeling with a full PDR model revealed that the dense phase has a filling factor of around 0.6 in both objects. The shape of the $^{12}$CO ladder was consistent with these components with heating dominated by grain photoelectric heating. An extra excitation component for the highest J lines (J > 20) is required for S 106.

The manner in which gas accretes and orbits within circumnuclear rings has direct implications for the star formation process. In particular, gas may be compressed and shocked at the inflow points, resulting in bursts of star formation at these locations. Afterwards the gas and young stars move together through the ring. In addition, star formation may occur throughout the ring, if and when the gas reaches sufficient density to collapse under gravity. These two scenarios for star formation in rings are often referred to as the `pearls on a string' and `popcorn' paradigms. In this paper, we use new Herschel PACS observations, obtained as part of the KINGFISH Open Time Key Program, along with archival Spitzer and ground-based observations from the SINGS Legacy project, to investigate the heating and cooling of the interstellar medium in the nearby star-forming ring galaxy, NGC4736. By comparing spatially resolved estimates of the stellar FUV flux available for heating, with the gas and dust cooling derived from the FIR continuum and line emission, we show that while star formation is indeed dominant at the inflow points in NGC 4736, additional star formation is needed to balance the gas heating and cooling throughout the ring. This additional component most likely arises from the general increase in gas density in the ring over its lifetime. Our data provide strong evidence, therefore, for a combination of the two paradigms for star formation in the ring in NGC4736.

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta_H) and molecular hydrogen fraction, f(H2). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042+-0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H2 fractions. We also infer zeta_H throughout our sample, and find a log-normal distribution with mean log(zeta_H)=-15.75, (zeta_H=1.78x10^-16 s^-1), and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H3+ observations. Ionization rates in the Galactic center tend to be 10--100 times larger than found in the Galactic disk, also in accord with prior studies.

Recent studies have confirmed the presence of buckminsterfullerene (C$_{60}$) in different interstellar and circumstellar environments. However, several aspects regarding C$_{60}$ in space are not well understood yet, such as the formation and excitation processes, and the connection between C$_{60}$ and other carbonaceous compounds in the interstellar medium, in particular polycyclic aromatic hydrocarbons (PAHs). In this paper we study several photodissociation regions (PDRs) where C$_{60}$ and PAHs are detected and the local physical conditions are reasonably well constrained, to provide observational insights into these questions. C$_{60}$ is found to emit in PDRs where the dust is cool ($T_d = 20-40$ K) and even in PDRs with cool stars. These results exclude the possibility for C$_{60}$ to be locked in grains at thermal equilibrium in these environments. We observe that PAH and C$_{60}$ emission are spatially uncorrelated and that C$_{60}$ is present in PDRs where the physical conditions (in terms of radiation field and hydrogen density) allow for full dehydrogenation of PAHs, with the exception of Ced 201. We also find trends indicative of an increase in C$_{60}$ abundance within individual PDRs, but these trends are not universal. These results support models where the dehydrogenation of carbonaceous species is the first step towards C$_{60}$ formation. However, this is not the only parameter involved and C$_{60}$ formation is likely affected by shocks and PDR age.

The [CII] 157.74 $\mu$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $\Sigma_{\rm [CII]}-\Sigma_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $\Sigma_{\rm[CII]}-\Sigma_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $\Sigma_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies.

We describe the assignment of a previously unidentified interstellar absorption line to ArH$^+$ and discuss its relevance in the context of hydride absorption in diffuse gas with a low H$_2$ fraction. The column densities along several lines of sight are determined and discussd in the framework of chemical models. The column densities of ArH$^+$ are compared to those of other species, tracing interstellar medium (ISM) components with different H$_2$ abundances. Chemical models are constructed, taking UV radiation and cosmic ray ionization into account. Due to the detection of two isotopologues, $^{36}$ArH$^+$ and $^{38}$ArH$^+$, we are confident about the carrier assignment to ArH$^+$. NeH$^+$ is not detected with a limit of [NeH$^+$]/[ArH$^+$] $\le$ 0.1. The derived column densities agree well with the predictions of chemical models. ArH$^+$ is a unique tracer of gas with a fractional H$_2$ abundance of $10^{-4}- 10^{-3}$ and shows little correlation with H$_2$O$^+$, which traces gas with a fractional H$_2$ abundance of $\approx $0.1. A careful analysis of variations in the ArH$^+$, OH$^+$, H$_2$O$^+$ and HF column densities promises to be a faithful tracer of the distribution of the H$_2$ fractional abundance, providing unique information on a poorly known phase in the cycle of interstellar matter, its transition from atomic diffuse gas to dense molecular gas traced by CO emission. Abundances of these species put strong observational constraints upon magnetohydrodynamical (MHD) simulations of the interstellar medium, and potentially could evolve into a tool to characterize the ISM. Paradoxically, the ArH$^+$ molecule is a better tracer of \new{almost} purely atomic hydrogen gas than H{\sc i} itself, since H{\sc i} can also be present in gas with a significant molecular content, but ArH$^+$ singles out gas that is $>99.9$\% atomic.

We derive the CO-to-H2 conversion factor, X_CO = N(H2)/I_CO, across the Perseus molecular cloud on sub-parsec scales by combining the dust-based N(H2) data with the I_CO data from the COMPLETE Survey. We estimate an average X_CO ~ 3 x 10^19 cm^-2 K^-1 km^-1 s and find a factor of ~3 variations in X_CO between the five sub-regions in Perseus. Within the individual regions, X_CO varies by a factor of ~100, suggesting that X_CO strongly depends on local conditions in the interstellar medium. We find that X_CO sharply decreases at Av < 3 mag but gradually increases at Av > 3 mag, with the transition occurring at Av where I_CO becomes optically thick. We compare the N(HI), N(H2), I_CO, and X_CO distributions with two models of the formation of molecular gas, a one-dimensional photodissociation region (PDR) model and a three-dimensional magnetohydrodynamic (MHD) model tracking both the dynamical and chemical evolution of gas. The PDR model based on the steady state and equilibrium chemistry reproduces our data very well but requires a diffuse halo to match the observed N(HI) and I_CO distributions. The MHD model generally matches our data well, suggesting that time-dependent effects on H2 and CO formation are insignificant for an evolved molecular cloud like Perseus. However, we find interesting discrepancies, including a broader range of N(HI), likely underestimated I_CO, and a large scatter of I_CO at small Av. These discrepancies likely result from strong compressions/rarefactions and density fluctuations in the MHD model.

We present spectral line images of [CI] 809 GHz, CO J=1-0 115 GHz and HI 1.4 GHz line emission, and calculate the corresponding C, CO and H column densities, for a sinuous, quiescent Giant Molecular Cloud about 5 kpc distant along the l=328{\deg} sightline (hereafter G328) in our Galaxy. The [CI] data comes from the High Elevation Antarctic Terahertz (HEAT) telescope, a new facility on the summit of the Antarctic plateau where the precipitable water vapor falls to the lowest values found on the surface of the Earth. The CO and HI datasets come from the Mopra and Parkes/ATCA telescopes, respectively. We identify a filamentary molecular cloud, ~75 x 5 pc long with mass ~4 x 10E4 Msun and a narrow velocity emission range of just 4 km/s. The morphology and kinematics of this filament are similar in CO, [CI] and HI, though in the latter appears as self-absorption. We calculate line fluxes and column densities for the three emitting species, which are broadly consistent with a PDR model for a GMC exposed to the average interstellar radiation field. The [C/CO] abundance ratio averaged through the filament is found to be approximately unity. The G328 filament is constrained to be cold (Tdust < 20K) by the lack of far-IR emission, to show no clear signs of star formation, and to only be mildly turbulent from the narrow line width. We suggest that it may represent a GMC shortly after formation, or perhaps still be in the process of formation.

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows.

We observed H2 line emission with Spitzer-IRS toward M17-SW and modeled the data with our PDR code. Derived gas density values of up to few times 10^7 cm^-3 indicate that H2 emission originates in high-density clumps. We discover that the PDR code can be utilized to map the amount of intervening extinction obscuring the H2 emission layers, and thus we obtain the radial profile of A_V relative to the central ionizing cluster NGC 6618. The extinction has a positive radial gradient, varying between 15--47 mag over the projected distance of 0.9--2.5 pc from the primary ionizer, CEN 1. These high extinction values are in good agreement with previous studies of A_V toward stellar targets in M17-SW. The ratio of data to PDR model values is used to infer the global line-of-sight structure of the PDR surface, which is revealed to resemble a concave surface relative to NGC 6618. Such a configuration confirms that this PDR can be described as a bowl-shaped boundary of the central H II region in M17. The derived structure and physical conditions are important for interpreting the fine-structure and rotational line emission from the PDR.

We present kiloparsec (kpc) spatial resolution maps of the CO-to-H2 conversion factor (alpha_co) and dust-to-gas ratio (DGR) in 26 nearby, star-forming galaxies. We have simultaneously solved for alpha_co and DGR by assuming that the DGR is approximately constant on kpc scales. With this assumption, we can combine maps of dust mass surface density, CO integrated intensity and HI column density to solve for both alpha_co and DGR with no assumptions about their value or dependence on metallicity or other parameters. Such a study has just become possible with the availability of high resolution far-IR maps from the Herschel key program KINGFISH, 12CO J=(2-1) maps from the IRAM 30m large program HERACLES and HI 21-cm line maps from THINGS. We use a fixed ratio between the (2-1) and (1-0) lines to present our alpha_co results on the more typically used 12CO J=(1-0) scale and show using literature measurements that variations in the line ratio do not effect our results. In total, we derive 782 individual solutions for alpha_co and DGR. On average, alpha_co = 3.1 Msun pc^-2 (K km s^-1)^-1 for our sample with a standard deviation of 0.3 dex. Within galaxies we observe a generally flat profile of alpha_co as a function of galactocentric radius. However, most galaxies exhibit a lower alpha_co in the central kpc---a factor of ~2 below the galaxy mean, on average. In some cases, the central alpha_co value can be factors of 5 to 10 below the standard Milky Way (MW) value of alpha_co,MW =4.4 Msun pc^-2 (K km s^-1)^-1. While for alpha_co we find only weak correlations with metallicity, DGR is well-correlated with metallicity, with an approximately linear slope. Finally, we present several recommendations for choosing an appropriate alpha_co for studies of nearby galaxies.

We present the first results from a new carbon monoxide (CO) survey of the southern Galactic plane being conducted with the Mopra radio telescope in Australia. The 12CO, 13CO and C18O J=1-0 lines are being mapped over the l = 305-345 deg, b = +/- 0.5 deg portion of the 4th quadrant of the Galaxy, at 35" spatial and 0.1 km/s spectral resolution. The survey is being undertaken with two principal science objectives: (i) to determine where and how molecular clouds are forming in the Galaxy and (ii) to probe the connection between molecular clouds and the "missing" gas inferred from gamma-ray observations. We describe the motivation for the survey, the instrumentation and observing techniques being applied, and the data reduction and analysis methodology. In this paper we present the data from the first degree surveyed, l = 323-324 deg, b = +/- 0.5 deg. We compare the data to the previous CO survey of this region and present metrics quantifying the performance being achieved; the rms sensitivity per 0.1 km/s velocity channel is ~1.5K for 12CO and ~0.7K for the other lines. We also present some results from the region surveyed, including line fluxes, column densities, molecular masses, 12CO/13CO line ratios and 12CO optical depths. We also examine how these quantities vary as a function of distance from the Sun when averaged over the 1 square degree survey area. Approximately 2 x 10E6 MSun of molecular gas is found along the G323 sightline, with an average H2 number density of nH2 ~ 1 cm-3 within the Solar circle. The CO data cubes will be made publicly available as they are published.

CO line emission represents the most accessible and widely used tracer of the molecular interstellar medium. This renders the translation of observed CO intensity into total H2 gas mass critical to understand star formation and the interstellar medium in our Galaxy and beyond. We review the theoretical underpinning, techniques, and results of efforts to estimate this CO-to-H2 "conversion factor," Xco, in different environments. In the Milky Way disk, we recommend a conversion factor Xco = 2x10^{20} cm^-2/(K km/s)^-1 with +/-30% uncertainty. Studies of other "normal galaxies" return similar values in Milky Way-like disks, but with greater scatter and systematic uncertainty. Departures from this Galactic conversion factor are both observed and expected. Dust-based determinations, theoretical arguments, and scaling relations all suggest that Xco increases with decreasing metallicity, turning up sharply below metallicity ~1/3-1/2 solar in a manner consistent with model predictions that identify shielding as a key parameter. Based on spectral line modeling and dust observations, Xco appears to drop in the central, bright regions of some but not all galaxies, often coincident with regions of bright CO emission and high stellar surface density. This lower Xco is also present in the overwhelmingly molecular interstellar medium of starburst galaxies, where several lines of evidence point to a lower CO-to-H2 conversion factor. At high redshift, direct evidence regarding the conversion factor remains scarce; we review what is known based on dynamical modeling and other arguments.

The first complete submillimetre spectrum (190-670um) of the Seyfert 2 galaxy NGC1068 has been observed with the SPIRE Fourier Transform Spectrometer onboard the {\it Herschel} Space Observatory. The sequence of CO lines (Jup=4-13), lines from water, the fundamental rotational transition of HF, two o-H_2O+ lines and one line each from CH+ and OH+ have been detected, together with the two [CI] lines and the [NII]205um line. The observations in both single pointing mode with sparse image sampling and in mapping mode with full image sampling allow us to disentangle two molecular emission components, one due to the compact circum-nuclear disk (CND) and one from the extended region encompassing the star forming ring (SF-ring). Radiative transfer models show that the two CO components are characterized by density of n(H_2)=10^4.5 and 10^2.9 cm^-3 and temperature of T=100K and 127K, respectively. The comparison of the CO line intensities with photodissociation region (PDR) and X-ray dominated region (XDR) models, together with other observational constraints, such as the observed CO surface brightness and the radiation field, indicate that the best explanation for the CO excitation of the CND is an XDR with density of n(H_2) 10^4 cm^-3 and X-ray flux of 9 erg s^-1 cm^-2, consistent with illumination by the active galactic nucleus, while the CO lines in the SF-ring are better modeled by a PDR. The detected water transitions, together with those observed with the \her \sim PACS Spectrometer, can be modeled by an LVG model with low temperature (T_kin \sim 40K) and high density (n(H_2) in the range 10^6.7-10^7.9 cm^-3).

We characterize the dust in NGC628 and NGC6946, two nearby spiral galaxies in the KINGFISH sample. With data from 3.6um to 500um, dust models are strongly constrained. Using the Draine & Li (2007) dust model, (amorphous silicate and carbonaceous grains), for each pixel in each galaxy we estimate (1) dust mass surface density, (2) dust mass fraction contributed by polycyclic aromatic hydrocarbons (PAH)s, (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in regions with high starlight intensity. We obtain maps for the dust properties, which trace the spiral structure of the galaxies. The dust models successfully reproduce the observed global and resolved spectral energy distributions (SEDs). The overall dust/H mass ratio is estimated to be 0.0082+/-0.0017 for NGC628, and 0.0063+/-0.0009 for NGC6946, consistent with what is expected for galaxies of near-solar metallicity. Our derived dust masses are larger (by up to a factor 3) than estimates based on single-temperature modified blackbody fits. We show that the SED fits are significantly improved if the starlight intensity distribution includes a (single intensity) "delta function" component. We find no evidence for significant masses of cold dust T<12K. Discrepancies between PACS and MIPS photometry in both low and high surface brightness areas result in large uncertainties when the modeling is done at PACS resolutions, in which case SPIRE, MIPS70 and MIPS160 data cannot be used. We recommend against attempting to model dust at the angular resolution of PACS.

Taking advantage of the sensitivity and angular resolution of the Herschel Space Observatory at far-infrared and submm wavelengths, we aim to characterize the physical properties of cold dust within nearby galaxies and study the robustness of the parameters we derive using different modified blackbody models. For a pilot subsample of the KINGFISH program, we perform 2 temperature fits of the Spitzer and Herschel photometric data (24 to 500um), with a warm and a cold component, globally and in each resolution element.At global scales, we observe ranges of values for beta_c(0.8 to 2.5) and Tc(19.1 to 25.1K).We compute maps of our parameters with beta fixed or free to test the robustness of the temperature and dust surface density maps we deduce. When the emissivity is fixed, we observe temperature gradients as a function of radius.When the emissivity is fitted as a free parameter, barred galaxies tend to have uniform fitted emissivities.Gathering resolved elements in a Tc-beta_c diagram underlines an anti-correlation between the two parameters.It remains difficult to assess whether the dominant effect is the physics of dust grains, noise, or mixing along the line of sight and in the beam. We finally observe in both cases that the dust column density peaks in central regions of galaxies and bar ends (coinciding with molecular gas density enhancements usually found in these locations).We also quantify how the total dust mass varies with our assumptions about the emissivity index as well as the influence of the wavelength coverage used in the fits. We show that modified blackbody fits using a shallow emissivity (beta_c < 2.0) lead to significantly lower dust masses compared to the beta_c < 2.0 case, with dust masses lower by up to 50% if beta_c=1.5 for instance.The working resolution affects our total dust mass estimates: masses increase from global fits to spatially-resolved fits.

We model the production of OH+, H2O+, and H3O+ in interstellar clouds, using a steady state photodissociation region code that treats the freeze-out of gas species, grain surface chemistry, and desorption of ices from grains. The code includes PAHs, which have important effects on the chemistry. All three ions generally have two peaks in abundance as a function of depth into the cloud, one at A_V<~1 and one at A_V~3-8, the exact values depending on the ratio of incident ultraviolet flux to gas density. For relatively low values of the incident far ultraviolet flux on the cloud ({\chi}<~ 1000; {\chi}= 1= local interstellar value), the columns of OH+ and H2O+ scale roughly as the cosmic ray primary ionization rate {\zeta}(crp) divided by the hydrogen nucleus density n. The H3O+ column is dominated by the second peak, and we show that if PAHs are present, N(H3O+) ~ 4x10^{13} cm^{-2} independent of {\zeta}(crp) or n. If there are no PAHs or very small grains at the second peak, N(H3O+) can attain such columns only if low ionization potential metals are heavily depleted. We also model diffuse and translucent clouds in the interstellar medium, and show how observations of N(OH+)/N(H) and N(OH+)/N(H2O+) can be used to estimate {\zeta}(crp)/n, {\chi}/n and A_V in them. We compare our models to Herschel observations of these two ions, and estimate {\zeta}(crp) ~ 4-6 x 10^-16 (n/100 cm^-3) s^-1 and \chi/n = 0.03 cm^3 for diffuse foreground clouds towards W49N.

NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared (FIR) [CII]158um and [OI]63um lines and mid-infrared H2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([CII]158um+[OI]63um)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7um PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [CII]158{\mu}m/PAH(5.5 - 14um) is found. PAHs in the ring are responsible for a factor of two more [CII]158um and [OI]63um emission per unit mass than PAHs in the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high intensity photodissociation regions (PDRs), in which case G0 ~ 10^2.3 and nH ~ 10^3.5 cm^-3 in the ring. For these values of G0 and nH PDR models cannot reproduce the observed H2 emission. Much of the the H2 emission in the starburst ring could come from warm regions in the diffuse ISM that are heated by turbulent dissipation or shocks.

Using the Herschel Space Observatory's Heterodyne Instrument for the Far-Infrared (HIFI), we have observed para-chloronium (H2Cl+) toward six sources in the Galaxy. We detected interstellar chloronium absorption in foreground molecular clouds along the sight-lines to the bright submillimeter continuum sources Sgr A (+50 km/s cloud) and W31C. Both the para-H2-35Cl+ and para-H2-37Cl+ isotopologues were detected, through observations of their 1(11)-0(00) transitions at rest frequencies of 485.42 and 484.23 GHz, respectively. For an assumed ortho-to-para ratio of 3, the observed optical depths imply that chloronium accounts for ~ 4 - 12% of chlorine nuclei in the gas phase. We detected interstellar chloronium emission from two sources in the Orion Molecular Cloud 1: the Orion Bar photodissociation region and the Orion South condensation. For an assumed ortho-to-para ratio of 3 for chloronium, the observed emission line fluxes imply total beam-averaged column densities of ~ 2.0E+13 cm-2 and ~ 1.2E+13 cm-2, respectively, for chloronium in these two sources. We obtained upper limits on the para-H2-35Cl+ line strengths toward H2 Peak 1 in the Orion Molecular cloud and toward the massive young star AFGL 2591. The chloronium abundances inferred in this study are typically at least a factor ~10 larger than the predictions of steady-state theoretical models for the chemistry of interstellar molecules containing chlorine. Several explanations for this discrepancy were investigated, but none has proven satisfactory, and thus the large observed abundances of chloronium remain puzzling.

The physical state of interstellar gas and dust is dependent on the processes which heat and cool this medium. To probe heating and cooling of the ISM over a large range of infrared surface brightness, on sub-kiloparsec scales, we employ line maps of [C \ii] 158 $\mu$m, [O \one] 63 $\mu$m, and [N \ii] 122 $\mu$m in NGC 1097 and NGC 4559, obtained with the PACS spectrometer onboard {\it Herschel}. We matched new observations to existing Spitzer-IRS data that trace the total emission of polycyclic aromatic hydrocarbons (PAHs). We confirm at small scales in these galaxies that the canonical measure of photoelectric heating efficiency, ([C \ii] + [O \one])/TIR, decreases as the far-infrared color, $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m), increases. In contrast, the ratio of far-infrared (far-IR) cooling to total PAH emission, ([C \ii] + [O \one])/PAH, is a near constant $\sim$6% over a wide range of far-infrared color, 0.5 \textless\ $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m) $\lesssim$ 0.95. In the warmest regions, where $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m) $\gtrsim$ 0.95, the ratio ([C \ii] + [O \one])/PAH drops rapidly to 4%. We derived representative values of the local UV radiation density, $G_0$, and the gas density, $n_H$, by comparing our observations to models of photodissociation regions. The ratio $G_0/n_H$, derived from fine-structure lines, is found to correlate with the mean dust-weighted starlight intensity, $<U>$ derived from models of the IR SED. Emission from regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small grains are susceptible to ionization effects. We note that there is a shift in the 7.7 / 11.3 $\mu$m PAH ratio in regions that exhibit a deficit in ([C \ii] + [O \one])/PAH, suggesting that small grains are ionized in these environments.

Aims. Young stars interact vigorously with their surroundings, as evident from the highly rotationally excited CO (up to Eup=4000 K) and H2O emission (up to 600 K) detected by the Herschel Space Observatory in embedded low-mass protostars. Our aim is to construct a model that reproduces the observations quantitatively, to investigate the origin of the emission, and to use the lines as probes of the various heating mechanisms. Methods. The model consists of a spherical envelope with a bipolar outflow cavity. Three heating mechanisms are considered: passive heating by the protostellar luminosity, UV irradiation of the outflow cavity walls, and C-type shocks along the cavity walls. Line fluxes are calculated for CO and H2O and compared to Herschel data and complementary ground-based data for the protostars NGC1333 IRAS2A, HH 46 and DK Cha. The three sources are selected to span a range of evolutionary phases and physical characteristics. Results. The passively heated gas in the envelope accounts for 3-10% of the CO luminosity summed over all rotational lines up to J=40-39; it is best probed by low-J CO isotopologue lines such as C18O 2-1 and 3-2. The UV-heated gas and the C-type shocks, probed by 12CO 10-9 and higher-J lines, contribute 20-80% each. The model fits show a tentative evolutionary trend: the CO emission is dominated by shocks in the youngest source and by UV-heated gas in the oldest one. This trend is mainly driven by the lower envelope density in more evolved sources. The total H2O line luminosity in all cases is dominated by shocks (>99%). The exact percentages for both species are uncertain by at least a factor of 2 due to uncertainties in the gas temperature as function of the incident UV flux. However, on a qualitative level, both UV-heated gas and C-type shocks are needed to reproduce the emission in far-infrared rotational lines of CO and H2O.

We present 15-20 um spectral maps towards the reflection nebula NGC2023 obtained with the Infrared Spectrograph in short-wavelength, high-resolution mode on board the Spitzer Space Telescope. These spectra reveal emission from PAHs, C60, and H2 superposed on a dust continuum. These emission components exhibit distinct spatial distributions: with increasing distance from the illuminating star, we observe the PAH emission followed by the dust continuum emission and the H2 emission. The C60 emission is located closest to the illuminating star in the south while in the north, it seems to be associated with the H/H2 transition. Emission from PAHs and PAH-related species produce features at 15.8, 16.4, 17.4, and 17.8 um and the 15-18 um plateau. These different PAH features show distinct spatial distributions. The 15.8 um band and 15-18 um plateau correlate with the 11.2 um PAH band and with each other, and are attributed to large, neutral PAHs. Conversely, the 16.4 um feature correlates with the 12.7 um PAH band, suggesting that both arise from species that are favored by the same conditions that favor PAH cations. The PAH contribution to the 17.4 um band is displaced towards the illuminating star with respect to the 11.2 and 12.7 um emission and is assigned to doubly ionized PAHs and/or a subset of cationic PAHs. The spatial distribution of the 17.8 um band suggests it arises from both neutral and cationic PAHs. In contrast to their intensities, the profiles of the PAH bands and the 15-18 um plateau do not vary spatially. Consequently, we conclude that the carrier of the 15-18 um plateau is distinct from that of the PAH bands.

New far-infrared and sub-millimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially-integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500um emission shows evidence for a sub-millimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photo-dissociation regions is found to be (21+/-4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.

The KINGFISH project (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) is an imaging and spectroscopic survey of 61 nearby (d < 30 Mpc) galaxies, chosen to cover a wide range of galaxy properties and local interstellar medium (ISM) environments found in the nearby Universe. Its broad goals are to characterize the ISM of present-day galaxies, the heating and cooling of their gaseous and dust components, and to better understand the physical processes linking star formation and the ISM. KINGFISH is a direct descendant of the Spitzer Infrared Nearby Galaxies Survey (SINGS), which produced complete Spitzer imaging and spectroscopic mapping and a comprehensive set of multi-wavelength ancillary observations for the sample. The Herschel imaging consists of complete maps for the galaxies at 70, 100, 160, 250, 350, and 500 microns. The spectal line imaging of the principal atomic ISM cooling lines ([OI]63um, [OIII]88um, [NII]122,205um, and [CII]158um) covers the subregions in the centers and disks that already have been mapped in the mid-infrared with Spitzer. The KINGFISH and SINGS multi-wavelength datasets combined provide panchromatic mapping of the galaxies sufficient to resolve individual star-forming regions, and tracing the important heating and cooling channels of the ISM, across a wide range of local extragalactic ISM environments. This paper summarizes the scientific strategy for KINGFISH, the properties of the galaxy sample, the observing strategy, and data processing and products. It also presents a combined Spitzer and Herschel image atlas for the KINGFISH galaxies, covering the wavelength range 3.6 -- 500 microns. All imaging and spectroscopy data products will be released to the Herschel user generated product archives.

We use the IRS on Spitzer to observe the southern part of the reflection nebula NGC 2023, including the Southern Ridge, which is a photodissociation region par excellence excited by HD 37903. Five pure-rotational H2 emission lines are detected and mapped over and around the Southern Ridge in order to compare with predicted level column densities from theoretical PDR models. We find very good agreement between PDR model predictions and emission line intensities and ratios measured with Spitzer, leading us to conclude that grain photoelectric heating sufficiently warms the gas to produce the observed H2 line emission via collisional excitation. On the Southern Ridge, we infer a hydrogen nucleus density n_H approx 2x10^5 cm^-3 and radiation field strength \chi\ approx 10^4 relative to the local Galactic interstellar radiation field. This high value for \chi\ independently predicts a distance toward HD 37903 of 300 pc, and is consistent with the most recent HIPPARCOS results. Over the map we find that both n_H and \chi\ vary by a factor of sim 3. Such 2-D variations provide clues about the underlying 3-D structure of the Southern Ridge field, which appears to be the tip of a molecular cloud. We also map variations in excitation temperature and the ortho-to-para ratio, the latter attaining values of sim 1.5---2.0 on the Southern Ridge, and find that PDR modeling can readily reproduce observed ortho-to-para ratios that are < 3 for rotational excitation dominated by collisional processes. Last, the stars Sellgren C and G are discovered to be resolved on archival HST images into two point sources each, with separations of lesssim 0.5".

We present Herschel SPIRE-FTS observations of Arp~220, a nearby ULIRG. The FTS continuously covers 190 -- 670 microns, providing a good measurement of the continuum and detection of several molecular and atomic species. We detect luminous CO (J = 4-3 to 13-12) and water ladders with comparable total luminosity; very high-J HCN absorption; OH+, H2O+, and HF in absorption; and CI and NII. Modeling of the continuum yields warm dust, with T = 66 K, and an unusually large optical depth of ~5 at 100 microns. Non-LTE modeling of the CO shows two temperature components: cold molecular gas at T ~ 50 K and warm molecular gas at T ~1350 K. The mass of the warm gas is 10% of the cold gas, but dominates the luminosity of the CO ladder. The temperature of the warm gas is in excellent agreement with H2 rotational lines. At 1350 K, H2 dominates the cooling (~20 L_sun/M_sun) in the ISM compared to CO (~0.4 L_sun/M_sun). We found that only a non-ionizing source such as the mechanical energy from supernovae and stellar winds can excite the warm gas and satisfy the energy budget of ~20 L_sun/M_sun. We detect a massive molecular outflow in Arp 220 from the analysis of strong P-Cygni line profiles observed in OH+, H2O+, and H2O. The outflow has a mass > 10^{7} M_sun and is bound to the nuclei with velocity < 250 km/s. The large column densities observed for these molecular ions strongly favor the existence of an X-ray luminous AGN (10^{44} ergs/s) in Arp 220.

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars.

We present the detection of extended (~30 kpc^2) dust emission in the tidal \hi\ arm near NGC 3077 (member of the M\,81 triplet) using SPIRE on board Herschel. Dust emission in the tidal arm is typically detected where the \hi\ column densities are >10^21 cm^-2. The SPIRE band ratios show that the dust in the tidal arm is significantly colder (~13 K) than in NGC 3077 itself (~31 K), consistent with the lower radiation field in the tidal arm. The total dust mass in the tidal arm is ~1.8 x 10^6 M_sun (assuming beta=2), i.e. substantially larger than the dust mass associated with NGC 3077 (~2 x 10^5 M_sun). Where dust is detected, the dust-to-gas ratio is 6+/-3 x 10^-3, consistent within the uncertainties with what is found in NGC 3077 and nearby spiral galaxies with Galactic metallicities. The faint HII regions in the tidal arm can not be responsible for the detected enriched material and are not the main source of the dust heating in the tidal arm. We conclude that the interstellar medium (atomic HI, molecules and dust) in this tidal feature was pre-enriched and stripped off NGC 3077 during its recent interaction (~3 x 10^8 yr ago) with M 82 and M 81. This implies that interaction can efficiently remove heavy elements and enriched material (dust, molecular gas) from galaxies. As interactions were more frequent at large lookback times, it is conceivable that they could substantially contribute (along with galactic outflows) to the enrichment of the intergalactic medium.

Using photometry of NGC 1097 from the Herschel PACS (Photodetector Array Camera and Spectrometer) instrument, we study the resolved properties of thermal dust continuum emission from a circumnuclear starburst ring with a radius ~ 900 pc. These observations are the first to resolve the structure of a circumnuclear ring at wavelengths that probe the peak (i.e. lambda ~ 100 micron) of the dust spectral energy distribution. The ring dominates the far-infrared (far-IR) emission from the galaxy - the high angular resolution of PACS allows us to isolate the ring's contribution and we find it is responsible for 75, 60 and 55% of the total flux of NGC 1097 at 70, 100 and 160 micron, respectively. We compare the far-IR structure of the ring to what is seen at other wavelengths and identify a sequence of far-IR bright knots that correspond to those seen in radio and mid-IR images. The mid- and far-IR band ratios in the ring vary by less than +/- 20% azimuthally, indicating modest variation in the radiation field heating the dust on ~ 600 pc scales. We explore various explanations for the azimuthal uniformity in the far-IR colors of the ring including a lack of well-defined age gradients in the young stellar cluster population, a dominant contribution to the far-IR emission from dust heated by older (> 10 Myr) stars and/or a quick smoothing of local enhancements in dust temperature due to the short orbital period of the ring. Finally, we improve previous limits on the far-IR flux from the inner ~ 600 pc of NGC 1097 by an order of magnitude, providing a better estimate of the total bolometric emission arising from the active galactic nucleus and its associated central starburst.

Stellar density and bar strength should affect the temperatures of the cool (T ~ 20-30 K) dust component in the inner regions of galaxies, which implies that the ratio of temperatures in the circumnuclear regions to the disk should depend on Hubble type. We investigate the differences between cool dust temperatures in the central 3 kpc and disk of 13 nearby galaxies by fitting models to measurements between 70 and 500 microns. We attempt to quantify temperature trends in nearby disk galaxies, with archival data from Spitzer/MIPS and new observations with Herschel/SPIRE, which were acquired during the first phases of the Herschel observations for the KINGFISH (key insights in nearby galaxies: a far-infrared survey with Herschel) sample. We fit single-temperature modified blackbodies to far-infrared and submillimeter measurements of the central and disk regions of galaxies to determine the temperature of the component(s) emitting at those wavelengths. We present the ratio of central-region-to-disk-temperatures of the cool dust component of 13 nearby galaxies as a function of morphological type. We find a significant temperature gradient in the cool dust component in all galaxies, with a mean center-to-disk temperature ratio of 1.15 +/- 0.03. The cool dust temperatures in the central ~3 kpc of nearby galaxies are 23(+/-3)% hotter for morphological types earlier than Sc, and only 9(+/-3)% hotter for later types. The temperature ratio is also correlated with bar strength, with only strongly barred galaxies having a ratio over 1.2. The strong radiation field in the high stellar density of a galactic bulge tends to heat the cool dust component to higher temperatures, at least in early-type spirals with relatively large bulges, especially when paired with a strong bar.

NGC 1097 is a nearby SBb galaxy with a Seyfert nucleus and a bright starburst ring. We study the physical properties of the interstellar medium (ISM) in the ring using spatially resolved far-infrared spectral maps of the circumnuclear starburst ring of NGC 1097, obtained with the PACS spectrometer on board the Herschel Space Telescope. In particular, we map the important ISM cooling and diagnostic emission lines of [OI] 63 $\mu$m, [OIII] 88 $\mu$m, [NII] 122 $\mu$m, [CII] 158 $\mu$m and [NII] 205 $\mu$m. We observe that in the [OI] 63 $\mu$m, [OIII] 88 $\mu$m, and [NII] 122 $\mu$m line maps, the emission is enhanced in clumps along the NE part of the ring. We observe evidence of rapid rotation in the circumnuclear ring, with a rotation velocity of ~220$ km s$^{-1}$ (inclination uncorrected) measured in all lines. The [OI] 63 $\mu$m/[CII] 158 $\mu$m ratio varies smoothly throughout the central region, and is enhanced on the northeastern part of the ring, which may indicate a stronger radiation field. This enhancement coincides with peaks in the [OI] 63 $\mu$m and [OIII] 88 $\mu$m maps. Variations of the [NII] 122 $\mu$m/[NII] 205 $\mu$m ratio correspond to a range in the ionized gas density between 150 and 400 cm$^{-3}$.

The mass of molecular gas in an interstellar cloud is often measured using line emission from low rotational levels of CO, which are sensitive to the CO mass, and then scaling to the assumed molecular hydrogen H_2 mass. However, a significant H_2 mass may lie outside the CO region, in the outer regions of the molecular cloud where the gas phase carbon resides in C or C+. Here, H_2 self-shields or is shielded by dust from UV photodissociation, where as CO is photodissociated. This H_2 gas is "dark" in molecular transitions because of the absence of CO and other trace molecules, and because H_2 emits so weakly at temperatures 10 K < T < 100 K typical of this molecular component. This component has been indirectly observed through other tracers of mass such as gamma rays produced in cosmic ray collisions with the gas and far-infrared/submillimeter wavelength dust continuum radiation. In this paper we theoretically model this dark mass and find that the fraction of the molecular mass in this dark component is remarkably constant (~ 0.3 for average visual extinction through the cloud with mean A_V ~ 8) and insensitive to the incident ultraviolet radiation field strength, the internal density distribution, and the mass of the molecular cloud as long as mean A_V, or equivalently, the product of the average hydrogen nucleus column and the metallicity through the cloud, is constant. We also find that the dark mass fraction increases with decreasing mean A_V, since relatively more molecular H_2 material lies outside the CO region in this case.

Although they are only minor constituents of the interstellar medium, halogen-containing molecules are of special interest because of their unique thermochemistry. Here, we present a theoretical study of the chemistry of interstellar molecules containing the halogen elements chlorine and fluorine. We have modeled both diffuse and dense molecular clouds, making use of updated estimates for the rates of several key chemical processes. We present predictions for the abundances of the three halogen molecules that have been detected to date in the interstellar medium: HF, CF+ and HCl. As in our previous study of fluorine-bearing interstellar molecules, we predict HF to be the dominant gas-phase reservoir of fluorine within both diffuse and dense molecular clouds; we expect the Herschel Space Observatory to detect widespread absorption in the HF J=1-0 transition. Our updated model now overpredicts the CF+ abundance by a factor ~10 relative to observations of the Orion Bar; this discrepancy has widened because we now adopt a laboratory measurement of the CF+ dissociative recombination rate that is smaller than the estimate we adopted previously. This disagreement suggests that the reaction of C+ with HF proceeds more slowly than the capture rate assumed in our model; a laboratory measurement of this reaction rate would be very desirable. Our model predicts diffuse cloud HCl abundances that are similar to those predicted previously and detected tentatively toward zeta Oph. Two additional species are potentially detectable from photodissociation regions: the H2Cl+ and HCl+ molecular ions. Ortho-H2Cl+ has its lowest-lying transition in the millimeter spectral region observable from the ground, and the lowest rotational transition of HCl+ is observable with Herschel's HIFI instrument.

We present a mid-infrared spectroscopic data cube of the central part of 30 Doradus, observed with Spitzer's IRS and MIPS/SED mode. Aromatic dust emission features and emission lines from molecular and atomic hydrogen are detected but not particularly strong. The dominant spectral features are emission lines from moderately ionized species of argon, neon, and sulphur, which are used to determine the physical conditions in the ionized gas. The ionized gas excitation shows strong variations on parsec scales, some of which can plausibly be associated with individual hot stars. We fit the ionic line strengths with photoionization and shock models, and find that photoionization dominates in the region. The ionization parameter U traces the rim of the central bubble, as well as highlighting isolated sources of ionization, and at least one quiescent clump. The hardness of the ionizing radiation field T_rad reveals several "hot spots" that are either the result of individual very hot stars or trace the propagation of the diffuse ionizing field through the surrounding neutral cloud. Consistent with other measurements of giant molecular hydrogen regions, log(U) ranges between -3 and -0.75, and T_rad between 30000 and 85000K.

We use observations of the CI, CII, HI, and H_2 column densities along lines of sight in the Galactic plane to determine the formation rate of H_2 on grains and to determine chemical reaction rates with Polycyclic Aromatic Hydrocarbons. Photodissociation region models are used to find the best fit parameters to the observed columns. We find the H_2 formation rate on grains has a low rate (R ~ 1 x 10^(-17) cm^(3) s^(-1)) along lines of sight with low column density (A_V < 0.25) and low molecular fraction (f_(H_2) < 10^(-4)). At higher column densities (0.25 < A_V <2.13), we find a rate of R ~ 3.5x10^(-17) cm^(3) s^(-1). The lower rate at low column densities could be the result of grain processing by interstellar shocks which may deplete the grain surface area or process the sites of H +H formation, thereby inhibiting H_2 production. Alternatively, the formation rate may be normal, and the low molecular fraction may be the result of lines of sight which graze larger clouds. Such lines of sight would have a reduced H_2 self-shielding compared to the line-of-sight column. We find the reaction C^+ +PAH^- --> C + PAH^0 is best fit with a rate 2.4 x 10^(-7) \Phi_PAH T_2^(-0.5) cm^(3) s^(-1) with T_2= T/100 K and the reaction C^+ + PAH^0 --> C + PAH^+ is best fit with a rate 8.8x 10^(-9)\Phi_PAH cm^(3) s^(-1). In high column density gas we find \Phi_PAH ~ 0.4. In low column density gas, \Phi_PAH is less well constrained with \Phi_PAH ~ 0.2 - 0.4.

We have mapped the warm molecular gas traced by the H_2 S(0) - H_2 S(5) pure rotational mid-infrared emission lines over a radial strip across the nucleus and disk of M51 (NGC 5194) using the Infrared Spectrograph (IRS) on the Spitzer Space Telescope. The six H_2 lines have markedly different emission distributions. We obtained the H_2 temperature and surface density distributions by assuming a two temperature model: a warm (T = 100 - 300 K) phase traced by the low J (S(0) - S(2)) lines and a hot phase (T = 400 - 1000 K) traced by the high J (S(2) - S(5)) lines. The lowest molecular gas temperatures are found within the spiral arms (T ~ 155 K), while the highest temperatures are found in the inter-arm regions (T > 700 K). The warm gas surface density reaches a maximum of 11 M_sun/pc^2 in the northwestern spiral arm, whereas the hot gas surface density peaks at 0.24 M_sun/pc^2 at the nucleus. The spatial offset between the peaks in the warm and hot phases and the differences in the distributions of the H_2 line emission suggest that the warm phase is mostly produced by UV photons in star forming regions while the hot phase is mostly produced by shocks or X-rays associated with nuclear activity. The warm H_2 is found in the dust lanes of M51, spatially offset from the brightest HII regions. The warm H_2 is generally spatially coincident with the cold molecular gas traced by CO (J = 1 - 0) emission, consistent with excitation of the warm phase in dense photodissociation regions (PDRs). In contrast, the hot H_2 is most prominent in the nuclear region. Here, over a 0.5 kpc radius around the nucleus of M51, the hot H_2 coincides with [O IV](25.89 micron) and X-ray emission indicating that shocks and/or X-rays are responsible for exciting this phase.

We present new Spitzer photometry of the Eagle Nebula (M16, containing the optical cluster NGC 6611) combined with near-infrared photometry from 2MASS. We use dust radiative transfer models, mid-infrared and near-infrared color-color analysis, and mid-infrared spectral indices to analyze point source spectral energy distributions, select candidate young stellar objects (YSOs), and constrain their mass and evolutionary state. Comparison of the different protostellar selection methods shows that mid-infrared methods are consistent, but as has been known for some time, near-infrared-only analysis misses some young objects. We reveal more than 400 protostellar candidates, including one massive young stellar object (YSO) that has not been previously highlighted. The YSO distribution supports a picture of distributed low-level star formation, with no strong evidence of triggered star formation in the ``pillars''. We confirm the youth of NGC 6611 by a large fraction of infrared-excess sources, and reveal a younger cluster of YSOs in the nearby molecular cloud. Analysis of the YSO clustering properties shows a possible imprint of the molecular cloud's Jeans length. Multiwavelength mid-IR imaging thus allows us to analyze the protostellar population, to measure the dust temperature and column density, and to relate these in a consistent picture of star formation in M16.

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.

We have identified 230 Tycho-2 Spectral Catalog stars that exhibit 8 micron mid-infrared extraphotospheric excesses in the MidCourse Space Experiment (MSX) and Spitzer Space Telescope Galactic Legacy MidPlane Survey Extraordinaire (GLIMPSE) surveys. Of these, 183 are either OB stars earlier than B8 in which the excess plausibly arises from a thermal bremsstrahlung component or evolved stars in which the excess may be explained by an atmospheric dust component. The remaining 47 stars have spectral classifications B8 or later and appear to be main sequence or late pre-main-sequence objects harboring circumstellar disks. Six of the 47 stars exhibit multiple signatures characteristic of pre-main-sequence circumstellar disks, including emission lines, near-infrared K-band excesses, and X-ray emission. Approximately one-third of the remaining 41 sources have emission lines suggesting relative youth. Of the 25 GLIMPSE stars with SST data at >24 microns, 20 also show an excess at 24 microns. Three additional objects have 24 micron upper limits consistent with possible excesses, and two objects have photospheric measurements at 24 microns. Six MSX sources had a measurement at wavelengths >8 microns. We modeled the excesses in 26 stars having two or more measurements in excess of the expected photospheres as single-component blackbodies. We determine probable disk temperatures and fractional infrared luminosities in the range 191 < T < 787 and 3.9x10^-4 < L_IR/L_* < 2.7x10^-1. We estimate a lower limit on the fraction of Tycho-2 Spectral Catalog main-sequence stars having mid-IR, but not near-IR, excesses to be 1.0+-0.3%.

We investigate the diffuse absolute calibration of the InfraRed Array Camera on the Spitzer Space Telescope at 8.0microns using a sample of 43 HII regions with a wide range of morphologies near GLON=312deg. For each region we carefully measure sky-subtracted,point-source- subtracted, areally-integrated IRAC 8.0-micron fluxes and compare these with Midcourse Space eXperiment (MSX) 8.3-micron images at two different spatial resolutions, and with radio continuum maps. We determine an accurate median ratio of IRAC 8.0-micron/MSX\8.3-micron fluxes, of 1.55+/-0.15. From robust spectral energy distributions of these regions we conclude that the present 8.0-micron diffuse calibration of the SST is 36% too high compared with the MSX validated calibration, perhaps due to scattered light inside the camera. This is an independent confirmation of the result derived for the diffuse calibration of IRAC by the Spitzer Science Center (SSC). From regression analyses we find that 843-MHz radio fluxes of HII regions and mid-infrared (MIR) fluxes are linearly related for MSX at 8.3-microns and Spitzer at 8.0 microns, confirming the earlier MSX result by Cohen & Green. The median ratio of MIR/843-MHz diffuse continuum fluxes is 600 times smaller in nonthermal than thermal regions, making it a sharp discriminant. The ratios are largely independent of morphology up to a size of ~24 arcsec. We provide homogeneous radio and MIR morphologies for all sources. MIR morphology is not uniquely related to radio structure. Compact regions may have MIR filaments and/or diffuse haloes, perhaps infrared counter- parts to weakly ionized radio haloes found around compact HII regions. We offer two IRAC colour-colour plots as quantitative diagnostics of diffuse HII regions.

We discuss the first astronomical detection of the CF+ (fluoromethylidynium) ion, obtained by observations of the J=1-0 (102.6 GHz), J=2-1 (205.2 GHz) and J=3-2 (307.7 GHz) rotational transitions toward the Orion Bar region. Our search for CF+, carried out using the IRAM 30m and APEX 12m telescopes, was motivated by recent theoretical models that predict CF+ abundances of a few times 1.E-10 in UV-irradiated molecular regions where C+ is present. The CF+ ion is produced by exothermic reactions of C+ with HF. Because fluorine atoms can react exothermically with H2, HF is predicted to be the dominant reservoir of fluorine, not only in well-shielded regions but also in the surface layers of molecular clouds where the C+ abundance is large. The observed CF+ line intensities imply the presence of CF+ column densities of at least 1.E+12 cm-2 over a region of size at least ~ 1 arcmin, in good agreement with theoretical predictions. They provide support for our current theories of interstellar fluorine chemistry, which suggest that hydrogen fluoride should be ubiquitous in interstellar gas clouds and widely detectable in absorption by future satellite and airborne observatories.

Recent Chandra and XMM observations of distant quasars have shown strong local (z~0) X-ray absorption lines from highly ionized gas, primarily He-like oxygen. The nature of these X-ray absorbers, i.e., whether they are part of the hot gas associated with the Milky Way or part of the intragroup medium in the Local Group, remains a puzzle due to the uncertainties in the distance. We present in this paper a survey of 20 AGNs with Chandra and XMM archival data. About 40% of the targets show local OVII He-alpha absorption with column densities around 10^16 cm^-2; in particular, OVII absorption is present in all the high quality spectra. We estimate that the sky covering fraction of this OVII-absorbing gas is at least 63%, at 90% confidence, and likely to be unity given enough high-quality spectra. Based on (1) the expected number of absorbers along sight lines toward distant AGNs, (2) joint analysis with X-ray emission measurements, and (3) mass estimation, we argue that the observed X-ray absorbers are part of the hot gas associated with our Galaxy. Future observations will significantly improve our understanding of the covering fraction and provide robust test of this result.

We report the first astronomical detection of the CF+ (fluoromethylidynium) ion obtained by recent observations of its J = 1 - 0 (102.6 GHz), J = 2 - 1 (205.2 GHz), and J = 3 - 2 (307.7 GHz) pure rotational emissions toward the Orion Bar. Our search for CF+, carried out using the IRAM 30m and APEX 12m telescopes, was motivated by recent theoretical models that predict CF+ abundances of a few x E-10 in UV-irradiated molecular regions where C+ is present. The measurements confirm the predictions. They provide support for our current theories of interstellar fluorine chemistry, which suggest that hydrogen fluoride should be ubiquitous in interstellar gas clouds.

The GLIMPSE (Galactic Legacy Mid-Plane Survey Extraordinaire) Point Source Catalog of ~ 30 million mid-infrared sources towards the inner Galaxy, 10 < |l| < 65 degrees and |b| < 1 degree, was used to determine the distribution of stars in Galactic longitude, latitude, and apparent magnitude. The counts versus longitude can be approximated by the modified Bessel function N=N_0*(l/l_0)*K_1(l/l_0), where l_0 is insensitive to limiting magnitude, band choice, and side of Galactic center: l_0= 17-30 degrees with a best fit value in the the 4.5 micron band of l_0=24 +/- 4 degrees. Modeling the source distribution as an exponential disk yields a radial scale length of H= 3.9 +/- 0.6 kpc. There is a pronounced north-south asymmetry in source counts for |l| < 30 degrees, with ~ 25% more stars in the north. For l=10-30 degrees, there is a strong enhancement of stars of m= 11.5-13.5 mag. A linear bar passing through the Galactic center with half-length R_bar=4.4 +/- 0.5 kpc, tilted by phi=44 +/- 10 degrees to the Sun-Galactic Center line, provides the simplest interpretation of this data. We examine the possibility that enhanced source counts at l=26-28 degrees, 31.5-34 degrees, and 306-309 degrees are related to Galactic spiral structure. Total source counts are depressed in regions where the counts of red objects (m_K-m_[8.0] >3) peak. In these areas, the counts are reduced by extinction due to molecular gas and/or high diffuse backgrounds associated with star formation.

Spitzer IRAC 3.6-8 micron photometry obtained as part of the GLIMPSE survey has revealed mid-infrared excesses for 33 field stars with known spectral types in a 1.2 sq. degree field centered on the southern Galactic HII region RCW49. These stars comprise a subset of 184 stars with known spectral classification, most of which were pre-selected to have unusually red IR colors. We propose that the mid-IR excesses are caused by circumstellar dust disks that are either very late remnants of stellar formation or debris disks generated by planet formation. Of these 33 stars, 29 appear to be main-sequence stars based on optical spectral classifications. Five of the 29 main-sequence stars are O or B stars with excesses that can be plausibly explained by thermal bremsstrahlung emission, and four are post main-sequence stars. The lone O star is an O4V((f)) at a spectrophotometric distance of 3233+ 540- 535 pc and may be the earliest member of the Westerlund 2 cluster. Of the remaining 24 main-sequence stars, 18 have SEDs that are consistent with hot dusty debris disks, a possible signature of planet formation. Modeling the excesses as blackbodies demonstrates that the blackbody components have fractional bolometric disk-to-star luminosity ratios, LIR/L*, ranging from 10^-3 to 10^-2 with temperatures ranging from 220 to 820 K. The inferred temperatures are more consistent with asteroid belts rather than the cooler temperatures expected for Kuiper belts. Mid-IR excesses are found in all spectral types from late B to early K.

We present a theoretical investigation of the chemistry of fluorine-bearing molecules in diffuse and dense interstellar gas clouds. The chemistry of interstellar fluorine is qualitatively different from that of any other element, because - unlike the neutral atoms of any other element found in diffuse or dense molecular clouds - atomic fluorine undergoes an exothermic reaction with molecular hydrogen. Over a wide range of conditions attained within interstellar gas clouds, the product of that reaction - hydrogen fluoride - is predicted to be the dominant gas-phase reservoir of interstellar fluorine nuclei. Our model predicts HF column densities ~ 1.E+13 cm-2 in dark clouds and column densities as large as 1.E-11 cm-2 in diffuse interstellar gas clouds with total visual extinctions as small as 0.1 mag. Such diffuse clouds will be detectable by means of absorption line spectroscopy of the J = 1 - 0 transition at 243.2 micron using the Stratospheric Observatory for Infrared Astronomy (SOFIA) and the Herschel Space Observatory (HSO). The CF+ ion is predicted to be the second most abundant fluorine-bearing molecule, with typical column densities a factor ~ 100 below those of HF; with its lowest two rotational transitions in the millimeter-wave spectral region, CF+ may be detectable from ground-based observatories. HF absorption in quasar spectra is a potential probe of molecular gas at high redshift, providing a possible bridge between the UV/optical observations capable of probing H2 in low column density systems and the radio/millimeter-wavelength observations that probe intervening molecular clouds of high extinction and large molecular fraction; at redshifts beyond ~ 0.3, HF is potentially detectable from ground-based submillimeter observatories in several atmospheric transmission windows.

We present a new planetary nebula, first identified in images from the Australia Telescope Compact Array, although not recognized at that time. Recent observations with the Spitzer Space Telescope during the GLIMPSE Legacy program have rediscovered the object. The high-resolution radio and infrared images enable the identification of the central star or its wind, the recognition of the radio emission as thermal, and the probable presence of polycylic aromatic hydrocarbons in and around the source. These lead to the conclusion that G313.3+00.3 is a planetary nebula. This object is of particular interest because it was discovered solely through radio and mid-infrared imaging, without any optical (or near-infrared) confirmation, and acts as a proof of concept for the discovery of many more highly extinguished planetary nebulae. G313.3+00.3 is well-resolved by both the instruments with which it was identified, and suffers extreme reddening due to its location in the Scutum-Crux spiral arm.

We determine and tabulate A(lambda)/A(K), the wavelength dependence of interstellar extinction, in the Galactic plane for 1.25um<lambda<8.0um along two lines of sight: l~42 and l~284. The first is a relatively quiescent and unremarkable region; the second contains the giant HII region RCW 49 as well as a ``field'' region unrelated to the cluster and nebulosity. Areas near these Galactic longitudes were imaged at J, H, and K bands by 2MASS and at 3--8um by Spitzer for the GLIMPSE Legacy program. We measure the mean values of the color excess ratios (A(lambda)-A(K))/(A(J)-A(K)) directly from the color distributions of observed stars. The extinction ratio between two of the filters, e.g. A(J)/A(K), is required to calculate A(lambda)/A(K) from those measured ratios. We use the apparent JHK magnitudes of giant stars along our two sightlines, and fit the reddening as a function of magnitude (distance) to determine A(J)/kpc, A(K)/kpc, and A(J)/A(K). Our values of A(lambda)/A(K) show a flattening across the 3--8um wavelength range, roughly consistent with the Lutz 1996 extinction measurements derived for the sightline toward the Galactic center.

Spitzer Space Telescope imaging as part of the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) reveals a previously unidentified low-latitude rich star cluster near l=31.3 degrees, b=-0.1 degrees. Near-infrared JHK' photometry from the Wyoming Infrared Observatory indicates an extinction of A_V ~ 15+/-3 mag for cluster members. Analysis of 13CO features along the same sightline suggests a probable kinematic distance of 3.1 -- 5.2 kpc. The new cluster has an angular diameter of ~1-2 pc, a total magnitude m_{K_0}=2.1 corrected for extinction, and a luminosity of M_K ~ -10.3 at 3.1 kpc. In contrast to young massive Galactic clusters with ages <100 Myr, the new cluster has no significant radio emission. Comparison to theoretical K-band luminosity functions indicates an age of at least several Gyr and a mass of at least 10^5 solar masses. Unlike known old open clusters, this new cluster lies in the inner Galaxy at R_{GC} ~ 6.1 kpc. We designate this object ``GLIMPSE-C01'' and classify it as a Milky Way globular cluster passing through the Galactic disk. We also identify a region of star formation and fan-shaped outflows from young stellar objects in the same field as the cluster. The cluster's passage through the Galactic molecular layer may have triggered this star formation activity.

The nearby spiral galaxy NGC 7331 was spectrally mapped from 5-38um using all modules of Spitzer's IRS spectrograph. A strong new dust emission feature, presumed due to PAHs, was discovered at 17.1um. The feature's intensity is nearly half that of the ubiquitous 11.3um band. The 7-14um spectral maps revealed significant variation in the 7.7 and 11.3um PAH features between the stellar ring and nucleus. Weak [OIV] 25.9um line emission was found to be centrally concentrated in the nucleus, with an observed strength over 10% of the combined neon line flux, indicating an AGN or unusually active massive star photo-ionization. Two [SIII] lines fix the characteristic electron density in the HII regions at n_e < ~200 cm^-3. Three detected H_2 rotational lines, tracing warm molecular gas, together with the observed IR continuum, are difficult to match with standard PDR models. Either additional PDR heating or shocks are required to simultaneously match lines and continuum.

GLIMPSE imaging using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope indicates that star formation is ongoing in the RCW 49 giant H II region. A photometric comparison of the sources in RCW 49 to a similar area to its north finds that at least 300 stars brighter than 13th magnitude in band [3.6] have infrared excesses inconsistent with reddening due to foreground extinction. These are likely young stellar objects (YSOs) more massive than 2.5 M_sun suggesting that thousands more low-mass stars are forming in this cloud. Some of the YSOs are massive (B stars) and therefore very young, suggesting that a new generation of star formation is occurring, possibly triggered by stellar winds and shocks generated by the older (2-3 Myr) central massive cluster. The Spitzer IRAC camera has proven to be ideally suited for distinguishing young stars from field stars, and the GLIMPSE survey of the Galactic Plane will likely find thousands of new star formation regions.

GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire), a SIRTF Legacy Science Program, will be a fully sampled, confusion-limited infrared survey of the inner two-thirds of the Galactic disk with a pixel resolution of \~1.2" using the Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 microns. The survey will cover Galactic latitudes |b| <1 degree and longitudes |l|=10 to 65 degrees (both sides of the Galactic center). The survey area contains the outer ends of the Galactic bar, the Galactic molecular ring, and the inner spiral arms. The GLIMPSE team will process these data to produce a point source catalog, a point source data archive, and a set of mosaicked images. We summarize our observing strategy, give details of our data products, and summarize some of the principal science questions that will be addressed using GLIMPSE data. Up-to-date documentation, survey progress, and information on complementary datasets are available on the GLIMPSE web site: www.astro.wisc.edu/glimpse.

We present a study of the [CII] 157.74 um fine-structure line in a sample of 15 ultraluminous infrared (IR) galaxies (L_IR>10^12 Lsun; ULIRGs) using the Long Wavelength Spectrometer (LWS) on the Infrared Space Observatory (ISO). We confirm the observed order of magnitude deficit (compared to normal and starburst galaxies) in the strength of the [CII] line relative to the far-IR dust continuum emission found in our initial report (Luhman et al. 1998), but here with a sample that is twice as large. This result suggests that the deficit is a general phenomenon affecting 4/5 ULIRGs. We present an analysis using observations of generally acknowledged photodissociation region (PDR) tracers ([CII], [OI] 63 and 145 um, and FIR continuum emission), which suggests that a high UV flux G_o incident on a moderate density n PDR could explain the deficit. However, comparisons with other ULIRG observations, including CO (1-0), [CI] (1-0), and 6.2 um polycyclic aromatic hydrocarbon (PAH) emission, suggest that high G_o/n PDRs alone cannot produce a self-consistent solution that is compatible with all of the observations. We propose that non-PDR contributions to the FIR continuum can explain the apparent [CII] deficiency. Here, unusually high G_o and/or n physical conditions in ULIRGs as compared to those in normal and starburst galaxies are not required to explain the [CII] deficit. Dust-bounded photoionization regions, which generate much of the FIR emission but do not contribute significant [CII] emission, offer one possible physical origin for this additional non-PDR component. Such environments may also contribute to the observed suppression of FIR fine-structure emission from ionized gas and PAHs, as well as the warmer FIR colors found in ULIRGs. The implications for observations at higher redshifts are also revisited.

Much of the interstellar medium in disk galaxies is in the form of neutral atomic hydrogen, H I. This gas can be in thermal equilibrium at relatively low temperatures, T < 300 K (the cold neutral medium, or CNM) or at temperatures somewhat less than 10^4 K (the warm neutral medium, or WNM). These two phases can coexist over a narrow range of pressures, P_min < P < P_max. We determine P_min and P_max in the plane of the Galaxy as a function of Galactocentric radius R using recent determinations of the gas heating rate and the gas phase abundances of interstellar gas. We provide an analytic approximation for P_min as a function of metallicity, far-ultraviolet radiation field, and the ionization rate of atomic hydrogen. Over most of the disk of the Galaxy, the H I must be in two phases: the weight of the H I in the gravitational potential of the Galaxy is large enough to generate thermal pressures exceeding P_min, so that turbulent pressure fluctuations can produce cold gas that is thermally stable; and the mean density of the H I is too low for the gas to be all CNM. Our models predict the presence of CNM gas to R = 16-18 kpc, somewhat farther than previous estimates. We also examine the potential impact of turbulent heating on our results and provide expressions for the heating rate as a function of Galactic radius.

We calculate the coupled hydrostatic and ionization structures of pressure-supported gas clouds that are confined by gravitationally dominant dark-matter (DM) mini-halos and by an external bounding pressure provided by a hot medium. We focus on clouds that are photoionized and heated by the present-day background metagalactic field and determine the conditions for the formation of warm (WNM), and multi-phased (CNM/WNM) neutral atomic hydrogen (HI) cores in the DM-dominated clouds. We consider LCDM dark-matter halos, and we compute models for a wide range of halo masses, total cloud gas masses, and external bounding pressures. We present models for the pressure-supported HI structures observed in the Local Group dwarf galaxies Leo A and Sag DIG. We then construct minihalo models for the multi-phased (and low-metallicity) compact high-velocity HI clouds (CHVCs). If the CHVCs are drawn from the same family of halos that successfully reproduce the dwarf galaxy observations, then the CHVCs must be "circumgalactic objects" with characteristic distances of 150 kpc. For such systems we find that multi-phased behavior occurs for peak WNM HI column densities between 2e19 and 1e20 cm^-2, consistent with observations. If the large population of CHVCs represent "missing low-mass satellites" of the Galaxy, then these clouds must be pressure-confined to keep the gas neutral. For an implied CHVC minihalo scale velocity of v_s=12 km s^-1, the confining pressure must exceed ~50 cm^-3 K. A hot (2e6 K) Galactic corona could provide the required pressure at 150 kpc.

We present the first full FIR spectrum of Centaurus A (NGC 5128) from 43 - 196.7 um. The data was obtained with the ISO Long Wavelength Spectrometer (LWS). We conclude that the FIR emission in a 70 arcsec beam centred on the nucleus is dominated by star formation rather than AGN activity. The flux in the far-infrared lines is ~ 1 % of the total FIR: the [CII] line flux is ~ 0.4 % FIR and the [OI] line is ~ 0.2 %, with the remainder arising from [OIII], [NII] and [NIII] lines. These are typical values for starburst galaxies. The ratio of the [NIII]/[NII] line intensities from the HII regions in the dust lane can be modelled as a ~ 6 million year old starburst. This suggests that the galaxy underwent either a recent merger or a merger which triggered a series of bursts. We estimate that < 5 % of the observed [CII] arises in the cold neutral medium (CNM) and that ~ 10 % arises in the warm ionized medium (WIM). The main contributors to the [CII] emission are the PDRs, which are located throughout the dust lane and in regions beyond where the bulk of the molecular material lies. On scales of ~ 1 kpc the average physical properties of the PDRs are modelled with a gas density, n ~ 1000 cm^-3, an incident far-UV field, G ~ 100 times the local Galactic field, and a gas temperature of ~ 250 K.

Photodissociation Region (PDR) models are computed over a wide range of physical conditions, from those appropriate to giant molecular clouds illuminated by the interstellar radiation field to the conditions experienced by circumstellar disks very close to hot massive stars. These models use the most up-to-date values of atomic and molecular data, the most current chemical rate coefficients, and the newest grain photoelectric heating rates which include treatments of small grains and large molecules. In addition, we examine the effects of metallicity and cloud extinction on the predicted line intensities. Results are presented for PDR models with densities over the range n=10^1-10^7 cm^-3 and for incident far-ultraviolet radiation fields over the range G_0=10^-0.5-10^6.5, for metallicities Z=1 and 0.1 times the local Galactic value, and for a range of PDR cloud sizes. We present line strength and/or line ratio plots for a variety of useful PDR diagnostics: [C II] 158 micron, [O I] 63 and 145 micron, [C I] 370 and 609 micron, CO J=1-0, J=2-1, J=3-2, J=6-5 and J=15-14, as well as the strength of the far-infrared continuum. These plots will be useful for the interpretation of Galactic and extragalactic far infrared and submillimeter spectra observable with ISO, SOFIA, SWAS, FIRST and other orbital and suborbital platforms. As examples, we apply our results to ISO and ground based observations of M82, NGC 278, and the Large Magellenic Cloud.