New ALMA CO, HCO+, and 100 GHz continuum observations and Gemini NIFS K-band spectra are combined with previous radio to X-ray data to study Feature i, a highly luminous starburst clump on an outer arm of the interacting galaxy NGC 2207. This clump has an optically-opaque dust cone extending out of its 170 pc core. The measured CO accounts for the dark cone extinction if almost all the gas and dust there is in front of the star clusters. An associated approaching CO outflow has v_z ~ 16 km/s, an estimated molecular mass 8 x 10^6 Msun, and rises to heights ~ 0.9 kpc. A receding CO outflow on the far side with v_z ~ 28 km/s is less extensive. The observed star formation in the core over 10 Myr can supply the dark cone kinetic energy via supernovae and stellar winds. Other signs of intense activity are variable radio continuum suggesting an embedded supernova or other outburst, X-ray emission, and Br-gamma and He I lines with 82 km/s line widths. According to previous models, the retrograde encounter suffered by NGC 2207 caused loss of angular momentum. This compressed its outer disk. We suggest that the resulting inward crashing gas stream hit a massive HI cloud as it was being compressed by a spiral density wave and triggered the observed starburst.

The depth-dependent abundance of both gas-phase and solid-state water within dense, quiescent, molecular clouds is important to both the cloud chemistry and gas cooling. Where water is in the gas phase, it's free to participate in the network of ion-neutral reactions that lead to a host of oxygen-bearing molecules, and its many ortho and para energy levels make it an effective coolant for gas temperatures greater than 20K. Where water is abundant as ice on grain surfaces, and unavailable to cool the gas, significant amounts of oxygen are removed from the gas phase, suppressing the gas-phase chemical reactions that lead to a number of oxygen-bearing species, including O2. Models of FUV-illuminated clouds predict that the gas-phase water abundance peaks in the range Av ~3 and 8mag of the cloud surface, depending on the gas density and FUV field strength. Deeper within such clouds, water is predicted to exist mainly as ice on grain surfaces. More broadly, these models are used to analyze a variety of other regions, including outflow cavities associated with young stellar objects and the surface layers of protoplanetary disks. In this paper, we report the results of observational tests of FUV-illuminated cloud models toward the Orion Molecular Ridge and Cepheus B using data obtained from the Herschel Space Observatory and the Five College Radio Astronomy Observatory. Toward Orion, 2220 spatial positions were observed along the face-on Orion Ridge in the H2O 110-101 557GHz and NH3 J,K=1,0-0,0 572GHz lines. Toward Cepheus B, two strip scans were made in the same lines across the edge-on ionization front. These new observations demonstrate that gas-phase water exists primarily within a few magnitudes of dense cloud surfaces, strengthening the conclusions of an earlier study based on a much smaller data set, and indirectly supports the prediction that water ice is quite abundant in dense clouds.

The current paradigm of Galactic Center (GC) gas motions and star formation envisions sequential star formation in streams of gas as they pass near the supermassive black hole, Sgr A*. This is based on the relative positions of dense molecular clouds, the very young star-forming region Sgr B2, the much older region Sgr C, and the several Myr old Arches and Quintuplet Clusters. Because Sgr B1 is found with Sgr B2 in a common envelope of molecular gas and far-infrared emission, the two sources are thought to be physically related, even though there are indicators of a significantly greater age for Sgr B1. To clarify the status of Sgr B1, we have mapped it with the FIFI-LS spectrometer on SOFIA in the far-infrared lines of [O III] 52 and 88 micron. From the ratios of these lines and lines measured with the Spitzer Infrared Spectrograph, we find that there are at least eight separate sub-regions that must contain the stars that excite the gas. We infer spectral energy distributions (SEDs) of the ionizing sources from models and find they are in agreement only with SEDs of late O stars augmented at the highest frequencies with interstellar X-rays from fast shocks. We suggest that although the gas, from its velocity structure, must be part of the very young Sgr B2 complex, the stars that are ionizing the gas were not formed there but are the remnants of a previous generation of star formation in the GC.

Recent observations from Herschel allow the identification of important mechanisms responsible for the heating of gas surrounding low-mass protostars and its subsequent cooling in the far-infrared (FIR). Shocks are routinely invoked to reproduce some properties of the far-IR spectra, but standard models fail to reproduce the emission from key molecules, e.g. H$_2$O. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars (Class 0/I). The Herschel-PACS spectral maps covering $\sim55-210$ $\mu$m with a field-of-view of $\sim$50'' are used to quantify the gas excitation conditions and spatial extent using rotational transitions of H$_{2}$O, high-$J$ CO, and OH, as well as [O I] and [C II]. We confirm that a warm ($\sim$300 K) CO reservoir is ubiquitous and that a hotter component ($760\pm170$ K) is frequently detected around protostars. The line emission is extended beyond $\sim$1000 AU spatial scales in 40/90 objects, typically in molecular tracers in Class 0 and atomic tracers in Class I objects. High-velocity emission ($\gtrsim90$ km s$^{-1}$) is detected in only 10 sources in the [O I] line, suggesting that the bulk of [O I] arises from gas that is moving slower than typical jets. Line flux ratios show an excellent agreement with models of $C$-shocks illuminated by UV photons for pre-shock densities of $\sim$$10^5$ cm$^{-3}$ and UV fields 0.1-10 times the interstellar value. The far-IR molecular and atomic lines are a unique diagnostic of feedback from UV emission and shocks in envelopes of deeply embedded protostars.

Far-infrared spectroscopy reveals gas cooling and its underlying heating due to physical processes taking place in the surroundings of protostars. These processes are reflected in both the chemistry and excitation of abundant molecular species. Here, we present the Herschel-PACS far-IR spectroscopy of 90 embedded low-mass protostars from the WISH (van Dishoeck et al. 2011), DIGIT (Green et al. 2013), and WILL surveys (Mottram et al. 2017). The $5\times5$ spectra covering the $\sim50''\times50''$ field-of-view include rotational transitions of CO, H$_2$O, and OH lines, as well as fine-structure [O I] and [C II] in the $\sim$50-200 $\mu$m range. The CO rotational temperatures (for $J_\mathrm{u}\geq14)$ are typically $\sim$300 K, with some sources showing additional components with temperatures as high as $\sim$1000 K. The H$_2$O / CO and H$_2$O / OH flux ratios are low compared to stationary shock models, suggesting that UV photons may dissociate some H$_2$O and decrease its abundance. Comparison to C shock models illuminated by UV photons shows good agreement between the line emission and the models for pre-shock densities of $10^5$ cm$^{-3}$ and UV fields 0.1-10 times the interstellar value. The far-infrared molecular and atomic lines are the unique diagnostic of shocks and UV fields in deeply-embedded sources.

We present an analysis of archival ISO observations of pure-rotational lines of H2 in three pointings in the central 3 parsecs of the Galaxy: toward the Southwest region and Northeast region of the Galactic center Circumnuclear Disk, and toward the supermassive black hole Sgr A*. We detect pure rotational lines from 0-0 S(0) to S(13), as well as a number of rovibrationally excited transitions. From the pure rotational lines, we are able to describe the molecular gas with three discrete temperature components: a `hot' component between 500-600 K, a `hotter' component at 1250-1350 K, and a `hottest' component at > 2600 K. Toward Sgr A*, likely due to a combination of poorer baselines and weaker emission, we only detect a single hot component, at 1100 K. The observed excitation is consistent with heating via C-shocks. We also fit a continuous temperature distribution to the the S(1) through S(7) lines by assuming a power-law distribution of temperatures. We measure power law indices of n = 3.22 for the Northeast region and n = 2.83 for the Southwest region, with a smaller index indicating a higher fraction of warm gas. These indices are lower than those measured for other galaxies or other Galactic center clouds, which are measured to have n= 4-6. If we extrapolate this temperature distribution down to a cutoff temperature of 50 K, then the total molecular gas mass that we would measure for the Southwest (Northeast) region is 32 (140) % of the total molecular gas mass inferred from the dust emission, and 26 (125) % of the total molecular gas mass inferred from the CO emission for these regions. Ultimately, we find that the disagreement in the amount of mass recovered for these two regions means that this method cannot yet be verified to yield a reliable and independent estimate of the mass in the local region of the Circumnuclear Disk.

