We present results from Atacama Large Millimeter/submillimeter Array (ALMA) observations of CS from the nearby galaxy NGC 1068 ($\sim14$ Mpc). This Seyfert 2 barred galaxy possesses a circumnuclear disc (CND, $r\sim200$ pc) and a starburst ring (SB ring, $r\sim1.3$ kpc). These high-resolution maps ($\sim0.5$", $\sim35$ pc) allow us to analyse specific sub-regions in the galaxy and investigate differences in line intensity ratios and physical conditions, particularly those between the CND and SB ring. Local thermodynamic equilibrium (LTE) analysis of the gas is used to calculate CS densities in each sub-region, followed by non-LTE analysis conducted using the radiative transfer code RADEX to fit observations and constrain gas temperature, CS column density and hydrogen density. Finally, the chemical code UCLCHEM is used to reconstruct the gas, allowing an insight into its origin and chemical history. The density of hydrogen in the CND is found to be $\geq10^5$ cm$^{-2}$, although exact values vary, reaching $10^6$ cm$^{-2}$ at the AGN. The conditions in the two arms of the SB ring appear similar to one another, though the density found ($\sim10^4$ cm$^{-2}$) is lower than in the CND. The temperature in the CND increases from east to west, and is also overall greater than found in the SB ring. These modelling methods indicate the requirement for multi-phase gas components in order to fit the observed emission over the galaxy. A larger number of high resolution transitions across the SLED may allow for further constraining of the conditions, particularly in the SB ring.

A large effort has been made to detect warm gas in the planet formation zone of circumstellar discs using space and ground-based near infrared facilities. GV Tau N, the most obscured component of the GV Tau system, is an outstanding source, being one of the first targets detected in HCN and the only one detected in CH$_4$ so far. Although near infrared observations have shed light on its chemical content, the physical structure and kinematics of the circumstellar matter remained unknown. We use interferometric images of the HCN 3-2 and $^{13}$CO 3-2 lines, and far-IR observations of $^{13}$CO, HCN, CN and H$_2$O transitions to discern the morphology, kinematics, and chemistry of the dense gas close to the star. These observations constitute the first detection of H$_2$O towards GV Tau N. Moreover, ALMA high spatial resolution (~ 7 au) images of the continuum at 1.1 mm and the HCN 3-2 line resolve different gas components towards GV Tau N, a gaseous disc with R~25 au, an ionized jet, and one (or two) molecular outflows. The asymmetric morphology of the gaseous disc shows that it has been eroded by the jet. All observations can be explained if GV Tau N is binary, and the primary component has a highly inclined individual disc relative to the circumbinary disc. We discuss the origin of the water and the other molecules emission according to this scenario. In particular, we propose that the water emission would come from the disrupted gaseous disc and the molecular outflows.

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

CO$_2$ ice has a phase transition at 35 K when its structure changes from amorphous to crystalline. Using Reflection absorption Infrared Spectroscopy (RAIRS), \"Oberg et al. observed that the photodesorption yield of CO$_2$ ice deposited at 60 K and irradiated at 18 K is 40% lower than that of CO$_2$ ice deposited and irradiated at 18 K. In this work, CO$_2$ ices were deposited at 16-60 K and UV-irradiated at 16 K to rule out the temperature effect and figure out the relationship between photodesorption yield and ice structure. IR spectroscopy is a common method used for measurement of the photodesorption yield in ices. We found that undetectable C atoms produced in irradiated CO$_2$ ice can account for 33% of the amount of depleted CO$_2$ molecules in the ice. A quantitative calibration of QMS was therefore performed to convert the measured ion current into photodesorption yield. During various irradiation periods, the dominant photodesorbing species were CO, O$_2$, and CO$_2$, and their photodesorption yields in CO$_2$ ices deposited at different temperature configurations were almost the same, indicating that ice morphology has no effect on the photodesorption yield of CO$_2$ ice. In addition, we found that the lower desorption yield reported by Mart\'in-Dom\'enech et al. is due to a linear relationship between the photodesorption yield and the combination of energy distribution of Microwave-Discharge Hydrogen-flow Lamp (MDHL) and UV absorption cross section of ices.

We investigate the fueling and the feedback of nuclear activity in the Seyfert 2 galaxy NGC1068, by studying the distribution and kinematics of molecular gas in the torus and its connections. We use ALMA to image the emission of a set of molecular gas tracers in the circumnuclear disk (CND) and the torus of the galaxy using the CO(2-1), CO(3-2) and HCO+(4-3) lines with spatial resolutions ~0.03"-0.09"(2-6pc). ALMA resolves the CND as an asymmetric ringed disk of D~400pc-size and mass of ~1.4x10^8 Msun. The inner edge of the ring is associated with edge-brightened arcs of NIR polarized emission identified with the working surface of the AGN ionized wind. ALMA proves the existence of a molecular torus of M_torus ~ 3x10^5Msun, which extends over a large range of spatial scales D=10-30pc around the central engine. The new observations evidence the density radial stratification of the torus: the HCO+(4-3) torus, with a full size D=11pc, is a factor of 2-3 smaller than its CO(2-1) and CO(3-2) counterparts, which have full-sizes D=26pc and D=28pc, respectively. The torus is connected to the CND through a network of gas streamers. The kinematics of molecular gas show strong departures from circular motions in the torus, the gas streamers, and the CND. These velocity distortions are interconnected and are part of a 3D outflow that reflects the effects of AGN feedback on the kinematics of molecular gas across a wide range of spatial scales. We conclude that a wide-angle AGN wind launched from the accretion disk is impacting a sizeable fraction of the gas inside the torus (~0.4-0.6 x M_torus). However, a large gas reservoir (~1.2-1.8 x 10^5Msun) close to the equatorial plane of the torus remains unaffected by the AGN wind and can continue fueling the AGN for ~1-4Myr. AGN fueling seems nevertheless thwarted on intermediate scales (15pc < r < 50pc).

High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2, harbors one of the closest such systems, making it an excellent target for case studies. We investigate the morphology, stability and dynamical properties of the hub-filament system on basis of 13CO and C18O observations obtained with the IRAM-30m telescope and H2 column density maps derived from Herschel dust emission observations. We identified the filamentary network and characterized the individual filaments as either main (converging into the hub) or secondary (converging to a main filament) filaments. The main filaments have line masses of 30-100 Msun/pc and show signs of fragmentation. The secondary filaments have line masses of 12-60 Msun/pc and show fragmentation only sporadically. In the context of Ostriker's hydrostatic filament model, the main filaments are thermally super-critical. If non-thermal motions are included, most of them are trans-critical. Most of the secondary filaments are roughly trans-critical regardless of whether non-thermal motions are included or not. From the main filaments, we estimate a mass accretion rate of 10(-4)-10(-3) Msun/pc into the hub. The secondary filaments accrete into the main filaments with a rate of 0.1-0.4x10(-4) Msun/pc. The main filaments extend into the hub. Their velocity gradients increase towards the hub, suggesting acceleration of the gas. We estimate that with the observed infall velocity, the mass-doubling time of the hub is ~2.5 Myr, ten times larger than the free-fall time, suggesting a dynamically old region. These timescales are comparable with the chemical age of the HII region. Inside the hub, the main filaments show a ring- or a spiral-like morphology that exhibits rotation and infall motions. One possible explanation for the morphology is that gas is falling into the central cluster following a spiral-like pattern.

AB Aur is a Herbig Ae star hosting a well-known transitional disk. Because of its proximity and low inclination angle, it is an excellent object to study planet formation. Our goal is to investigate the chemistry and dynamics of the molecular gas component in the AB Aur disk, and its relation with the prominent horseshoe shape observed in continuum mm emission. We used the NOEMA interferometer to map with high angular resolution the J = 3-2 lines of HCO+ and HCN. By combining both, we can gain insight into the AB Aur disk structure. Chemical segregation is observed in the AB Aur disk: HCO+ shows intense emission toward the star position, at least one bright molecular bridge within the dust cavity, and ring-like emission at larger radii, while HCN is only detected in an annular ring that is coincident with the dust ring and presents an intense peak close to the dust trap. We use HCO+ to investigate the gas dynamics inside the cavity. The observed bright HCO+ bridge connects the compact central source with the outer dusty ring. This bridge can be interpreted as an accretion flow from the outer ring to the inner disk/jet system proving gas accretion through the cavity.

Aims. Our goal is to complete the inventory of S-bearing molecules and their abundances in the prototypical photodissociation region (PDR) the Horsehead nebula to gain insight into sulphur chemistry in UV irradiated regions. Based on the WHISPER millimeter (mm) line survey, our goal is to provide an improved and more accurate description of sulphur species and their abundances towards the core and PDR positions in the Horsehead. Methods. The Monte Carlo Markov chain (MCMC) methodology and the molecular excitation and radiative transfer code RADEX were used to explore the parameter space and determine physical conditions and beam-averaged molecular abundances. Results. A total of 13 S-bearing species (CS, SO, SO2, OCS, H2CS - both ortho and para - HDCS, C2S, HCS+, SO+, H2S, S2H, NS and NS+) have been detected in the two targeted positions. This is the first detection of SO+ in the Horsehead and the first detection of NS+ in any PDR. We find a differentiated chemical behaviour between C-S and O-S bearing species within the nebula. The C-S bearing species C2S and o-H2CS present fractional abundances a factor grater than two higher in the core than in the PDR. In contrast, the O-S bearing molecules SO, SO2, and OCS present similar abundances towards both positions. A few molecules, SO+, NS, and NS+, are more abundant towards the PDR than towards the core, and could be considered as PDR tracers. Conclusions. This is the first complete study of S-bearing species towards a PDR. Our study shows that CS, SO, and H2S are the most abundant S-bearing molecules in the PDR with abundances of a few 1E-9. We recall that SH, SH+, S, and S+ are not observable at the wavelengths covered by the WHISPER survey. At the spatial scale of our observations, the total abundance of S atoms locked in the detected species is < 1E-8, only ~0.1% of the cosmic sulphur abundance.

GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) to determine the fractional abundances of CO, HCO+, HCN, CS, SO, HCS+, and N2H+ in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from Av~3 to ~20mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of (1-5)x10$^3$ cm$^{-3}$; and ii) the dense phase, located at Av>10mag, with molecular hydrogen densities >10$^4$ cm$^{-3}$. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone where C/H~8x10$^{-5}$ and C/O~0.8-1, until the beginning of the dense phase at Av~10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H~(0.4 - 2.2) x10$^{-6}$, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H~8x10$^{-8}$). Based on our chemical modeling, we constrain the value of $\zeta_{\rm H_2}$ to ~(0.5 - 1.8) x10$^{-16}$ s$^{-1}$ in the translucent cloud.

We have developed the first gas-grain chemical model for oxygen fractionation (also including sulphur fractionation) in dense molecular clouds, demonstrating that gas-phase chemistry generates variable oxygen fractionation levels, with a particularly strong effect for NO, SO, O2, and SO2. This large effect is due to the efficiency of the neutral 18O + NO, 18O + SO, and 18O + O2 exchange reactions. The modeling results were compared to new and existing observed isotopic ratios in a selection of cold cores. The good agreement between model and observations requires that the gas-phase abundance of neutral oxygen atoms is large in the observed regions. The S16O/S18O ratio is predicted to vary substantially over time showing that it can be used as a sensitive chemical proxy for matter evolution in dense molecular clouds.

To our knowledge, R Mon is the only B0 star in which a gaseous Keplerian disk has been detected. However, there is some controversy about the spectral type of R Mon. Some authors propose that it could be a later B8e star, where disks are more common. We have re-evaluated the spectral type of R Mon using the available continuum data and UVES emission lines. We used a power-law disk model to fit previous 12 CO 1-0 and 2-1 interferometric observations and the PACS CO data to investigate the disk structure. Interferometric detections of 13 CO J=1-0, HCO+ 1-0, and CN 1-0 lines using the IRAM Plateau de Bure Interferometer (PdBI) are presented. The HCN 1-0 line was not detected. Our analysis confirms that R Mon is a B0 star. The disk model compatible with the 12 CO 1-0 and 2-1 interferometric observations falls short of predicting the observed fluxes of the 14<Ju<31 PACS lines; this is consistent with the scenario in which some contribution to these lines is coming from a warm envelope and/or UV-illuminated outflow walls. More interestingly, the upper limits to the fluxes of the Ju>31 CO lines suggest the existence of a region empty of CO at R<20 au in the proto-planetary disk. The intense emission of the HCO+ and CN lines shows the strong influence of UV photons on gas chemistry. The observations gathered in this paper are consistent with the presence around R Mon of a transition disk with a cavity about 20 au. This size is similar to the photoevaporation radius that supports the interpretation that UV photoevaporation is main disk dispersal mechanism in massive stars

The Barnard 1b core shows signatures of being at the earliest stages of low-mass star formation, with two extremely young and deeply embedded protostellar objects. Hence, this core is an ideal target to study the structure and chemistry of the first objects formed in the collapse of prestellar cores. We present ALMA Band 6 spectral line observations at ~0.6'' of angular resolution towards Barnard 1b. We have extracted the spectra towards both protostars, and used a Local Thermodynamic Equilibrium (LTE) model to reproduce the observed line profiles. B1b-S shows rich and complex spectra, with emission from high energy transitions of complex molecules, such as CH3OCOH and CH3CHO, including vibrational level transitions. We have tentatively detected for the first time in this source emission from NH2CN, NH2CHO, CH3CH2OH, CH2OHCHO, CH3CH2OCOH and both aGg' and gGg' conformers of (CH2OH)2. This is the first detection of ethyl formate (CH3CH2OCOH) towards a low-mass star forming region. On the other hand, the spectra of the FHSC candidate B1b-N are free of COMs emission. In order to fit the observed line profiles in B1b-S, we used a source model with two components: an inner hot and compact component (200 K, 0.35'') and an outer and colder one (60 K, 0.6''). The resulting COM abundances in B1b-S range from 1e-13 for NH2CN and NH2CHO, up to 1e-9 for CH3OCOH. Our ALMA Band 6 observations reveal the presence of a compact and hot component in B1b-S, with moderate abundances of complex organics. These results indicate that a hot corino is being formed in this very young Class 0 source.

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.

Evidence is mounting that the small bodies of our Solar System, such as comets and asteroids, have at least partially inherited their chemical composition from the first phases of the Solar System formation. It then appears that the molecular complexity of these small bodies is most likely related to the earliest stages of star formation. It is therefore important to characterize and to understand how the chemical evolution changes with solar-type protostellar evolution. We present here the Large Program "Astrochemical Surveys At IRAM" (ASAI). Its goal is to carry out unbiased millimeter line surveys between 80 and 272 GHz of a sample of ten template sources, which fully cover the first stages of the formation process of solar-type stars, from prestellar cores to the late protostellar phase. In this article, we present an overview of the surveys and results obtained from the analysis of the 3 mm band observations. The number of detected main isotopic species barely varies with the evolutionary stage and is found to be very similar to that of massive star-forming regions. The molecular content in O- and C- bearing species allows us to define two chemical classes of envelopes, whose composition is dominated by either a) a rich content in O-rich complex organic molecules, associated with hot corino sources, or b) a rich content in hydrocarbons, typical of Warm Carbon Chain Chemistry sources. Overall, a high chemical richness is found to be present already in the initial phases of solar-type star formation.

