Aims. 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. Methods. A number of benchmark models have been created, covering low and high gas densities n = 10 3 , 10 5.5 cm −3 and far ultraviolet intensities χ = 10, 10 5 in units of the Draine field (FUV: 6 < h ν < 13.6 eV). 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 self consistently by solving the thermal balance, thus testing the modeling of the heating and cooling mechanisms accounting for the detailed energy balance throughout the clouds. Results. 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 such as the effect of PAHs on the electron density or the temperature dependence of the dissociation of CO by cosmic ray induced secondary photons, 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. Nevertheless, the remaining spread in the computed observables such as the atomic fine-structure line intensities serves as a warning that there is still a considerable uncertainty when interpreting astronomical data with our models.

We present a comprehensive analysis of a broad band spectral line survey of the Orion Kleinmann-Low nebula (Orion KL), one of the most chemically rich regions in the Galaxy, using the HIFI instrument on board the Herschel Space Observatory. This survey spans a frequency range from 480 to 1907 GHz at a resolution of 1.1 MHz. These observations thus encompass the largest spectral coverage ever obtained toward this high-mass star-forming region in the sub-mm with high spectral resolution, and include frequencies > 1 THz where the Earth's atmosphere prevents observations from the ground. In all, we detect emission from 39 molecules (79 isotopologues). Combining this dataset with ground based mm spectroscopy obtained with the IRAM 30 m telescope, we model the molecular emission from the mm to the far-IR using the XCLASS program which assumes local thermodynamic equilibrium (LTE). Several molecules are also modeled with the MADEX non-LTE code. Because of the wide frequency coverage, our models are constrained by transitions over an unprecedented range

We report observations of three rotational transitions of molecular oxygen (O 2 ) in emission from the H 2 Peak

Photodissociation regions (PDRs) define the transition zone between an ionized and a dark molecular region. They consist of neutral gas which interacts with far-ultraviolet radiation and are characterized by strong infrared line emission. Various numerical codes treating onedimensional PDRs have been developed in the past, simulating the complexity of chemical reactions occurring and providing a better understanding of the structure of a PDR. In this paper we present the three-dimensional code, 3D-PDR, which can treat PDRs of arbitrary density distribution. The code solves the chemistry and the thermal balance self-consistently within a given three-dimensional cloud. It calculates the total heating and cooling functions at any point in a given PDR by adopting an escape probability method. It uses a HEALPIxbased ray tracing scheme to evaluate the attenuation of the far-ultraviolet radiation in the PDR and the propagation of the far-infrared/submm line emission out of the PDR. We present benchmarking results and apply 3D-PDR to (i) a uniform-density spherical cloud interacting with a plane-parallel external radiation field, (ii) a uniform-density spherical cloud interacting with a two-component external radiation field and (iii) a cometary globule interacting with a planeparallel external radiation field. We find that the code is able to reproduce the benchmarking results of various other one-dimensional numerical codes treating PDRs. We also find that the accurate treatment of the radiation field in the fully three-dimensional treatment of PDRs can in some cases leads to different results when compared to a standard one-dimensional treatment.

We present SHARC-2 350 µm data on 20 luminous z ∼ 2 starbursts with S 1.2mm >2 mJy from the Spitzerselected samples of Lonsdale et al. and Fiolet et al. All the sources were detected, with S 350µm >25 mJy for 18 of them. With the data, we determine precise dust temperatures and luminosities for these galaxies using both single-temperature fits and models with power-law mass-temperature distributions. We derive appropriate formulae to use when optical depths are non-negligible. Our models provide an excellent fit to the 6 µm-2 mm measurements of local starbursts. We find characteristic single-component temperatures T 1 ≃35.5±2.2 K and integrated infrared (IR) luminosities around 10 12.9±0.1 L ⊙ for the SWIRE-selected sources. Molecular gas masses are estimated at ≃4×10 10 M ⊙ , assuming κ 850µm =0.15 m 2 kg −1 and a submillimeter-selected galaxy (SMG)-like gas-to-dust mass ratio. The best-fit models imply 2 kpc emission scales. We also note a tight correlation between rest-frame 1.4 GHz radio and IR luminosities confirming star formation as the predominant power source. The far-IR properties of our sample are indistinguishable from the purely submillimeter-selected populations from current surveys. We therefore conclude that our original selection criteria, based on mid-IR colors and 24 µm flux densities, provides an effective means for the study of SMGs at z∼1.5-2.5.