We report the discovery of water maser emission at frequencies above 1 THz. Using the GREAT instrument on SOFIA, we have detected emission in the 1.296411 THz 8(27)-7(34) transition of water toward three oxygen-rich evolved stars: W Hya, U Her, and VY CMa. An upper limit on the 1.296 THz line flux was obtained toward R Aql. Near-simultaneous observations of the 22.23508 GHz 6(16)-5(23) water maser transition were carried out towards all four sources using the Effelsberg 100m telescope. The measured line fluxes imply 22 GHz / 1.296 THz photon luminosity ratios of 0.012, 0.12, and 0.83 respectively for W Hya, U Her, and VY CMa, values that confirm the 22 GHz maser transition to be unsaturated in W Hya and U Her. We also detected the 1.884888 THz 8(45)-7(53) transition toward W Hya and VY CMa, and the 1.278266 THz 7(43)-6(52) transition toward VY CMa. Like the 22 GHz maser transition, all three of the THz emission lines detected here originate from the ortho-H2O spin isomer. Based upon a model for the circumstellar envelope of W Hya, we estimate that stimulated emission is responsible for ~ 85% of the observed 1.296 THz line emission, and thus that this transition may be properly described as a terahertz-frequency maser. In the case of the 1.885 THz transition, by contrast, our W Hya model indicates that the observed emission is dominated by spontaneous radiative decay, even though a population inversion exists.

ALMA observations of CO(1-0) emission in the interacting galaxies IC 2163 and NGC 2207 are used to determine the properties of molecular clouds and their association with star-forming regions observed with the Hubble Space Telescope. Half of the CO mass is in 249 clouds each more massive than 4.0x10^5Mo. The mass distribution functions for the CO clouds and star complexes in a galactic-scale shock front in IC 2163 both have a slope on a log-log plot of -0.7, similar to what is observed in Milky Way clouds. The molecular cloud mass function is steeper in NGC 2207. The CO distribution in NGC 2207 also includes a nuclear ring, a mini-bar, and a mini-starburst region that dominates the 24micron, radio, and Halpha emission in both galaxies. The ratio of the sum of the masses of star complexes younger than 30 Myr to the associated molecular cloud masses is ~4%. The maximum age of star complexes in the galactic-scale shock front in IC 2163 is about 200 Myr, the same as the interaction time of the two galaxies, suggesting the destruction of older complexes in the eyelids.

ALMA observations in the CO 1 - 0 line of the interacting galaxies IC 2163 and NGC 2207 at 2" x 1.5" resolution reveal how the encounter drives gas to pile up in narrow, ~ 1 kpc wide, "eyelids" in IC 2163. IC 2163 and NGC 2207 are involved in a grazing encounter, which has led to development in IC 2163 of an eye-shaped (ocular) structure at mid-radius and two tidal arms. The CO data show that there are large velocity gradients across the width of each eyelid, with a mixture of radial and azimuthal streaming of gas at the outer edge of the eyelid relative to its inner edge. The sense of the radial streaming in the eyelids is consistent with the idea that gas from the outer part of IC 2163 flows inward until its radial streaming slows down abruptly and the gas piles up in the eyelids. The radial compression at the eyelids causes an increase in the gas column density by direct radial impact and also leads to a high rate of shear. We find a strong correlation between the molecular column densities and the magnitude of dv/dR across the width of the eyelid at fixed values of azimuth. Substantial portions of the eyelids have high velocity dispersion in CO, indicative of elevated turbulence there.

Giant molecular clouds contain supersonic turbulence that can locally heat small fractions of gas to over 100 K. We run shock models for low-velocity, C-type shocks propagating into gas with densities between 10^3 and 10^5 cm^(-3) and find that CO lines are the most important cooling lines. Comparison to photodissociation region (PDR) models indicates that mid-J CO lines (J = 8-7 and higher) should be dominated by emission from shocked gas. In Papers I and II we presented CO J = 3-2, 8-7, and 9-8 observations toward four primarily quiescent clumps within infrared dark clouds. Here we fit PDR models to the combined spectral line energy distributions and show that the PDR models that best fit the low-J CO emission underpredict the mid-J CO emission by orders of magnitude, strongly hinting at a hot gas component within these clumps. The low-J CO data clearly show that the integrated intensities of both the CO J = 8-7 and 9-8 lines are anomalously high, such that the line ratio can be used to characterize the hot gas component. Shock models are reasonably consistent with the observed mid-J CO emission, with models with densities near 10^(4.5) cm^(-3) providing the best agreement. Where this mid-J CO is detected, the mean volume filling factor of the hot gas is 0.1%. Much of the observed mid-J CO emission, however, is also associated with known protostars and may be due to protostellar feedback.

CO observations of the interacting galaxies IC 2163 and NGC 2207 are combined with HI, Halpha and 24 microns to study the star formation rate (SFR) surface density as a function of the gas surface density. More than half of the high SFR regions are HI dominated. When compared to other galaxies, these HI-dominated regions have excess SFRs relative to their molecular gas surface densities but normal SFRs relative to their total gas surface densities. The HI-dominated regions are mostly located in the outer part of NGC 2207, where the HI velocity dispersion is high, 40 - 50 km/s. We suggest that the star-forming clouds in these regions have envelopes at lower densities than normal, making them predominantly atomic, and cores at higher densities than normal because of the high turbulent Mach numbers. This is consistent with theoretical predictions of a flattening in the density probability distribution function for compressive, high Mach number turbulence.

Infrared dark clouds are kinematically complex molecular structures in the interstellar medium that can host sites of massive star formation. We present 4 square arcminute maps of the 12CO, 13CO, and C18O J = 3 to 2 lines from selected locations within the C and F (G028.37+00.07 and G034.43+00.24) infrared dark clouds (IRDCs), as well as single pointing observations of the 13CO and C18O J = 2 to 1 lines towards three cores within these clouds. We derive CO gas temperatures throughout the maps and find that CO is significantly frozen out within these IRDCs. We find that the CO depletion tends to be the highest near column density peaks, with maximum depletion factors between 5 and 9 in IRDC F and between 16 and 31 in IRDC C. We also detect multiple velocity components and complex kinematic structure in both IRDCs. Therefore, the kinematics of IRDCs seem to point to dynamically evolving structures yielding dense cores with considerable depletion factors.

We performed Herschel HIFI, PACS and SPIRE observations towards the molecular cloud interacting supernova remnant G349.7+0.2. An extremely broad emission line was detected at 557 GHz from the ground state transition 1_{10}-1_{01} of ortho-water. This water line can be separated into three velocity components with widths of 144, 27 and 4 km/s. The 144 km/s component is the broadest water line detected to date in the literature. This extremely broad line width shows importance of probing shock dynamics. PACS observations revealed 3 additional ortho-water lines, as well as numerous high-J carbon monoxide (CO) lines. No para-water lines were detected. The extremely broad water line is indicative of a high velocity shock, which is supported by the observed CO rotational diagram that was reproduced with a J-shock model with a density of 10^4 cm^{-3} and a shock velocity of 80 km/s. Two far-infrared fine-structure lines, [O~I] at 145 micron and [C~II] line at 157 micron, are also consistent with the high velocity J-shock model. The extremely broad water line could be simply from short-lived molecules that have not been destroyed in high velocity J-shocks; however, it may be from more complicated geometry such as high-velocity water bullets or a shell expanding in high velocity. We estimate the CO and H2O densities, column densities, and temperatures by comparison with RADEX and detailed shock models. Detection of Extremely Broad Water Emission from the molecular cloud interacting Supernova Remnant G349.7+0.2

We present VLA H I and 6 cm radio continuum observations of the spiral NGC 3145 and H I observations of its two companions, NGC 3143 and PGC 029578. In optical images NGC 3145 has stellar arms that appear to cross, forming "X"-features. Our radio continuum observations rule out shock fronts at 3 of the 4 "X"-features. In the middle-to-outer disk, the H I line-profiles of NGC 3145 are skewed. Relative to the disk, the gas in the skewed wing of the line-profiles has z-motions away from us on the approaching side of the galaxy and z-motions of about the same magnitude (about 40 km/s) towards us on the receding side. These warping motions imply that there has been a perturbation with a sizeable component perpendicular to the disk over large spatial scales. Two features in NGC 3145 have velocities indicating that they are out-of-plane tidal arms. One is an apparent branch of a main spiral arm; the velocity of the branch is 150 km/s greater than the spiral arm where they appear to intersect in projection. The other is an arm that forms 3 of the "X"-features. It differs in velocity by 56 km/s from the disk at the same projected location. Based on its SFR and H I properties, NGC 3143 is the more likely of the two companions to have interacted with NGC 3145 recently. A simple analytic model demonstrates that an encounter between NGC 3143 and NGC 3145 is a plausible explanation for the observed warping motions in NGC 3145.

