M33 is a gas rich spiral galaxy of the Local Group. We investigate the relationship between the two major gas cooling lines and the total infrared (TIR) dust continuum. We mapped the emission of gas and dust in M33 using the far-infrared lines of [CII] and [OI](63um) and the TIR. The line maps were observed with Herschel/PACS. These maps have 50pc resolution and form a ~370pc wide stripe along its major axis covering the sites of bright HII regions, but also more quiescent arm and inter-arm regions from the southern arm at 2kpc galacto-centric distance to the south out to 5.7kpc distance to the north. Full-galaxy maps of the continuum emission at 24um from Spitzer/MIPS, and at 70um, 100um, and 160um from PACS were combined to obtain a map of the TIR. TIR and [CII] intensities are correlated over more than two orders of magnitude. The range of TIR translates to a range of far ultraviolet (FUV) emission of G0,obs~2 to 200 in units of the average Galactic radiation field. The binned [CII]/TIR ratio drops with rising TIR, with large, but decreasing scatter. Fits of modified black bodies (MBBs) to the continuum emission were used to estimate dust mass surface densities and total gas column densities. A correction for possible foreground absorption by cold gas was applied to the [OI] data before comparing it with models of photon dominated regions (PDRs). Most of the ratios of [CII]/[OI] and ([CII]+[OI])/TIR are consistent with two model solutions. The median ratios are consistent with one solution at n~2x10^2 cm-3, G0~60, and and a second low-FUV solution at n~10^4 cm-3, G0~1.5. The bulk of the gas along the lines-of-sight is represented by a low-density, high-FUV phase with low beam filling factors ~1. A fraction of the gas may, however, be represented by the second solution.