We present subarcsecond 1.3 mm continuum ALMA observations towards the Orion Molecular Cloud 1 South (OMC-1S) region, down to a spatial resolution of 74 AU, which reveal a total of 31 continuum sources. We also present subarcsecond 7 mm continuum VLA observations of the same region, which allow to further study fragmentation down to a spatial resolution of 40 AU. By applying a Mean Surface Density of Companions method we find a characteristic spatial scale at ~560 AU, and we use this spatial scale to define the boundary of 19 `cores' in OMC-1S as groupings of millimeter sources. We find an additional characteristic spatial scale at ~2900 AU, which is the typical scale of the filaments in OMC-1S, suggesting a two-level fragmentation process. We measured the fragmentation level within each core and find a higher fragmentation towards the southern filament. In addition, the cores of the southern filament are also the densest (within 1100 AU) cores in OMC-1S. This is fully consistent with previous studies of fragmentation at spatial scales one order of magnitude larger, and suggests that fragmentation down to 40 AU seems to be governed by thermal Jeans processes in OMC-1S.

We present the first detection of gas phase S2H in the Horsehead, a moderately UV-irradiated nebula. This confirms the presence of doubly sulfuretted species in the interstellar medium and opens a new challenge for sulfur chemistry. The observed S2H abundance is ~5x10$^{-11}$, only a factor 4-6 lower than that of the widespread H2S molecule. H2S and S2H are efficiently formed on the UV-irradiated icy grain mantles. We performed ice irradiation experiments to determine the H2S and S2H photodesorption yields. The obtained values are ~1.2x10$^{-3}$ and <1x10$^{-5}$ molecules per incident photon for H2S and S2H, respectively. Our upper limit to the S2H photodesorption yield suggests that photo-desorption is not a competitive mechanism to release the S2H molecules to the gas phase. Other desorption mechanisms such as chemical desorption, cosmic-ray desorption and grain shattering can increase the gaseous S2H abundance to some extent. Alternatively, S2H can be formed via gas phase reactions involving gaseous H2S and the abundant ions S+ and SH+. The detection of S2H in this nebula could be therefore the result of the coexistence of an active grain surface chemistry and gaseous photo-chemistry.

We present far-infrared observations of Monoceros R2 (a giant molecular cloud at approximately 830 pc distance, containing several sites of active star formation), as observed at 70 {\mu}m, 160 {\mu}m, 250 {\mu}m, 350 {\mu}m, and 500 {\mu}m by the Photodetector Array Camera and Spectrometer (PACS) and Spectral and Photometric Imaging Receiver (SPIRE) instruments on the Herschel Space Observatory as part of the Herschel imaging survey of OB young stellar objects (HOBYS) Key programme. The Herschel data are complemented by SCUBA-2 data in the submillimetre range, and WISE and Spitzer data in the mid-infrared. In addition, C18O data from the IRAM 30-m Telescope are presented, and used for kinematic information. Sources were extracted from the maps with getsources, and from the fluxes measured, spectral energy distributions were constructed, allowing measurements of source mass and dust temperature. Of 177 Herschel sources robustly detected in the region (a detection with high signal-to-noise and low axis ratio at multiple wavelengths), including protostars and starless cores, 29 are found in a filamentary hub at the centre of the region (a little over 1% of the observed area). These objects are on average smaller, more massive, and more luminous than those in the surrounding regions (which together suggest that they are at a later stage of evolution), a result that cannot be explained entirely by selection effects. These results suggest a picture in which the hub may have begun star formation at a point significantly earlier than the outer regions, possibly forming as a result of feedback from earlier star formation. Furthermore, the hub may be sustaining its star formation by accreting material from the surrounding filaments.

The formation epoch of protostellar disks is debated because of the competing roles of rotation, turbulence, and magnetic fields in the early stages of low-mass star formation. Magnetohydrodynamics simulations of collapsing cores predict that rotationally supported disks may form in strongly magnetized cores through ambipolar diffusion or misalignment between the rotation axis and the magnetic field orientation. Detailed studies of individual sources are needed to cross check the theoretical predictions. We present 0.06-0.1" resolution images at 350 GHz toward B1b-N and B1b-S, which are young class 0 protostars, possibly first hydrostatic cores. The images have been obtained with ALMA, and we compare these data with magnetohydrodynamics simulations of a collapsing turbulent and magnetized core. The submillimeter continuum emission is spatially resolved by ALMA. Compact structures with optically thick 350 GHz emission are detected toward both B1b-N and B1b-S, with 0.2 and 0.35" radii (46 and 80 au at the Perseus distance of 230 pc), within a more extended envelope. The flux ratio between the compact structure and the envelope is lower in B1b-N than in B1b-S, in agreement with its earlier evolutionary status. The size and orientation of the compact structure are consistent with 0.2" resolution 32 GHz observations obtained with the Very Large Array as a part of the VANDAM survey, suggesting that grains have grown through coagulation. The morphology, temperature, and densities of the compact structures are consistent with those of disks formed in numerical simulations of collapsing cores. Moreover, the properties of B1b-N are consistent with those of a very young protostar, possibly a first hydrostatic core. These observations provide support for the early formation of disks around low-mass protostars.

We study the feedback of star formation and nuclear activity on the chemistry of molecular gas in NGC1068, a nearby (D=14Mpc) Seyfert 2 barred galaxy, by analyzing if the abundances of key molecular species like ethynyl (C2H), a classical tracer of PDR, change in the different environments of the disk of the galaxy. We have used ALMA to map the emission of the hyperfine multiplet of C2H(N=1-0) and its underlying continuum emission in the central r~35"(2.5kpc)-region of the disk of NGC1068 with a spatial resolution 1.0"x0.7"(50-70pc). We have developed a set of time-dependent chemical models to determine the origin of the C2H gas. A sizeable fraction of the total C2H line emission is detected from the r~1.3kpc starburst (SB) ring. However, the brightest C2H emission originates from a r~200pc off-centered circumnuclear disk (CND), where evidence of a molecular outflow has been previously found in other molecular tracers imaged by ALMA. We also detect significant emission that connects the CND with the outer disk. We derived the fractional abundances of C2H (X(C2H)) assuming LTE conditions. Our estimates range from X(C2H)~a few 10^-8 in the SB ring up to X(C2H)~ a few 10^-7 in the outflow region. PDR models that incorporate gas-grain chemistry are able to account for X(C2H) in the SB ring for moderately dense (n(H2)>10^4 cm^-3) and moderately UV-irradiated gas (UV-field<10xDraine field) in a steady-state regime. However, the high fractional abundances estimated for C2H in the outflow region can only be reached at very early times (T< 10^2-10^3 yr) in models of UV/X-ray irradiated dense gas (n(H2)>10^4-10^5) cm^-3). We interpret that the transient conditions required to fit the high values of X(C2H) in the outflow are likely due to UV/X-ray irradiated non-dissociative shocks associated with the highly turbulent interface between the outflow and the molecular gas in NGC1068.

The extremely young Class 0 object B1b-S and the first hydrostatic core (FSHC) candidate, B1b-N, provide a unique opportunity to study the chemical changes produced in the elusive transition from the prestellar core to the protostellar phase. We present 40"x70" images of Barnard 1b in the 13CO 1->0, C18O 1->0, NH2D 1_{1,1}a->1_{0,1}s, and SO 3_2->2_1 lines obtained with the NOEMA interferometer. The observed chemical segregation allows us to unveil the physical structure of this young protostellar system down to scales of ~500au. The two protostellar objects are embedded in an elongated condensation, with a velocity gradient of ~0.2-0.4 m s^{-1} au^{-1} in the east-west direction, reminiscent of an axial collapse. The NH2D data reveal cold and dense pseudo-disks (R~500-1000 au) around each protostar. Moreover, we observe evidence of pseudo-disk rotation around B1b-S. We do not see any signature of the bipolar outflows associated with B1b-N and B1b-S, which were previously detected in H2CO and CH3OH, in any of the imaged species. The non-detection of SO constrains the SO/CH3OH abundance ratio in the high-velocity gas.

One serious challenge for planet formation is the rapid inward drift of pebble-sized dust particles in protoplanetary disks. Dust trapping at local maxima in the disk gas pressure has received much theoretical attention but still lacks observational support. The cold dust emission in the AB Aur disk forms an asymmetric ring at a radius of about 120 au, which is suggestive of dust trapping in a gas vortex. We present high spatial resolution (0".58x0".78 ~ 80x110 au) NOEMA observations of the 1.12 mm and 2.22 mm dust continuum emission from the AB Aur disk. Significant azimuthal variations of the flux ratio at both wavelengths indicate a size segregation of the large dust particles along the ring. Our continuum images also show that the intensity variations along the ring are smaller at 2.22 mm than at 1.12 mm, contrary to what dust trapping models with a gas vortex have predicted. Our two-fluid (gas+dust) hydrodynamical simulations demonstrate that this feature is well explained if the gas vortex has started to decay due to turbulent diffusion, and dust particles are thus losing the azimuthal trapping on different timescales depending on their size. The comparison between our observations and simulations allows us to constrain the size distribution and the total mass of solid particles in the ring, which we find to be of the order of 30 Earth masses, enough to form future rocky planets.

The early stages of low-mass star formation are likely to be subject to intense ionization by protostellar energetic MeV particles. As a result, the surrounding gas is enriched in molecular ions, such as HCO$^{+}$ and N$_{2}$H$^{+}$. Nonetheless, this phenomenon remains poorly understood for Class 0 objects. Recently, based on Herschel observations taken as part of the key program Chemical HErschel Surveys of Star forming regions (CHESS), a very low HCO$^{+}$/N$_{2}$H$^{+}$ abundance ratio of about 3-4, has been reported toward the protocluster OMC-2 FIR4. This finding suggests a cosmic-ray ionization rate in excess of 10$^{-14}$ s$^{-1}$, much higher than the canonical value of $\zeta$ = 3$\times$10$^{-17}$ s$^{-1}$ (value expected in quiescent dense clouds). To assess the specificity of OMC-2 FIR4, we have extended this study to a sample of sources in low- and intermediate mass. More specifically, we seek to measure the HCO$^{+}$/N$_2$H$^{+}$ abundance ratio from high energy lines (J $\ge$ 6) toward this source sample in order to infer the flux of energetic particles in the warm and dense gas surrounding the protostars. We use observations performed with the Heterodyne Instrument for the FarInfrared spectrometer on board the Herschel Space Observatory toward a sample of 9 protostars. We report HCO$^{+}$/N$_2$H$^{+}$ abundance ratios in the range of 5 up to 73 toward our source sample. The large error bars do not allow us to conclude whether OMC-2~FIR4 is a peculiar source. Nonetheless, an important result is that the measured HCO$^{+}$/N$_2$H$^{+}$ ratio does not vary with the source luminosity. At the present time, OMC-2 FIR4 remains the only source where a high flux of energetic particles is clearly evident. More sensitive and higher angular resolution observations are required to further investigate this process.

Far-UV photons (E<13.6 eV) from hot massive stars regulate, or at least influence, the heating, ionization, and chemistry of most of the neutral interstellar medium. Investigating the interaction between FUV radiation and interstellar matter thus plays an important role in astrochemistry. We have used ALMA to mosaic a small field of the Orion Bar where the critical transition from atomic to molecular gas takes place. These observations provide an unprecedented sharp view of this transition layer (~1" resolution or ~414 AU). The resulting images (so far in the rotational emission of CO, HCO+, H13CO+, SO+, SO, and reactive ions SH+ and HOC+) show the small-scale structure in gas density and temperature, and the steep abundance gradients. The images reveal a pattern of high-density substructures, photo-ablative gas flows and instabilities at the edge of the molecular cloud. We have also used the IRAM 30m telescope to carry out a line-survey of the illuminated edge of the Bar in the mm domain. Our observations reveal the presence of complex organic molecules that were not expected in such a harsh environment. In particular, we have reported the first detection of the unstable cis conformer of formic acid (HCOOH) in the ISM. The energy barrier to internal rotation (the conversion from trans to cis) is approximately 4827 cm-1. Hence, this detection is surprising. The low inferred trans-to-cis abundance ratio of 2.8+/-1.0 supports a photoswitching mechanism: a given conformer absorbs a FUV stellar photon that radiatively excites the molecule to electronic states above the interconversion barrier. Subsequent fluorescent decay leaves the molecule in a different conformer form. This mechanism, which we have specifically studied with ab initio quantum calculations, was not considered so far in astrochemistry although it can affect the structure of a variety of molecules in PDRs. (abridged)

We investigate the presence of COMs in strongly UV-irradiated interstellar molecular gas. We have carried out a complete millimetre line survey using the IRAM30m telescope towards the edge of the Orion Bar photodissociation region (PDR), close to the H2 dissociation front, a position irradiated by a very intense far-UV (FUV) radiation field. These observations have been complemented with 8.5 arcsec resolution maps of the H2CO 5(1,5)-4(1,4) and C18O 3-2 emission at 0.9 mm. Despite being a harsh environment, we detect more than 250 lines from COMs and related precursors: H2CO, CH3OH, HCO, H2CCO, CH3CHO, H2CS, HCOOH, CH3CN, CH2NH, HNCO, H13-2CO, and HC3N (in decreasing order of abundance). For each species, the large number of detected lines allowed us to accurately constrain their rotational temperatures (Trot) and column densities (N). Owing to subthermal excitation and intricate spectroscopy of some COMs (symmetric- and asymmetric-top molecules such as CH3CN and H2CO, respectively), a correct determination of N and Trot requires building rotational population diagrams of their rotational ladders separately. We also provide accurate upper limit abundances for chemically related molecules that might have been expected, but are not conclusively detected at the edge of the PDR (HDCO, CH3O, CH3NC, CH3CCH, CH3OCH3, HCOOCH3, CH3CH2OH, CH3CH2CN, and CH2CHCN). A non-LTE LVG excitation analysis for molecules with known collisional rate coefficients, suggests that some COMs arise from different PDR layers but we cannot resolve them spatially. In particular, H2CO and CH3CN survive in the extended gas directly exposed to the strong FUV flux (Tk = 150-250 K and Td > 60 K), whereas CH3OH only arises from denser and cooler gas clumps in the more shielded PDR interior (Tk = 40-50 K). We find a HCO/H2CO/CH3OH = 1/5/3 abundance ratio. These ratios are different from those inferred in hot cores and shocks.

We report high angular resolution (4.9" x 3.0") images of reactive ions SH+, HOC+, and SO+ toward the Orion Bar photodissociation region (PDR). We used ALMA-ACA to map several rotational lines at 0.8 mm, complemented with multi-line observations obtained with the IRAM 30m telescope. The SH+ and HOC+ emission is restricted to a narrow layer of 2"- to 10"-width (~800 to 4000 AU depending on the assumed PDR geometry) that follows the vibrationally excited H2^* emission. Both ions efficiently form very close to the H/H2 transition zone, at a depth of A_V < 1 mag into the neutral cloud, where abundant C+, S+, and H2^* coexist. SO+ peaks slightly deeper into the cloud. The observed ions have low rotational temperatures (T_rot~10-30 K << T_k) and narrow line-widths (~2-3 km/s), a factor of ~2 narrower that those of the lighter reactive ion CH+. This is consistent with the higher reactivity and faster radiative pumping rates of CH+ compared to the heavier ions, which are driven relatively faster toward smaller velocity dispersion by elastic collisions and toward lower T_rot by inelastic collisions. We estimate column densities and average physical conditions from a non-LTE excitation model (n(H2)~10^5-10^6 cm^-3, n(e^-)~10 cm^-3, and T_k~200 K). Regardless of the excitation details, SH+ and HOC+ clearly trace the most exposed layers of the UV-irradiated molecular cloud surface, whereas SO+ arises from slightly more shielded layers.

Extragalactic observations of water emission can provide valuable insights into the excitation of the interstellar medium. In addition, extragalactic megamasers are powerful probes of kinematics close to active nuclei. Therefore, it is paramount to determine the true origin of the water emission, whether it is excited by processes close to an AGN or in star-forming regions. We use ALMA Band 5 science verification observations to analyse the emission of the 183 GHz water line in Arp 220 on sub-arcsecond scales, in conjunction with new ALMA Band 7 data at 325 GHz. Specifically, the nature of the process leading to the excitation of emission at these water lines is studied in this context. Supplementary 22 GHz VLA observations are used to better constrain the parameter space in the excitation modelling of the water lines. We detect 183 GHz H2O and 325 GHz water emission towards the two compact nuclei at the center of Arp 220, being brighter in Arp 220 West. The emission at these two frequencies is compared to previous single-dish data and does not show evidence of variability. The 183 and 325 GHz lines show similar spectra and kinematics, but the 22 GHz profile is significantly different in both nuclei due to a blend with an NH3 absorption line. Our findings suggest that the most likely scenario to cause the observed water emission in Arp 220 is a large number of independent masers originating from numerous star-forming regions.

We report on subarcsecond observations of complex organic molecules (COMs) in the high-mass protostar IRAS20126+4104 with the Plateau de Bure Interferometer in its most extended configurations. In addition to the simple molecules SO, HNCO and H2-13CO, we detect emission from CH3CN, CH3OH, HCOOH, HCOOCH3, CH3OCH3, CH3CH2CN, CH3COCH3, NH2CN, and (CH2OH)2. SO and HNCO present a X-shaped morphology consistent with tracing the outflow cavity walls. Most of the COMs have their peak emission at the putative position of the protostar, but also show an extension towards the south(east), coinciding with an H2 knot from the jet at about 800-1000 au from the protostar. This is especially clear in the case of H2-13CO and CH3OCH3. We fitted the spectra at representative positions for the disc and the outflow, and found that the abundances of most COMs are comparable at both positions, suggesting that COMs are enhanced in shocks as a result of the passage of the outflow. By coupling a parametric shock model to a large gas-grain chemical network including COMs, we find that the observed COMs should survive in the gas phase for about 2000 yr, comparable to the shock lifetime estimated from the water masers at the outflow position. Overall, our data indicate that COMs in IRAS20126+4104 may arise not only from the disc, but also from dense and hot regions associated with the outflow.

SiO(3-2) and HNCO(6-5) emission has been imaged in NGC 1068 with the Plateau de Bure Interferometer (PdBI). We perform an LTE and RADEX analysis to determine the column densities and physical characteristics of the gas emitting these two lines. We then use a chemical model to determine the origin of the emission. There is a strong SiO peak to the East of the AGN, with weak detections to the West. This distribution contrasts that of HNCO, which is detected more strongly to the West. The SiO emission peak in the East is similar to the peak of the molecular gas mass traced by CO. HNCO emission is offset from this peak by as much as 80 pc ( 1"). We compare velocity integrated line ratios in the East and West. We confirm that SiO emission strongly dominates in the East, while the reverse is true in the West. We use RADEX to analyse the possible gas conditions that could produce such emission. We find that, in both East and West, we cannot constrain a single temperature for the gas. We run a grid of chemical models of potential shock processes in the CND and find that SiO is significantly enhanced during a fast (60 km/s) shock but not during a slow (20 km/s) shock, nor in a gas not subjected to shocks at all. We find the inverse for HNCO, whose abundance increases during slow shocks and in warm non-shocked gas. High SiO and low HNCO indicated a fast shock, while high HNCO and low SiO indicates either a slow shock or warm, dense, non-shocked gas. The East Knot is therefore likely to contain gas that is heavily shocked. From chemical modelling, gas in the West Knot may be non-shocked, or maybe undergoing a much milder shock event. When taking into account RADEX results, the milder shock event is the more likely of the two scenarios.

The Orion Bar is the archetypal edge-on molecular cloud surface illuminated by strong ultraviolet radiation from nearby massive stars. Owing to the close distance to Orion (about 1,350 light-year), the effects of stellar feedback on the parental cloud can be studied in detail. Visible-light observations of the Bar(1) show that the transition between the hot ionised gas and the warm neutral atomic gas (the ionisation front) is spatially well separated from the transition from atomic to molecular gas (the dissociation front): about 15 arcseconds or 6,200 astronomical units (one astronomical unit is the Earth-Sun distance). Static equilibrium models(2,3) used to interpret previous far-infrared and radio observations of the neutral gas in the Bar(4,5,6) (typically at 10-20 arcsecond resolution) predict an inhomogeneous cloud structure consisting of dense clumps embedded in a lower density extended gas component. Here we report one-arcsecond-resolution millimetre-wave images that allow us to resolve the molecular cloud surface. In contrast to stationary model predictions(7,8,9), there is no appreciable offset between the peak of the H2 vibrational emission (delineating the H/H2 transition) and the edge of the observed CO and HCO+ emission. This implies that the H/H2 and C+/C/CO transition zones are very close. These observations reveal a fragmented ridge of high-density substructures, photoablative gas flows and instabilities at the molecular cloud surface. The results suggest that the cloud edge has been compressed by a high-pressure wave that currently moves into the molecular cloud. The images demonstrate that dynamical and nonequilibrium effects are important for the cloud evolution.

The CO+ reactive ion is thought to be a tracer of the boundary between a HII region and the hot molecular gas. In this study, we present the spatial distribution of the CO+ rotational emission toward the Mon R2 star-forming region. The CO+ emission presents a clumpy ring-like morphology, arising from a narrow dense layer around the HII region. We compare the CO+ distribution with other species present in photon-dominated regions (PDR), such as [CII] 158 mm, H2 S(3) rotational line at 9.3 mm, polycyclic aromatic hydrocarbons (PAHs) and HCO+. We find that the CO+ emission is spatially coincident with the PAHs and [CII] emission. This confirms that the CO+ emission arises from a narrow dense layer of the HI/H2 interface. We have determined the CO+ fractional abundance, relative to C+ toward three positions. The abundances range from 0.1 to 1.9x10^(-10) and are in good agreement with previous chemical model, which predicts that the production of CO+ in PDRs only occurs in dense regions with high UV fields. The CO+ linewidth is larger than those found in molecular gas tracers, and their central velocity are blue-shifted with respect to the molecular gas velocity. We interpret this as a hint that the CO+ is probing photo-evaporating clump surfaces.

In a previous study of the L1157 B1 shocked cavity, a comparison between NH$_3$(1$_0$-$0_0$) and H$_2$O(1$_{\rm 10}$--1$_{\rm 01}$) transitions showed a striking difference in the profiles, with H$_2$O emitting at definitely higher velocities. This behaviour was explained as a result of the high-temperature gas-phase chemistry occurring in the postshock gas in the B1 cavity of this outflow. If the differences in behaviour between ammonia and water are indeed a consequence of the high gas temperatures reached during the passage of a shock, then one should find such differences to be ubiquitous among chemically rich outflows. In order to determine whether the difference in profiles observed between NH$_3$ and H$_2$O is unique to L1157 or a common characteristic of chemically rich outflows, we have performed Herschel-HIFI observations of the NH$_3$(1$_0$-0$_0$) line at 572.5 GHz in a sample of 8 bright low-mass outflow spots already observed in the H$_2$O(1$_{\rm 10}$--1$_{\rm 01}$) line within the WISH KP. We detected the ammonia emission at high-velocities at most of the outflows positions. In all cases, the water emission reaches higher velocities than NH$_3$, proving that this behaviour is not exclusive of the L1157-B1 position. Comparisons with a gas-grain chemical and shock model confirms, for this larger sample, that the behaviour of ammonia is determined principally by the temperature of the gas.

Barnard B1b has revealed as one of the most interesting globules from the chemical and dynamical point of view. It presents a rich molecular chemistry characterized by large abundances of deuterated and complex molecules. Furthermore, it hosts an extremely young Class 0 object and one candidate to First Hydrostatic Core (FHSC). Our aim was to determine the cosmic ray ionization rate and the depletion factors in this extremely young star forming region. We carried out a spectral survey towards Barnard 1b as part of the IRAM Large program ASAI using the IRAM 30-m telescope at Pico Veleta (Spain). This provided a very complete inventory of neutral and ionic C-, N- and S- bearing species with, up to our knowledge, the first secure detections of the deuterated ions DCS+ and DOCO+. We used a state-of-the-art pseudo-time-dependent gas-phase chemical model to determine the value of the cosmic ray ionization rate and the depletion factors. The observational data were well fitted with $\zeta_{H_2}$ between 3E-17 s$^{-1}$ and 1E-16 s$^{-1}$. Elemental depletions were estimated to be ~10 for C and O, ~1 for N and ~25 for S. Barnard B1b presents similar depletions of C and O than those measured in pre-stellar cores. The depletion of sulfur is higher than that of C and O but not as extreme as in cold cores. In fact, it is similar to the values found in some bipolar outflows, hot cores and photon-dominated regions. Several scenarios are discussed to account for these peculiar abundances. We propose that it is the consequence of the initial conditions (important outflows and enhanced UV fields in the surroundings) and a rapid collapse (~0.1 Myr) that permits to maintain most S- and N-bearing species in gas phase to great optical depths. The interaction of the compact outflow associated with B1b-S with the surrounding material could enhance the abundances of S-bearing molecules, as well.

We have used the Atacama Large Millimeter Array (ALMA) to map the emission of the CO(6-5) molecular line and the 432 {\mu}m continuum emission from the 300 pc-sized circumnuclear disk (CND) of the nearby Seyfert 2 galaxy NGC 1068 with a spatial resolution of ~4 pc. These observations spatially resolve the CND and, for the first time, image the dust emission, the molecular gas distribution, and the kinematics from a 7-10 pc-diameter disk that represents the submillimeter counterpart of the putative torus of NGC 1068. We fitted the nuclear spectral energy distribution of the torus using ALMA and near and mid-infrared (NIR/MIR) data with CLUMPY models. The mass and radius of the best-fit solution for the torus are both consistent with the values derived from the ALMA data alone: Mgas_torus=(1+-0.3)x10^5 Msun and Rtorus=3.5+-0.5 pc. The dynamics of the molecular gas in the torus show non-circular motions and enhanced turbulence superposed on the rotating pattern of the disk. The kinematic major axis of the CO torus is tilted relative to its morphological major axis. By contrast with the nearly edge-on orientation of the H2O megamaser disk, we have found evidence suggesting that the molecular torus is less inclined (i=34deg-66deg) at larger radii. The lopsided morphology and complex kinematics of the torus could be the signature of the Papaloizou-Pringle instability, long predicted to likely drive the dynamical evolution of active galactic nuclei (AGN) tori.

Transitional disks are structures of dust and gas around young stars with large inner cavities in which planet formation may occur. Lopsided dust distributions are observed in the dust continuum emission at millimeter wavelengths. These asymmetrical structures can be explained as the result of an enhanced gas density vortex where the dust is trapped potentially promoting the rapid growth to the planetesimal scale. AB Aur hosts a transitional disk with a clear horseshoe morphology which strongly suggests the presence of a dust trap. Our goal is to investigate its formation and the possible effects on the gas chemistry. We used the NOEMA interferometer to image the 1mm continuum dust emission and the 13CO J=2->1, C18O J=2->1, SO J=56->45 and H2CO J=303->202 rotational lines. Line integrated intensity ratio images are built to investigate the chemical changes within the disk. We have used a single point (n,T) chemical model to investigate the lifetime of gaseous CO, H2CO and SO in the dust trap. Our model shows that for densities >10^7 cm^-3, the SO molecules are depleted (directly frozen or converted into SO2 and then frozen out) in less than 0.1 Myr. The lower SO abundance towards the dust trap could indicate that a larger fraction of the gas is in a high density environment. Gas dynamics, grain growth and gas chemistry are coupled in the planet formation process. Because of the strong dependence of SO abundance on the gas density, the sulfur chemistry can be used as a chemical diagnostic to detect the birthsites of future planets. However, the large uncertainties inherent to chemical models and the limited knowledge of the disk physical structure and initial conditions are important drawbacks.

We present the first observations of H$^{13}$CN$(1-0)$, H$^{13}$CO$^+(1-0)$ and SiO$(2-1)$ in NGC\,6240, obtained with the IRAM PdBI. Combining a Markov Chain Monte Carlo (MCMC) code with Large Velocity Gradient (LVG) modelling, and with additional data from the literature, we simultaneously fit three gas phases and six molecular species to constrain the physical condition of the molecular gas, including mass$-$luminosity conversion factors. We find $\sim10^{10}M_\odot$ of dense molecular gas in cold, dense clouds ($T_{\rm k}\sim10$\,K, $n_{{\rm H}_2}\sim10^6$\,cm$^{-3}$) with a volume filling factor $<0.002$, embedded in a shock heated molecular medium ($T_{\rm k}\sim2000$\,K, $n_{{\rm H}_2}\sim10^{3.6}$\,cm$^{-3}$), both surrounded by an extended diffuse phase ($T_{\rm k}\sim200$\,K, $n_{{\rm H}_2}\sim10^{2.5}$\,cm$^{-3}$). We derive a global $\alpha_{\rm CO}=1.5^{7.1}_{1.1}$ with gas masses $\log_{10}\left(M / [M_\odot]\right)=10.1_{10.0}^{10.8}$, dominated by the dense gas. We also find $\alpha_{\rm HCN} = 32^{89}_{13}$, which traces the cold, dense gas. The [$^{12}$C]/[$^{13}$C] ratio is only slightly elevated ($98^{230}_{65}$), contrary to the very high [CO]/[$^{13}$CO] ratio (300-500) reported in the literature. However, we find very high [HCN]/[H$^{13}$CN] and [HCO$^+$]/[H$^{13}$CO$^+$] abundance ratios $(300^{500}_{200})$ which we attribute to isotope fractionation in the cold, dense clouds.

We present the first 7.5'x11.5' velocity-resolved map of the [CII]158um line toward the Orion molecular cloud-1 (OMC-1) taken with the Herschel/HIFI instrument. In combination with far-infrared (FIR) photometric images and velocity-resolved maps of the H41alpha hydrogen recombination and CO J=2-1 lines, this data set provides an unprecedented view of the intricate small-scale kinematics of the ionized/PDR/molecular gas interfaces and of the radiative feedback from massive stars. The main contribution to the [CII] luminosity (~85%) is from the extended, FUV-illuminated face of the cloud G_0>500, n_H>5x10^3 cm^-3) and from dense PDRs (G_0~10^4, n_H~10^5 cm^-3) at the interface between OMC-1 and the HII region surrounding the Trapezium cluster. Around 15% of the [CII] emission arises from a different gas component without CO counterpart. The [CII] excitation, PDR gas turbulence, line opacity (from [13CII]) and role of the geometry of the illuminating stars with respect to the cloud are investigated. We construct maps of the [CII]/FIR and FIR/M_Gas ratios and show that [CII]/FIR decreases from the extended cloud component (10^-2-10^-3) to the more opaque star-forming cores (10^-3-10^-4). The lowest values are reminiscent of the "[CII] deficit" seen in local ultra-luminous IR galaxies hosting vigorous star formation. Spatial correlation analysis shows that the decreasing [CII]/FIR ratio correlates better with the column density of dust through the molecular cloud than with FIR/M_Gas. We conclude that the [CII] emitting column relative to the total dust column along each line of sight is responsible for the observed [CII]/FIR variations through the cloud.

Cyanogen (NCCN) is the simplest member of the series of dicyanopolyynes. It has been hypothesized that this family of molecules can be important constituents of interstellar and circumstellar media, although the lack of a permanent electric dipole moment prevents its detection through radioastronomical techniques. Here we present the first solid evidence of the presence of cyanogen in interstellar clouds through the detection of its protonated form toward the cold dark clouds TMC-1 and L483. Protonated cyanogen (NCCNH+) has been identified through the J=5-4 and J=10-9 rotational transitions using the 40m radiotelescope of Yebes and the IRAM 30m telescope. We derive beam averaged column densities for NCCNH+ of (8.6+/-4.4)e10 cm-2 in TMC-1 and (3.9+/-1.8)e10 cm-2 in L483, which translate to fairly low fractional abundances relative to H2, in the range (1-10)e-12. The chemistry of protonated molecules in dark clouds is discussed, and it is found that, in general terms, the abundance ratio between the protonated and non protonated forms of a molecule increases with increasing proton affinity. Our chemical model predicts an abundance ratio NCCNH+/NCCN of 1e-4, which implies that the abundance of cyanogen in dark clouds could be as high as (1-10)e-8 relative to H2, i.e., comparable to that of other abundant nitriles such as HCN, HNC, and HC3N.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

Planets are thought to form in the gas and dust disks around young stars. In particular, it has been proposed that giant planets can form through the gravitational instability of massive extended disks around intermediate-mass stars. However, we still lack direct observations to constrain this mechanism. We have spatially resolved the 8.6 and 11.2 $\mu$m emission of a massive protoplanetary disk seen edge on around an A star, Gomez's Hamburger (GoHam), using VISIR at the Very Large Telescope. A compact region situated at a projected distance of $350\pm50$ AU south of the central star is found to have a reduced emission.This asymmetry is fully consistent with the presence of a cold density structure, or clump, identified in earlier CO observations, and we derive physical characteristics consistent with those observations: a mass of 0.8-11.4 Jupiter masses (for a dust-to-gas mass ratio of 0.01), a radius of about 10$^2$ astronomical units, and a local density of about $10^{7}$ cm$^{-3}$. Based on this evidence, we argue that this clump, which we call GoHam b, is a promising candidate for a young protoplanet formed by gravitational instability that might be representative of the precursors of massive planets observed around A stars, such as HR 8799 or Beta pictoris. More detailed studies at high angular resolution are needed to better constrain the physical properties of this object to confirm this proposal.

Aims. Our goal is to determine the molecular composition of the circumstellar disk around AB Aurigae (hereafter, AB Aur). AB Aur is a prototypical Herbig Ae star and the understanding of its disk chemistry is of paramount importance to understand the chemical evolution of the gas in warm disks. Methods. We used the IRAM 30-m telescope to perform a sensitive search for molecular lines in AB Aur as part of the IRAM Large program ASAI (A Chemical Survey of Sun-like Star-forming Regions). These data were complemented with interferometric observations of the HCO+ 1-0 and C17O 1-0 lines using the IRAM Plateau de Bure Interferometer (PdBI). Single-dish and interferometric data were used to constrain chemical models. Results. Throughout the survey, several lines of CO and its isotopologues, HCO+, H2CO, HCN, CN and CS, were detected. In addition, we detected the SO 54-33 and 56-45 lines, confirming the previous tentative detection. Comparing to other T Tauri's and Herbig Ae disks, AB Aur presents low HCN 3-2/HCO+ 3-2 and CN 2-1/HCN 3-2 line intensity ratios, similar to other transition disks. AB Aur is the only protoplanetary disk detected in SO thus far. Conclusions. We modeled the line profiles using a chemical model and a radiative transfer 3D code. Our model assumes a flared disk in hydrostatic equilibrium. The best agreement with observations was obtained for a disk with a mass of 0.01 Msun , Rin=110 AU, Rout=550 AU, a surface density radial index of 1.5 and an inclination of 27 deg. The intensities and line profiles were reproduced within a factor of 2 for most lines. This agreement is reasonable taking into account the simplicity of our model that neglects any structure within the disk. However, the HCN 3-2 and CN 2-1 line intensities were predicted more intense by a factor of >10. We discuss several scenarios to explain this discrepancy.

Intermediate-mass young stellar objects (YSOs) provide a link to understand how feedback from shocks and UV radiation scales from low to high-mass star forming regions. Aims: Our aim is to analyze excitation of CO and H$_2$O in deeply-embedded intermediate-mass YSOs and compare with low-mass and high-mass YSOs. Methods: Herschel/PACS spectral maps are analyzed for 6 YSOs with bolometric luminosities of $L_\mathrm{bol}\sim10^2 - 10^3$ $L_\odot$. The maps cover spatial scales of $\sim 10^4$ AU in several CO and H$_2$O lines located in the $\sim55-210$ $\mu$m range. Results: Rotational diagrams of CO show two temperature components at $T_\mathrm{rot}\sim320$ K and $T_\mathrm{rot}\sim700-800$ K, comparable to low- and high-mass protostars probed at similar spatial scales. The diagrams for H$_2$O show a single component at $T_\mathrm{rot}\sim130$ K, as seen in low-mass protostars, and about $100$ K lower than in high-mass protostars. Since the uncertainties in $T_\mathrm{rot}$ are of the same order as the difference between the intermediate and high-mass protostars, we cannot conclude whether the change in rotational temperature occurs at a specific luminosity, or whether the change is more gradual from low- to high-mass YSOs. Conclusions: Molecular excitation in intermediate-mass protostars is comparable to the central $10^{3}$ AU of low-mass protostars and consistent within the uncertainties with the high-mass protostars probed at $3\cdot10^{3}$ AU scales, suggesting similar shock conditions in all those sources.

In the theory of star formation, the first hydrostatic core (FHSC) phase is a critical step in which a condensed object emerges from a prestellar core. This step lasts about one thousand years, a very short time compared with the lifetime of prestellar cores, and therefore is hard to detect unambiguously. We present IRAM Plateau de Bure observations of the Barnard 1b dense molecular core, combining detections of H2CO and CH3OH spectral lines and dust continuum at 2.3" resolution (~ 500 AU). The two compact cores B1b-N and B1b-S are detected in the dust continuum at 2mm, with fluxes that agree with their spectral energy distribution. Molecular outflows associated with both cores are detected. They are inclined relative to the direction of the magnetic field, in agreement with predictions of collapse in turbulent and magnetized gas with a ratio of mass to magnetic flux somewhat higher than the critical value, \mu ~ 2 - 7. The outflow associated with B1b-S presents sharp spatial structures, with ejection velocities of up to ~ 7 kms from the mean velocity. Its dynamical age is estimated to be ~2000 yrs. The B1b-N outflow is smaller and slower, with a short dynamical age of ~1000 yrs. The B1b-N outflow mass, mass-loss rate, and mechanical luminosity agree well with theoretical predictions of FHSC. These observations confirm the early evolutionary stage of B1b-N and the slightly more evolved stage of B1b-S.

In dense parts of interstellar clouds (> 10^5 cm^-3), dust & gas are expected to be in thermal equilibrium, being coupled via collisions. However, previous studies have shown that the temperatures of the dust & gas may remain decoupled even at higher densities. We study in detail the temperatures of dust & gas in the photon-dominated region S 140, especially around the deeply embedded infrared sources IRS 1-3 and at the ionization front. We derive the dust temperature and column density by combining Herschel PACS continuum observations with SOFIA observations at 37 $\mu$m and SCUBA at 450 $\mu$m. We model these observations using greybody fits and the DUSTY radiative transfer code. For the gas part we use RADEX to model the CO 1-0, CO 2-1, 13CO 1-0 and C18O 1-0 emission lines mapped with the IRAM-30m over a 4' field. Around IRS 1-3, we use HIFI observations of single-points and cuts in CO 9-8, 13CO 10-9 and C18O 9-8 to constrain the amount of warm gas, using the best fitting dust model derived with DUSTY as input to the non-local radiative transfer model RATRAN. We find that the gas temperature around the infrared sources varies between 35 and 55K and that the gas is systematically warmer than the dust by ~5-15K despite the high gas density. In addition we observe an increase of the gas temperature from 30-35K in the surrounding up to 40-45K towards the ionization front, most likely due to the UV radiation from the external star. Furthermore, detailed models of the temperature structure close to IRS 1 show that the gas is warmer and/or denser than what we model. Finally, modelling of the dust emission from the sub-mm peak SMM 1 constrains its luminosity to a few ~10^2 Lo. We conclude that the gas heating in the S 140 region is very efficient even at high densities, most likely due to the deep UV penetration from the embedded sources in a clumpy medium and/or oblique shocks.

We present interferometric observations of the CN 1-0 (113.491 GHz), N2H+ 1-0 (93.173 GHz), H(41)a (92.034 GHz), CH3CN (91.987 GHz), CS 3-2 (146.969 GHz), c-C3H2 3-2 (145.089 GHz), H2CO 2-1 (145.603 GHz) and HC3N 16-15 (145.601 GHz) lines towards M82, carried out with the IRAM Plateau de Bure Interferometer (PdBI). PDR chemical modelling is used to interpret these observations. Our results show that the abundances of N2H+, CS and H13 CO+ remain quite constant across the galaxy confirming that these species are excellent tracers of the dense molecular gas. On the contrary, the abundance of CN increases by a factor of 3 in the inner x2 bar orbits. The [CN]/[N2 H+ ] ratio is well correlated with the H(41)a emission at all spatial scales down to 100 pc. Chemical modelling shows that the variations in the [CN]/[N2H+] ratio can be explained as the consequence of differences in the local intestellar UV field and in the average cloud sizes within the nucleus of the galaxy. Our high-spatial resolution imaging of the starburst galaxy M 82 shows that the star formation activity has a strong impact on the chemistry of the molecular gas. In particular, the entire nucleus behaves as a giant photon-dominated region (PDR) whose chemistry is determined by the local UV flux. The detection of N2H+ shows the existence of a population of clouds with Av >20 mag all across the galaxy plane. These clouds constitute the molecular gas reservoir for the formation of new stars and, although distributed all along the nucleus, the highest concentration occurs in the outer x1 bar orbits (R = 280 pc).

We study the spatial distribution and chemistry of small hydrocarbons in the Orion Bar PDR. We used the IRAM-30m telescope to carry out a millimetre line survey towards the Orion Bar edge, complemented with ~2'x2' maps of the C2H and c-C3H2 emission. We analyse the excitation of the detected hydrocarbons and constrain the physical conditions of the emitting regions with non-LTE radiative transfer models. We compare the inferred column densities with updated gas-phase photochemical models including 13CCH and C13CH isotopomer fractionation. ~40% of the lines in the survey arise from hydrocarbons (C2H, C4H, c-C3H2, c-C3H, C13CH, 13CCH, l-C3H and l-H2C3). We detect new lines from l-C3H+ and improve its rotational spectroscopic constants. Anions or deuterated hydrocarbons are not detected: [C2D]/[C2H]<0.2%, [C2H-]/[C2H]<0.007% and [C4H-]/[C4H]<0.05%. Our gas-phase models can reasonably match the observed column densities of most hydrocarbons (within factors <3). Since the observed spatial distribution of the C2H and c-C3H2 emission is similar but does not follow the PAH emission, we conclude that, in high UV-flux PDRs, photodestruction of PAHs is not a necessary requirement to explain the observed abundances of the smallest hydrocarbons. Instead, gas-phase endothermic reactions (or with barriers) between C+, radicals and H2 enhance the formation of simple hydrocarbons. Observations and models suggest that the [C2H]/[c-C3H2] ratio (~32 at the PDR edge) decreases with the UV field attenuation. The observed low cyclic-to-linear C3H column density ratio (<3) is consistent with a high electron abundance (Xe) PDR environment. In fact, the poorly constrained Xe gradient influences much of the hydrocarbon chemistry in the more UV-shielded gas. We propose that reactions of C2H isotopologues with 13C+ and H atoms can explain the observed [C13CH]/[13CCH]=1.4(0.1) fractionation level.

We present the results of interferometric spectral line observations of Arp 220 at 3.5mm and 1.2mm from the Plateau de Bure Interferometer (PdBI), imaging the two nuclear disks in H$^{13}$CN$(1 - 0)$ and $(3 - 2)$, H$^{13}$CO$^+(1 - 0)$ and $(3 - 2)$, and HN$^{13}$C$(3 - 2)$ as well as SiO$(2 - 1)$ and $(6 - 5)$, HC$^{15}$N$(3 - 2)$, and SO$(6_6 - 5_5)$. The gas traced by SiO$(6 - 5)$ has a complex and extended kinematic signature including a prominent P Cygni profile, almost identical to previous observations of HCO$^+(3 - 2)$. Spatial offsets $0.1''$ north and south of the continuum centre in the emission and absorption of the SiO$(6 - 5)$ P Cygni profile in the western nucleus (WN) imply a bipolar outflow, delineating the northern and southern edges of its disk and suggesting a disk radius of $\sim40$ pc, consistent with that found by ALMA observations of Arp 220. We address the blending of SiO$(6 - 5)$ and H$^{13}$CO$^+(3 - 2)$ by considering two limiting cases with regards to the H$^{13}$CO$^+$ emission throughout our analysis. Large velocity gradient (LVG) modelling is used to constrain the physical conditions of the gas and to infer abundance ratios in the two nuclei. Our most conservative lower limit on the [H$^{13}$CN]/[H$^{13}$CO$^+$] abundance ratio is 11 in the WN, cf. 0.10 in the eastern nucleus (EN). Comparing these ratios to the literature we argue on chemical grounds for an energetically significant AGN in the WN driving either X-ray or shock chemistry, and a dominant starburst in the EN.

The determination of the physical conditions in molecular clouds is a key step towards our understanding of their formation and evolution of associated star formation. We investigate the density, temperature, and column density of both dust and gas in the photodissociation regions (PDRs) located at the interface between the atomic and cold molecular gas of the NGC 7023 reflection nebula. We study how young stars affect the gas and dust in their environment. Our approach combining both dust and gas delivers strong constraints on the physical conditions of the PDRs. We find dense and warm molecular gas of high column density in the PDRs.

The enormous radiative and mechanical luminosities of massive stars impact a vast range of scales and processes, from the reionization of the universe, to the evolution of galaxies, to the regulation of the interstellar medium, to the formation of star clusters, and even to the formation of planets around stars in such clusters. Two main classes of massive star formation theory are under active study, Core Accretion and Competitive Accretion. In Core Accretion, the initial conditions are self-gravitating, centrally concentrated cores that condense with a range of masses from the surrounding, fragmenting clump environment. They then undergo relatively ordered collapse via a central disk to form a single star or a small-N multiple. In this case, the pre-stellar core mass function has a similar form to the stellar initial mass function. In Competitive Accretion, the material that forms a massive star is drawn more chaotically from a wider region of the clump without passing through a phase of being in a massive, coherent core. In this case, massive star formation must proceed hand in hand with star cluster formation. If stellar densities become very high near the cluster center, then collisions between stars may also help to form the most massive stars. We review recent theoretical and observational progress towards understanding massive star formation, considering physical and chemical processes, comparisons with low and intermediate-mass stars, and connections to star cluster formation.

We report the tentative detection in space of the nitrosylium ion, NO$^+$. The observations were performed towards the cold dense core Barnard 1-b. The identification of the NO$^+$ $J$=2--1 line is supported by new laboratory measurements of NO$^+$ rotational lines up to the $J$=8--7 transition (953207.189\,MHz), which leads to an improved set of molecular constants: $B_0 = 59597.1379(62)$\,MHz, $D_0 = 169.428(65)$\,kHz, and $eQq_0(\textrm{N}) = -6.72(15)$\,MHz. The profile of the feature assigned to NO$^+$ exhibits two velocity components at 6.5 and 7.5 km s$^{-1}$, with column densities of $1.5 \times 10^{12}$ and $6.5\times10^{11}$ cm$^{-2}$, respectively. New observations of NO and HNO, also reported here, allow to estimate the following abundance ratios: $X$(NO)/$X$(NO$^+$)$\simeq511$, and $X$(HNO)/$X$(NO$^+$)$\simeq1$. This latter value provides important constraints on the formation and destruction processes of HNO. The chemistry of NO$^+$ and other related nitrogen-bearing species is investigated by the means of a time-dependent gas phase model which includes an updated chemical network according to recent experimental studies. The predicted abundance for NO$^+$ and NO is found to be consistent with the observations. However, that of HNO relative to NO is too high. No satisfactory chemical paths have been found to explain the observed low abundance of HNO. HSCN and HNCS are also reported here with an abundance ratio of $\simeq1$. Finally, we have searched for NNO, NO$_2$, HNNO$^+$, and NNOH$^+$, but only upper limits have been obtained for their column density, except for the latter for which we report a tentative 3-$\sigma$ detection.

We present a detailed analysis of ALMA Bands 7 and 9 data of CO, HCO+, HCN and CS, augmented with Plateau de Bure Interferometer (PdBI) data of the ~ 200 pc circumnuclear disk (CND) and the ~ 1.3 kpc starburst ring (SB ring) of NGC~1068, a nearby (D = 14 Mpc) Seyfert 2 barred galaxy. We aim at determining the physical characteristics of the dense gas present in the CND and whether the different line intensity ratios we find within the CND as well as between the CND and the SB ring are due to excitation effects (gas density and temperature differences) or to a different chemistry. We estimate the column densities of each species in Local Thermodynamic Equilibrium (LTE). We then compute large one-dimensional non-LTE radiative transfer grids (using RADEX) by using first only the CO transitions, and then all the available molecules in order to constrain the densities, temperatures and column densities within the CND. We finally present a preliminary set of chemical models to determine the origin of the gas. We find that in general the gas in the CND is very dense (> 10^5 cm^-3) and hot (T> 150K), with differences especially in the temperature across the CND. The AGN position has the lowest CO/HCO+, CO/HCN and CO/CS column density ratios. RADEX analyses seem to indicate that there is chemical differentiation across the CND. We also find differences between the chemistry of the SB ring and some regions of the CND; the SB ring is also much colder and less dense than the CND. Chemical modelling does not succeed in reproducing all the molecular ratios with one model per region, suggesting the presence of multi-gas phase components. The LTE, RADEX and chemical analyses all indicate that more than one gas-phase component is necessary to uniquely fit all the available molecular ratios within the CND.

We investigate the fueling and the feedback of star formation and nuclear activity in NGC1068, a nearby (D=14Mpc) Seyfert 2 barred galaxy, by analyzing the distribution and kinematics of the molecular gas in the disk. We have used ALMA to map the emission of a set of dense molecular gas tracers (CO(3-2), CO(6-5), HCN(4-3), HCO+(4-3) and CS(7-6)) and their underlying continuum emission in the central r ~ 2kpc of NGC1068 with spatial resolutions ~ 0.3"-0.5" (~ 20-35pc). Molecular line and dust continuum emissions are detected from a r ~ 200pc off-centered circumnuclear disk (CND), from the 2.6kpc-diameter bar region, and from the r ~ 1.3kpc starburst (SB) ring. Most of the emission in HCO+, HCN and CS stems from the CND. Molecular line ratios show dramatic order-of-magnitude changes inside the CND that are correlated with the UV/X-ray illumination by the AGN, betraying ongoing feedback. The gas kinematics from r ~ 50pc out to r ~ 400pc reveal a massive (M_mol ~ 2.7 (+0.9, -1.2) x 10^7 Msun) outflow in all molecular tracers. The tight correlation between the ionized gas outflow, the radio jet and the occurrence of outward motions in the disk suggests that the outflow is AGN-driven. The outflow rate estimated in the CND, dM/dt ~ 63 (+21, -37) Msun yr^-1, is an order of magnitude higher than the star formation rate at these radii, confirming that the outflow is AGN-driven. The power of the AGN is able to account for the estimated momentum and kinetic luminosity of the outflow. The CND mass load rate of the CND outflow implies a very short gas depletion time scale of <=1 Myr.

The massive star-forming region Mon R2 hosts the closest ultra-compact HII region that can be spatially resolved with current single-dish telescopes. We used the IRAM-30m telescope to carry out an unbiased spectral survey toward two important positions (namely IF and MP2), in order to studying the chemistry of deuterated molecules toward Mon R2. We found a rich chemistry of deuterated species at both positions, with detections of C2D, DCN, DNC, DCO+, D2CO, HDCO, NH2D, and N2D+ and their corresponding hydrogenated species and isotopologs. Our high spectral resolution observations allowed us to resolve three velocity components: the component at 10 km/s is detected at both positions and seems associated with the layer most exposed to the UV radiation from IRS 1; the component at 12 km/s is found toward the IF position and seems related to the molecular gas; finally, a component at 8.5 km/s is only detected toward the MP2 position, most likely related to a low-UV irradiated PDR. We derived the column density of all the species, and determined the deuterium fractions (Dfrac). The values of Dfrac are around 0.01 for all the observed species, except for HCO+ and N2H+ which have values 10 times lower. The values found in Mon R2 are well explained with pseudo-time-dependent gas-phase model in which deuteration occurs mainly via ion-molecule reactions with H2D+, CH2D+ and C2HD+. Finally, the [H13CN]/[HN13C] ratio is very high (~11) for the 10 km/s component, which also agree with our model predictions for an age of ~0.01-0.1 Myr. The deuterium chemistry is a good tool for studying star-forming regions. The low-mass star-forming regions seem well characterized with Dfrac(N2H+) or Dfrac(HCO+), but it is required a complete chemical modeling to date massive star-forming regions, because the higher gas temperature together with the rapid evolution of massive protostars.

NGC 7129 FIRS 2 (hereafter FIRS 2) is an intermediate-mass (2 to 8 Msun) protostar located at a distance of 1250 pc. High spatial resolution observations are required to resolve the hot core at its center. We present a molecular survey from 218200 MHz to 221800 MHz carried out with the IRAM Plateau de Bure Interferometer. These observations were complemented with a long integration single-dish spectrum taken with the IRAM 30m telescope. We used a Local Thermodynamic Equilibrium (LTE) single temperature code to model the whole dataset. The interferometric spectrum is crowded with a total of ~300 lines from which a few dozens remain unidentified yet. The spectrum has been modeled with a total of 20 species and their isomers, isotopologues and deuterated compounds. Complex molecules like methyl formate (CH3OCHO), ethanol (CH3CH2OH), glycolaldehyde (CH2OHCHO), acetone (CH3COCH3), dimethyl ether (CH3OCH3), ethyl cyanide (CH3CH2CN) and the aGg' conformer of ethylene glycol (aGg'-(CH2OH)_2) are among the detected species. The detection of vibrationally excited lines of CH3CN, CH3OCHO, CH3OH, OCS, HC3N and CH3CHO proves the existence of gas and dust at high temperatures. In fact, the gas kinetic temperature estimated from the vibrational lines of CH3CN, ~405 K, is similar to that measured in massive hot cores. Our data allow an extensive comparison of the chemistry in FIRS~2 and the Orion hot core. We find a quite similar chemistry in FIRS 2 and Orion. Most of the studied fractional molecular abundances agree within a factor of 5. Larger differences are only found for the deuterated compounds D2CO and CH2DOH and a few molecules (CH3CH2CN, SO2, HNCO and CH3CHO). Since the physical conditions are similar in both hot cores, only different initial conditions (warmer pre-collapse phase in the case of Orion) and/or different crossing time of the gas in the hot core can explain this behavior.

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 micron (or 1.2 mm in two cases) and the Spectral Energy Distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations.

Context. Monoceros R2 (Mon R2), at a distance of 830 pc, is the only ultra-compact H ii region (UC H ii) where its associated photon-dominated region (PDR) can be resolved with the Herschel Space Observatory. Aims. Our aim is to investigate observationally the kinematical patterns in the interface regions (i.e., the transition from atomic to molecular gas) associated with Mon R2. Methods. We used the HIFI instrument onboard Herschel to observe the line profiles of the reactive ions CH+, OH+ and H2O+ toward different positions in Mon R2. We derive the column density of these molecules and compare them with gas-phase chemistry models. Results. The reactive ion CH+ is detected both in emission (at central and red-shifted velocities) and in absorption (at blue-shifted velocities). OH+ is detected in absorption at both blue- and red-shifted velocities, with similar column densities. H2O+ is not detected at any of the positions, down to a rms of 40 mK toward the molecular peak. At this position, we find that the OH+ absorption originates in a mainly atomic medium, and therefore is associated with the most exposed layers of the PDR. These results are consistent with the predictions from photo-chemical models. The line profiles are consistent with the atomic gas being entrained in the ionized gas flow along the walls of the cavity of the H ii region. Based on this evidence, we are able to propose a new geometrical model for this region. Conclusions. The kinematical patterns of the OH+ and CH+ absorption indicate the existence of a layer of mainly atomic gas for which we have derived, for the first time, some physical parameters and its dynamics.

The evolution of deeply embedded young stellar objects (YSOs) has traditionally been traced through dust continuum spectral energy distributions (SEDs), but the use of CO excitation as an evolutionary probe has not yet been explored due to lack of high-J CO observations. The aim is to constrain the physical characteristics (excitation, kinematics, column density) of the warm gas toward 26 low-mass Class 0 and I YSOs using spectrally-resolved Herschel Space Observatory data of high-J lines of CO. Data are complemented by ground-based observations from APEX and the JCMT to compare those with the colder gas traced by lower-J CO lines. This is the first large spectrally resolved high-J CO survey conducted for these types of sources. The median excitation temperatures for 12CO, 13CO and C18O derived from single-temperature fits to the J_u=2-10 integrated intensities are ~70K, 48K and 37K, respectively, with no significant difference between Class 0 and I sources and no trend with M_env or L_bol. Thus, in contrast with the continuum SEDs, the spectral line energy distributions (SLEDs) do not show an evolution during the embedded stage. In contrast, the integrated line intensities of all CO isotopologs show a clear decrease with evolutionary stage as the envelope is dispersed. The H2O 110-101/CO 10-9 intensity ratio does not change significantly with velocity, in contrast with the H2O/CO 3-2 ratio, indicating that CO 10-9 is the lowest transition for which the line wings probe the same warm shocked gas as H2O. Modeling of the full suite of C18O lines indicates an abundance profile for Class 0 sources that is consistent with a freeze-out zone below 25 K and evaporation at higher temperatures, with but some fraction of the CO transformed into other species in the cold phase. In contrast, observations for two Class I sources in Ophiuchus are consistent with a constant high CO abundance profile.

The high resolution spectrum of the v4 band of NH3D+ has been measured by difference frequency IR laser spectroscopy in a multipass hollow cathode discharge cell. From the set of molecular constants obtained from the analysis of the spectrum, a value of 262817(6) MHz (3sigma) has been derived for the frequency of the 1(0)-0(0) rotational transition. This value supports the assignment to NH3D+ of lines at 262816.7 MHz recorded in radio astronomy observations in Orion-IRc2 and the cold prestellar core B1-bS.

We report on the detection of a narrow feature at 262816.73 MHz towards Orion and the cold prestellar core B1-bS, that we attribute to the 1(0)-0(0) line of the deuterated Ammonium ion, NH3D+. The observations were performed with the IRAM 30m radio telescope. The carrier has to be a light molecular species as it is the only feature detected over 3.6 GHz of bandwidth. The hyperfine structure is not resolved indicating a very low value for the electric quadrupolar coupling constant of Nitrogen which is expected for NH3D+ as the electric field over the N nucleus is practically zero. Moreover, the feature is right at the predicted frequency for the 1(0)-0(0) transition of the Ammonium ion, 262817(6) MHz (3sigma), using rotational constants derived from new infrared data obtained in our laboratory in Madrid. The estimated column density is 1.1(0.2)e12 cm-2. Assuming a deuterium enhancement similar to that of NH2D, we derive N(NH4+) sim 2.6e13 cm-2, i.e., an abundance for Ammonium of a few 1e(-11).

Broadband spectral surveys of protostars offer a rich view of the physical, chemical and dynamical structure and evolution of star-forming regions. The Herschel Space Observatory opened up the terahertz regime to such surveys, giving access to the fundamental transitions of many hydrides and to the high-energy transitions of many other species. A comparative analysis of the chemical inventories and physical processes and properties of protostars of various masses and evolutionary states is the goal of the Herschel CHEmical Surveys of Star forming regions (CHESS) key program. This paper focusses on the intermediate-mass protostar, OMC-2 FIR 4. We obtained a spectrum of OMC-2 FIR 4 in the 480 to 1902 GHz range with the HIFI spectrometer onboard Herschel and carried out the reduction, line identification, and a broad analysis of the line profile components, excitation, and cooling. We detect 719 spectral lines from 40 species and isotopologs. The line flux is dominated by CO, H2O, and CH3OH. The line profiles are complex and vary with species and upper level energy, but clearly contain signatures from quiescent gas, a broad component likely due to an outflow, and a foreground cloud. We find abundant evidence for warm, dense gas, as well as for an outflow in the field of view. Line flux represents 2% of the 7 L_Sol luminosity detected with HIFI in the 480 to 1250 GHz range. Of the total line flux, 60% is from CO, 13% from H2O and 9% from CH3OH. A comparison with similar HIFI spectra of other sources is set to provide much new insight into star formation regions, a case in point being a difference of two orders of magnitude in the relative contribution of sulphur oxides to the line cooling of Orion KL and OMC-2 FIR 4.

OMC-2 FIR 4 is one of the closest known young intermediate-mass protoclusters, located at a distance of 420 pc in Orion. This region is one of the few where the complete 500-2000 GHz spectrum has been observed with the heterodyne spectrometer HIFI on board the Herschel satellite, and unbiased spectral surveys at 0.8, 1, 2 and 3 mm have been obtained with the JCMT and IRAM 30-m telescopes. In order to investigate the morphology of this region, we used the IRAM Plateau de Bure Interferometer to image OMC-2 FIR 4 in the 2-mm continuum emission, as well as in DCO+(2-1), DCN(2-1), C34S(3-2), and several CH3OH lines. In addition, we analysed observations of the NH3(1,1) and (2,2) inversion transitions made with the Very Large Array of the NRAO. The resulting maps have an angular resolution which allows us to resolve structures of 5", equivalent to 2000 AU. Our observations reveal three spatially resolved sources within OMC-2 FIR 4, of one or several solar masses each, with hints of further unresolved substructure within them. Two of these sources have elongated shapes and are associated with dust continuum emission peaks, thus likely containing at least one molecular core each. One of them also displays radio continuum emission, which may be attributed to a young B3-B4 star that dominates the overall luminosity output of the region. The third source identified displays a DCO+(2-1) emission peak, and weak dust continuum emission. Its higher abundance of DCO+ relative to the other two regions suggests a lower temperature and therefore its possible association with either a younger low-mass protostar or a starless core. It may alternatively be part of the colder envelope of OMC-2 FIR 4. Our interferometric observations evidence the complexity of this region, where multiple cores, chemical differentiation and an ionised region all coexist within an area of only 10000 AU.

We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope and Herschel. We determine the column densities of these species, and compare their spatial distributions with that of polycyclic aromatic hydrocarbon (PAH). We compare the observational results with different chemical models to explore the relative importance of gas-phase, grain-surface and time-dependent chemistry in these environments. The emission of the small hydrocarbons show different patterns. The CCH emission is extended while CH and c-C$_3$H$_2$ are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of c-C$_3$H$_2$ and CCH shows spatial variations up to a factor of a few, increasing from $N(c-C$_3$H$_2$)/N(CCH)\approx0.004$ in the envelope to a maximum of $\sim0.015-0.029$ towards the 8$\mu$m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite constant over a wide range of G$_0$, whereas the abundance of c-C$_3$H$_2$ is higher in low-UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce c-C$_3$H$_2$. In the molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons abundances. Our study shows that CCH and c-C$_3$H$_2$ present a complex chemistry in which UV photons, grain-surface chemistry and time dependent effects contribute to determine their abundances.

We observationally investigate the relation between the photoelectric heating efficiency in PDRs and the charge of PAHs, which are considered to play a key role in photoelectric heating. Using PACS onboard Herschel, we observed six PDRs spanning a wide range of FUV radiation fields (G_0=100-10^5). To measure the photoelectric heating efficiency, we obtained the intensities of the main cooling lines, i.e., the [OI]63um, 145um, and [CII]158um, as well as the FIR continuum intensity. We used Spitzer/IRS spectroscopic mapping observations to investigate the MIR PAH features in the same regions. We decomposed the MIR PAH emission into that of neutral (PAH^0) and positively ionized (PAH^+) species to derive the fraction of the positively charged PAHs, and compare it to the photoelectric heating efficiency. The heating efficiency traced by ([OI]63um+[OI]145um+[CII]158um) / TIR, ranges between 0.1% and 0.9% in different sources, and the fraction of PAH^+ relative to (PAH^0 + PAH^+) spans from 0(+11)% to 87(+/-10)%. All positions with a high PAH^+ fraction show a low heating efficiency, and all positions with a high heating efficiency have a low PAH^+ fraction, supporting the scenario in which a positive grain charge results in a decreased heating efficiency. Theoretical estimates of the photoelectric heating efficiency show a stronger dependence on the charging parameter gamma=G_0 T^{1/2}/n_e than the observed efficiency reported in this study, and the discrepancy is significant at low gamma. The photoelectric heating efficiency on PAHs, traced by ([OI]63um+[OI]145um+[CII]158um) / (PAH+[OI]63um+[OI]145um+[CII]158um), shows a much better match between the observations and the theoretical estimates. The good agreement of the photoelectric heating efficiency on PAHs with a theoretical model indicates the dominant contribution of PAHs to the photoelectric heating. (abridged for arXiv)

As part of the Herschel Guaranteed Time Key Programme CHESS, we present the discovery of a diffuse gas component in the foreground of the intermediate-mass protostar OMC-2 FIR 4, located in the Orion A region. Making use of the full HIFI spectrum of OMC-2 FIR 4 obtained in CHESS, we detected several ground-state lines from OH+, H2O+, HF, and CH+, all of them seen in absorption against the dust continuum emission of the protostar's envelope. The lines peak at a velocity of 9 km/s, which is blue-shifted by 2 km/s with respect to the systemic velocity of OMC-2 FIR 4 (V = 11.4 km/s). We derived column densities for each species, as well as an upper limit to the column density of the undetected H3O+. In order to model and characterise the foreground cloud, we used the Meudon PDR code to run a homogeneous grid of models that spans a reasonable range of densities, visual extinctions, cosmic ray ionisation rates and far-ultraviolet (FUV) radiation fields. The results of our modelling indicate that the foreground cloud is composed of predominantly neutral diffuse gas (n(H) = 100 cm-3) and is heavily irradiated by an external source of FUV that most likely arises from the nearby Trapezium OB association. The cloud is 6 pc thick and bears many similarities with the so-called C+ interface between Orion-KL and the Trapezium cluster, 2 pc south of OMC-2 FIR 4. We conclude that the foreground cloud we detected is an extension of the C+ interface seen in the direction of Orion KL, and interpret it to be the remains of the parental cloud of OMC-1, which extends from OMC-1 up to OMC-2.

Recently, we introduced detailed isotopic chemistry into the KOSMA-tau model for photon-dominated regions (PDRs) to give theoretical predictions for the abundance of the carbon isotopologues as a function of PDR parameters. Combined with radiative transfer computations for specific geometries, we estimated the possible intensity ratio of the [CII]/[13CII] lines. Here, we compare these predictions with new observations. We performed Herschel/HIFI observations of the [CII] 158micron line in a number of PDRs. In all sources we observed at least two hyperfine components of the [13CII] transition allowing to determine the [CII]/[13CII] intensity ratio, after some revision of the intrinsic hyperfine ratios. Comparing the intensity ratios with the results from the updated KOSMA-tau model, we identify cases dominated by chemical fractionation and cases dominated by the optical depth of the main isotopic line. An observable enhancement of the [CII]/[13CII] intensity ratio due to chemical fractionation depends mostly on geometry and velocity structure, and less on the gas density and radiation field. In our observations the [CII]/[13CII] ratio for the integrated line intensity was always dominated by the optical depth of the main isotopic line. However, an enhanced intensity ratio is found for particular velocity components in a few sources: the red-shifted material in the ultracompact HII region Mon R2, the wings of the turbulent profile in the Orion Bar, and possibly a blue wing in NGC7023. The mapping of the [13CII] lines in the Orion Bar allows to derive a C+ column density map confirming the temperature stratification of the C+ layer, in agreement with the chemical stratification of the Bar. The C+ column densities for all sources show that at the position of the [CII] peak emission, a dominant fraction of the gas-phase carbon is in the form of C+.

In order to study the fragmentation of massive dense cores, which constitute the cluster cradles, we observed with the PdBI in the most extended configuration the continuum at 1.3 mm and the CO(2-1) emission of four massive cores. We detect dust condensations down to ~0.3 Msun and separate millimeter sources down to 0.4" or ~1000 AU, comparable to the sensitivities and separations reached in optical/infrared studies of clusters. The CO(2-1) high angular resolution images reveal high-velocity knots usually aligned with previously known outflow directions. This, in combination with additional cores from the literature observed at similar mass sensitivity and spatial resolution, allowed us to build a sample of 18 protoclusters with luminosities spanning 3 orders of magnitude. Among the 18 regions, ~30% show no signs of fragmentation, while 50% split up into ~4 millimeter sources. We compiled a list of properties for the 18 massive dense cores, such as bolometric luminosity, total mass, and mean density, and found no correlation of any of these parameters with the fragmentation level. In order to investigate the combined effects of magnetic field, radiative feedback and turbulence in the fragmentation process, we compared our observations to radiation magneto-hydrodynamic simulations, and found that the low-fragmented regions are well reproduced in the magnetized core case, while the highly-fragmented regions are consistent with cores where turbulence dominates over the magnetic field. Overall, our study suggests that the fragmentation in massive dense cores could be determined by the initial magnetic field/turbulence balance in each particular core.

Aims: We study the emission of molecular gas in 3C236, a FR II radio source at z~0.1, and search for the footprints of AGN feedback. 3C236 shows signs of a reactivation of its AGN triggered by a recent minor merger episode. Observations have also previously identified an extreme HI outflow in this source. Methods: The IRAM PdBI has been used to study the distribution and kinematics of molecular gas in 3C236 by imaging with high spatial resolution the emission of the 12CO(2-1) line in the nucleus of the galaxy. We have searched for outflow signatures in the CO map. We have also derived the SFR in 3C236 using data available from the literature at UV, optical and IR wavelengths, to determine the star-formation efficiency of molecular gas. Results: The CO emission in 3C236 comes from a spatially resolved 2.6 kpc disk with a regular rotating pattern. Within the limits imposed by the sensitivity and velocity coverage of the CO data, we do not detect any outflow signatures in the cold molecular gas. The disk has a cold gas mass M(H2)~2.1x10^9 Msun. We determine a new value for the redshift of the source zCO=0.09927. The similarity between the CO and HI profiles indicates that the deep HI absorption in 3C236 can be accounted for by a rotating HI structure, restricting the evidence of HI outflow to the most extreme velocities. In the light of the new redshift, the analysis of the ionized gas kinematics reveals a 1000 km/s outflow. As for the CO emitting gas, outflow signatures are nevertheless absent in the warm molecular gas emission traced by infrared H2 lines. The star-formation efficiency in 3C236 is consistent with the value measured in normal galaxies, which follow the canonical KS-law. This result, confirmed to hold in other young radio sources examined in this work, is in stark contrast with the factor of 10-50 lower SFE that seems to characterize evolved powerful radio galaxies.

Large spatial-spectral surveys are more and more common in astronomy. This calls for the need of new methods to analyze such mega- to giga-pixel data-cubes. In this paper we present a method to decompose such observations into a limited and comprehensive set of components. The original data can then be interpreted in terms of linear combinations of these components. The method uses non-negative matrix factorization (NMF) to extract latent spectral end-members in the data. The number of needed end-members is estimated based on the level of noise in the data. A Monte-Carlo scheme is adopted to estimate the optimal end-members, and their standard deviations. Finally, the maps of linear coefficients are reconstructed using non-negative least squares. We apply this method to a set of hyperspectral data of the NGC 7023 nebula, obtained recently with the HIFI instrument onboard the Herschel space observatory, and provide a first interpretation of the results in terms of 3-dimensional dynamical structure of the region.

Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new observations obtained with HIFI and the IRAM-30m telescope. Using a large velocity gradient approach, we model the line intensities and derive an average abundance of H2O and NH3 across the region. Finally, we model the line profiles with a non-local radiative transfer model and compare these results with the abundance predicted by the Meudon PDR code. Results. The variations of the line profiles and intensities indicate complex geometrical and kinematical patterns. The H2O lines present a strong absorption at the ambient velocity and emission in high velocity wings towards the HII region. The spatial distribution of the o-H2^18O line shows that the its emission arises in the PDR surrounding the HII region. By modeling the o-H2^18O emission we derive a mean abundance of o-H2O of ~10^-8 relative to H2. The ortho-H2O abundance is however larger, ~1x10^-7, in the high velocity wings. Possible explanations for this larger abundance include an expanding hot PDR and/or an outflow. Ammonia seems to be present only in the envelope with an average abundance of ~2x10^-9 relative to H2. Conclusions. The Meudon PDR code can account for the measured water abundance in the high velocity gas as long as we assume that it originates from a <1 mag hot expanding layer of the PDR, i.e. that the outflow has only a minor contribution to this emission. To explain the abundances in the rest of the cloud the molecular freeze out and grain surface chemistry would need to be included.

The Plateau de Bure Interferometer has been used to map the continuum emission at 3.4 mm and 1.1 mm together with the J=1->0 and J=3->2 lines of HCN and HCO+ towards the binary star GV Tau. The 3.4 mm observations did not resolve the binary components and the HCN J=1->0 and HCO+ J=1->0 line emissions trace the circumbinary disk and the flattened envelope. However, the 1.1 mm observations resolved the individual disks of GV Tau N and GV Tau S and allowed us to study their chemistry. We detected the HCN 3->2 line only towards the individual disk of GV Tau N, and the emission of the HCO+ 3->2 line towards GV Tau S. Simple calculations indicate that the 3->2 line of HCN is formed in the inner R<12 AU of the disk around GV Tau N where the HCN/HCO+ abundance ratio is >300. On the contrary, this ratio is <1.6 in the disk around GV Tau S. The high HCN abundance measured in GV Tau N is well explained by photo-chemical processes in the warm (>400K) and dense disk surface.

Ultracompact (UC) HII regions constitute one of the earliest phases in the formation of a massive star and are characterized by extreme physical conditions (Go>10^5 Habing field and n>10^6 cm^-3). The UC HII Mon R2 is the closest one and therefore an excellent target to study the chemistry in these complex regions. We carried out a 3mm and 1mm spectral survey using the IRAM 30-m telescope towards three positions that represent different physical environments in Mon R2: (i) the ionization front (IF) at (0",0"); two peaks in the molecular cloud (ii) MP1 at the offset (+15",-15") and (iii) MP2 at the farther offset (0",40"). In addition, we carried out extensive modeling to explain the chemical differences between the three observed regions. We detected more than thirty different species. We detected SO+ and C4H suggesting that UV radiation plays an important role in the molecular chemistry of this region. We detected the typical PDR molecules CN, HCN, HCO, C2H, and c-C3H2. While the IF and the MP1 have a chemistry similar to that found in high UV field and dense PDRs like the Orion Bar, the MP2 is more similar to lower UV/density PDRs like the Horsehead nebula. We also detected complex molecules that are not usually found in PDRs (CH3CN, H2CO, HC3N, CH3OH and CH3C2H). Sulfur compounds CS, HCS+, C2S, H2CS, SO and SO2 and the deuterated species DCN and C2D were also identified. [DCN]/[HCN]=0.03 and [C2D]/[C2H]=0.05, are among the highest in warm regions. Our results show that the high UV/dense PDRs present a different chemistry from that of the low UV case. Abundance ratios like [CO+]/[HCO+] or [HCO]/[HCO+] are good diagnostics to differentiate between them. In Mon R2 we have the two classes of PDRs, a high UV PDR towards the IF and the adjacent molecular bar and a low-UV PDR which extends towards the north-west following the border of the cloud.

Intermediate mass protostars, the bridge between the very common solar-like protostars and the more massive, but rarer, O and B stars, can only be studied at high physical spatial resolutions in a handful of clouds. In this paper we present and analyze the continuum results from an observing campaign at the Submillimeter Array targeting two well-studied intermediate mass protostars in Orion, NGC 2071 and L1641 S3 MMS 1. The extended SMA (eSMA) probes structure at angular resolutions up to 0.2", revealing protostellar disks on scales of 200 AU. Continuum flux measurements on these scales indicate that a significant amount of mass, a few tens of M{\odot}, are present. Envelope, stellar, and disk masses are derived using both compact, extended and eSMA configurations and compared against SED-fitting models. We hypothesize that fragmentation into three components occurred within NGC 2071 at an early time, when the envelopes were less than 10% of their current masses, e.g. < 0.5 M{\odot}. No fragmentation occurred for L1641 S3 MMS 1. For NGC 2071 evidence is given that the bulk of the envelope material currently around each source was accreted after the initial fragmentation. In addition, about 30% of the total core mass is not yet associated to one of the three sources. A global accretion model is favored and a potential accretion history of NGC 2071 is presented. It is shown that the relatively low level of fragmentation in NGC 2071 was stifled compared to the expected fragmentation from a Jean's argument.

NGC 7129 FIRS 2 is a young intermediate-mass (IM) protostar, which is associated with two energetic bipolar outflows and displays clear signs of the presence of a hot core. It has been extensively observed with ground based telescopes and within the WISH Guaranteed Time Herschel Key Program. We present new observations of the C18O 3-2 and the HDO 3_{12}-2_{21} lines towards NGC 7129 FIRS 2. Combining these observations with Herschel data and modeling their emissions, we constrain the C18O and HDO abundance profiles across the protostellar envelope. In particular, we derive the abundance of C18O and HDO in the hot core. The intensities of the C18O lines are well reproduced assuming that the C18O abundance decreases through the protostellar envelope from the outer edge towards the centre until the point where the gas and dust reach the CO evaporation temperature (~20-25 K) where the C18O is released back to the gas phase. Once the C18O is released to the gas phase, the modelled C18O abundance is found to be ~1.6x10^{-8}, which is a factor of 10 lower than the reference abundance. This result is supported by the non-detection of C18O 9-8, which proves that even in the hot core (T_k>100 K) the CO abundance must be 10 times lower than the reference value. Several scenarios are discussed to explain this C18O deficiency. One possible explanation is that during the pre-stellar and protostellar phase, the CO is removed from the grain mantles by reactions to form more complex molecules. Our HDO modeling shows that the emission of HDO 3_{12}-2_{21} line is maser and comes from the hot core (T_k>100 K). Assuming the physical structure derived by Crimier et al. (2010), we determine a HDO abundance of ~0.4 - 1x10^{-7} in the hot core of this IM protostar, similar to that found in the hot corinos NGC 1333 IRAS 2A and IRAS 16293-2422.

Observations with the Plateau de Bure Interferometer in the most extended configuration toward two intermediate-mass star-forming regions, IRAS22198+6336 and AFGL5142, reveal the presence of several complex organic molecules at ~500 AU scales, confirming the presence of hot cores in both regions. The hot cores are not rich in CN-bearing molecules, as often seen in massive hot cores, and are mainly traced by CH3CH2OH, (CH2OH)2, CH3COCH3, and CH3OH, with additionally CH3CHO, CH3OD and HCOOD for IRAS22198+6336, and C6H, and O13CS for AFGL5142. The emission of complex molecules is resolved down to sizes of ~300 and ~600 AU, for IRAS22198+6336 and AFGL5142, respectively, and most likely is tracing protostellar disks rather than flattened envelopes or toroids as usually found. This is specially clear for the case of IRAS22198+6336, where we detect a velocity gradient for all the mapped molecules perpendicular to the most chemically rich outflow of the region, yielding a dynamic mass >4 Msun. As for AFGL5142, the hot core emission is resolved into two elongated cores separated 1800 AU. A detailed comparison of the complex molecule peaks to the new CO(2-1) data and H2O maser data from literature suggests that also for AFGL5142 the complex molecules are mainly associated with disks, except for a faint and extended molecular emission found to the west, which is possibly produced in the interface between one of the outflows and the dense surrounding gas.

Herbig Ae stars (HAe) are the precursors of Vega-type systems and, therefore, crucial objects in planet formation studies. Thus far, only a few disks associated with HAe stars have been studied using millimetre interferometers. Our aim is to determine the dust evolution and the lifetime of the disks associated with Herbig Ae stars. We imaged the continuum emission at 3 mm and 1.3 mm of the Herbig Ae/Be stars BD+61154, RR Tau, VY Mon and LkHa 198 using the Plateau de Bure Interferometer (PdBI). These stars are in the upper end of the stellar mass range of the Herbig Ae stars (stellar mass greater than 3 solar masses). Our measurements were used to complete the Spectral Energy Distribution (SED). The modelling of the SED, in particular the FIR-mm part, allow us to determine the masses and dust properties of these disks. We detected the disks associated with BD+61154, RR Tau and VY Mon with disk masses of 0.35 Msun, 0.05 Msun and 0.40 Msun respectively. The disk around LkHa 198 was not detected with an upper limit to the disk mass of 0.004 Msun. We detected, however, the disks associated with the younger stellar objects LkHa 198--IR and LkHa 198-mm that are located in the vicinity of LkHa 198. The fitting of the mm part of the SED reveal that the grains in the mid-plane of the disks around BD+61154, RR Tau and VY Mon have sizes of 1--1000 microns. Therefore, grains have not grown to centimetre sizes in these disks yet. These massive (M>3 Msun) and young (about 1 Myr) HAe stars are surrounded by massive (>0.04 Msun) disks with grains of micron-millimetre sizes. Although grain growth is proceeding in these disks, their evolutionary stage is prior to the formation of planetesimals. These disks are less evolved than those detected around T Tauri and Herbig Be stars.

Context: L1448-mm is the prototype of a low-mass Class 0 protostar driving a high-velocity jet. Given its bright H2O spectra observed with ISO, L1448-mm is an ideal laboratory to observe heavy water (HDO) emission. Aims: Our aim is to image the HDO emission in the protostar surroundings, the possible occurrence of HDO emission also investigating off L1448-mm, towards the molecular outflow. Methods: We carried out observations of L1448-mm in the HDO(1_10-1_11) line at 80.6 GHz, an excellent tracer of HDO column density, with the IRAM Plateau de Bure Interferometer. Results: We image for the first time HDO emission around L1448-mm. The HDO structure reveals a main clump at velocities close to the ambient one towards the the continuum peak that is caused by the dust heated by the protostar. In addition, the HDO map shows tentative weaker emission at about 2000 AU from the protostar towards the south, which is possibly associated with the walls of the outflow cavity opened by the protostellar wind. Conclusions: Using an LVG code, modelling the density and temperature profile of the hot-corino, and adopting a gas temperature of 100 K and a density of 1.5 10^8 cm^-3, we derive a beam diluted HDO column density of about 7 10^13 cm^-2, corresponding to a HDO abundance of about 4 10^-7. In addition, the present map supports the scenario where HDO can be efficiently produced in shocked regions and not uniquely in hot corinos heated by the newly born star.

We present Herschel observations of the water molecule in the massive star-forming region W3 IRS5. The o-H17O 110-101, p-H18O 111-000, p-H2O 22 202-111, p-H2O 111-000, o-H2O 221-212, and o-H2O 212-101 lines, covering a frequency range from 552 up to 1669 GHz, have been detected at high spectral resolution with HIFI. The water lines in W3 IRS5 show well-defined high-velocity wings that indicate a clear contribution by outflows. Moreover, the systematically blue-shifted absorption in the H2O lines suggests expansion, presumably driven by the outflow. No infall signatures are detected. The p-H2O 111-000 and o-H2O 212-101 lines show absorption from the cold material (T ~ 10 K) in which the high-mass protostellar envelope is embedded. One-dimensional radiative transfer models are used to estimate water abundances and to further study the kinematics of the region. We show that the emission in the rare isotopologues comes directly from the inner parts of the envelope (T > 100 K) where water ices in the dust mantles evaporate and the gas-phase abundance increases. The resulting jump in the water abundance (with a constant inner abundance of 10^{-4}) is needed to reproduce the o-H17O 110-101 and p-H18O 111-000 spectra in our models. We estimate water abundances of 10^{-8} to 10^{-9} in the outer parts of the envelope (T < 100 K). The possibility of two protostellar objects contributing to the emission is discussed.

Very few molecular species have been detected in circumstellar disks surrounding young stellar objects. We are carrying out an observational study of the chemistry of circumstellar disks surrounding T Tauri and Herbig Ae stars. First results of this study are presented in this note. We used the EMIR receivers recently installed at the IRAM 30m telescope to carry a sensitive search for molecular lines in the disks surrounding AB Aur, DM Tau, and LkCa 15. We detected lines of the molecules HCO+, CN, H2CO, SO, CS, and HCN toward AB Aur. In addition, we tentatively detected DCO+ and H2S lines. The line profiles suggest that the CN, HCN, H2CO, CS and SO lines arise in the disk. This makes it the first detection of SO in a circumstellar disk. We have unsuccessfully searched for SO toward DM Tau and LkCa 15, and for c-C3H2 toward AB Aur, DM Tau, and LkCa 15. Our upper limits show that contrary to all the molecular species observed so far, SO is not as abundant in DM Tau as it is in AB Aur. Our results demonstrate that the disk associated with AB Aur is rich in molecular species. Our chemical model shows that the detection of SO is consistent with that expected from a very young disk where the molecular adsorption onto grains does not yet dominate the chemistry.

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. The spectrally-resolved HIFI data are complemented by ground-based observations of lower-J CO and isotopologue lines. The 12CO 10-9 profiles are dominated by broad (FWHM 25-30 km s^-1) emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 100-200 K. Several CO and 13CO line profiles also reveal a medium-broad component (FWHM 5-10 km s^-1), seen prominently in H2O lines. Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The narrow C18O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars.

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species.

HERSCHEL-HIFI observations of water from the intermediate mass protostar NGC7129 FIRS 2 provide a powerful diagnostic of the physical conditions in this star formation environment. Six spectral settings, covering four H216O and two H218O lines, were observed and all but one H218O line were detected. The four H2 16 O lines discussed here share a similar morphology: a narrower, \approx 6 km/s, component centered slightly redward of the systemic velocity of NGC7129 FIRS 2 and a much broader, \approx 25 km/s component centered blueward and likely associated with powerful outflows. The narrower components are consistent with emission from water arising in the envelope around the intermediate mass protostar, and the abundance of H2O is constrained to \approx 10-7 for the outer envelope. Additionally, the presence of a narrow self-absorption component for the lowest energy lines is likely due to self-absorption from colder water in the outer envelope. The broader component, where the H2O/CO relative abundance is found to be \approx 0.2, appears to be tracing the same energetic region that produces strong CO emission at high J.

Using Herschel's HIFI instrument we have observed [C II] along a cut through S140 and high-J transitions of CO and HCO+ at two positions on the cut, corresponding to the externally irradiated ionization front and the embedded massive star forming core IRS1. The HIFI data were combined with available ground-based observations and modeled using the KOSMA-tau model for photon dominated regions. Here we derive the physical conditions in S140 and in particular the origin of [C II] emission around IRS1. We identify three distinct regions of [C II] emission from the cut, one close to the embedded source IRS1, one associated with the ionization front and one further into the cloud. The line emission can be understood in terms of a clumpy model of photon-dominated regions. At the position of IRS1, we identify at least two distinct components contributing to the [C II] emission, one of them a small, hot component, which can possibly be identified with the irradiated outflow walls. This is consistent with the fact that the [C II] peak at IRS1 coincides with shocked H2 emission at the edges of the outflow cavity. We note that previously available observations of IRS1 can be well reproduced by a single-component KOSMA-tau model. Thus it is HIFI's unprecedented spatial and spectral resolution, as well as its sensitivity which has allowed us to uncover an additional hot gas component in the S140 region.

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

Early results from the Herschel Space Observatory revealed the water cation H2O+ to be an abundant ingredient of the interstellar medium. Here we present new observations of the H2O and H2O+ lines at 1113.3 and 1115.2 GHz using the Herschel Space Observatory toward a sample of high-mass star-forming regions to observationally study the relation between H2O and H2O+ . Nine out of ten sources show absorption from H2O+ in a range of environments: the molecular clumps surrounding the forming and newly formed massive stars, bright high-velocity outflows associated with the massive protostars, and unrelated low-density clouds along the line of sight. Column densities per velocity component of H2 O+ are found in the range of 10^12 to a few 10^13 cm-2 . The highest N(H2O+) column densities are found in the outflows of the sources. The ratios of H2O+/H2O are determined in a range from 0.01 to a few and are found to differ strongly between the observed environments with much lower ratios in the massive (proto)cluster envelopes (0.01-0.1) than in outflows and diffuse clouds. Remarkably, even for source components detected in H2O in emission, H2O+ is still seen in absorption.

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

The Heterodyne Instrument for the Far Infrared (HIFI) onboard the Herschel Space Observatory allows the first observations of light diatomic molecules at high spectral resolution and in multiple transitions. Here, we report deep integrations using HIFI in different lines of hydrides towards the high-mass star forming region AFGL 2591. Detected are CH, CH+, NH, OH+, H2O+, while NH+ and SH+ have not been detected. All molecules except for CH and CH+ are seen in absorption with low excitation temperatures and at velocities different from the systemic velocity of the protostellar envelope. Surprisingly, the CH(JF,P = 3/2_2,- - 1/2_1,+) and CH+(J = 1 - 0, J = 2 - 1) lines are detected in emission at the systemic velocity. We can assign the absorption features to a foreground cloud and an outflow lobe, while the CH and CH+ emission stems from the envelope. The observed abundance and excitation of CH and CH+ can be explained in the scenario of FUV irradiated outflow walls, where a cavity etched out by the outflow allows protostellar FUV photons to irradiate and heat the envelope at larger distances driving the chemical reactions that produce these molecules.

'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.

Water in Star-forming regions with Herschel (WISH) is a Herschel Key Program investigating the water chemistry in young stellar objects (YSOs) during protostellar evolution. Hydroxyl (OH) is one of the reactants in the chemical network most closely linked to the formation and destruction of H2O. High-temperature chemistry connects OH and H2O through the OH + H2 <-> H2O + H reactions. Formation of H2O from OH is efficient in the high-temperature regime found in shocks and the innermost part of protostellar envelopes. Moreover, in the presence of UV photons, OH can be produced from the photo-dissociation of H2O. High-resolution spectroscopy of the OH 163.12 micron triplet towards HH 46 and NGC 1333 IRAS 2A was carried out with the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel. The low- and intermediate-mass YSOs HH 46, TMR 1, IRAS 15398-3359, DK Cha, NGC 7129 FIRS 2, and NGC 1333 IRAS 2A were observed with the Photodetector Array Camera and Spectrometer (PACS) in four transitions of OH and two [OI] lines. The OH transitions at 79, 84, 119, and 163 micron and [OI] emission at 63 and 145 micron were detected with PACS towards the class I low-mass YSOs as well as the intermediate-mass and class I Herbig Ae sources. No OH emission was detected from the class 0 YSO NGC 1333 IRAS 2A, though the 119 micron was detected in absorption. With HIFI, the 163.12 micron was not detected from HH 46 and only tentatively detected from NGC 1333 IRAS 2A. The combination of the PACS and HIFI results for HH 46 constrains the line width (FWHM > 11 km/s) and indicates that the OH emission likely originates from shocked gas. This scenario is supported by trends of the OH flux increasing with the [OI] flux and the bolometric luminosity. Similar OH line ratios for most sources suggest that OH has comparable excitation temperatures despite the different physical properties of the sources.

We performed a sensitive search for the ground-state emission lines of ortho- and para-water vapor in the DM Tau protoplanetary disk using the Herschel/HIFI instrument. No strong lines are detected down to 3sigma levels in 0.5 km/s channels of 4.2 mK for the 1_{10}--1_{01} line and 12.6 mK for the 1_{11}--0_{00} line. We report a very tentative detection, however, of the 1_{10}--1_{01} line in the Wide Band Spectrometer, with a strength of T_{mb}=2.7 mK, a width of 5.6 km/s and an integrated intensity of 16.0 mK km/s. The latter constitutes a 6sigma detection. Regardless of the reality of this tentative detection, model calculations indicate that our sensitive limits on the line strengths preclude efficient desorption of water in the UV illuminated regions of the disk. We hypothesize that more than 95-99% of the water ice is locked up in coagulated grains that have settled to the midplane.

Mon R2, at a distance of 830 pc, is the only ultracompact HII region (UC HII) where the photon-dominated region (PDR) between the ionized gas and the molecular cloud can be resolved with Herschel. HIFI observations of the abundant compounds 13CO, C18O, o-H2-18O, HCO+, CS, CH, and NH have been used to derive the physical and chemical conditions in the PDR, in particular the water abundance. The 13CO, C18O, o-H2-18O, HCO+ and CS observations are well described assuming that the emission is coming from a dense (n=5E6 cm-3, N(H2)>1E22 cm-2) layer of molecular gas around the UC HII. Based on our o-H2-18O observations, we estimate an o-H2O abundance of ~2E-8. This is the average ortho-water abundance in the PDR. Additional H2-18O and/or water lines are required to derive the water abundance profile. A lower density envelope (n~1E5 cm-3, N(H2)=2-5E22 cm-2) is responsible for the absorption in the NH 1_1-0_2 line. The emission of the CH ground state triplet is coming from both regions with a complex and self-absorbed profile in the main component. The radiative transfer modeling shows that the 13CO and HCO+ line profiles are consistent with an expansion of the molecular gas with a velocity law, v_e =0.5 x (r/Rout)^{-1} km/s, although the expansion velocity is poorly constrained by the observations presented here.

We report the first detection of the ground-state rotational transition of the methylidyne cation CH+ towards the massive star-forming region DR21 with the HIFI instrument onboard the Herschel satellite. The line profile exhibits a broad emission line, in addition to two deep and broad absorption features associated with the DR21 molecular ridge and foreground gas. These observations allow us to determine a CH+ J=1-0 line frequency of 835137 +/- 3 MHz, in good agreement with a recent experimental determination. We estimate the CH+ column density to be a few 1e13 cm^-2 in the gas seen in emission, and > 1e14 cm^-2 in the components responsible for the absorption, which is indicative of a high line of sight average abundance [CH+]/[H] > 1.2x10^-8. We show that the CH+ column densities agree well with the predictions of state-of-the-art C-shock models in dense UV-illuminated gas for the emission line, and with those of turbulent dissipation models in diffuse gas for the absorption lines.

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.

SWAS and Odin provided stringent upper limits on the gas phase water abundance of dark clouds (x(H2O) < 7x10^-9). We investigate the chemistry of water vapor in starless cores beyond the previous upper limits using the highly improved angular resolution and sensitivity of Herschel and measure the abundance of water vapor during evolutionary stages just preceding star formation. High spectral resolution observations of the fundamental ortho water (o-H2O) transition (557 GHz) were carried out with Herschel HIFI toward two starless cores: B68, a Bok globule, and L1544, a prestellar core embedded in the Taurus molecular cloud complex. The rms in the brightness temperature measured for the B68 and L1544 spectra is 2.0 and 2.2 mK, respectively, in a velocity bin of 0.59 km s^-1. The continuum level is 3.5+/-0.2 mK in B68 and 11.4+/-0.4 mK in L1544. No significant feature is detected in B68 and the 3 sigma upper limit is consistent with a column density of o-H2O N(o-H2O) < 2.5x10^13 cm^-2, or a fractional abundance x(o-H2O) < 1.3x10^-9, more than an order of magnitude lower than the SWAS upper limit on this source. The L1544 spectrum shows an absorption feature at a 5 sigma level from which we obtain the first value of the o-H2O column density ever measured in dark clouds: N(o-H2O) = (8+/-4)x10^12 cm^-2. The corresponding fractional abundance is x(o-H2O) ~ 5x10^-9 at radii > 7000 AU and ~2x10^-10 toward the center. The radiative transfer analysis shows that this is consistent with a x(o-H2O) profile peaking at ~10^-8, 0.1 pc away from the core center, where both freeze-out and photodissociation are negligible. Herschel has provided the first measurement of water vapor in dark regions. Prestellar cores such as L1544 (with their high central densities, strong continuum, and large envelopes) are very promising tools to finally shed light on the solid/vapor balance of water in molecular clouds.

We are carrying out a physical and chemical study of the protostellar envelopes in a representative sample of IM Class 0 protostars. In our first paper (Crimier et al. 2010), we determined the physical structure (density-temperature radial profiles) of the protostellar envelopes. Here, we study the CO depletion and N2H+ deuteration. We observed the millimeter lines of C18O, C17O, N2H+ and N2D+ toward the protostars using the IRAM 30m telescope. Based on these observations, we derived the C18O, N2H+ and N2D+ radial abundance profiles across their envelopes using a radiative transfer code. In addition, we modeled the chemistry of the protostellar envelopes. All the C18O 1-0 maps are well fit assuming that the C18O abundance decreases inwards within the protostellar envelope until the gas and dust reach the CO evaporation temperature, 20-25K, where the CO is released back to the gas phase. The N2H+ deuterium fractionation in Class 0 IMs is [N2D+]/[N2H+]=0.005-0.014, two orders of magnitude higher than the elemental [D/H] value in the interstellar medium, but a factor of 10 lower than in pre-stellar clumps. Chemical models account for the C18O and N2H+ observations if we assume the CO abundance is 2 times lower than the canonical value in the inner envelope. This could be the consequence of the CO being converted into CH3OH on the grain surfaces prior to the evaporation and/or the photodissociation of CO by the stellar UV radiation. The deuterium fractionation is not fitted by chemical models. This discrepancy is very likely caused by the simplicity of our model that assumes spherical geometry and neglects important phenomena like the effect of bipolar outflows and UV radiation from the star. More important, the deuterium fractionation is dependent on the ortho-to-para H2 ratio, which is not likely to reach the steady-state value in the dynamical time scales of these protostars.

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

Aims: We present preliminary results of the first Herschel spectroscopic observations of NGC7129 FIRS2, an intermediate mass star-forming region. We attempt to interpret the observations in the framework of an in-falling spherical envelope. Methods: The PACS instrument was used in line spectroscopy mode (R=1000-5000) with 15 spectral bands between 63 and 185 microns. This provided good detections of 26 spectral lines seen in emission, including lines of H2O, CO, OH, O I, and C II. Results: Most of the detected lines, particularly those of H2O and CO, are substantially stronger than predicted by the spherical envelope models, typically by several orders of magnitude. In this paper we focus on what can be learned from the detected CO emission lines. Conclusions: It is unlikely that the much stronger than expected line emission arises in the (spherical) envelope of the YSO. The region hot enough to produce such high excitation lines within such an envelope is too small to produce the amount of emission observed. Virtually all of this high excitation emission must arise in structures such as as along the walls of the outflow cavity with the emission produced by a combination of UV photon heating and/or non-dissociative shocks.

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.

The molecular gas in the DR21 massive star formation region is known to be affected by the strong UV field from the central star cluster and by a fast outflow creating a bright shock. The relative contribution of both heating mechanisms is the matter of a long debate. By better sampling the excitation ladder of various tracers we provide a quantitative distinction between the different heating mechanisms. HIFI observations of mid-J transitions of CO and HCO+ isotopes allow us to bridge the gap in excitation energies between observations from the ground, characterizing the cooler gas, and existing ISO LWS spectra, constraining the properties of the hot gas. Comparing the detailed line profiles allows to identify the physical structure of the different components. In spite of the known shock-excitation of H2 and the clearly visible strong outflow, we find that the emission of all lines up to > 2 THz can be explained by purely radiative heating of the material. However, the new Herschel/HIFI observations reveal two types of excitation conditions. We find hot and dense clumps close to the central cluster, probably dynamically affected by the outflow, and a more widespread distribution of cooler, but nevertheless dense, molecular clumps.

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158\mu\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

'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.

This paper is part of a multi-species survey of line emission from the molecular gas in the circum-nuclear disk (CND) of the Seyfert 2 galaxy NGC1068. Single-dish observations have provided evidence that the abundance of silicon monoxide(SiO) in the CND of NGC1068 is enhanced by 3-4 orders of magnitude with respect to the values typically measured in quiescent molecular gas in the Galaxy. We aim at unveiling the mechanism(s) underlying the SiO enhancement. We have imaged with the IRAM Plateau de Bure interferometer the emission of the SiO(2-1) and CN(2--1) lines in NGC1068 at 150pc and 60pc spatial resolution, respectively. We have also obtained complementary IRAM 30m observations of HNCO and methanol (CH3OH) lines. SiO is detected in a disk of 400pc size around the AGN. SiO abundances in the CND of (1-5)xE-09 are about 1-2 orders of magnitude above those measured in the starburst ring. The overall abundance of CN in the CND is high: (0.2-1)xE-07. The abundances of SiO and CN are enhanced at the extreme velocities of gas associated with non-circular motions close to the AGN (r<70pc). Abundances measured for CN and SiO, and the correlation of CN/CO and SiO/CO ratios with hard X-ray irradiation, suggest that the CND of NGC1068 has become a giant X-ray dominated region (XDR). The extreme properties of molecular gas in the circum-nuclear molecular disk of NGC1068 result from the interplay between different processes directly linked to nuclear activity. Whereas XDR chemistry offers a simple explanation for CN and SiO in NGC1068, the relevance of shocks deserves further scrutiny. The inclusion of dust grain chemistry would help solve the controversy regarding the abundances of other molecular species, like HCN, which are under-predicted by XDR models.

Context: Intermediate mass protostars provide a bridge between low- and high-mass protostars. Furthermore, they are an important component of the UV interstellar radiation field. Despite their relevance, little is known about their formation process. Aims: We present a systematic study of the physical structure of five intermediate mass, candidate Class 0 protostars. Our two goals are to shed light on the first phase of intermediate mass star formation and to compare these protostars with low- and high-mass sources. Methods: We derived the dust and gas temperature and density profiles of the sample. We analysed all existing continuum data on each source and modelled the resulting SED with the 1D radiative transfer code DUSTY. The gas temperature was then predicted by means of a modified version of the code CHT96. Results: We found that the density profiles of five out of six studied intermediate mass envelopes are consistent with the predictions of the "inside-out" collapse theory.We compared several physical parameters, like the power law index of the density profile, the size, the mass, the average density, the density at 1000 AU and the density at 10 K of the envelopes of low-, intermediate, and high-mass protostars. When considering these various physical parameters, the transition between the three groups appears smooth, suggesting that the formation processes and triggers do not substantially differ.

Mid-infrared (IR) observations of polycyclic aromatic hydrocarbons (PAHs) and molecular hydrogen emission are a potentially powerful tool to derive physical properties of dense environments irradiated by intense UV fields. We present new, spatially resolved, \emph{Spitzer} mid-IR spectroscopy of the high UV-field and dense photodissocation region (PDR) around Monoceros R2, the closest ultracompact \hII region, revealing the spatial structure of ionized gas, PAHs and H$_2$ emissions. Using a PDR model and PAH emission feature fitting algorithm, we build a comprehensive picture of the physical conditions prevailing in the region. We show that the combination of the measurement of PAH ionization fraction and of the ratio between the H$_2$ 0-0 S(3) and S(2) line intensities, respectively at 9.7 and 12.3 $\mu$m, allows to derive the fundamental parameters driving the PDR: temperature, density and UV radiation field when they fall in the ranges $T = 250-1500 $K, $n_H=10^4-10^6$cm$^{-3}$, $G_0=10^3-10^5$ respectively. These mid-IR spectral tracers thus provide a tool to probe the similar but unresolved UV-illuminated surface of protoplanetary disks or the nuclei of starburst galaxies.

We have carried out high-angular resolution (1.4") observations in the continuum at 3.1mm and in the N2H+ 1-0, CH3CN 5_k-4_k and 13CS 2-1 lines using the Plateau de Bure Interferometer (PdBI) towards the intermediate mass (IM) protostar IRAS21391+5802 (IC1396N). In addition, we have merged the PdBI images with previous BIMA (continuum data at 1.2mm and 3.1mm) and single-dish (N2H+ 1-0) data to have a comprehensive description of the region. The combination of our data with BIMA and 30m data show that the bipolar outflow associated has completely eroded the initial molecular globule. The 1.2mm and 3.1mm continuum emissions are extended along the outflow axis tracing the warm walls of the biconical cavity. Most of the molecular gas, however, is located in an elongated feature in the direction perpendicular to the outflow. A strong chemical differentiation is detected across the molecular toroid, with the N2H+ 1-0 emission absent in the inner region.This chemical differentiation can be understood in terms of the different gas kinetic temperature. The [CH3CN]/[N2H+] ratio increases by 5 orders of magnitude with gas temperature, for temperatures between 20K and 100K. The CH3CN abundance towards IRAM 2A, the most massive protostellar core, is similar to that found in hot corinos and lower than that expected towards IM and high mass hot cores. This could indicate that IRAM 2A is a low mass or at most Herbig Ae star (IRAM 2A) instead of the precursor of a massive Be star. Alternatively, the low CH3CN abundance could also be the consequence of IRAM 2A being a Class 0/I transition object which has already formed a small photodissociation region (PDR).

Water is a key species in many astrophysical environments, but it is particularly important in proto-planetary disks. So far,observations of water in these objects have been scarce, but the situation should soon change thanks to the Herschel satellite. We report here a theoretical study of the water line spectrum of a proto-planetary disk surrounding Ae stars. We show that several lines will be observable with the HIFI instrument onboard the Herschel Space Observatory. We predict that some maser lines could also be observable with ground telescopes and we discuss how the predictions depend not only on the adopted physical and chemical model but also on the set of collisional coefficients used and on the H2 ortho to para ratio through its effect on collisional excitation. This makes the water lines observations a powerful, but dangerous -if misused- diagnostic tool.

With the incorporation of high-J molecular lines, we aim to constrain the physical conditions of the dense gas in the central region of the Seyfert 2 galaxy NGC 1068 and to determine signatures of the AGN or the starburst contribution. We used the James Clerk Maxwell Telescope to observe the J=4-3 transition of HCN, HNC, and HCO+, as well as the CN N_J=2_{5/2}-1_{3/2} and N_J=3_{5/2}-2_{5/2}, in NGC 1068. We estimate the excitation conditions of HCN, HNC, and CN, based on the line intensity ratios and radiative transfer models. We find that the bulk emission of HCN, HNC, CN, and the high-J HCO+ emerge from dense gas n(H_2)>=10^5 cm^-3). However, the low-J HCO+ lines (dominating the HCO+ column density) trace less dense (n(H_2)<10^5 cm^-3) and colder (T_K<=20 K) gas, whereas the high-J HCO+ emerges from warmer (>30 K) gas than the other molecules. The HCO+ J=4-3 line intensity, compared with the lower transition lines and with the HCN J=4-3 line, support the influence of a local XDR environment. The estimated N(CN)/N(HCN)~1-4 column density ratios are indicative of an XDR/AGN environment with a possible contribution of grain-surface chemistry induced by X-rays or shocks.

We completed a search for circumstellar disks around Herbig Be stars using the NRAO Very Large Array (VLA) and the IRAM Plateau de Bure (PdB) interferometers. We present our new VLA and PdBI data for the three objects MWC 297, Z CMa, and LKHa 215. We constructed the SED from near-IR to centimeter wavelengths by adding our millimeter and centimeter data to the available data at other wavelengths, mainly Spitzer images. The entire SED was fitted using a disk+envelope model. In addition, we compiled all the disk millimeter observations in the literature and completed a statistical analysis of all the data. We show that the disk mass is usually only a small percentage (less than 10%) of the mass of the entire envelope in HBe stars. For the disks, there are large source-to-source variations. Two disks in our sample, R Mon and Z CMa, have similar sizes and masses to those found in T Tauri and Herbig Ae stars. The disks around MWC 1080 and MWC 297 are, however, smaller (rout<100 AU). We did not detect the disks towards MWC 137 and LkHa 215 at millimeter wavelengths, which limits the mass and the size of the possible circumstellar disks. A comparison between our data and previous results for T Tauri and Herbig Ae stars indicates that although massive disks (0.1 Msun) are found in young objects (10^4 yr), the masses of the disks around Herbig Be stars are usually 5-10 times lower than those around lower mass stars. We propose that disk photoevaporation is responsible for this behavior. In Herbig Be stars, the UV radiation disperses the gas in the outer disk on a timescale of a few 10^5 yr. Once the outer part of the disk has vanished, the entire gaseous disk is photoevaporated on a very short timescale (10^5 yr) and only a small, dusty disk consisting of large grains remains.

In this Paper we analyze the mid-infrared (mid-IR) emission of very small dust particles in a sample of 12 protoplanetary disks to see how they are connected to interstellar dust particles and to investigate the possibility that their emission can be used as a probe of the physical conditions and evolution of the disk. We define a basis made of three mid-IR template spectra PAH$^0$, PAH$^+$ and VSGs that were derived from the analysis of reflection nebulae, and an additional PAH$^x$ spectrum that was introduced by Joblin et al. (2008) for the analysis of the spectra of planetary nebulae. From the optimization of the fit of 12 star+disk spectra, using a linear combination of the 4 template spectra, we found that an additional small grain component with a broad feature at 8.3 $\mu$m is needed. We find that the fraction of VSG emission in disks decreases with increasing stellar temperature. VSGs appear to be destroyed by UV photons at the surface of disks, thus releasing free PAH molecules, which are eventually ionized as it is observed in photodissociation regions. On the opposite, we observe that the fraction of PAH$^x$ increases with increasing star temperature except in the case of B stars where they are absent. We argue that this is compatible with the identification of PAH$^x$ as large ionized PAHs, most likely emitting in regions of the disk that are close to the star. Finally, we provide a UV-dependant scheme to explain the evolution of PAHs and VSGs in protoplanetary disks. We show that A stars modify the size spectrum of PAHs and VSGs in favor of large PAHs while B stars destroy even the largest PAHs up to large radii in the disk. These results allow us to put new constrains on the properties of two sources: IRS 48 and "Gomez's Hamburger" which are poorly characterized.

We present the first results of a high-resolution study of the distribution and kinematics of molecular gas in two nearby radio galaxies, 4C31.04 and 3C293, representative of two different stages of evolution in radio-loud active galactic nuclei (AGN). These observations, conducted with the IRAM Plateau de Bure Interferometer (PdBI), map with unprecedented spatial resolution (0.5"-1") and sensitivity the emission and absorption of key molecular species such as CO, HCN and HCO+. We report on the detection of a kinematically disturbed and massive (Mgas~10^10Msun) molecular/dusty disk of ~1.4kpc-size fueling the central engine of the compact symmetric object (CSO) 4C31.04. We also report on the detection of a massive (Mgas~10^10Msun) regularly rotating ~7kpc-size disk in the FRII radio galaxy 3C293. A complex system of molecular line absorptions is detected against the mm-continuum source of this galaxy (AGN and jet). We compare the properties of the molecular disks in the two sources and discuss them in the light of the different theories describing the evolution of radio galaxies.

The molecular gas composition in the inner 1 kpc disk of the starburst galaxy M82 resembles that of Galactic Photon Dominated Regions (PDRs). In particular, large abundances of the reactive ions HOC+ and CO+ have been measured in the nucleus of this galaxy. To investigate the origin of the large abundances of reactive ions in M82, we have completed our previous 30m HOC+ J=1-0 observations with the higher excitation HCO+ and HOC+ J=4-3 and 3-2 rotational lines. In addition, we have obtained with the IRAM Plateau de Bure Interferometer (PdBI) a 4" resolution map of the HOC+ 1-0 emission, the first ever obtained in a Galactic or extragalactic source. Our HOC+ interferometric image shows that the emission of the HOC+ 1-0 line is mainly restricted to the nuclear disk, with the maxima towards the E. and W. molecular peaks. In addition, line excitation calculations imply that the HOC+ emission arises in dense gas. Therefore, the HOC+ emission is arising in the dense PDRs embedded in the M82 nuclear disk, rather than in the intercloud phase and/or wind. We have improved our previous chemical model of M82 by (i) using the new version of the Meudon PDR code, (ii) updating the chemical network, and (iii) considering two different types of clouds (with different thickness). Most molecular observations (HCO+, HOC+, CO+, CN, HCN, H3O+) are well explained assuming that ~ 87 % of the mass of the molecular gas is forming small clouds (Av=5 mag) while only ~ 13 % of the mass is in large molecular clouds (Av=50 mag). Such small number of large molecular clouds suggests that M82 is an old starburst, where star formation has almost exhausted the molecular gas reservoir.

We have carried out observations at millimeter and centimeter wavelengths towards VV Ser using the Plateau de Bure Interferometer and the Very Large Array. This allows us to compute the SED from near infrared to centimeter wavelengths. The modeling of the full SED has provided insight into the dust properties and a more accurate value of the disk mass. The mass of dust in the disk around VV Ser is found to be about 4 10^(-5) Msun, i.e. 400 times larger than previous estimates. Moreoever, the SED can only be accounted for assuming dust stratification in the vertical direction across the disk. The existence of small grains (0.25--1 micron) in the disk surface is required to explain the emission at near- and mid-infrared wavelengths. The fluxes measured at millimeter wavelengths imply that the dust grains in the midplane have grown up to very large sizes, at least to some centimeters.

The use of specific tracers of the dense molecular gas phase can help to explore the feedback of activity on the interstellar medium (ISM) in galaxies. This information is a key to any quantitative assessment of the efficiency of the star formation process in galaxies. We present the results of a survey devoted to probe the feedback of activity through the study of the excitation and chemistry of the dense molecular gas in a sample of local universe starbursts and active galactic nuclei (AGNs). Our sample includes also 17 luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs). From the analysis of the LIRGs/ULIRGs subsample, published in Gracia-Carpio et al.(2007) we find the first clear observational evidence that the star formation efficiency of the dense gas, measured by the L_FIR/L_HCN ratio, is significantly higher in LIRGs and ULIRGs than in normal galaxies. Mounting evidence of overabundant HCN in active environments would even reinforce the reported trend, pointing to a significant turn upward in the Kennicutt-Schmidt law around L_FIR=10^11 L_sun. This result has major implications for the use of HCN as a tracer of the dense gas in local and high-redshift luminous infrared galaxies.

We present new observations made with the IRAM 30m telescope of the J=1-0 and 3-2 lines of HCN and HCO^+ used to probe the dense molecular gas content in a sample of 17 local luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs). These observations have allowed us to derive an updated version of the power law describing the correlation between the FIR luminosity (L_FIR) and the HCN(1-0) luminosity (L'_HCN(1-0)) of local and high-redshift galaxies. We present the first clear observational evidence that the star formation efficiency of the dense gas (SFE_dense), measured as the L_FIR/L'_HCN(1-0) ratio, is significantly higher in LIRGs and ULIRGs than in normal galaxies, a result that has also been found recently in high-redshift galaxies. This may imply a statistically significant turn upward in the Kennicutt-Schmidt law derived for the dense gas at L_FIR >= 10^11 L_sun. We have used a one-phase Large Velocity Gradient (LVG) radiative transfer code to fit the three independent line ratios derived from our observations. The results of this analysis indicate that the [HCN]/[HCO^+] abundance ratios could be up to one order of magnitude higher than normal in a significant number of LIRGs and ULIRGs of our sample. An overabundance of HCN at high L_FIR implies that the reported trend in the L_FIR/L'_HCN ratio as a function of L_FIR would be underestimating a potentially more dramatic change of the SFE_dense. Results obtained with two-phase LVG models corroborate that the L'_HCN(1-0)-to-M_dense conversion factor must be lowered at high L_FIR. We discuss the implications of these findings for the use of HCN as a tracer of the dense molecular gas in local and high-redshift luminous infrared galaxies.

In this Paper we present high angular resolution imaging of 3 intermediate-mass (IM) stars using the Plateau de Bure Interferometer (PdBI). In particular we present the chemical study we have carried out towards the IM hot core NGC 7129--FIRS 2. This is the first chemical study in an IM hot core and provides important hints to understand the dependence of the hot core chemistry on the stellar luminosity. We also present our high angular resolution (0.3") images of the borderline Class 0-Class I object IC1396 N. These images trace the warm region of this IM protostar with unprecedent detail (0.3"\sim200 AU at the distance of IC1396 N) and provide the first detection of a cluster of IM hot cores. Finally, we present our interferometric continuum and spectroscopic images of the disk around the Herbig Be star R Mon. We have determined the kinematics and physical structure of the disk associated with this B0 star. The low spectral index derived from the dust emission as well as the flat geometry of the disk suggest a more rapid evolution of the disks associated with massive stars. In the Discussion, we dare to propose a possible evolutionary sequence for the warm circumstellar material around IM stars.

We investigate the mm-morphology of IC1396N with unprecedented spatial resolution to analyze its dust and molecular gas properties, and draw comparisons with objects of similar mass. We have carried out sensitive observations in the most extended configurations of the IRAM Plateau de Bure interferometer, to map the thermal dust emission at 3.3 and 1.3mm, and the emission from the $J$=13$_k\to12_k$ hyperfine transitions of methyl cyanide (CH$_3$CN). We unveil the existence of a sub-cluster of hot cores in IC1396N, distributed in a direction perpendicular to the emanating outflow. The cores are embedded in a common envelope of extended and diffuse dust emission. We find striking differences in the dust properties of the cores ($\beta\simeq$ 0) and the surrounding envelope ($\beta\simeq$ 1), very likely testifying to differences in the formation and processing of dust material. The CH$_3$CN emission peaks towards the most massive hot core and is marginally extended in the outflow direction.

We report the detection of a massive (M(gas) > 5x10^9 Msun) molecular/dusty disk of 1.4kpc-size fueling the central engine of the Compact Symmetric Object (CSO) 4C31.04 based on high-resolution (0.5''--1.2'') observations done with the IRAM Plateau de Bure interferometer (PdBI). These observations allow for the first time to detect and map the continuum emission from dust at 218GHz in the disk of a CSO. The case for a massive disk is confirmed by the detection of strong HCO+(1--0) line emission and absorption. The molecular gas mass of 4C31.04 is in the range 0.5x10^10--5x10^10Msun. While the distribution and kinematics of the gas correspond roughly to those of a rotating disk, we find evidence of distortions and non-circular motions suggesting that the disk is not in a dynamically relaxed state. We discuss the implications of these results for the general understanding of the evolution of radio galaxies.

The transition between the low density groups of T Tauri stars and the high density clusters around massive stars occurs in the intermediate-mass (IM) range (M$_*$$\sim$2--8 M$_\odot$). High spatial resolution studies of IM young stellar objects (YSO) can provide important clues to understand the clustering in massive star forming regions. Aims: Our aim is to search for clustering in IM Class 0 protostars. The high spatial resolution and sensitivity provided by the new A configuration of the Plateau de Bure Interferometer (PdBI) allow us to study the clustering in these nearby objects. Methods: We have imaged three IM Class 0 protostars (Serpens-FIRS 1, IC 1396 N, CB 3) in the continuum at 3.3 and 1.3mm using the PdBI. The sources have been selected with different luminosity to investigate the dependence of the clustering process on the luminosity of the source. Results: Only one millimeter (mm) source is detected towards the low luminosity source Serpens--FIRS 1. Towards CB 3 and IC1396 N, we detect two compact sources separated by $\sim$0.05 pc. The 1.3mm image of IC 1396 N, which provides the highest spatial resolution, reveal that one of these cores is splitted in, at least, three individual sources.

We present new high angular resolution observations of the 13CO 1-0 rotational line towards the HBe star R Mon, obtained with the IRAM Bure Interferometer, along with the previous results in the observed transitions 12CO 1-0 and 2-1. We have used a flat disk model to fit the 12CO 1-0 and 2-1, and 13CO 1-0 emission, in a strip perpendicular to the outflow axis. The model assumes standard abundances (X(12CO) = 8x10^-5, X(13CO) = 9x10^7), radial potential temperature and density laws (T (K) = T0 r^-q, rho (cm^-3) = rho0 r^-p), and Local Thermodynamic Equilibrium. The 13CO and 12CO emission is consistent with a flat disk at an inclination angle of 20 degrees in Keplerian rotation around the star. The gaseous disk is fitted with a mass of 0.014 Msun, an outer radius of 1500 AU, a temperature of 4500 K at the inner radius (1 AU), and values of q=0.62 and p=1.3 for the indexes of the temperature and density laws. This new 13CO observations allow us to conclude that the disk around R Mon is flat. Our result confirms previous works suggesting a predominant flat geometry in disks around early Be stars.

We present high-angular resolution observations of the circumstellar disk around the massive Herbig Be star R Mon (M~8 Msun) in the continuum at 2.7mm and 1.3mm and the CO 1->0 and 2->1 rotational lines. Based on the new 1.3mm continuum image we estimate a disk mass (gas+dust) of 0.007 Msun and an outer radius of <150 AU. Our CO images are consistent with the existence of a Keplerian rotating gaseous disk around this star. Up to our knowledge, this is the most clear evidence for the existence of Keplerian disks around massive stars reported thus far. The mass and physical characteristics of this disk are similar to thoseof the more evolved T Tauri stars and indicate a shorter timescale for the evolution and dispersal of circumstellar disks around massive stars which lose most of their mass before the star becomes visible.

Observational constraints provided by high resolution and high sensitivity observations of external galaxies made in the millimeter and submillimeter range have started to put on a firm ground the study of extragalactic chemistry of molecular gas. In particular, the availability of multi-species and multi-line surveys of nearby galaxies is central to the interpretation of existent and forthcoming millimeter observations of the high redshift universe. Probing the physical and chemical status of molecular gas in starbursts and active galaxies (AGN) requires the use of specific tracers of the relevant energetic phenomena that are known to be at play in these galaxies: large-scale shocks, strong UV fields, cosmic rays and X-rays. We present below the first results of an ongoing survey, allying the IRAM 30m telescope with the Plateau de Bure interferometer(PdBI), devoted to study the chemistry of molecular gas in a sample of starbursts and AGN of the local universe. These observations highlight the existence of a strong chemical differentiation in the molecular disks of starbursts and AGN.

We present the detection of the reactive ion CO+ towards the prototypical starburst galaxy M82. This is the first secure detection of this short-lived ion in an external galaxy. Values of [CO+]/[HCO+]>0.04 are measured across the inner 650pc of the nuclear disk of M82. Such high values of the [CO+]/[HCO+] ratio had only been previously measured towards the atomic peak in the reflection nebula NGC7023. This detection corroborates that the molecular gas reservoir in the M82 disk is heavily affected by the UV radiation from the recently formed stars. Comparing the column densities measured in M82 with those found in prototypical Galactic photon-dominated regions (PDRs), we need \~20 clouds along the line of sight to explain our observations. We have completed our model of the molecular gas chemistry in the M82 nucleus. Our PDR chemical model successfully explains the [CO+]/[HCO+] ratios measured in the M~82 nucleus but fails by one order of magnitude to explain the large measured CO+ column densities (~1--4x10^{13} cm^{-2}). We explore possible routes to reconcile the chemical model and the observations.

(Abridged) We present the first high-resolution (5.6''x5.1'') images of the emission of silicon monoxide (SiO) in the nearby spiral IC342, obtained with the IRAM Plateau de Bure Interferometer. We have mapped the emission of the SiO(2-1) and H13CO+(1-0) lines in a 0.9kpcx1.3kpc region around the nucleus of IC342. The bulk of the line emission comes from a 290pc spiral arm and a r~80pc nuclear ring. The distribution of SiO clouds is mostly anticorrelated with the continuum emission at 87GHz, dominated by thermal free-free bremsstrahlung. The SiO/H13CO+ intensity ratio increases by an order of magnitude from the nuclear ring to the spiral arm. Furthermore the gas kinematics show significant differences between SiO and H13CO+ over the spiral arm, where the linewidths of SiO are 2 times larger than that of H13CO+. The average abundance of SiO in the inner r~320pc, X(SiO)>2x10^{-10}, evidences that shock chemistry is at work in IC342. We have compared the emission of SiO with another tracer of molecular shocks: the emission of methanol (CH3OH). The significant difference of the X(SiO) measured between the spiral arm (~a few 10^{-9}) and the nuclear ring (~10^{-10}) is not echoed by a comparable variation in the SiO/CH3OH intensity ratio. This implies that the typical shock velocities should be similar in both regions. In contrast, the fraction of shocked molecular gas should be 5-7 times larger in the spiral arm (up to ~10% of the available molecular gas mass over the arm region) compared to the nuclear ring. We conclude that the large-scale shocks revealed by the SiO map of IC342 are mostly unrelated with star formation and arise instead in a pre-starburst phase. Shocks are driven by cloud-cloud collisions along the potential well of the IC342 bar. The general implications for the current understanding of galaxy evolution are discussed.

Aiming to learn about the chemistry of the dense PDR around the ultracompact (UC) HII region in Mon R2, we have observed a series of mm-wavelength transitions of C3H2 and C2H. In addition, we have traced the distribution of other molecules, such as H13CO+, SiO, HCO, and HC3N. These data, together with the reactive ions recently detected, have been considered to determine the physical conditions and to model the PDR chemistry. We then identified two kind of molecules. The first group, formed by the reactive ions (CO+, HOC+) and small hydrocarbons (C2H, C3H2), traces the surface layers of the PDR and is presumably exposed to a high UV field (hence we called it as "high UV", or HUV). HUV species is expected to dominate for visual absorptions 2 < Av < 5 mag. A second group (less exposed to the UV field, and hence called "low UV", or LUV) includes HCO and SiO, and is mainly present at the edges of the PDR (Av > 5 mag). While the abundances of the HUV molecules can be explained by gas phase models, this is not the case for the studied LUV ones. Although some efficient gas-phase reactions might be lacking, grain chemistry sounds like a probable mechanism able to explain the observed enhancement of HCO and SiO. Within this scenario, the interaction of UV photons with grains produces an important effect on the molecular gas chemistry and constitutes the first evidence of an ionization front created by the UC HII region carving its host molecular cloud. The physical conditions and kinematics of the gas layer which surrounds the UC HII region were derived from the HUV molecules. Molecular hydrogen densities > 4 10^6 cm^(-3) are required to reproduce the observations. Such high densities suggest that the HII region could be pressure-confined by the surrounding high density molecular gas.

We report high angular resolution (~0.6''x0.5'' at 1.3mm) observations of the Class 0 intermediate-mass (IM) protostar NGC 7129--FIRS 2 using the PdBI. Our observations show the existence of a hot core at the position of the Class 0 object. This is, up to our knowledge, the first IM hot core detected so far. Interferometric maps of the region in the the continuum at 3mm and 1.3mm and in the CH3CN 5_k-->4_k, CH3OH 5_kk'-->4_kk', and D2CO 4_04-->3_03 lines are presented in this paper. Enhanced CH3CN and CH3OH abundances are measured towards the hot core (X(CH3CN)~7.010^{-9}, X(CH3OH)~3 10^{-8}-- a few 10^{-7}). While intense D2CO emission is detected towards the hot core, the N2D+ line has not been detected in our interferometric observations. The different behaviors of D2CO and N2D+ emissions suggest different formation mechanisms for the two species and different deuteration processes for H2CO and N2H+ (surface and gas-phase chemistry, respectively). Finally, the spectrum of the large bandwidth correlator show a forest of lines at the hot core position revealing that this object is extraordinarily rich in complex molecules.

The nucleus of M82 has been mapped in several 3mm and 1mm lines of CN, HCN, C2H, c-C3H2, CH3C2H, HC3N and HOC+ using the IRAM 30m telescope. These species have been purposely selected as good tracers of photon-dominated chemistry. We have measured [CN]/[HCN] ~ 5 in the inner 650 pc galaxy disk. Furthermore, we have detected the HOC+ 1--0 line with an intensity similar to that of the H13CO+ 1--0 line. This implies a [HCO+]/[HOC+] ratio of ~40. These results corroborate the existence of a giant photo-dissociation region (PDR) in the nucleus of M82. In fact, the low [HCO+]/[HOC+] ratio can only be explained if the nucleus of M82 is formed by small (r<0.02-0.2 pc) and dense (n ~ a few 10^4--10^5 cm^{-3}) clouds immersed in an intense UV field (G_0 ~ 10^4 in units of the Habing field). The detection of the hydrocarbons c-C3H2 and CH3C2H in the nucleus of M82 suggests that a complex carbon chemistry is developing in this giant PDR.

We have carried out a molecular survey of the Class 0 IM protostar NGC 7129 -- FIRS 2 (hereafter FIRS 2) and the Herbig Be star LkH$\alpha$ 234 with the aim of studying the chemical evolution of the envelopes of intermediate-mass (IM) young stellar objects (YSOs). Both objects have similar luminosities (~500 Lsun) and are located in the same molecular cloud which minimizes the chemical differences due to different stellar masses or initial cloud conditions. Moreover, since they are located at the same distance, we have the same spatial resolution in both objects. A total of 17 molecular species (including rarer isotopes) have been observed in both objects and the structure of their envelopes and outflows is determined with unprecedent detail. Our results show that the protostellar envelopes are dispersed and warmed up during the evolution to become a pre-main sequence star. In fact, the envelope mass decreases by a factor >5 from FIRS 2 to LkH$\alpha$234, while the kinetic temperature increases from ~13K to 28K. On the other hand, there is no molecular outflow associated with LkH$\alpha$234. The molecular outflow seems to stop before the star becomes visible. These physical changes strongly affect the chemistry of their envelopes. Based on our results in FIRS2 and LkH$\alpha$ 234, we propose some abundance ratios that can be used as chemical clocks for the envelopes of IM YSOs. The SiO/CS, CN/N2H+, HCN/N2H+, DCO+/HCO+ and D2CO/DCO+ ratios are good diagnostics of the protostellar evolutionary stage.

The 500 central pc of the Galaxy (hereafter GC) exhibit a widespread gas component with a kinetic temperature of 100-200 K. The bulk of this gas is not associated to the well-known thermal radio continuum or far infrared sources like Sgr A or Sgr B. How this gas is heated has been a longstanding problem. With the aim of studying the thermal balance of the neutral gas and dust in the GC, we have observed 18 molecular clouds located at projected distances far from thermal continuum sources with the Infrared Space Observatory (ISO). In this paper we present observations of several fine structure lines and the full continuum spectra of the dust between 40 and 190 microns. A warm dust component with a temperature between 27 and 42 K is needed to fit the spectra. We have compared the gas and the dust emission with the predictions from J-type and C-type shocks and photodissociation region (PDRs) models. We conclude that the dust and the fine structure lines observations are best explained by a PDR with a density of 10$^3$ cm^-3 and an incident far-ultraviolet field 10$^3$ times higher than the local interstellar radiation field. PDRs can naturally explain the discrepancy between the gas and the dust temperatures. However, these PDRs can only account for 10-30% of the total H2 column density with a temperature of ~ 150 K. We discuss other possible heating mechanisms (short version).

There is observational evidence that nuclear winds and X-rays can heavily influence the physical conditions and chemical abundances of molecular gas in the circumnuclear disks (CND) of Active Galactic Nuclei (AGN). In this paper we probe the chemical status of molecular gas in the CND of NGC1068, a prototypical Seyfert 2 galaxy. Precedent claims that the chemistry of molecular gas in the nucleus of NGC1068 is abnormal by galactic standards were based on the high HCN/CO luminosity ratio measured in the CND. Results from new observations obtained in this survey have served to derive abundances of molecular species such as SiO, CN, HCO+, HOC+, H13CO+ and HCO. These estimates are complemented by a re-evaluation of molecular abundances for HCN, CS and CO, based on previously published single-dish and interferometer observations of NGC1068. We report on the first detection of SiO emission in the CND of NGC1068. The estimated large abundance of SiO in the CND, X(SiO)~5-10x10^{-9}, cannot be attributed to shocks related to star formation, as there is little evidence of a recent starburst in the nucleus of NGC1068. Alternatively, we propose that silicon chemistry is driven by intense X-ray processing of molecular gas. We also report on the first extragalactic detection of the reactive ion HOC+. Most remarkably, the estimated HCO+/HOC+ abundance ratio in the nucleus of NGC1068, ~30-80, is the smallest ever measured in molecular gas. The abundances derived for all molecules that have been the subject of this survey are compared with the predictions of models invoking either oxygen-depletion or X-ray chemistry in molecular gas. Our conclusions favour an overall scenario where the CND of NGC1068 has become a giant X-ray Dominated Region (XDR).

We have studied the feedback influence that the central engine of the Seyfert 2 galaxy NGC 1068 may have on the chemistry of the 200 pc circumnuclear gas disk (CND). With this purpose, we have conducted a multi-species/multi-transition survey of molecular gas in the CND of NGC 1068 using the IRAM 30m telescope. Abundances of several molecular species have been estimated, including HCN, CN, CS, HCO+, SiO and HOC+. We report on the detection of significant SiO emission in this galaxy, as well as on the first extragalactic detection of the active radical HOC+. We conclude that the chemistry of the molecular gas reservoir in the CND can be best explained in the framework of X-rays Dominated Regions (XDR) models.

We present the main results of an extragalactic survey aiming to study the chemistry of molecular gas in a limited sample of starburst galaxies (SB) and AGN hosts. Observations have been carried out with the IRAM 30m telescope and the Plateau de Bure Interferometer (PdBI). The high resolution/sensitivity of the PdBI has made possible to obtain high quality images of the galaxies using specific molecular gas tracers of Shock Chemistry, Photon Dominated Regions (PDR) and X-ray Dominated Regions (XDR). The occurrence of large-scale shocks and the propagation of PDR chemistry in starbursts can be studied. We also discuss the onset of XDR chemistry in AGN.

We have carried out a high-sensitivity search for circumstellar disks around Herbig Be stars in the continuum at 1.4mm and 2.7mm using the IRAM interferometer at the Plateau de Bure (PdBI) . In this letter, we report data on three well studied B0 stars, MWC 1080, MWC 137 and R Mon. The two latter have also been observed in the continuum at 0.7 cm and 1.3 cm using the NRAO Very Large Array (VLA) . We report the detection of circumstellar disks around MWC 1080 and R Mon with masses of Md ~ 0.003 and 0.01 Msun, respectively, while for MWC 137 we estimate a disk mass upper limit of 0.007 Msun. Our results show that the ratio Md/M* is at least an order of magnitude lower in Herbig Be stars than in Herbig Ae and T Tauri stars.

We report the first detection of the reactive ions CO+ and HOC+ towards ultracompact (UC) HII regions, particularly in Mon R2 and G29.96-0.02. We have observed two positions in Mon R2, namely the peak of the UC HII region and a position in the high density molecular cloud which bounds it. CO+ and HOC+ were detected at the UC HII region but not at the molecular cloud, as expected if the CO$^+$ and HOC$^+$ emissions arise in the PDR surrounding the \uch. The measured CO$^+$ and HOC$^+$ column densities are of the order of 10$^{11}$ cm$^{-2}$ in both sources, which yields a strikingly low [HCO$^+$]/[HOC$^+$] abundance ratio of 460 in Mon R2. These values are similar to those found in some other well-known PDRs, like NGC 7023 and the Orion Bar. We briefly discuss the chemical implications of these results.

We present infrared and millimeter observations of molecular gas, dust and ionized gas towards a sample of clouds distributed along the 500 central pc of the Galaxy. The clouds were selected to investigate the physical state, in particular the high gas temperatures, of the Galactic center region (GCr) clouds located far from far-infrared of thermal radio continuum sources. We have found that there is ionized gas associated with the molecular gas. The ionizing radiation is hard (~35000 K) but diluted due to the inhomogeneity of the medium. We estimate that ~30 % of the warm molecular gas observed in the GCr clouds is heated by ultra-violet radiation in photo-dissociation regions.

We present a high-resolution (~ 5'') image of the nucleus of M82 showing the presence of widespread emission of the formyl radical (HCO). The HCO map, the first obtained in an external galaxy, reveals the existence of a structured disk of ~ 650 pc full diameter. The HCO distribution in the plane mimics the ring morphology displayed by other molecular/ionized gas tracers in M82. More precisely, rings traced by HCO, CO and HII regions are nested, with the HCO ring lying in the outer edge of the molecular torus. Observations of HCO in galactic clouds indicate that the abundance of HCO is strongly enhanced in the interfaces between the ionized and molecular gas. The surprisingly high overall abundance of HCO measured in M82 (X(HCO) ~ 4x10^{-10}) indicates that its nuclear disk can be viewed as a giant Photon Dominated Region (PDR) of ~ 650 pc size. The existence of various nested gas rings, with the highest HCO abundance occurring at the outer ring (X(HCO) ~ 0.8x10^{-9}), suggests that PDR chemistry is propagating in the disk. We discuss the inferred large abundances of HCO in M82 in the context of a starburst evolutionary scenario, picturing the M82 nucleus as an evolved starburst.

We present some recent studies on the physical conditions and the chemistry of the interstellar medium (ISM) in the Galactic center region (GCR). Large scale shocks in the context of a bar potential can play a role in explaining the observed properties, at least in some particular regions. Alternatively, we propose that the observed properties (inhomogeneity, high temperatures of the molecular gas, rich chemistry,...) of the GCR ISM can be the consequence of the recent past ($10^6-10^7$ years ago) of massive star formation in the GCR.

In this Letter we present the first images of the emission of SiO and H13CO+ in the nucleus of the starburst galaxy M82. Contrary to other molecular species that mainly trace the distribution of the star-forming molecular gas within the disk, the SiO emission extends noticeably out of the galaxy plane. The bulk of the SiO emission is restricted to two major features. The first feature, referred to as the SiO supershell, is an open shell of 150 pc diameter, located 120 pc west from the galaxy center. The SiO supershell represents the inner front of a molecular shell expanding at 40 km/s, produced by mass ejection around a supercluster of young stars containing supernova remnant SNR 41.95+57.5. The second feature is a vertical filament, referred to as the SiO chimney, emanating from the disk at 200 pc east from the galaxy center. The SiO chimney reaches a 500 pc vertical height, and it is associated with the most prominent chimney identified in radio continuum maps. The kinematics, morphology, and fractional abundances of the SiO gas features in M82 can be explained in the framework of shocked chemistry driven by local episodes of gas ejection from the starburst disk.

We present molecular line observations of the star-forming cloud around RNO6 along with a newly discovered nearby molecular cloud that we name RNO6NW. Both clouds display striking similarities in their cometary structures and overall kinematics. By using 13CO line observations, we estimate that these clouds have similar sizes (~4.5 pc) and masses (~200 solar masses). Both molecular clouds RNO6 and RNO6NW are active in star formation. From new high resolution near-IR narrowband images, we confirm that RNO6 hosts an embedded IR cluster that includes a Herbig Be star. A conspicuous H2 filament is found to delineate the dense cometary head of the globule. RNO6NW hosts at least two IR sources and a bipolar molecular outflow of ~0.9 pc of length and ~0.5 solar masses. We show that the cometary structure of both clouds has been created by the UV radiation from numerous OB stars lying ~1.5 degree to the north. Such OB stars are associated with the double cluster h and chi Persei, and are probably members of the PerOB1 association. Thus star formation inside these clouds has been very likely triggered by the Radiation Driven Implosion (RDI) mechanism. From comparison to RDI theoretical models, we find that the similar kinematics and morphology of both clouds is well explained if they are at a re-expansion phase. Triggered sequential star formation also explains the observed spatial distribution of the members of the near-IR cluster inside the RNO6 cloud, and the morphology of the H2 filament. We conclude that the RNO6 and RNO6NW clouds are high-mass counterparts to the cometary globules of smaller masses which have been studied up to now. Thus our observations demonstrate that the RDI mechanism can produce, not only low mass stars in small globules, but also intermediate mass stars and clusters in massive clouds.

We present HHT observations in the N=3-2 rotational transition of the CN radical toward selected positions of the Trapezium region and of the molecular Ridge in the Orion molecular cloud. Two of the positions in the Ridge were also observed in the N=2-1 line of CN and 13CN. The N=3-2 CN lines have been combined with observations of the N=2-1 and N=1-0 transitions of CN, and of the N=2-1 of 13CN to estimate the physical conditions and CN abundances in the molecular gas. We analyze in detail the excitation of the CN lines and find that the hyperfine ratios of the N=3-2 line are always close to the Local Thermodynamic Equilibrium (LTE) optically thin values even in the case of optically thick emission. This is due to different excitation temperatures for the different hyperfine lines. From the line intensity ratios between the different CN transitions we derive H_2 densities of ~10^5cm^-3 for the molecular Ridge and of ~3x10^6cm^-3 for the Trapezium region. The CN column densities are one order of magnitude larger in the Ridge than in the Trapezium region, but the CN to H_2 ratio is similar both in the Trapezium and in the Ridge. The combination of the low CN column densities, high H_2 densities and relatively high CN abundances toward the Trapezium region requires that the CN emission arises from a thin layer with a depth along the line of sight of only \~5x10^15cm. This high density thin layer of molecular gas seems to be related with material that confines the rear side of the HII region Orion A. However the molecular layer is not moving as expected from the expansion of the HII region, but it is ``static'' with respect to the gas in the molecular cloud. We discuss the implication of a high density ``static'' layer in the evolution of an HII region.

We present observations of C2H5OH toward molecular clouds in Sgr A, Sgr B2 and associated with thermal and non-thermal features in the Galactic center. C2H5OH emission in Sgr A and Sgr B2 is widespread, but not uniform. C2H5OH emission is much weaker or it is not detected in some molecular clouds in both complexes, in particular those with radial velocities between 70 and 120 km/s. While most of the clouds associated with the thermal features do not show C2H5OH emission, that associated with the Non-Thermal Radio Arc shows emission. The fractional abundance of C2H5OH in most of the clouds with radial velocities between 0 and 70 km/s in Sgr A and Sgr B2 is relatively high, of few 1e-8. The C2H5OH abundance decreases by more than one order of magnitude (aprox. 1e-9) in the clouds associated with the thermal features. The large abundance of C2H5OH in the gas-phase indicates that C2H5OH has formed in grains and released to gas-phase by shocks in the last aprox. 1e5 years.

We present ISO observations of several H2 pure-rotational lines (from S(0) to S(5)) towards a sample of 16 molecular clouds distributed along the central ~ 500 pc of the Galaxy. We also present C18O and 13CO J=1->0 and J=2->1 observations of these sources made with the IRAM-30m telescope. With the CO data we derive H2 densities of 10e(3.5-4.0) cm-3 and H2 column densities of a few 10e22 cm-2. We have corrected the H2 data for ~ 30 magnitudes of visual extinction using a self-consistent method. In every source, we find that the H2 emission exhibits a large temperature gradient. The S(0) and S(1) lines trace temperatures (T) of ~150 K while the S(4) and S(5) lines indicate temperatures of ~ 600K. The warm H2 column density is typically ~1-2 x 10e22 cm-2, and is predominantly gas with T=150 K. This is the first direct estimate of the total column density of the warm molecular gas in the Galactic center region. These warm H2 column densities represent a fraction of ~ 30 % of the gas traced by the CO isotopes emission. The cooling by H2 in the warm component is comparable to that by CO. Comparing our H2 and CO data with available ammonia NH3 observations from literature one obtains relatively high NH3 abundances of a few 10e(-7) in both the warm and the cold gas. A single shock or Photo-Dissociation Region (PDR) cannot explain all the observed H2 lines. Alternatives for the heating mechanisms are discussed.

In this letter we present a map of the GC in the 1-0 line of SiO covering the same region than the ASCA observations of the Fe 6.4 keV line. The SiO emission is restricted to molecular clouds with radial velocity between 10 and 60 km/s in the Sgr A and Sgr B complexes. We find a correlation between the SiO morphology and the spatial distribution of the Fe 6.4 keV line, on the large scale and also within Sgr A and Sgr B. The SiO abundance increases by a factor of 20 in the regions with strong Fe 6.4 keV line. This indicates that the Fe 6.4 keV line mainly arises from molecular clouds with large gas phase abundance of refractory elements. We discuss the implications of the correlation on the origin of the hard X-rays, and the heating and the chemistry of the molecular clouds in the GC

We present ISO observations of the S(0), S(1), S(2), and S(3) rotational lines of molecular hydrogen from two molecular clouds near the Galactic Center (GC). We have also measured continuum dust emission at infrared wavelengths with ISO and the rotational radio lines J=1-0 of 13CO and C18O and J=2-1 of C18O with the IRAM-30m telescope. Using the dust continuum spectra and the CO lines we derive a total visual extinction of 15-20 magnitudes toward these GC clouds. After correcting the H2 data for extinction, the gas temperatures are \~250 K and the column densities of warm gas are ~2e21 cm-2. This is the first direct measure of the H2 column densities of the warm component; with this, we estimate an NH3 abundance in the warm gas of ~ 2e-7. The column density of warm gas is, at least, a factor of 100 larger than the corresponding column densities derived from the warm dust. The observed ortho-to-para ratio (OTPR) is ~ 1, clearly below the local thermodynamical equilibrium (LTE) OTPR for gas at 250 K of ~3. Low velocity shocks (~ 10 km/s) are the most likely explanation for the column densities of warm gas and dust and the non-LTE H2 OTPR.

This paper presents the first high-resolution SiO map made in an external galaxy. The nucleus of the nearby barred spiral NGC253 has been observed simultaneously in the v=0, J=2-1 line of SiO and in the J=1-0 line of H13CO+ with the IRAM interferometer, with a resolution of 7.5"x2.6". Emission from SiO and H13CO+ is extended in the nucleus of NGC253. The bulk of the SiO/H13CO+ emission arises from a 600pcx250pc circumnuclear disk (CND) with a double ringed structure. The inner ring, of radius r=60pc, viewed edge-on along PA=51deg, hosts the nuclear starburst; the outer pseudo-ring opens out as a spiral-like arc up to r=300pc. The kinematics of the gaseous disk, characterized by strong non-circular motions, is interpreted in terms of the resonant response of the gas to the barred potential. The inner ring would correspond to the inner Inner Lindblad Resonance(iILR), whereas the outer region is linked to the onset of a trailing spiral wave across the outer ILR (oILR). Most notably, we report the detection of a molecular gas counterpart of the giant outflow of hot gas, previously seen in X-ray and optical lines, and tentatively identified as a dust chimney in the the 450mu continuum band.The SiO shows a high average fractional abundance in the CND of X(SiO)=1.5x10E-10. This is more than an order of magnitude above the predicted value of a PDR. Moreover, X(SiO) varies at least by an order of magnitude between the inner starburst region, which dominates the global emission, where we derive X(SiO)=1--2x10E-10, and the outer region, where X(SiO) reaches a few 10E-9. SiO abundance is also significantly enhanced in the outflow (X(SiO)=3--5x10E-10). Different mechanisms are explored to explain the unlike chemical processing of molecular gas within the nucleus.

We present SWS and LWS ISO observations towards a strip across the photodissociation region (PDR) in the reflection nebula NGC 7023. SWS02 and LWS01 spectra have been taken towards the star and the brightest infrared filaments located NW and SW from the star (hereafter NW and SW PDRs). In addition, SWS02 spectra have been taken towards two intermediate positions. While all the PDR tracers (the [CII] 157.7$\mic$, [OI] 63.2 and 145.6$\mic$, [HI] 21cm and the H$_2$ rotational lines) present a ring-like morphology with the peaks toward the NW and SW PDRs and a minimum around the star, the SiII emission is filling the hole of this ring with the peak towards the star. This morphology can only be explained if the SiII emission arise in the lowest extinction layers of the PDR (A$_v$ $<$ 2 mag) and the HII region. At least 20% -- 30% of the Si must be in gas phase in these layers. For A$_v$ $\geq$ 2 mag, the Si is mainly in solid form ($\delta$ Si = -1.3). In both, the NW and SW PDRs, the intensities of the H$_2$ rotational lines can only be fitted by assuming an ortho-to-para-H$_2$ ratio lower than 3 in gas with rotation temperatures from 400 to 700 K. Therefore, there is a non-equilibrium OTP ratio in the region. The comparison between the OTP ratio derived from H$_2$ vibrational lines and the pure H$_2$ rotational lines, shows that the OTP ratio increases from $\sim$ 1.5 to 3 across the photodissociation region with larger values in the less shielded gas (A$_v$$<$ 0.7 mag). This behavior is interpreted as a consequence of an advancing photodissociation front.

We have observed the S(0), S(1), S(2), S(3), S(4) and S(5) rotational lines of molecular hydrogen (H2) towards the peak of the photodissociation region (PDR) associated with the reflection nebula NGC 7023. The observed H2 line ratios show that they arise in warm gas with kinetic temperatures ~300 - 700 K. However, the data cannot be fitted by an ortho- to para- (OTP) ratio of 3. An OTP ratio in the range ~1.5 - 2 is necessary to explain our observations. This is the first detection of a non-equilibrium OTP ratio measured from the H2 pure-rotational lines in a PDR. The existence of a dynamical PDR is discussed as the most likely explanation for this low OTP ratio.

We have mapped the J=1-0 line of SiO in a 1 degree x 12 arcmin. (lxb) region around the Galactic center (GC) with an angular resolution of 2 arcmin. (approx. 4 pc). In contrast to the spatial distribution of other high dipole moment molecules like CS, whose emission is nearly uniform, the SiO emission is very fragmented and it is only associated with some molecular clouds. In particular, it is remarkable that the SiO emission closely follows the non-thermal radio arc in the GC. The SiO clouds are more extended than the beam with typical sizes between 4 and 20 pc. High angular resolution (26 arcsec.) mapping in the J=1-0 line of SiO toward the molecular clouds in Sgr B2 and Sgr A shows that the SiO emission is relatively smooth with structures of typically 2 pc. From the line intensities of the J=1-0, J=3-2 and J=5-4 transitions of SiO we derive H2 densities for these clouds of a few 1e4 cm-3. The SiO fractional abundances are approx. 1e-9 for the SiO clouds and 1e-10 for the other molecular clouds in the GC. The characteristics (size and H2 densities) of the SiO emission in the GC are completely different from those observed in the Galactic disk, where the SiO emission arises from much smaller regions with larger H2 densities. We briefly discuss the implications of the SiO emission in the molecular clouds of the GC. We conclude that the particular chemistry in these clouds is probably related to large scale fast shocks occurring in the Galactic center region.

We have detected CO+ toward the photon-dominated region (PDR) associated with the reflection nebula NGC 7023. This is the first detection of CO+ in the vicinity of a Be star. A CO+ column density of ~ 3E11 cm-2 has been derived toward the PDR peak. We have, however, not detected CO+ in a well shielded clump of the adjacent molecular cloud, where the CO+/HCO+ abundance ratio is at least 100 times lower than in the PDR. Our results show, for the first time, that CO^+ column densities as large as ~ 3E11 cm-2 can be produced in regions with incident UV fields of just a few 1E3 (in units of Habing field) and densities of ~ 1E5 cm-3. Furthermore, since the ionization potential of CO is larger than 13.6 eV, our data rule out the direct photoionization of CO as a significant CO+ formation mechanism.