Molecular oxygen, O_2, has been the target of ground-based and space-borne searches for decades. Of the thousands of lines of sight surveyed, only those toward Rho Oph and Orion H_2 Peak 1 have yielded detections of any statistical significance. The detection of the O_2 N_J =3_3 -1_2 and 5_4 - 3_4 lines at 487.249 GHz and 773.840 GHz, respectively, toward Rho Ophiuchus has been attributed to a short-lived peak in the time-dependent, cold-cloud O_2 abundance, while the detection of the O_2 N_J =3_3 - 1_2, 5_4 - 3_4 lines, plus the 7_6 - 5_6 line at 1120.715 GHz, toward Orion has been ascribed to time-dependent preshock physical and chemical evolution and low-velocity (12 km/s) non-dissociative C-type shocks, both of which are fully shielded from far-ultraviolet (FUV) radiation, plus a postshock region that is exposed to a FUV field. We report a re-interpretation of the Orion O_2 detection based on new C-type shock models that fully incorporate the significant effects the presence of even a weak FUV field can have on the preshock gas, shock structure and postshock chemistry. In particular, we show that a family of solutions exists, depending on the FUV intensity, that reproduces both the observed O_2 intensities and O_2 line ratios. The solution in closest agreement with the shock parameters inferred for H_2 Peak 1 from other gas tracers assumes a 23 km/s shock impacting gas with a preshock density of 8x10^4 cm^-3 and G_0 =1, substantially different from that inferred for the fully-shielded shock case. As pointed out previously, the similarity between the LSR velocity of all three O_2 lines (~11 km/s) and recently measured H_2O 5_32 - 4_41 maser emission at 620.701 GHz toward H_2 Peak 1 suggests that the O_2 emission arises behind the same shocks responsible for the maser emission, though the O_2 emission is almost certainly more extended than the localized high density maser spots

Infrared dark clouds (IRDCs) are dense, molecular structures in the interstellar medium that can harbour sites of high-mass star formation. IRDCs contain supersonic turbulence, which is expected to generate shocks that locally heat pockets of gas within the clouds. We present observations of the CO J = 8-7, 9-8, and 10-9 transitions, taken with the Herschel Space Observatory, towards four dense, starless clumps within IRDCs (C1 in G028.37+00.07, F1 and F2 in G034.43+0007, and G2 in G034.77-0.55). We detect the CO J = 8-7 and 9-8 transitions towards three of the clumps (C1, F1, and F2) at intensity levels greater than expected from photodissociation region (PDR) models. The average ratio of the 8-7 to 9-8 lines is also found to be between 1.6 and 2.6 in the three clumps with detections, significantly smaller than expected from PDR models. These low line ratios and large line intensities strongly suggest that the C1, F1, and F2 clumps contain a hot gas component not accounted for by standard PDR models. Such a hot gas component could be generated by turbulence dissipating in low velocity shocks.

Giant molecular clouds contain supersonic turbulence and magnetohydrodynamic simulations predict that this turbulence should decay rapidly. Such turbulent dissipation has the potential to create a warm (T ~100 K) gas component within a molecular cloud. We present observations of the CO J = 5-4 and 6-5 transitions, taken with the Herschel Space Observatory, towards the Perseus B1-East 5 region. We combine these new observations with archival measurements of lower rotational transitions and fit photodissociation region models to the data. We show that Perseus B1-E5 has an anomalously large CO J = 6-5 integrated intensity, consistent with a warm gas component existing within the region. This excess emission is consistent with predictions for shock heating due to the dissipation of turbulence in low velocity shocks with the shocks having a volume filling factor of 0.15 per cent. We find that B1-E has a turbulent energy dissipation rate of 3.5 x 10$^{32}$ erg / s and a dissipation time-scale that is only a factor of 3 larger than the flow crossing time-scale.

The energy to desorb atomic oxygen from an interstellar dust grain surface, $E_{\rm des}$, is an important controlling parameter in gas-grain models; its value impacts the temperature range over which oxygen resides on a dust grain. However, no prior measurement has been done of the desorption energy. We report the first direct measurement of $E_{\rm des}$ for atomic oxygen from dust grain analogs. The values of $E_{\rm des}$ are $1660\pm 60$~K and $1850\pm 90$~K for porous amorphous water ice and for a bare amorphous silicate film, respectively, or about twice the value previously adopted in simulations of the chemical evolution of a cloud. We use the new values to study oxygen chemistry as a function of depth in a molecular cloud. For $n=10^4$ cm$^{-3}$ and $G_0$=10$^2$ ($G_0$=1 is the average local interstellar radiation field), the main result of the adoption of the higher oxygen binding energy is that H$_2$O can form on grains at lower visual extinction $A_{\rm V}$, closer to the cloud surface. A higher binding energy of O results in more formation of OH and H$_2$O on grains, which are subsequently desorbed by FUV radiation, with consequences for gas-phase chemistry. For higher values of $n$ and $G_0$, the higher binding energy can lead to a large increase in the column of H$_2$O but a decrease in the column of O$_2$.