The central area (40"x40") of the bipolar nebula S106 was mapped in the OI line at 63.2 micron with high angular (6") and spectral resolution, using GREAT on board SOFIA. The OI emission distribution is compared to the CO 16-15, CII 158 micron, and CO 11-10 lines, mm-molecular lines, and continuum. It is composed of several velocity components in the range from -30 km/s to 25 km/s. The high-velocity blue- and redshifted emission can be explained as arising from accelerated photodissociated (PDR) gas associated with a dark lane close to the massive binary system S106 IR, and from shocks caused by the stellar wind and/or a disk--envelope interaction. At velocities from -9 to -4 km/s and 0.5 to 8 km/s line wings are observed that we attribute to cooling in PDRs created by the ionizing radiation impinging on the cavity walls. The bulk velocity range is dominated by PDR emission from the clumpy molecular cloud. Modelling the emission in the different velocity ranges with the KOSMA-tau code constrains a radiation field chi of a few times 10^4 and densities n of a few times 10^4 cm^-3. Considering self-absorption of the OI line results in higher densities (up to 10^6 cm^-3) only for the gas component seen at high blue- and red velocities. The dark lane has a mass of 275 Msun and shows a velocity difference of 1.4 km/s along its projected length of 1 pc, determined from H13CO+ 1-0 mapping. It can be interpreted as a massive accretion flow, or the remains of it, linked to S106 IR/FIR. The most likely explanation is that the binary system is at a stage of its evolution where gas accretion is counteracted by the stellar winds and radiation, leading to the very complex observed spatial and kinematic emission distribution of the various tracers.

In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J lines that trace warm molecular gas at PDR interfaces. We observed the 12CO and 13CO ladders in two prototypical PDRs, the Orion Bar and NGC 7023 NW using the instruments onboard Herschel. We also considered line emission from key species in the gas cooling of PDRs (C+, O, H2) and other tracers of PDR edges such as OH and CH+. All the intensities are collected from Herschel observations, the literature and the Spitzer archive and are analyzed using the Meudon PDR code. A grid of models was run to explore the parameter space of only two parameters: thermal gas pressure and a global scaling factor that corrects for approximations in the assumed geometry. We conclude that the emission in the high-J CO lines, which were observed up to Jup=23 in the Orion Bar (Jup=19 in NGC7023), can only originate from small structures of typical thickness of a few 1e-3 pc and at high thermal pressures (Pth~1e8 K cm-3). Compiling data from the literature, we found that the gas thermal pressure increases with the intensity of the UV radiation field given by G0, following a trend in line with recent simulations of the photoevaporation of illuminated edges of molecular clouds. This relation can help rationalising the analysis of high-J CO emission in massive star formation and provides an observational constraint for models that study stellar feedback on molecular clouds.

Photon Dominated Regions (PDRs) are interfaces between the mainly ionized and mainly molecular material around young massive stars. Analysis of the physical and chemical structure of such regions traces the impact of far-ultraviolet radiation of young massive stars on their environment. We present results on the physical and chemical structure of the prototypical high UV-illumination edge-on Orion Bar PDR from an unbiased spectral line survey with a wide spectral coverage. A spectral scan from 480-1250 GHz and 1410-1910 GHz at 1.1 MHz resolution was obtained by the HIFI instrument onboard the Herschel Space Observatory. For molecules with multiple transitions we used rotational diagrams to obtain excitation temperatures and column densities. For species with a single detected transition we used an optically thin LTE approximation. In case of species with available collisional rates, we also performed a non-LTE analysis to obtain kinetic temperatures, H2 volume densities, and column densities. About 120 lines corresponding to 29 molecules (including isotopologues) have been detected in the Herschel/HIFI line survey, including 11 transitions of CO, 7 transitions of 13CO, 6 transitions of C18O, 10 transitions of H2CO, and 6 transitions of H2O. Most species trace kinetic temperatures in the range between 100 and 150 K and H2 volume densities in the range between 10^5 and 10^6 cm^-3. The species with temperatures and / or densities outside of this range include the H2CO transitions tracing a very high temperature (315 K) and density (1.4x10^6 cm^-3) component and SO corresponding to the lowest temperature (56 K) measured as a part of this line survey. The observed lines/species reveal a range of physical conditions (gas density /temperature) involving structures at high density / high pressure, obsoleting the traditional 'clump / interclump' picture of the Orion Bar.

The radiative feedback of massive stars on molecular clouds creates pillars, globules and other features at the interface between the HII region and molecular cloud. We present here Herschel observations between 70 and 500 micron of the immediate environment of the Cygnus OB2 association, performed within the HOBYS program. All structures were detected based on their appearance at 70 micron, and have been classified as pillars, globules, evaporating gasous globules (EGGs), proplyd-like objects, and condensations. From the 70 and 160 micron flux maps, we derive the local FUV field on the PDR surfaces. In parallel, we use a census of the O-stars to estimate the overall FUV-field, that is 10^3-10^4 G_0 close to the central OB cluster (within 10 pc) and decreases down to a few tens G_0, in a distance of 50 pc. From a SED fit to the four longest Herschel wavelengths, we determine column density and temperature maps and derive masses, volume densities and surface densities for these structures. We find that the morphological classification corresponds to distinct physical properties. Pillars and globules have the longest estimated photoevaporation lifetimes, a few 10^6 yr, while all other features should survive less than that. These lifetimes are consistent with that found in simulations of turbulent, UV-illuminated clouds. We propose a tentative evolutionary scheme in which pillars can evolve into globules, which in turn then evolve into EGGs, condensations and proplyd-like objects.

We aim to understand the contribution of the ionized, atomic and molecular phases of the ISM to the [CII] emission from clouds near the dynamical center and the BCLMP302 HII region in the north of the nearby galaxy M33 at a spatial resolution of 50pc. We combine high resolution [CII] spectra taken with the HIFI spectrometer onboard the Herschel satellite with [CII] Herschel-PACS maps and ground-based observations of CO(2-1) and HI. All data are at a common spatial resolution of 50pc. Typically, the [CII] lines have widths intermediate between the narrower CO(2-1) and broader HI line profiles. We decomposed the [CII] spectra in terms of contribution from molecular and atomic gas detected in CO(2-1) and HI, respectively. We find that the relative contribution of molecular and atomic gas traced by CO(2-1) and HI varies depends mostly on the local physical conditions and geometry. We estimate that 11-60% and 5-34% of the [CII] intensities in the center and in BCLMP302, respectively, arise at velocities showing no CO(2-1) or HI emission and could arise in CO-dark molecular gas. The deduced strong variation in the [CII] emission not associated with CO and HI cannot be explained in terms of differences in A_v, far-ultraviolet radiation field, and metallicity between the two studied regions. Hence the relative amounts of diffuse (CO-dark) and dense molecular gas possibly vary on spatial scales smaller than 50pc. Based on the emission measure observed at radio wavelengths we estimate the contribution of ionized gas at a few positions to lie between 10-25%. The correlations between the intensities of tracers corresponding to the same velocity range as [CII], differ from the correlation derived from PACS data. The results in this paper emphasize the need for velocity-resolved observations to discern the contribution of different components of the ISM to [CII] emission. (abridged)

We study the wavelength dependence of the dust emission as a function of position and environment across the disk of M33 at a linear resolution of 160 pc using Spitzer and Herschel photometric data. Expressing the emissivity of the dust as a power law, the power-law exponent (beta) is estimated from two independent approaches designed to properly treat the degeneracy between beta and the dust temperature. Both beta and the dust temperature are higher in the inner disk than in the outer disk, contrary to reported beta-T anti-correlations found in other sources. In the cold + warm dust model, the warm component and the ionized gas (Halpha) have a very similar distribution across the galaxy, demonstrating that the model separates the components in an appropriate fashion. The flocculent spiral arms and the dust lanes are evident in the map of the cold component. Both cold and warm dust column densities are high in star forming regions and reach their maxima toward the giant star forming complexes NGC604 and NGC595. beta declines from close to 2 in the center to about 1.3 in the outer disk. beta is positively correlated with star formation and with molecular gas column, as traced by Halpha and CO emission. The lower dust emissivity index in the outer parts of M33 is likely related to the reduced metallicity (different grain composition) and possibly different size distribution. It is not due to the decrease in stellar radiation field or temperature in a simple way because the FIR-bright regions in the outer disk also have a low beta. Like most spirals, M33 has a (decreasing) radial gradient in star formation and molecular-to-atomic gas ratio such that the regions bright in Halpha or CO tend to trace the inner disk, making it difficult to distinguish between their effects on the dust.

We aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis of M33, in combination with Herschel PACS and SPIRE continuum maps, IRAM 30m CO 2-1 and VLA HI data to study the variation of velocity integrated intensities. The ratio of [CII] emission over the far-infrared continuum is used as a proxy for the heating efficiency, and models of photon-dominated regions are used to study the local physical densities, FUV radiation fields, and average column densities of the molecular clouds. The heating efficiency stays constant at 0.8% in the inner 4.5kpc radius of the galaxy where it starts to increase to reach values of ~3% in the outskirts at about 6kpc radial distance. The rise of efficiency is explained in the framework of PDR models by lowered volume densities and FUV fields, for optical extinctions of only a few magnitudes at constant metallicity. In view of the significant fraction of HI emission stemming from PDRs, and for typical pressures found in the Galactic cold neutral medium (CNM) traced by HI emission, the CNM contributes ~15% to the observed [CII] emission in the inner 2kpc radius of M33. The CNM contribution remains largely undetermined in the south, while positions between 2 and 7.3kpc radial distance in the north of M33 show a contribution of ~40%+-20%.

The abundances of interstellar CH+ and SH+ are not well understood as their most likely formation channels are highly endothermic. Using data from Herschel, we study the formation of CH+ and SH+ in a typical high UV-illumination photon-dominated region (PDR), the Orion Bar. Herschel/HIFI provides velocity-resolved data of CH+ 1-0 and 2-1 and three hyperfine transitions of SH+. Herschel/PACS provides information on the excitation and spatial distribution of CH+ (up to J=6-5). The widths of the CH+ 2-1 and 1-0 transitions are of ~5 km/s, significantly broader than the typical width of dense gas tracers in the Orion Bar (2-3 km/s) and are comparable to the width of tracers of the interclump medium such as C+ and HF. The detected SH+ transitions are narrower compared to CH+ and have line widths of 3 km/s, indicating that SH+ emission mainly originates in denser condensations. Non-LTE radiative transfer models show that electron collisions affect the excitation of CH+ and SH+, and that reactive collisions need to be taken into account to calculate the excitation of CH+. Comparison to PDR models shows that CH+ and SH+ are tracers of the warm surface region (AV<1.5) of the PDR with temperatures between 500-1000 K. We have also detected the 5-4 transition of CF+ (FWHM=1.9 km/s) with an intensity that is consistent with previous observations of lower-J CF+ transitions toward the Orion Bar. A comparison to PDR models indicate that the internal vibrational energy of H2 can explain the formation of CH+ for typical physical conditions in the Orion Bar near the ionization front. H2 vibrational excitation is the most likely explanation of SH+ formation as well. The abundance ratios of CH+ and SH+ trace the destruction paths of these ions, and through that, indirectly, the ratios of H, H2 and electron abundances as a function of depth into the cloud.

(Abridged) The 30 Dor region in the Large Magellanic Cloud (LMC) is the most vigorous star-forming region in the Local Group. Star formation in this region is taking place in low-metallicity molecular gas that is exposed to an extreme far--ultraviolet (FUV) radiation field powered by the massive compact star cluster R136. We used the NANTEN2 telescope to obtain high-angular resolution observations of the 12CO 4-3, 7-6, and 13CO 4-3 rotational lines and [CI] 3P1-3P0 and 3P2-3P1 fine-structure submillimeter transitions in 30Dor-10, the brightest CO and FIR-emitting cloud at the center of the 30Dor region. We derived the properties of the low-metallicity molecular gas using an excitation/radiative transfer code and found a self-consistent solution of the chemistry and thermal balance of the gas in the framework of a clumpy cloud PDR model. We compared the derived properties with those in the N159W region, which is exposed to a more moderate far-ultraviolet radiation field compared with 30Dor-10, but has similar metallicity. We also combined our CO detections with previously observed low-J CO transitions to derive the CO spectral-line energy distribution in 30Dor-10 and N159W. The separate excitation analysis of the submm CO lines and the neutral carbon fine structure lines shows that the mid-J CO and [CI]-emitting gas in the 30Dor-10 region has a temperature of about 160 K and a H2 density of about 10^4 cm^-3. We find that the molecular gas in 30Dor-10 is warmer and has a lower beam filling factor compared to that of N159W, which might be a result of the effect of a strong FUV radiation field heating and disrupting the low--metallicity molecular gas. We use a clumpy PDR model (including the [CII] line intensity reported in the literature) to constrain the FUV intensity to about chi_0 ~ 3100 and an average total H density of the clump ensemble of about 10^5 cm^-3 in 30Dor-10.

Power spectra of de-projected images of late-type galaxies in gas and/or dust emission are very useful diagnostics of the dynamics and stability of their interstellar medium. Previous studies have shown that the power spectra can be approximated as two power-laws, a shallow one at large scales (larger than 500 pc) and a steeper one at small scales, with the break between the two corresponding to the line-of-sight thickness of the galaxy disk. We present a thorough analysis of the power spectra of the dust and gas emission at several wavelengths in the nearby galaxy M33. In particular, we use the recently obtained images at five wavelengths by PACS and SPIRE onboard Herschel. The large dynamical range (2-3 dex in scale) of most images allow us to determine clearly the change in slopes from -1.5 to -4, with some variations with wavelength. The break scale is increasing with wavelength, from 100 pc at 24 and 100micron to 350 pc at 500micron, suggesting that the cool dust lies in a thicker disk than the warm dust, may be due to star formation more confined to the plane. The slope at small scale tends to be steeper at longer wavelength, meaning that the warmer dust is more concentrated in clumps. Numerical simulations of an isolated late-type galaxy, rich in gas and with no bulge, like M33, are carried out, in order to better interpret these observed results. Varying the star formation and feedback parameters, it is possible to obtain a range of power-spectra, with two power-law slopes and breaks, which nicely bracket the data. The small-scale power-law is indeed reflecting the 3D behaviour of the gas layer, steepening strongly while the feedback smoothes the structures, by increasing the gas turbulence. M33 appears to correspond to a fiducial model with an SFR of $\sim$ 0.7 Mo/yr, with 10% supernovae energy coupled to the gas kinematics.

The clumpy density structure of photon-dominated regions is well established, but the physical properties of the clumps and of the surrounding interclump medium are only approximately known. The aim of this paper is to constrain the physical and chemical conditions in the Orion Bar, a prototypical nearby photon-dominated region. We present observations of the HF J=1-0 line, which appears in emission toward the Orion Bar, and compare the brightness of the line to non-LTE radiative transfer calculations. The large width of the HF line suggests an origin of the emission in the interclump gas, but collisional excitation by H2 in the interclump gas underpredicts the observed line intensity by factors of 3-5. In contrast, an origin of the line in the dense clumps requires a density of ~10^9 cm^-3, 10-100 times higher than previous estimates, which is unlikely. However, electron impact excitation reproduces our observations for T = 100 K and n(e) = 10 cm^-3, as expected for the interclump gas. We conclude that HF emission is a signpost of molecular gas with a high electron density. Similar conditions may apply to active galactic nuclei where HF also appears in emission.

Context: The emission line of [CII] at 158 micron is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. Aims: Disentangling the relative contributions of the different ISM phases to [CII] emission, is a major topic of the HerM33es program, a Herschel key project to study the ISM in the nearby spiral galaxy M33. Methods: Using PACS, we have mapped the emission of [CII] 158 micron, [OI] 63 micron, and other FIR lines in a 2'x2' region of the northern spiral arm of M33, centered on the HII region BCLMP302. At the peak of H-alpha emission, we have observed in addition a velocity resolved [CII] spectrum using HIFI. We use scatterplots to compare these data with PACS 160 micron continuum maps, and with maps of CO and HI data, at a common resolution of 12 arcsec or 50 pc. Maps of H-alpha and 24 micron emission observed with Spitzer are used to estimate the SFR. We have created maps of the [CII] and [OI] 63 micron emission and detected [NII] 122 micron and NIII 57 micron at individual positions. Results: The [CII] line observed with HIFI is significantly broader than that of CO, and slightly blue-shifted. In addition, there is little spatial correlation between [CII] observed with PACS and CO over the mapped region. There is even less spatial correlation between [CII] and the atomic gas traced by HI. Detailed comparison of the observed intensities towards the HII region with models of photo ionization and photon dominated regions, confirms that a significant fraction, 20--30%, of the observed [CII] emission stems from the ionized gas and not from the molecular cloud. The gas heating efficiency, using the ratio between [CII] and the TIR as a proxy, varies between 0.07 and 1.5%, with the largest variations found outside the HII region.

We present observations of the rotational ortho-water ground transition, the two lowest para-water transitions, and the ground transition of ionised ortho-water in the archetypal starburst galaxy M82, performed with the HIFI instrument on the Herschel Space Observatory. These observations are the first detections of the para-H2O(111-000) (1113\,GHz) and ortho-H2O+(111-000) (1115\,GHz) lines in an extragalactic source. All three water lines show different spectral line profiles, underlining the need for high spectral resolution in interpreting line formation processes. Using the line shape of the para-H2O(111-000) and ortho-H2O+(111-000) absorption profile in conjunction with high spatial resolution CO observations, we show that the (ionised) water absorption arises from a ~2000 pc^2 region within the HIFI beam located about ~50 pc east of the dynamical centre of the galaxy. This region does not coincide with any of the known line emission peaks that have been identified in other molecular tracers, with the exception of HCO. Our data suggest that water and ionised water within this region have high (up to 75%) area-covering factors of the underlying continuum. This indicates that water is not associated with small, dense cores within the ISM of M82 but arises from a more widespread diffuse gas component.

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations were carried out using the HERA at the 30m IRAM telescope in the CO(2-1) line achieving a resolution of 12"x2.6 km/s, enabling individual GMCs to be resolved. The observed region mainly along the major axis out to a radius of 8.5 kpc, and covers the strip observed with HIFI/PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK at 2.6 km/s velocity resolution. The CO(2-1) luminosity of the observed region is 1.7\pm0.1x10^7 Kkm/s pc^2, corresponding to H2 masses of 1.9x10^8 Msun (including He), calculated with a NH2/ICO twice the Galactic value due to the half-solar metallicity of M33. HI 21 cm VLA archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in CASA, yielding a series of datacubes with resolutions ranging from 5" to 25". The HI mass within a radius of 8.5 kpc is estimated to be 1.4x10^9 Msun. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9\pm0.1 kpc whereas the atomic gas surface density is constant at Sigma_HI=6\pm2 Msun/pc^2 deprojected to face-on. The central kiloparsec H_2 surface density is Sigma_H2=8.5\pm0.2 Msun/pc^2. The star formation rate per unit molecular gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to Halpha and FIR brightness, is constant with radius. The SFE appears 2-4 times greater than of large spiral galaxies. A morphological comparison of molecular and atomic gas with tracers of star formation shows good agreement between these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated more than 8 kpc from the galaxy center, are presented.

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $\sigma$ using SPIRE and recent CO and \HI\ observations; a variation in $\sigma$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33.

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

Almost all properties of a photodissociation region (PDR) depend on its metallicity. The heating and cooling efficiencies that determine the temperature of the gas and dust, the dust composition, as well as the elemental abundances that influence the chemical structure of the PDR are just three examples that demonstrate the importance of metallicity effects in PDRs. PDRs are often associated with sites of star formation. If we want to understand the star formation history of our own Galaxy and of distant low-metallicity objects we need to understanding how metallicity acts on PDR physics and chemistry.

Studying molecular gas in the central regions of the star burst galaxies NGC4945 and Circinus enables us to characterize the physical conditions and compare them to previous local and high-z studies. We estimate temperature, molecular density and column densities of CO and atomic carbon. Using model predictions we give a range of estimated CO/C abundance ratios. Using the new NANTEN2 4m sub-millimeter telescope in Pampa La Bola, Chile, we observed for the first time CO 4-3 and [CI] 3P1-3 P0 at the centers of both galaxies at linear scale of 682 pc and 732 pc respectively. We compute the cooling curves of 12CO and 13CO using radiative transfer models and estimate the physical conditions of CO and [CI]. The centers of NGC4945 and Circinus are very [CI] bright objects, exhibiting [CI] 3P1 - 3 P0 luminosities of 91 and 67Kkms-1kpc2, respectively. The [CI] 3P1-3 P0/CO 4-3 ratio of integrated intensities are large at 1.2 in NGC4945 and 2.8 in Circinus. Combining previous CO J= 1-0, 2-1 and 3-2 and 13CO J= 1-0, 2-1 studies with our new observations, the radiative transfer calculations give a range of densities, n(H2) = 10^3-3*104^cm-3, and a wide range of kinetic temperatures, Tkin = 20 - 100K, depending on the density. Future CO J= 7-6 and [CI] 2-1 observations will be important to resolve the ambiguity in the physical conditions and confirm the model predictions.

The Carina region is an excellent astrophysical laboratory for studying the feedback mechanisms of newly born, very massive stars within their natal giant molecular clouds (GMCs) at only 2.35 kpc distance. We use a clumpy PDR model to analyse the observed intensities of atomic carbon and CO and to derive the excitation conditions of the gas. The NANTEN2-4m submillimeter telescope was used to map the [CI] 3P1-3P0, 3P2-3P1 and CO 4-3, 7-6 lines in two 4'x4' regions of Carina where molecular material interfaces with radiation from the massive star clusters. One region is the northern molecular cloud near the compact OB cluster Tr14, and the second region is in the molecular cloud south of etaCar and Tr16. These data were combined with 13CO SEST spectra, HIRES/IRAS 60um and 100um maps of the FIR continuum, and maps of 8um IRAC/Spitzer and MSX emission. We used the HIRES far-infrared dust data to create a map of the FUV field heating the gas. The northern region shows an FUV field of a few 1000 in Draine units while the field of the southern region is about a factor 10 weaker. We constructed models consisting of an ensemble of small spherically symmetric PDR clumps within the 38" beam (0.43pc), which follow canonical power-law mass and mass-size distributions. We find that an average local clump density of 2x10**5 cm-3 is needed to reproduce the observed line emission at two selected interface positions. Stationary, clumpy PDR models reproduce the observed cooling lines of atomic carbon and CO at two positions in the Carina Nebula.

To study the properties of the interstellar medium in the prototypical merging system of the Antennae galaxies (NGC 4038 and NGC 4039), we have obtained CO(1-0), (2-1) and (3-2) line maps, as well as a map of the 870 micron continuum emission. Our results are analysed in conjunction with data from X-ray to radio wavelengths. In order to distinguish between exact coincidence and merely close correspondence of emission features, we compare the morphological structure of the different emission components at the highest available angular resolution. To constrain the physical state of the molecular gas, we apply models of photon dominated regions (PDRs) that allow us to fit CO and [CII] data, as well as other indicators of widespread PDRs in the Antennae system, particularly within the super giant molecular cloud (SGMC) complexes of the interaction region (IAR) between the two galaxies. The modeled clouds have cores with moderately high gas densities up to 4 10^4 / cm^3 and rather low kinetic temperatures <=25K). At present, all these clouds, including those near the galactic nuclei, show no signs of intense starburst activity. Thermal radio or mid-infrared emission are all observed to peak slightly offset from the molecular peaks. The total molecular gas mass of the Antennae system adds up to \~10^10 M_sun. In the vicinity of each galactic nucleus, the moleculargas mass, 1-2 10^9 M_sun, exceeds that of the Galactic centre region by a factorof almost 100. Furthermore, the gas does not seem to deviate much from the N_{H_2}/I_CO ratio typical of the disk of our Galaxy rather than our Galactic centre.

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.