Residual star formation at late times in early-type galaxies and their progenitors must be suppressed in order to explain the population of red, passively evolving systems we see today. Likewise, residual or newly accreted reservoirs of molecular gas that are fueling star formation must be destroyed. This suppression of star formation in early-type galaxies is now commonly attributed to active galactic nucleus (AGN) feedback wherein the reservoir of gas is heated and expelled during a phase of accretion onto the central supermassive black hole. However, direct observational evidence for a link between the destruction of this molecular gas and an AGN phase has been missing so far. We present new mm-wavelength observations from the IRAM 30 m telescope of a sample of lowredshift SDSS early-type galaxies currently undergoing this process of quenching of late-time star formation. Our observations show that the disappearance of the molecular gas coincides within less than 100 Myr with the onset of accretion onto the black hole and is too rapid to be due to star formation alone. Since our sample galaxies are not associated to powerful quasar activity or radio jets, we conclude that low-luminosity AGN episodes are sufficient to suppress residual star formation in early-type galaxies. This "suppression mode" of AGN feedback is very different from the "truncation mode" linked to powerful quasar activity during early phases of galaxy formation.

High resolution line spectra of star-forming regions are mines of information: they provide unique clues to reconstruct the chemical, dynamical, and physical structure of the observed source. We present the first results from the Herschel key project "Chemical HErschel Surveys of Star forming regions", CHESS. We report and discuss observations towards five CHESS targets, one outflow shock spot and four protostars with luminosities bewteen 20 and 2 × 10 5 L : L1157-B1, IRAS 16293-2422, OMC2-FIR4, AFGL 2591, and NGC 6334I. The observations were obtained with the heterodyne spectrometer HIFI on board Herschel, with a spectral resolution of 1 MHz. They cover the frequency range 555−636 GHz, a range largely unexplored before the launch of the Herschel satellite. A comparison of the five spectra highlights spectacular differences in the five sources, for example in the density of methanol lines, or the presence/absence of lines from S-bearing molecules or deuterated species. We discuss how these differences can be attributed to the different star-forming mass or evolutionary status.

We report the detection of absorption lines by the reactive ions OH + , H 2 O + and H 3 O + along the line of sight to the submillimeter continuum source G10.6−0.4 (W31C). We used the Herschel HIFI instrument in dual beam switch mode to observe the ground state rotational transitions of OH + at 971 GHz, H 2 O + at 1115 and 607 GHz, and H 3 O + at 984 GHz. The resultant spectra show deep absorption over a broad velocity range that originates in the interstellar matter along the line of sight to G10.6−0.4 as well as in the molecular gas directly associated with that source. The OH + spectrum reaches saturation over most velocities corresponding to the foreground gas, while the opacity of the H 2 O + lines remains lower than 1 in the same velocity range, and the H 3 O + line shows only weak absorption. For LSR velocities between 7 and 50 km s −1 we estimate total column densities of N(OH + ) ≥ 2.5 × 10 14 cm −2 , N(H 2 O + ) ∼6 × 10 13 cm −2 and N(H 3 O + ) ∼4.0 × 10 13 cm −2 . These detections confirm the role of O + and OH + in initiating the oxygen chemistry in diffuse molecular gas and strengthen our understanding of the gas phase production of water. The high ratio of the OH + by the H 2 O + column density implies that these species predominantly trace low-density gas with a small fraction of hydrogen in molecular form.