Protostars interact with their surroundings through jets and winds impacting on the envelope and creating shocks, but the nature of these shocks is still poorly understood. Our aim is to survey far-infrared molecular line emission from a uniform and significant sample of deeply-embedded low-mass young stellar objects in order to characterize shocks and the possible role of ultraviolet radiation in the immediate protostellar environment. Herschel/PACS spectral maps of 22 objects in the Perseus molecular cloud were obtained as part of the `William Herschel Line Legacy' survey. Line emission from H$_\mathrm{2}$O, CO, and OH is tested against shock models from the literature. Observed line ratios are remarkably similar and do not show variations with source physical parameters. Observations show good agreement with the shock models when line ratios of the same species are compared. Ratios of various H$_\mathrm{2}$O lines provide a particularly good diagnostic of pre-shock gas densities, $n_\mathrm{H}\sim10^{5}$ cm$^{-3}$, in agreement with typical densities obtained from observations of the post-shock gas. The corresponding shock velocities, obtained from comparison with CO line ratios, are above 20 km\,s$^{-1}$. However, the observations consistently show one-to-two orders of magnitude lower H$_\mathrm{2}$O-to-CO and H$_\mathrm{2}$O-to-OH line ratios than predicted by the existing shock models. The overestimated model H$_\mathrm{2}$O fluxes are most likely caused by an overabundance of H$_\mathrm{2}$O in the models since the excitation is well-reproduced. Illumination of the shocked material by ultraviolet photons produced either in the star-disk system or, more locally, in the shock, would decrease the H$_\mathrm{2}$O abundances and reconcile the models with observations. Detections of hot H$_\mathrm{2}$O and strong OH lines support this scenario.

We report observations of molecular oxygen (O$_2$) rotational transitions at 487 GHz, 774 GHz, and 1121 GHz toward Orion Peak A. The O2 lines at 487 GHz and 774 GHz are detected at velocities of 10-12 km/s with line widths 3 km/s; however, the transition at 1121 GHz is not detected. The observed line characteristics, combined with the results of earlier observations, suggest that the region responsible for the O$_2$ emission is 9" (6e16 cm) in size, and is located close to the H2 Peak 1position (where vibrationally-excited H$_2$ emission peaks), and not at Peak A, 23" away. The peak O2 column density is 1.1e18/cm2. The line velocity is close to that of 621 GHz water maser emission found in this portion of the Orion Molecular Cloud, and having a shock with velocity vector lying nearly in the plane of the sky is consistent with producing maximum maser gain along the line-of-sight. The enhanced O$_2$ abundance compared to that generally found in dense interstellar clouds can be explained by passage of a low-velocity C-shock through a clump with preshock density 2e4/cm3, if a reasonable flux of UV radiation is present. The postshock O$_2$ can explain the emission from the source if its line of sight dimension is ~10 times larger than its size on the plane of the sky. The special geometry and conditions required may explain why O$_2$ emission has not been detected in the cores of other massive star-forming molecular clouds.

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.

We report the results of a search for molecular oxygen (O2) toward the Orion Bar, a prominent photodissociation region at the southern edge of the HII region created by the luminous Trapezium stars. We observed the spectral region around the frequency of the O2 N_J = 3_3 - 1_2 transition at 487 GHz and the 5_4 - 3_4 transition at 774 GHz using the Heterodyne Instrument for the Far Infrared on the Herschel Space Observatory. Neither line was detected, but the 3sigma upper limits established here translate to a total line-of-sight O2 column density < 1.5 10^16 cm^-2 for an emitting region whose temperature is between 30K and 250 K, or < 1 10^16 cm^-2 if the O2 emitting region is primarily at a temperature of ~< 100 K. Because the Orion Bar is oriented nearly edge-on relative to our line of sight, the observed column density is enhanced by a factor estimated to be between 4 and 20 relative to the face-on value. Our upper limits imply that the face-on O2 column density is less than 4 10^15 cm^-2, a value that is below, and possibly well below, model predictions for gas with a density of 10^4 - 10^5 cm^-3 exposed to a far ultraviolet flux 10^4 times the local value, conditions inferred from previous observations of the Orion Bar. The discrepancy might be resolved if: (1) the adsorption energy of O atoms to ice is greater than 800 K; (2) the total face-on Av of the Bar is less than required for O2 to reach peak abundance; (3) the O2 emission arises within dense clumps with a small beam filling factor; or, (4) the face-on depth into the Bar where O2 reaches its peak abundance, which is density dependent, corresponds to a sky position different from that sampled by our Herschel beams.

Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a beam size of 10 arcmin. A larger telescope aperture such as that of the Herschel Space Observatory is required to resolve the O2 emission and to pinpoint its origin. We use the Heterodyne Instrument for the Far Infrared aboard Herschel to obtain high resolution O2 spectra toward selected positions in rho Oph A. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. The 487.2GHz line was clearly detected toward all three observed positions in rho Oph A. In addition, an oversampled map of the 773.8GHz transition revealed the detection of the line in only half of the observed area. Based on their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (> 50 K) adjacent to a much colder region, where temperatures are below 30 K. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on the derived gas temperatures, but fall into the range N(O2) = (3 to >6)e15/cm^2. Beam averaged O2 abundances are about 5e-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz of about 4 - 5 arcmin, encompassing the entire rho Oph A core.

Giant molecular clouds contain supersonic turbulence and simulations of magnetohydrodynamic turbulence show that these supersonic motions decay in roughly a crossing time, which is less than the estimated lifetimes of molecular clouds. Such a situation requires a significant release of energy. We run models of C-type shocks propagating into gas with densities around 10^3 cm^(-3) at velocities of a few km / s, appropriate for the ambient conditions inside of a molecular cloud, to determine which species and transitions dominate the cooling and radiative energy release associated with shock cooling of turbulent molecular clouds. We find that these shocks dissipate their energy primarily through CO rotational transitions and by compressing pre-existing magnetic fields. We present model spectra for these shocks and by combining these models with estimates for the rate of turbulent energy dissipation, we show that shock emission should dominate over emission from unshocked gas for mid to high rotational transitions (J >5) of CO. We also find that the turbulent energy dissipation rate is roughly equivalent to the cosmic-ray heating rate and that the ambipolar diffusion heating rate may be significant, especially in shocked gas.

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 report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2.

Radio continuum, Spitzer infrared, optical and XMM-Newton X-ray and UVM2 observations are used to study large-scale shock fronts, young star complexes, and the galactic nuclei in the interacting galaxies NGC 2207/IC 2163. There are two types of large-scale shock fronts in this galaxy pair. The shock front along the rim of the ocular oval in IC 2163 has produced vigorous star formation in a dusty environment. In the outer part of the companion side of NGC 2207, a large-scale front attributed to disk or halo scraping is particularly bright in the radio continuum but not in any tracers of recent star formation or in X-rays. This radio continuum front may be mainly in the halo on the back side of NGC 2207 between the two galaxies. Values of the flux density ratio S(8 um)/S(6 cm) of kpc-sized, Spitzer IRAC star-forming clumps in NGC 2207/IC 2163 are compared with those of giant H II regions in M81. We find evidence that in 2001 a radio supernova was present in the core of feature i, a mini-starburst on an outer arm of NGC 2207. X-ray emission is detected from the NGC 2207 nucleus and from nine discrete sources, one of which corresponds to SN 1999ec, and another may be a radio supernova or a background quasar. The X-ray luminosity and X-ray spectrum of the NGC 2207 nucleus suggests it is a highly absorbed, low luminosity AGN.

We report the results of a large-area study of water vapor along the Orion Molecular Cloud ridge, the purpose of which was to determine the depth-dependent distribution of gas-phase water in dense molecular clouds. We find that the water vapor measured toward 77 spatial positions along the face-on Orion ridge, excluding positions surrounding the outflow associated with BN/KL and IRc2, display integrated intensities that correlate strongly with known cloud surface tracers such as CN, C2H, 13CO J =5-4, and HCN, and less well with the volume tracer N2H+. Moreover, at total column densities corresponding to Av < 15 mag., the ratio of H2O to C18O integrated intensities shows a clear rise approaching the cloud surface. We show that this behavior cannot be accounted for by either optical depth or excitation effects, but suggests that gas-phase water abundances fall at large Av. These results are important as they affect measures of the true water-vapor abundance in molecular clouds by highlighting the limitations of comparing measured water vapor column densities with such traditional cloud tracers as 13CO or C18O. These results also support cloud models that incorporate freeze-out of molecules as a critical component in determining the depth-dependent abundance of water vapor.

'Water In Star-forming regions with Herschel' (WISH) is a key programme dedicated to studying the role of water and related species during the star-formation process and constraining the physical and chemical properties of young stellar objects. The Heterodyne Instrument for the Far-Infrared (HIFI) on the Herschel Space Observatory observed three deeply embedded protostars in the low-mass star-forming region NGC1333 in several H2-16O, H2-18O, and CO transitions. Line profiles are resolved for five H16O transitions in each source, revealing them to be surprisingly complex. The line profiles are decomposed into broad (>20 km/s), medium-broad (~5-10 km/s), and narrow (<5 km/s) components. The H2-18O emission is only detected in broad 1_10-1_01 lines (>20 km/s), indicating that its physical origin is the same as for the broad H2-16O component. In one of the sources, IRAS4A, an inverse P Cygni profile is observed, a clear sign of infall in the envelope. From the line profiles alone, it is clear that the bulk of emission arises from shocks, both on small (<1000 AU) and large scales along the outflow cavity walls (~10 000 AU). The H2O line profiles are compared to CO line profiles to constrain the H2O abundance as a function of velocity within these shocked regions. The H2O/CO abundance ratios are measured to be in the range of ~0.1-1, corresponding to H2O abundances of ~10-5-10-4 with respect to H2. Approximately 5-10% of the gas is hot enough for all oxygen to be driven into water in warm post-shock gas, mostly at high velocities.

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

In the framework of the Water in Star-forming regions with Herschel (WISH) key program, maps in water lines of several outflows from young stars are being obtained, to study the water production in shocks and its role in the outflow cooling. This paper reports the first results of this program, presenting a PACS map of the o-H2O 179 um transition obtained toward the young outflow L1157. The 179 um map is compared with those of other important shock tracers, and with previous single-pointing ISO, SWAS, and Odin water observations of the same source that allow us to constrain the water abundance and total cooling. Strong H2O peaks are localized on both shocked emission knots and the central source position. The H2O 179 um emission is spatially correlated with emission from H2 rotational lines, excited in shocks leading to a significant enhancement of the water abundance. Water emission peaks along the outflow also correlate with peaks of other shock-produced molecular species, such as SiO and NH3. A strong H2O peak is also observed at the location of the proto-star, where none of the other molecules have significant emission. The absolute 179 um intensity and its intensity ratio to the H2O 557 GHz line previously observed with Odin/SWAS indicate that the water emission originates in warm compact clumps, spatially unresolved by PACS, having a H2O abundance of the order of 10^-4. This testifies that the clumps have been heated for a time long enough to allow the conversion of almost all the available gas-phase oxygen into water. The total water cooling is ~10^-1 Lo, about 40% of the cooling due to H2 and 23% of the total energy released in shocks along the L1157 outflow.

'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CII] lines located between 63 and 186 um. The excitation and spatial distribution of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, [OI] emission is seen in low-velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km/s in both the red- and blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [OI] emission is attributed to J-type shocks in dense gas close to the protostar. In the scenario described here, the combined cooling by far-IR lines within the central spatial pixel is estimated to be 2 \times 10-2 L_sun, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope.

We report the results of spectroscopic mapping observations carried out toward protostellar outflows in the BHR71, L1157, L1448, NGC 2071, and VLA 1623 molecular regions using the Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations, covering the 5.2 - 37 micron spectral region, provide detailed maps of the 8 lowest pure rotational lines of molecular hydrogen and of the [SI] 25.25 micron and [FeII] 26.0 micron fine structure lines. The molecular hydrogen lines, believed to account for a large fraction of the radiative cooling from warm molecular gas that has been heated by a non-dissociative shock, allow the energetics of the outflows to be elucidated. Within the regions mapped towards these 5 outflow sources, total H2 luminosities ranging from 0.02 to 0.75 L(solar) were inferred for the sum of the 8 lowest pure rotational transitions. By contrast, the much weaker [FeII] 26.0 micron fine structure transition traces faster, dissociative shocks; here, only a small fraction of the fast shock luminosity emerges as line radiation that can be detected with Spitzer/IRS.

We model the temperature and chemical structure of molecular clouds as a function of depth into the cloud, assuming a cloud of constant density n illuminated by an external FUV (6 eV < E < 13.6 eV) flux G_0 (scaling factor in multiples of the local interstellar field). Extending previous photodissociation region models, we include the freezing of species, simple grain surface chemistry, and desorption (including FUV photodesorption) of ices. We also treat the opaque cloud interior with time-dependent chemistry. Here, under certain conditions, gas phase elemental oxygen freezes out as water ice and the elemental C/O abundance ratio can exceed unity, leading to complex carbon chemistry. Gas phase H2O and O2 peak in abundance at intermediate depth into the cloud, roughly A_V~3-8 from the surface, the depth proportional to ln(G_0/n). Closer to the surface, molecules are photodissociated. Deeper into the cloud, molecules freeze to grain surfaces. At intermediate depths photodissociation rates are attenuated by dust extinction, but photodesorption prevents total freezeout. For G_0 < 500, abundances of H2O and O2 peak at values ~10^(-7), producing columns ~10^(15) per cm^2, independent of G_0 and n. The peak abundances depend primarily on the product of the photodesorption yield of water ice and the grain surface area per H nucleus. At higher values of G_0, thermal desorption of O atoms from grains enhances the gas phase H2O peak abundance and column slightly, whereas the gas phase O2 peak abundance rises to ~10^(-5) and the column to ~2x10^(16) per cm^2. We present simple analytic equations for the abundances as a function of depth which clarify the dependence on parameters. The models are applied to observations of H2O, O2, and water ice in a number of sources, including B68, NGC 2024, and Rho Oph.

We report the results of spectroscopic mapping observations carried out toward a ~1 min x 1 min region within the northern lobe of the outflow from NGC 2071 using the Infrared Spectrograph (IRS) of the Spitzer Space Telescope. These observations covered the 5.2-37 um spectral region and have led to the detection of a number of ionic, atomic, and molecular lines, including fine-structure emission of Si+, Fe+, S++, S, the S(0)-S(7) pure rotational lines of H2, the R(3) and R(4) transitions of HD, and at least 11 transitions of H2O. In addition, the 6.2, 7.4, 7.6, 7.9, 8.6 and 11.3 um PAH emission bands were also observed and several transitions of OH were tentatively detected. Most of the detected line transitions were strong enough to map including, for the first time, three transitions of hot H2O. We find that: (1) the water emission is extended; (2) the extended emission is aligned with the outflow; and, (3) the spatial distribution of the water emission generally follows that observed for H2. Based on the measured line intensities, we derive an HD abundance relative to H2 of 1.1-1.8 10^-5 and an H2O number density of 12-2 cm^3. The H2 density in the water-emitting region is not well constrained by our observations, but is likely between 3 10^4 and 10^6 cm^3, yielding an H2O abundance relative to H2 of between 2 10^-5 and 6 10^-4. Future observations planned for the Herschel Space Observatory should greatly improve the density estimate, and thus our knowledge of the H2O abundance, for the water-emitting regions reported here. Finally, we note a possible departure from the H2O ortho-to-para ratio of 3:1 expected for water formed in hot post-shocked gas, suggesting that a significant fraction of the water vapor we detect may arise from H2O sputtered from cold dust grains.

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 present SWAS detections of the ground-state 1(10)-1(01) transition of o-H2O at 557 GHz in 18 molecular outflows. These results are combined with ground-based observations of the J=1-0 transitions of 12CO and 13CO obtained at the FCRAO and, for a subset of the outflows, data from ISO. Assuming the SWAS water line emission originates from the same gas traced by CO emission, we find that the outflowing gas in most outflows has an o-H2O abundance relative to H2 of between 10(-7) and 10(-6). Analysis of the water abundance as a function of outflow velocity reveals a strong dependence. The water abundance increases with velocity, and at the highest outflow velocities some outflows have relative o-H2O abundances of order 10(-4). However the mass of gas with such elevated water abundances represents less that 1% of the total outflow gas mass. The ISO LWS observations of high-J rotational lines of CO and the 179.5 micron transition of o-H2O provide evidence for a warmer outflow component than required to produce either the SWAS or FCRAO lines. The mass associated with the ISO emission is similar to that responsible for the highest velocity water emission detected by SWAS. The bulk of the outflowing gas has an abundance of o-H2O well below what would be expected if the gas has passed through a C-shock with shock velocities greater than 10 km/s. Gas-phase water can be depleted in the post-shock gas due to freeze-out onto grain mantles, however the rate of freeze-out is too slow to explain our results. Therefore we believe that only a small fraction of the outflowing molecular gas has passed through shocks strong enough to fully convert the gas-phase oxygen to water. This result has implications for the acceleration mechanism of the molecular gas in these outflows.

HST NICMOS narrowband images of the shocked molecular hydrogen emission in OMC-1 are analyzed to reveal new information on the BN/KL outflow. The outstanding morphological feature of this region is the array of molecular hydrogen ``fingers'' emanating from the general vicinity of IRc2 and the presence of several Herbig-Haro objects. The NICMOS images appear to resolve individual shock fronts. This work is a more quantitative and detailed analysis of our data from a previous paper (Schultz etal. 1999). Line strengths for the H_2 1--0 S(4) plus 2--1 S(6) lines at 1.89 micron are estimated from measurements with the Paschen_alpha continuum filter F190N at 1.90 micron, and continuum measurements at 1.66 and 2.15 micron. We compare the observed H_2 line strengths and ratios of the 1.89 micron and 2.12 micron 1--0 S(1) lines with models for molecular cloud shock waves. Most of the data cannot be fit by J-shocks, but are well matched by C-shocks with shock velocities in the range of 20--45 km/s and preshock densities of 10^{4} - 10^{6} cm^{-3}, similiar to values obtained in larger beam studies which averaged over many shocks. There is also some evidence that shocks with higher densities have lower velocities.

We present Spitzer IRS spectra (R ~600, 10 - 38 micron) of 38 positions in the Galactic Center (GC), all at the same Galactic longitude and spanning plus/minus 0.3 degrees in latitude. Our positions include the Arches Cluster, the Arched Filaments, regions near the Quintuplet Cluster, the ``Bubble'' lying along the same line-of-sight as the molecular cloud G0.11-0.11, and the diffuse interstellar gas along the line-of-sight at higher Galactic latitudes. From measurements of the [O IV], [Ne II], [Ne III], [Si II], [S III], [S IV], [Fe II], [Fe III], and H_2 S(0), S(1), and S(2) lines we determine the gas excitation and ionic abundance ratios. The Ne/H and S/H abundance ratios are ~ 1.6 times that of the Orion Nebula. The main source of excitation is photoionization, with the Arches Cluster ionizing the Arched Filaments and the Quintuplet Cluster ionizing the gas nearby and at lower Galactic latitudes including the far side of the Bubble. In addition, strong shocks ionize gas to O^{+3} and destroy dust grains, releasing iron into the gas phase (Fe/H ~ 1.3 times 10^{-6} in the Arched Filaments and Fe/H ~ 8.8 times 10^{-6} in the Bubble). The shock effects are particularly noticeable in the center of the Bubble, but O$^{+3}$ is present in all positions. We suggest that the shocks are due to the winds from the Quintuplet Cluster Wolf-Rayet stars. On the other hand, the H_2 line ratios can be explained with multi-component models of warm molecular gas in photodissociation regions without the need for H_2 production in shocks.

We report the results of spectroscopic mapping observations carried out toward small (1 x 1 arcmin) regions within the supernova remnants W44, W28, IC443, and 3C391 using the Infrared Spectrograph of the Spitzer Space Telescope. These observations, covering the 5.2 - 37 micron spectral region, have led to the detection of a total of 15 fine structure transitions of Ne+, Ne++, Si+, P+, S, S++, Cl+, Fe+, and Fe++; the S(0) - S(7) pure rotational lines of molecular hydrogen; and the R(3) and R(4) transitions of hydrogen deuteride. In addition to these 25 spectral lines, the 6.2, 7.7, 8.6, 11.3 and 12.6 micron PAH emission bands were also observed. Most of the detected line transitions have proven strong enough to map in several sources, providing a comprehensive picture of the relative distribution of the various line emissions observable in the Spitzer/IRS bandpass. A principal component analysis of the spectral line maps reveals that the observed emission lines fall into five distinct groups, each of which may exhibit a distinct spatial distribution: (1) lines of S and H2 (J > 2); (2) the H2 S(0) line; (3) lines of ions with appearance potentials less than 13.6 eV; (4) lines of ions with appearance potentials greater than 13.6 eV, not including S++; (5) lines of S++. Lines of group (1) likely originate in molecular material subject to a slow, nondissociative shock that is driven by the overpressure within the supernova remnant, and lines in groups (3) - (5) are associated primarily with dissociative shock fronts with a range of (larger) shock velocities. The H2 S(0) line shows a low-density diffuse emission component, and - in some sources - a shock-excited component.

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 report the detection of interstellar hydrogen deuteride (HD) toward the supernova remnant IC443, and the tentative detection of HD toward the Herbig Haro objects HH54 and HH7 and the star forming region GGD37 (Cepheus A West). Our detections are based upon spectral line mapping observations of the R(3) and R(4) rotational lines of HD, at rest wavelengths of 28.502 and 23.034 micron respectively, obtained using the Infrared Spectrograph onboard the Spitzer Space Telescope. The HD R(4)/R(3) line intensity ratio promises to be a valuable probe of the gas pressure in regions where it can be observed. The derived HD/H2 abundance ratios are 1.19(+0.35/-0.24)E-5, 1.80(+0.54/-0.32)E-5, and 1.41(+0.46/-0.33)E-5 respectively (68.3% confidence limits, based upon statistical errors alone) for IC443 (clump C), HH54, and HH7. If HD is the only significant reservoir of gas-phase deuterium in these sources, the inferred HD/H2 ratios are all consistent with a gas-phase elemental abundance [n(D)/n(H)](gas) ~ 7.5E-6, a factor 2 - 3 below the values obtained previously from observations of atomic deuterium in the local bubble and the Galactic halo. However, similarly low gas-phase deuterium abundances have been inferred previously for molecular gas clouds in the Orion region, and in atomic clouds along sight-lines within the Galactic disk to stars more distant than 500 pc from the Sun.

IC 2163 and NGC 2207 are interacting galaxies that have been well studied at optical and radio wavelengths and simulated in numerical models to reproduce the observed kinematics and morphological features. Spitzer IRAC and MIPS observations reported here show over 200 bright clumps from young star complexes. The brightest IR clump is a morphologically peculiar region of star formation in the western arm of NGC 2207. This clump, which dominates the Halpha and radio continuum emission from both galaxies, accounts for ~12% of the total 24mu m flux. Nearly half of the clumps are regularly spaced along some filamentary structure, whether in the starburst oval of IC 2163 or in the thin spiral arms of NGC 2207. This regularity appears to influence the clump luminosity function, making it peaked at a value nearly a factor of 10 above the completeness limit, particularly in the starburst oval. This is unlike the optical clusters inside the clumps, which have a luminosity function consistent with the usual power law form. The giant IR clumps presumably formed by gravitational instabilities in the compressed gas of the oval and the spiral arms, whereas the individual clusters formed by more chaotic processes, such as turbulence compression, inside these larger-scale structures.

Recent observations of the ground state transition of HDO at 464 GHz towards the protoplanetary disk of DM Tau have detected the presence of water vapor in the regions just above the outer disk midplane (Ceccarelli et al 2005). In the absence of non-thermal desorption processes, water should be almost entirely frozen onto the grain mantles and HDO undetectable. In this Letter we present a chemical model that explores the possibility that the icy mantles are photo-desorbed by FUV (6eV < h nu < 13.6eV) photons. We show that the average Interstellar FUV field is enough to create a layer of water vapor above the disk midplane over the entire disk. Assuming a photo-desorption yield of 10^{-3}, the water abundance in this layer is predicted to be ~ 3 x 10^{-7} and the average H2O column density is ~ 1.6x 10^{15} cm^{-2}. The predictions are very weakly dependent on the details of the model, like the incident FUV radiation field, and the gas density in the disk. Based on this model, we predict a gaseous HDO/H2O ratio in DM Tau of ~1%. In addition, we predict the ground state transition of water at 557 GHz to be undetectable with ODIN and/or HSO-HIFI.

We present numerical hydrodynamical models of the collision between the galaxies IC 2163 and NGC 2207. These models extend the results of earlier work where the galaxy discs were modeled one at a time. We confirm the general result that the collision is primarily planar, that is, at moderate inclination relative to the two discs, and prograde for IC 2163, but retrograde for NGC 2207. We list 34 specific morphological or kinematic features on a variety of scales, found with multi-waveband observations, which we use to constrain the models. The models are able to reproduce most of these features, with a relative orbit in which the companion (IC 2163) disc first side-swipes the primary (NGC 2207) disc on the west side, then moves around the edge of the primary disc to the north and to its current position on the east side. The models also provide evidence that the dark matter halo of NGC 2207 has only moderate extent. For IC 2163, the prolonged prograde disturbance in the model produces a tidal tail, and an oval or ocular waveform very much like the observed ones, including some fine structure. The retrograde disturbance in the model produces no strong waveforms within the primary galaxy. This suggests that the prominent spiral waves in NGC 2207 were present before the collision, and models with waves imposed in the initial conditions confirm that they would not be disrupted by the collision. With an initial central hole in the gas disc of the primary, and imposed spirals, the model also reproduces the broad ring seen in HI observations. We believe that this is now one of the best modeled systems of colliding galaxies. (Abridged)

We report the detection of the 2P_3/2 -> 2P_1/2 fine-structure line of C+ at 157.74 micron in SDSSJ114816.64+525150.3 (hereafter J1148+5251), the most distant known quasar, at z=6.42, using the IRAM 30-meter telescope. This is the first detection of the [CII] line at high redshift, and also the first detection in a Hyperluminous Infrared Galaxy (L_FIR > 10^13 Lsun). The [CII] line is detected at a significance level of 8 sigma and has a luminosity of 4.4 x 10^9 Lsun. The L_[CII]/L_FIR ratio is 2 x 10^-4, about an order of magnitude smaller than observed in local normal galaxies and similar to the ratio observed in local Ultraluminous Infrared Galaxies. The [CII] line luminosity indicates that the host galaxy of this quasar is undergoing an intense burst of star formation with an estimated rate of ~3000 Msun/yr. The detection of C+ in SDSS J1148+5251 suggests a significant enrichment of metals at z ~ 6 (age of the universe ~870 Myr), although the data are consistent with a reduced carbon to oxygen ratio as expected from chemical evolutionary models of the early phases of galaxy formation.

We observed Hubble V, the brightest HII region complex in the dwarf irregular galaxy NGC 6822, at near-infrared (1.8--2.4 um) wavelengths using the Cooled Grating Spectrometer 4 (CGS4) at the United Kingdom Infra-Red Telescope (UKIRT). The line emission maps of Hubble V show the typical structure of a photo-dissociation region (PDR) where an ionized core, traced by compact He I emission (2.0587 um) and Br_gamma emission (2.1661 um), is surrounded by an outer layer traced by molecular hydrogen (H_2) emission. The measured line ratios of H_2 2-1 S(1) (2.2477 um) / 1-0 S(1) (2.1218 um) from 0.2 to 0.6 and the un-shifted and un-resolved line profiles suggest that the H_2 emission originates purely from a PDR. By comparing the H_2 results with a PDR model, we conclude that Hubble V includes dense (10^{4.5} cm^{-3}) and warm PDRs. In this environment, most of the H_2 molecules are excited by far-UV photons (with a field strength of 10^{2-4} times that of the average interstellar field), although collisional processes de-excite H_2 and contribute significantly to the excitation of the first vibrational level. We expect that Hubble V is in the early stage of molecular cloud dissolution.

The gas at the surfaces of molecular clouds in galaxies is heated and dissociated by photons from young stars both near and far. HI resulting from the dissociation of molecular hydrogen H2 emits hyperfine line emission at 21 cm, and warmed CO emits dipole rotational lines such as the 2.6 mm line of CO(1-0). We use previously developed models for photodissociation regions (PDRs) to compute the intensities of these HI and CO(1-0) lines as a function of the total volume density n in the cloud and the far ultraviolet flux G0 incident upon it and present the results in units familiar to observers. The intensities of these two lines behave differently with changing physical conditions in the PDR, and, taken together, the two lines can provide a ground-based radio astronomy diagnostic for determining n and G0 separately in distant molecular clouds. This diagnostic is particularly useful in the range Gzero <~ 100, 10 cm^{-3} <~ n <~ 10^5 cm^{-3}, which applies to a large fraction of the volume of the interstellar medium in galaxies. If the molecular cloud is located near discrete sources of far-UV (FUV) emission, the PDR-generated HI and CO(1-0) emission on the cloud surface can be more easily identified, appearing as layered ``blankets'' or ``blisters'' on the side of the cloud nearest to the FUV source. As an illustration, we consider the Galactic object G216 -2.5, i.e. ``Maddalena's Cloud'', which has been previously identified as a large PDR in the Galaxy. We determine that this cloud has n ~ 200 cm^{-3}, G0 ~ 0.8, consistent with other data.

We have observed 6 clouds along the line-of-sight toward W49A using the Submillimeter Wave Astronomy Satellite (SWAS) and several ground-based observatories. The ortho-H2O 1-0 and OH (1665 and 1667 MHz) transitions are observed in absorption, whereas the low-J CO, 13CO, and C18O lines, as well as the [CI] 1-0 transition, are seen in emission. By using both the o-H218O and o-H2O absorption lines, we are able to constrain the column-averaged o-H_2O abundances in each line-of-sight cloud to within about an order of magnitude. Assuming the standard N(H2)/N(CO) ratio of 10^4, we find N(o-H2O)/N(H2) = 8.1 x 10^-8 - 4 x 10^-7 for three clouds with optically thin water lines. In three additional clouds, the H$_2$O lines are saturated so we have used observations of the H218O ground-state transition to find upper limits to the water abundance of 8.2x 10^-8 - 1.5x10^-6. We measure the OH abundance from the average of the 1665 and 1667 MHz observations and find N(OH)/N(H2) = 2.3x10^-7 - 1.1x10^-6. The o-H2O and OH abundances are similar to those determined for line-of-sight water absorption features towards W51 and Sgr B2 but are higher than those seen from water emission lines in molecular clouds. However, the clouds towards W49 have lower ratios of OH relative to H2O column densities than are predicted by simple models which assume that dissociative recombination is the primary formation pathway for OH and H2O. Building on the work of Neufeld et al. (2002), we present photo-chemistry models including additional chemical effects, which can also explain the observed OH and H2O column densities as well as the observed H2O/CO abundance ratios.

Using NASA's Submillimeter Wave Astronomy Satellite (SWAS) we have examined the production of water in quiescent and shocked molecular gas through a survey of the 556.936 GHz transition of ortho-H2O in the NGC 1333 molecular core. These observations reveal broad emission lines associated with the IRAS~2, IRAS~4, IRAS~7, and HH7-11 outflows. Towards 3 positions we detect narrow (~2-3 km/s) emission lines clearly associated with the ambient gas. The SWAS observations, with a resolution of 4', are supplemented with observations from the Infrared Space Observatory (ISO) and by an unbiased survey of a 17' x 15' area, with 50" resolution, in the low-J transitions of CO, 13CO, C18O, N2H+, CH3OH, and SiO. Using these combined data sets, with consistent assumptions, we find beam-averaged o-H2O abundances of > 10^{-6} relative to H2O for all four outflows. A comparison of SWAS and ISO water data is consistent with non-dissociative shock models, provided the majority of the 557 GHz emission arises from cool post-shock material with enhanced abundances. In the ambient gas the o-H2O abundance is found to lie between 0.1-1 x 10^{-7} relative to H2 and is enhanced when compared to cold pre-stellar molecular cores. A comparison of the water emission with tracers of dense condensations and shock chemistry finds no clear correlation. However, the water emission appears to be associated with the presence of luminous external heating sources which power the reflection nebula and the photodissociation (PDR) region. Simple PDR models are capable of reproducing the water and high-J 13CO emission, suggesting that a PDR may account for the excitation of water in low density undepleted gas as suggested by Spaans & van Dishoeck (2001).

Observations of W51 with the Submillimeter Wave Astronomy Satellite (SWAS) have yielded the first detection of water vapor in a diffuse molecular cloud. The water vapor lies in a foreground cloud that gives rise to an absorption feature at an LSR velocity of 6 km/s. The inferred H2O column density is 2.5E+13 cm-2. Observations with the Arecibo radio telescope of hydroxyl molecules at ten positions in W51 imply an OH column density of 8E+13 cm-2 in the same diffuse cloud. The observed H2O/OH ratio of ~ 0.3 is significantly larger than an upper limit derived previously from ultraviolet observations of the similar diffuse molecular cloud lying in front of HD 154368. The observed variation in H2O/OH likely points to the presence in one or both of these clouds of a warm (T > 400) gas component in which neutral-neutral reactions are important sources of OH and/or H2O.

(Abridged) We present the analysis of the main FIR fine structure lines emission in NGC1313 and NGC6946. We calculate that a component probably associated with the diffuse disks contributes <~40% in N6946 and ~30 % in N1313 to the total [CII] emission. The main PDR physical parameters responsible for the neutral atomic gas emission in N1313 and N6946 do not significantly differ from what Malhotra etal (2001) found by modelling the integrated emission of a sample of 60 normal galaxies,although there are evidences for a beam averaged contribution of a less active component inside NGC6946 higher than its contribution in the integrated emission of normal galaxies. CO and [CII] in N6946 are well correlated with a mean [CII]/CO ratio similar to that of the normal galaxies sample. In N1313 the [CII]/CO seems to systematically increase from the North to the south, along the S-shaped spiral arm, indicating much more inhomogeneous conditions than in N6946. HI and [CII] in N6946 are completely de-correlated, probably because they arise from different gas components. In N1313 we successfully detect two distinct gas components: a cirrus-like component where HI and [CII] are weakly correlated as observed in our Galaxy, and a component associated with dense PDRs completely de-correlated from HI as observed in N6946.Finally, we find that the HI residing in dense PDRs and presumably recently photo-dissociated, constitutes a few % of the total HI. In turn, this dense gas component produces most of the [CII] emission emitted by the atomic neutral medium.

BIMA CO 1-0 observations are presented of the spiral galaxies NGC 5394 and NGC 5395 that have undergone a recent, grazing encounter. In NGC 5394, 80% of the CO emission detected by BIMA is concentrated in the central 800 pc (FWHM) starburst region.In an encounter simulation that reproduces some of the main features of this galaxy pair, a considerable amount of gas in NGC 5394 falls into the central region early in the collision. The observed total gas distribution in the disk of NGC 5394 is lopsided, with more HI, CO, and H-alpha emission coming from the western or southwestern side. The innermost western arm of NGC 5394 is seen in CO and H-alpha emission, but the eastern inner-disk arm, which is very bright in the optical continuum, is not detected in CO or H-alpha emission. From a comparison of the radio continuum, H-alpha, 60 micron, and CO luminosities, we estimate that the average visual extinction of the starburst is 3 - 4 mag and the conversion factor N(H2)/I(CO) in the starburst is a factor of 3 - 4 below the standard value. Comparison of NGC 5394 with two other systems previously studied suggests that in prograde grazing encounters a central starburst may not develop until near the end of the ocular phase. Very little of the CO emission from NGC 5395 found in previous single-dish observations is detected by BIMA.

The most important cooling lines of the neutral interstellar medium (ISM) lie in the far-infrared (FIR). We present measurements by the Infrared Space Observatory Long Wavelength Spectrometer of seven lines from neutral and ionized ISM of 60 normal, star-forming galaxies. The galaxy sample spans a range in properties such as morphology, FIR colors (indicating dust temperature), and FIR/Blue ratios (indicating star-formation activity and optical depth). In two-thirds of the galaxies in this sample, the [CII] line is proportional to FIR dust continuum. The other one-third show a smooth decline in [CII]/FIR with increasing F60/F100 and FIR/B, spanning a range of a factor of more than 50. Two galaxies, at the warm and active extreme of the range have [CII]/FIR < 2 \times 10^{-4} (3-sigma upper limit). This is due to increased positive grain charge in the warmer and more active galaxies, which leads to less efficient heating by photoelectrons from dust grains. The ratio of the two principal photodissociation region (PDR) cooling lines [CII]/[OI] shows a tight correlation with F60/F100, indicating that both gas and dust temperatures increase together. We derive a theoretical scaling between [NII] and [CII] from ionized gas and use it to separate [CII] emission from neutral PDRs and ionized gas. Comparison of PDR models of Kaufman et al. (1999) with observed ratios of (a) [OI]/[CII] and ([CII]+[OI])/FIR and (b) [OI]/FIR and F60/F100 yields far-UV flux G0 and gas density n. The derived G0 scales as n to the power 1.4. We interpret this correlation as arising from Stromgren sphere scalings if much of the line and continuum luminosity arises near star-forming regions. The differences in G0 and n may be due to differences in the physical properties of the star-forming clouds.(Short abstract)

We present mid-infrared imaging and far-infrared (FIR) spectroscopy of 5 IBm galaxies observed by ISO as part of our larger study of the inter- stellar medium of galaxies. Most of the irregulars in our sample are very actively forming stars. The mid-infrared imaging was in a band centered at 6.75 microns that is dominated by polycyclic aromatic hydrocarbons (PAHs) and in a band centered at 15 microns that is dominated by small dust grains. The spectroscopy of 3 of the galaxies includes [CII]158 microns and [OI]63 microns, important coolants of photodissociation regions (PDRs), and [OIII]88 microns and [NII]122 microns, which come from ionized gas. [OI]145 microns and [OIII]52 microns were measured in one galaxy as well. These data are combined with PDR and HII region models to deduce properties of the interstellar medium of these galaxies.

We mapped a 13 by 22 pc region in emission from 492 GHz [C I] and, for the first time, 551 GHz ^{13}CO(5-4) in the giant molecular cloud M17SW, using the Submillimeter Wave Astronomy Satellite. The morphologies of the [C I] and ^{13}CO emission are strikingly similar. The extent and intensity of the [C I] and ^{13}CO(5-4) emission is explained as arising from photodissociation regions on the surfaces of embedded molecular clumps. Modeling of the ^{13}CO(5-4) emission in comparison to ^{13}CO(1-0) indicates a temperature gradient across the cloud, peaking to at least 63 K near the M17 ionization front and decreasing to at least 20 K at the western edge of the cloud. We see no correlation between gas density and column density. The beam-averaged column density of C I in the core is 1x10^{18} cm^-2, and the mean column density ratio N(C I)/N(CO) is about 0.4. The variations of N(C I)/N(CO) with position in M17SW indicate a similar clump size distribution throughout the cloud.

Hubble Space Telescope images of the galaxies NGC 2207 and IC 2163 show star formation and dust structures in a system that has experienced a recent grazing encounter. Tidal forces from NGC 2207 compressed and elongated the disk of IC 2163, forming an oval ridge of star formation. Gas flowing away from this ridge has thin parallel dust filaments transverse to the direction of motion. Numerical models suggest that the filaments come from flocculent spiral arms that were present before the interaction. A dust lane at the outer edge of the tidal tail is a shock front where the flow abruptly changes direction. A spiral arm of NGC 2207 that is backlit by IC 2163 is seen to contain several parallel, knotty filaments that are probably shock fronts in a density wave. Blue clusters of star formation inside these dust lanes show density wave triggering by local gravitational collapse. Spiral arms inside the oval of IC 2163 could be the result of ILR-related orbits in the tidal potential that formed the oval. Their presence suggests that tidal forces alone may initiate a temporary nuclear gas flow and eventual starburst without first forming a stellar bar. Several emission structures resembling jets 100-1000 pc long appear. There is a dense dark cloud with a conical shape 400 pc long and a bright compact cluster at the tip, and with a conical emission nebula of the same length that points away from the cluster in the other direction. This region coincides with a non-thermal radio continuum source that is 1000 times the luminosity of Cas A at 20 cm.

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.

Within Abell galaxy clusters containing wide-angle tailed radio sources, there is evidence of a ``prevailing wind'' which directs the WAT jets. We study the alignment of WAT jets and nearby clusters to test the idea that this wind may be a fossil of drainage along large-scale supercluster axes. We also test this idea with a study of the alignment of WAT jets and supercluster axes. Statistical test neighbours indicate no alignment of WAT jets towards nearest clusters, but do indicate approximately 98% confidence in alignment with the long axis of the supercluster in which the cluster lies. We find a preferred scale for such superclusters of order 25 Mpc $h^{-1}$.

Observations with the Hubble Space Telescope reveal an irregular network of dust spiral arms in the nuclear region of the interacting disk galaxy NGC 2207. The spirals extend from ~50 pc to ~300 pc in galactocentric radius, with a projected width of ~20 pc. Radiative transfer calculations determine the gas properties of the spirals and the inner disk, and imply a factor of ~4 local gas compression in the spirals. The gas is not strongly self-gravitating, nor is there a nuclear bar, so the spirals could not have formed by the usual mechanisms applied to main galaxy disks. Instead, they may result from acoustic instabilities that amplify at small galactic radii. Such instabilities may promote gas accretion into the nucleus.

Using the Long Wavelength Spectrometer (LWS) onboard the Infrared Space Observatory (ISO), we have observed thermal water vapor emission from a roughly circular field of view approximately 75 arc seconds in diameter centered on the Orion BN-KL region. The Fabry-Perot line strengths, line widths, and spectral line shifts observed in eight transitions between 71 and 125 microns show good agreement with models of thermal emission arising from a molecular cloud subjected to a magnetohydrodynamic C-type shock. Both the breadth and the relative strengths of the observed lines argue for emission from a shock rather than from warm quiescent gas in the Orion core. Though one of the eight transitions appears anomalously strong, and may be subject to the effects of radiative pumping, the other seven indicate an H2O/H2 abundance ratio of order 5E-4, and a corresponding gas-phase oxygen-to-hydrogen abundance ratio of order 4E-4. Given current estimates of the interstellar, gas-phase, oxygen and carbon abundances in the solar vicinity, this value is consistent with theoretical shock models that predict the conversion into water of all the gas-phase oxygen that is not bound as CO. The overall cooling provided by rotational transitions of H2O in this region appears to be comparable to the cooling through rotational lines of CO, but is an order of magnitude lower than cooling through H2 emission. However, the model that best fits our observations shows cooling by H2O and CO dominant in that portion of the post-shock region where temperatures are below ~ 800 K and neither vibrational nor rotational radiative cooling by H2 is appreciable.

We report about high spatial resolution observations made with the IRAM interferometer of the CO emission from a 1.1x1.1 kpc plane-of-sky field on a spiral arm of Messier 81. With a beam of 5''(90 pc), we identify 6 giant molecular cloud complexes with virial masses of about 10^6 Msun, including one associated with a giant HII region. The deduced N(H2)/I(CO) ratios are about 3 times larger on average than those measured near the solar neighborhood, suggesting that the complexes are not self-gravitationally bound except, possibly, for the complex associated with the giant HII region; they could be the average of several clouds of mass a few 10^5 Msun and diameter < 100 pc. The linewidths are very narrow with respect to the measured sizes, so that the size-linewidth relation for M 81 clouds is very different from that in the Milky Way. The narrow linewidths imply smaller virial masses than for Galactic complexes of the same size, and this is consistent with the weaker CO emission from the GMCs in M 81. We conclude from these observations that the molecular medium in M 81 differs from that in the Milky Way.