Stars form in the densest, coldest, most quiescent regions of molecular clouds. Molecules provide the only probes which can reveal the dynamics, physics, chemistry and evolution of these regions, but our understanding of the molecular inventory of sources and how this is related to their physical state and evolution is rudimentary and incomplete. The Spectral Legacy Survey (SLS) is one of seven surveys recently approved by the JCMT Board. Starting in 2007, the SLS will produce a spectral imaging survey of the content and distribution of all the molecules detected in the 345 GHz atmospheric window (between 332 GHz and 373 GHz) towards a sample of 5 sources. Our intended targets are: a low mass core (NGC1333 IRAS4), 3 high mass cores spanning a range of star forming environments and evolutionary states (W49, AFGL2591, and IRAS20126), and a PDR (the Orion Bar). The SLS will use the unique spectral imaging capabilities of HARP-B/ACSIS to study the molecular inventory and the physical structure of these objects, which span different evolutionary stages and physical environments, to probe their evolution during the star formation process. As its name suggests, the SLS will provide a lasting data legacy from the JCMT that is intended to benefit the entire astronomical community. As such, the entire data set (including calibrated spectral datacubes, maps of molecular emission, line identifications, and calculations of the gas temperature and column density) will be publicly available.

Aim: The aim of the paper is to understand the emission from the photon dominated regions in Cepheus B, estimate the column densities of neutral carbon in bulk of the gas in Cepheus B and to derive constraints on the factors which determine the abundance of neutral carbon relative to CO. Methods: This paper presents 15'x15' fully sampled maps of CI at 492 GHz and 12CO 4-3 observed with KOSMA at 1' resolution. The new observations have been combined with the FCRAO 12CO 1-0, IRAM-30m 13CO 2-1 and C18O 1-0 data, and far-infrared continuum data from HIRES/IRAS. The KOSMA-tau spherical PDR model has been used to understand the CI and CO emission from the PDRs in Cepheus B and to explain the observed variation of the relative abundances of both C^0 and CO. Results: The emission from the PDR associated with Cepheus B is primarily at V_LSR between -14 and -11 km s^-1. We estimate about 23% of the observed CII emission from the molecular hotspot is due to the ionized gas in the HII region. Over bulk of the material the C^0 column density does not change significantly, (2.0+-1.4)x10^17 cm^-2, although the CO column density changes by an order of magnitude. The observed \cbyco abundance ratio varies between 0.06 and 4 in Cepheus B. We find an anti-correlation of the observed C/CO abundance ratio with the observed hydrogen column density, which holds even when all previous observations providing C/CO ratios are included. Here we show that this observed variation of C/CO abundance with total column density can be explained only by clumpy PDRs consisting of an ensemble of clumps. At high H2 column densities high mass clumps, which exhibit low C/CO abundance, dominate, while at low column densities, low mass clumps with high C/CO abundance dominate.