We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

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.

We present observations of SO and $\rm SO_2$ lines toward the shocked regions along the L1157 chemically rich outflow, taken in the context of the Seeds Of Life In Space IRAM-NOrthern Extended Millimeter Array Large Program, and supported by data from Submillimeter Array and IRAM-30 m telescope at 1.1--3.6 mm wavelengths. We simultaneously analyze, for the first time, all of the brightest shocks in the blueshifted lobe, namely, B0, B1, and B2. We found the following. (1) SO and $\rm SO_2$ may trace different gas, given that the large(-scale) velocity gradient analysis indicates for $\rm SO_2$ a volume density ($\rm 10^5\text{--}10^6\,cm^{-3}$) denser than that of the gas emitting in SO by a factor up to an order of magnitude. (2) Investigating the 0.1 pc scale field of view, we note a tentative gradient along the path of the precessing jet. More specifically, $\rm \chi({SO/SO_2})$ decreases from the B0-B1 shocks to the older B2. (3) At a linear resolution of 500--1400 au, a tentative spatial displacement between the two emitting molecules is detected, with the SO peak closer (with respect to $\rm SO_2$) to the position where the recent jet is impinging on the B1 cavity wall. Our astrochemical modeling shows that the SO and $\rm SO_2$ abundances evolve on timescales less than about 1000 years. Furthermore, the modeling requires high abundances ($2\times10^{-6}$) of both $\rm H_2S/H$ and S/H injected in the gas phase due to the shock occurrence, so pre-frozen OCS only is not enough to reproduce our new observations.

We discuss the science drivers for ALMA Band 2 which spans the frequency range from 67 to 90 GHz. The key science in this frequency range are the study of the deuterated molecules in cold, dense, quiescent gas and the study of redshifted emission from galaxies in CO and other species. However, Band 2 has a range of other applications which are also presented. The science enabled by a single receiver system which would combine ALMA Bands 2 and 3 covering the frequency range 67 to 116 GHz, as well as the possible doubling of the IF bandwidth of ALMA to 16 GHz, are also considered.

Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze \nustar hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ($N_{\rm H}>10^{24}$ cm$^{-2}$) AGN to their bolometric luminosity. The X-ray flux upper limits are presently an order of magnitude below what XDR-driven chemistry model predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only $\sim3\%$ of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO$^+$ (for several $\mathrm{J}_\mathrm{upper}$ levels) and $L_\mathrm{AGN}/L_\mathrm{IR}$ for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO$^+$ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of whether SMBH accretion is ongoing.

We show the first interferometric maps of the 14N/15N ratio obtained with the Atacama Large Millimeter Array (ALMA) towards the Solar-like forming protocluster OMC-2 FIR4. We observed N2H+, 15NNH+, N15NH+ (1-0), and N2D+ (2-1), from which we derive the isotopic ratios 14N/15N and D/H. The target is one of the closest analogues of the environment in which our Sun may have formed. The ALMA images, having synthesised beam corresponding to ~600 au, show that the emission of the less abundant isotopologues is distributed in several cores of ~10" (i.e. ~0.02~pc or 4000 au) embedded in a more extended N2H+ emission. Overall, our results indicate that: (1) 14N/15N does not change across the region at core scales, and (2) 14N/15N does not depend on temperature variations. Our findings also suggest that the 14N/15N variations found in pristine Solar System objects are likely not inherited from the protocluster stage, and hence their reason has to be found elsewhere.

The SOLIS (Seeds Of Life In Space) IRAM/NOEMA Large Program aims at studying a set of crucial complex organic molecules in a sample of sources, with well-known physical structure, covering the various phases of Solar-type star formation. One representative object of the transition from the prestellar core to the protostar phases has been observed toward the Very Low Luminosity Object (VeLLO) called L1521F. This type of source is important to study to make the link between prestellar cores and Class 0 sources and also to constrain the chemical evolution during the process of star formation. Two frequency windows (81.6-82.6 GHz and 96.65-97.65 GHz) were used to observe the emission from several complex organics toward the L1521F VeLLO. Only 2 transitions of methanol (A+, E2) have been detected in the narrow window centered at 96.7 GHz (with an upper limit on E1) in a very compact emission blob (~7'' corresponding to ~1000au) toward the NE of the L1521F protostar. The CS 2-1 transition is also detected within the WideX bandwidth. Consistently, with what has been found in prestellar cores, the methanol emission appears ~1000au away from the dust peak. The location of the methanol blob coincides with one of the filaments previously reported in the literature. The Tex of the gas inferred from methanol is (10$\pm$2) K, while the H2 gas density (estimated from the detected CS 2-1 emission and previous CS 5-4 ALMA obs.) is a factor >25 higher than the density in the surrounding environment (n(H2) >10$^{7}$ cm$^{-3}$). From its compactness, low excitation temperature and high gas density, we suggest that the methanol emission detected with NOEMA is either a cold and dense shock-induced blob, recently formed ($\leq$ few hundred years) by infalling gas or a cold and dense fragment that may have just been formed as a result of the intense gas dynamics found within the L1521F VeLLO system.

Aims: The physical structure of a shock wave may take a form unique to its shock type, implying that the chemistry of each shock type is unique as well. We aim to investigate the different chemistries of J-type and C-type shocks in order to identify unique molecular tracers of both shock types. We apply these diagnostics to the protostellar outflow L1157 to establish whether the B2 clump could host shocks exhibiting type-specific behaviour. Of particular interest is the L1157-B2 clump, which has been shown to exhibit bright emission in S-bearing species and HNCO. Methods: We simulate, using a parameterised approach, a planar, steady-state J-type shock wave using UCLCHEM. We compute a grid of models using both C-type and J-type shock models to determine the chemical abundance of shock-tracing species as a function of distance through the shock and apply it to the L1157 outflow. We focus on known shock-tracing molecules such as H2O, HCN, and CH3 OH. Results: We find that a range of molecules including H2O and HCN have unique behaviour specific to a J-type shock, but that such differences in behaviour are only evident at low $v_{s}$ and low $n_{H}$. We find that CH3OH is enhanced by shocks and is a reliable probe of the pre-shock gas density. However, we find no difference between its gas-phase abundance in C-type and J-type shocks. Finally, from our application to L1157, we find that the fractional abundances within the B2 region are consistent with both C-type and J-type shock emission.

We study the r31=L'CO(3-2)/L'CO(1-0) luminosity line ratio in a sample of nearby (z < 0.05) galaxies: 25 star-forming galaxies (SFGs) from the xCOLD GASS survey, 36 hard X-ray selected AGN host galaxies from BASS and 37 infrared luminous galaxies from SLUGS. We find a trend for r31 to increase with star-formation efficiency (SFE). We model r31 using the UCL-PDR code and find that the gas density is the main parameter responsible for variation of r31, while the interstellar radiation field and cosmic ray ionization rate play only a minor role. We interpret these results to indicate a relation between SFE and gas density. We do not find a difference in the r31 value of SFGs and AGN host galaxies, when the galaxies are matched in SSFR (<r31>= 0.52 +/- 0.04 for SFGs and <r31> = 0.53 +/- 0.06 for AGN hosts). According to the results of UCL-PDR models, the X-rays can contribute to the enhancement of the CO line ratio, but only for strong X-ray fluxes and for high gas density (nH > 10$^4$ cm-3). We find a mild tightening of the Kennicutt-Schmidt relation when we use the molecular gas mass surface density traced by CO(3-2) (Pearson correlation coefficient R=0.83), instead of the molecular gas mass surface density traced by CO(1-0) (R=0.78), but the increase in correlation is not statistically significant (p-value=0.06). This suggests that the CO(3-2) line can be reliably used to study the relation between SFR and molecular gas for normal SFGs at high redshift, and to compare it with studies of low-redshift galaxies, as is common practice.

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

Magneto-hydrodynamic (MHD) shocks are violent events that inject large amounts of energy in the interstellar medium (ISM) dramatically modifying its physical properties and chemical composition. Indirect evidence for the presence of such shocks has been reported from the especial chemistry detected toward a variety of astrophysical shocked environments. However, the internal physical structure of these shocks remains unresolved since their expected spatial scales are too small to be measured with current instrumentation. Here we report the first detection of a fully spatially resolved, MHD shock toward the Infrared Dark Cloud (IRDC) G034.77-00.55. The shock, probed by Silicon Monoxide (SiO) and observed with the Atacama Large Millimetre/sub-millimetre Array (ALMA), is associated with the collision between the dense molecular gas of the cloud and a molecular gas flow pushed toward the IRDC by the nearby supernova remnant (SNR) W44. The interaction is occurring on sub-parsec spatial scales thanks to the enhanced magnetic field of the SNR, making the dissipation region of the MHD shock large enough to be resolved with ALMA. Our observations suggest that molecular flow-flow collisions can be triggered by stellar feedback, inducing shocked molecular gas densities compatible with those required for massive star formation.

In studies of the interstellar medium in galaxies, radiative transfer models of molecular emission are useful for relating molecular line observations back to the physical conditions of the gas they trace. However, doing this requires solving a highly degenerate inverse problem. In order to alleviate these degeneracies, the abundances derived from astrochemical models can be converted into column densities and fed into radiative transfer models. This enforces that the molecular gas composition used by the radiative transfer models be chemically realistic. However, because of the complexity and long running time of astrochemical models, it can be difficult to incorporate chemical models into the radiative transfer framework. In this paper, we introduce a statistical emulator of the UCLCHEM astrochemical model, built using neural networks. We then illustrate, through examples of parameter estimations, how such an emulator can be applied on real and synthetic observations.

We investigate differences in the molecular abundances between magnetically super- and sub-critical prestellar cores, performing three-dimensional non-ideal magnetohydrodynamical (MHD) simulations with varying densities and magnetic field strengths, and post-processing the results with a time-dependent gas-grain chemical code. Most molecular species show significantly more central depletion in subcritical models, due to the longer duration of collapse. However, the directly observable quantities - the molecule to hydrogen column density ratios - are generally too similar for observational data to discriminate between models. The profiles of N$_2$H$^+$ and HCO$^+$ show qualitative differences between supercritical and subcritical models on scales of $0.01 \, {\rm pc}$, which may allow the two cases to be distinguished. However, this requires knowledge of the hydrogen column density, which is not directly measureable, and predicted line intensity profiles from radiative transfer modelling are similar for these molecules. Other commonly observed species, such as HCN and CH$_3$OH, have line intensity profiles which differ more strongly between models, and so are more promising as tracers of the mechanism of cloud collapse.

We derive molecular-gas-phase $^{12}$C/$^{13}$C isotope ratios for the central few 100 pc of the three nearby starburst galaxies NGC 253, NGC 1068, and NGC 4945 making use of the $\lambda$ $\sim$ 3 mm $^{12}$CN and $^{13}$CN $N$ = 1--0 lines in the ALMA Band 3. The $^{12}$C/$^{13}$C isotopic ratios derived from the ratios of these lines range from 30 to 67 with an average of 41.6 $\pm$ 0.2 in NGC 253, from 24 to 62 with an average of 38.3 $\pm$ 0.4 in NGC 1068, and from 6 to 44 with an average of 16.9 $\pm$ 0.3 in NGC 4945. The highest $^{12}$C/$^{13}$C isotopic ratios are determined in some of the outskirts of the nuclear regions of the three starburst galaxies. The lowest ratios are associated with the northeastern and southwestern molecular peaks of NGC 253, the northeastern and southwestern edge of the mapped region in NGC 1068, and the very center of NGC 4945. In case of NGC 1068, the measured ratios suggest inflow from the outer part of NGC 1068 into the circum-nuclear disk through both the halo and the bar. Low $^{12}$C/$^{13}$C isotopic ratios in the central regions of these starburst galaxies indicate the presence of highly processed material.

High-resolution (0."03 - 0."09) (9 - 26 pc)) ALMA (100 - 350 GHz (3 - 0.8 mm)) and VLA 45 GHz measurements have been used to image continuum and spectral line emission from the inner region of the nearby infrared luminous galaxy IC 860. We detect compact (r=8-10 pc), luminous, mm continuum emission in the core of IC 860, with brightness temperatures $T_{\rm B}>$ 280 K. The 45 GHz continuum is equally compact, but fainter in flux. We suggest that the 3 to 0.8 mm continuum is opaque at mm-wavelengths, implying a very large H$_2$ column density $N$(H$_2$) of 1e26 cm-2, and that it emerges from hot dust with temperature $T_{\rm d} \sim 280$ K. Vibrationally excited lines of HCN $\nu_2$=1f J=4-3 and 3-2 (HCN-VIB) are seen in emission and resolved on scales of 40-50 pc. The HCN-VIB emission reveals a north-south nuclear velocity gradient with projected rotation velocities of v=100 kms$^{-1}$ at r=10 pc. The brightest HCN-VIB emission is oriented perpendicular to the velocity gradient. Ground state lines of HCN 3-2, 4-3, HC$^{15}$N 4-3, HCO$^+$ 3-2, 4-3 and CS 7-6, show complex line absorption and emission features. HCN and HCO$^+$ have red-shifted, reversed P-Cygni profiles consistent with gas inflows of $\sim$50 kms$^{-1}$. The absorption features can be traced from the north-east into the nucleus. In contrast, CS 7-6 shows blue-shifted line wings extending to -180 kms$^{-1}$. We suggest that a dense and slow outflow is hidden behind a foreground layer of inflowing gas. It appears that the centre of IC 860 is in a phase of rapid evolution where an inflow is building up the nuclear column density of gas. A slow, dense outflow may be signalling the onset of feedback. The inner, r=10 pc, IR luminosity can be powered by an accreting black hole and/or a compact starburst with a top-heavy initial mass function.

The main carrier of sulfur in dense clouds, where it is depleted from the gas phase, remains a mystery. Shock waves in young molecular outflows disrupt the ice mantles and allow us to directly probe the material that is ejected into the gas phase. A comprehensive study of sulfur-bearing species towards L1157-B1, a shocked region along a protostellar outflow, has been carried out as part of the IRAM-30m large program ASAI. The dataset contains over 100 lines of CCS, H$_2$CS, OCS, SO, SO$_2$ and isotopologues. The results of these observations are presented, complementing previous studies of sulfur-bearing species in the region. The column densities and fractional abundances of these species are measured and together these species account for 10\% of the cosmic sulfur abundance in the region. The gas properties derived from the observations are also presented, demonstrating that sulfur bearing species trace a wide range of different gas conditions in the region.

We report ALMA observations of CO(3-2) emission in the Seyfert galaxy NGC613, at a spatial resolution of 17pc, as part of our NUclei of GAlaxies sample. Our aim is to investigate the morphology and dynamics of the gas inside the central kpc, and to probe nuclear fueling and feedback phenomena. The morphology of CO(3-2) line emission reveals a 2-arm trailing nuclear spiral at $r\lesssim$100pc and a circumnuclear ring at ~350pc radius, that is coincident with the star-forming ring seen in the optical images. The molecular gas in the galaxy disk is in a remarkably regular rotation, however, the kinematics in the nuclear region is very skewed. The nuclear spectrum of CO and dense gas tracers HCN(4-3), HCO+(4-3), and CS(7-6) show broad wings up to $\pm$300km/s, associated with a molecular outflow emanating from the nucleus (r~25pc). We derive a molecular outflow mass $M_{out}$=2x10$^6$M$_\odot$ and a mass outflow rate of $\dot{M}_{out}=$27$\rm M_\odot yr^{-1}$. The molecular outflow energetics exceed the values predicted by AGN feedback models: the kinetic power of the outflow corresponds to 20%L_AGN and the momentum rate is $\dot{M}_{out}v\sim400L_{AGN}/c$. The outflow is mainly boosted by the AGN through entrainment by the radio jet, but given the weak nuclear activity of NGC613, we might be witnessing a fossil outflow, resulted from a strong past AGN that now has already faded. The nuclear trailing spiral observed in CO emission is inside the ILR ring of the bar. We compute the gravitational torques exerted in the gas to estimate the efficiency of the angular momentum exchange. The gravity torques are negative from 25 to 100pc and the gas loses its angular momentum in a rotation period, providing evidence of a highly efficient inflow towards the center. This phenomenon shows that the massive central black hole has a significant dynamical influence on the gas, triggering its fueling.

Many species of complex organic molecules (COMs) have been observed in several astrophysical environments but it is not clear how they are produced, particularly in cold, quiescent regions. One process that has been proposed as a means to enhance the chemical complexity of the gas phase in such regions is the explosion of the ice mantles of dust grains. In this process, a build up of chemical energy in the ice is released, sublimating the ices and producing a short lived phase of high density, high temperature gas. The gas-grain chemical code UCLCHEM has been modified to treat these explosions in order to model the observed abundances of COMs towards the TMC-1 region. It is found that, based on our current understanding of the explosion mechanism and chemical pathways, the inclusion of explosions in chemical models is not warranted at this time. Explosions are not shown to improve the model's match to the observed abundances of simple species in TMC-1. Further, neither the inclusion of surface diffusion chemistry, nor explosions, results in the production of COMs with observationally inferred abundances.

In star forming regions in our own Galaxy, the 14N/15N ratio is found to vary from $\sim$ 100 in meteorites, comets and protoplanetary disks up to $\sim$ 1000 in pre-stellar and star forming cores, while in external galaxies the very few single-dish large scale measurements of this ratio lead to values of 100-450. The extent of the contribution of isotopic fractionation to these variations is, to date, unknown. In this paper we present a theoretical chemical study of nitrogen fractionation in external galaxies in order to determine the physical conditions that may lead to a spread of the 14N/15N ratio from the solar value of $\sim$440 and hence evaluate the contribution of chemical reactions in the ISM to nitrogen fractionation. We find that the main cause of ISM enrichment of nitrogen fractionation is high gas densities, aided by high fluxes of cosmic rays.

Iram 30-m Observations towards eight protostellar outflow sources were taken in the 96-\SI{176}{\giga\hertz} range. Transitions of CH$_3$OH and CH$_3$CHO were detected in seven of them. The integrated emission of the transitions of each species that fell into the observed frequency range were measured and fit using RADEX and LTE models. Column densities and gas properties inferred from this fitting are presented. The ratio of the A and E-type isomers of CH$_3$OH indicate that the methanol observed in these outflows was formed on the grain surface. Both species demonstrate a reduction of terminal velocity in their line profiles in faster outflows, indicating destruction in the post-shock gas phase. This destruction, and a near constant ratio of the CH$_3$OH and CH$_3$CHO column densities imply it is most likely that CH$_3$CHO also forms on the grain surface.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

We report ALMA observations of CO(3-2) emission in a sample of 7 Seyfert/LINER galaxies at the unprecedented spatial resolution of 0.1'' = 4-9 pc. Our aim is to explore the close environment of AGN, and the dynamical structures leading to its fueling. The selected galaxies host low-luminosity AGN, and have a wide range of activity types, and barred or ringed morphologies. The observed maps reveal the existence of circum-nuclear disk structures, defined by their morphology and decoupled kinematics. We call these structures "molecular tori", even though they appear often as disks, without holes in the center. They have varied orientations, unaligned with their host galaxy. The radius of the tori ranges from 6 to 27 pc, and their mass from 0.7 10$^7$ to 3.9 10$^7$ Msun. At larger scale, the gas is always piled up in a few 100~pc scale resonant rings, that play the role of a reservoir to fuel the nucleus. In some cases, a trailing spiral is observed inside the ring, providing evidence for feeding processes. Most frequently, the torus and the AGN are slightly off-centered, with respect to the bar-resonant ring position, implying that the black hole is wandering by a few 10~pc amplitude around the center of mass of the galaxy. Our spatial resolution allows us to measure gas velocities inside the sphere of influence of the central black holes. By fitting the observations with different simulated cubes, varying the torus inclination and the black hole mass, it is possible to estimate the mass of the central black hole, which is in general difficult for such late-type galaxies, with only a pseudo-bulge.

The formyl radical HCO has been proposed as the basic precursor of many complex organic molecules such as methanol (CH$_3$OH) or glycolaldehyde (CH$_2$OHCHO). Using ALMA, we have mapped, for the first time at high angular resolution ($\sim$1$^{\prime\prime}$, $\sim$140 au), HCO towards the Solar-type protostellar binary IRAS 16293$-$2422, where numerous complex organic molecules have been previously detected. We also detected several lines of the chemically related species H$_2$CO, CH$_3$OH and CH$_2$OHCHO. The observations revealed compact HCO emission arising from the two protostars. The line profiles also show redshifted absorption produced by foreground material of the circumbinary envelope that is infalling towards the protostars. Additionally, IRAM 30m single-dish data revealed a more extended HCO component arising from the common circumbinary envelope. The comparison between the observed molecular abundances and our chemical model suggests that whereas the extended HCO from the envelope can be formed via gas-phase reactions during the cold collapse of the natal core, the HCO in the hot corinos surrounding the protostars is predominantly formed by the hydrogenation of CO on the surface of dust grains and subsequent thermal desorption during the protostellar phase. The derived abundance of HCO in the dust grains is high enough to produce efficiently more complex species such as H$_2$CO, CH$_3$OH, and CH$_2$OHCHO by surface chemistry. We found that the main formation route of CH$_2$OHCHO is the reaction between HCO and CH$_2$OH.

Grain surface chemistry and its treatment in gas-grain chemical models is an area of large uncertainty. Whilst laboratory experiments are making progress, there is still much that is unknown about grain surface chemistry. Further, the results and parameters produced by experiment are often not easily translated to the rate equation approach most commonly used in astrochemical modelling. It is possible that statistical methods can reduce the uncertainty in grain surface chemical networks. In this work, a simple model of grain surface chemistry in a molecular cloud is developed and a Bayesian inference of the reactions rates is performed through MCMC sampling. Using observational data of the solid state abundances of major chemical species in molecular clouds, the posterior distributions for the rates of seven reactions producing CO, CO$_2$, CH$_3$OH and H$_2$O are calculated, in a form that is suitable for rate equation models. This represents a vital first step in the development of a method to infer reaction rates from observations of chemical abundances in astrophysical environments.

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

For decades, the detection of phosphorus-bearing molecules in the interstellar medium was restricted to high-mass star-forming regions (as e.g. SgrB2 and Orion KL) and the circumstellar envelopes of evolved stars. However, recent higher-sensitivity observations have revealed that molecules such as PN and PO are present not only toward cold massive cores and low-mass star-forming regions with PO/PN ratios >1, but also toward the Giant Molecular Clouds in the Galactic Center known to be exposed to highly energetic phenomena such as intense UV radiation fields, shock waves and cosmic rays. In this paper, we carry out a comprehensive study of the chemistry of phosphorus-bearing molecules across different astrophysical environments which cover a range of physical conditions (cold molecular dark clouds, warm clouds, hot cores/hot corinos) and are exposed to different physical processes and energetic phenomena (proto-stellar heating, shock waves, intense UV radiation and cosmic rays). We show how the measured PO/PN ratio (either >1 as in e.g. hot molecular cores, or <1 as in UV strongly illuminated environments) can provide constraints on the physical conditions and energetic processing of the source. We propose that the reaction P + OH --> PO + H, not included in previous works, could be an efficient gas-phase PO formation route in shocks. Our modelling provides a template with which to study the detectability of P-bearing species not only in regions in our own Galaxy but also in extragalactic sources.

We develop analytic approximations to the density evolution of prestellar cores, based on the results of hydrodynamical simulations. We use these approximations as input for a time-dependent gas-grain chemical code to investigate the effects of differing modes of collapse on the molecular abundances in the core. We confirm that our method can provide reasonable agreement with an exact numerical solution of both the hydrodynamics and chemistry while being significantly less computationally expensive, allowing a large grid of models varying multiple input parameters to be run. We present results using this method to illustrate how the chemistry is affected not only by the collapse model adopted, but also by the large number of unknown physical and chemical parameters. Models which are initially gravitationally unstable predict similar abundances despite differing densities and collapse timescales, while ambipolar diffusion produces more extended inner depleted regions which are not seen in observations of prestellar cores. Molecular observations are capable of discriminating between modes of collapse despite the unknown values of various input parameters. We also investigate the evolution of the ambipolar diffusion timescale for a range of collapse modes, metallicities and cosmic ray ionization rates, finding that it remains comparable to or larger than the collapse timescale during the initial stages for all models we consider, but becomes smaller at later evolutionary stages. This confirms that ambipolar diffusion is an important process for diffuse gas, but becomes less significant as cores collapse to higher densities.

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.

NGC 4945 is one of the nearest (~3.8 Mpc; 1" ~ 19 pc) starburst galaxies. ALMA band 3 (3--4\,mm) observations of HCN, HCO+, CS, C3H2, SiO, HCO, and CH3C2H were carried out with ~2" resolution. The lines reveal a rotating nuclear disk of projected size 10" x 2" with position angle ~45 deg, inclination ~75 deg and an unresolved bright central core of size <2.5". The continuum source (mostly free-free radiation) is more compact than the nuclear disk by a linear factor of two but shows the same position angle and is centered 0.39" +_ 0.14" northeast of the nuclear accretion disk defined by H2O maser emission. Outside the nuclear disk, both HCN and CS delineate molecular arms on opposite sides of the dynamical center. These are connected by a (deprojected) 0.6 kpc sized molecular bridge, likely a dense gaseous bar seen almost ends-on, shifting gas from the front and back side into the nuclear disk. Modeling this nuclear disk located farther inside <100 pc) with tilted rings indicates a coplanar outflow reaching a characteristic deprojectd velocity of ~50 km/s. All our molecular lines, with the notable exception of CH3C2H, show significant absorption near the systemic velocity (~571 km/s), within a range of ~500-660 km/s. Apparently, only molecular transitions with low critical H2-density do not show absorption. The velocity field of the nuclear disk, derived from CH3C2H, provides evidence for rigid rotation in the inner few arcseconds and a dynamical mass of M = (2.1+_0.2) x 10^8 Mo inside a galactocentric radius of 2.45", with a significantly flattened rotation curve farther out. Velocity integrated line intensity maps with most pronounced absorption show molecular peak positions up to 1.5" southwest of the continuum peak, presumably due to absorption, which appears to be most severe slightly northeast of the nuclear maser disk.

In dense molecular clouds, interstellar grains are covered by mantles of iced molecules. The formation of the grain mantles has two important consequences: it removes species from the gas phase and promotes the synthesis of new molecules on the grain surfaces. The composition of the mantle is a strong function of the environment which the cloud belongs to. Therefore, clouds in high-zeta galaxies, where conditions -like temperature, metallicity and cosmic rays flux- are different from those in the Milky Way, will have different grain mantles. In the last years, several authors have suggested that silicate grains might grow by accretion of silicon bearing species on smaller seeds. This would occur simultaneously to the formation of the iced mantles and be greatly affected by its composition as a function of time. In this work, we present a numerical study of the grain mantle formation in high-zeta galaxies and we quantitatively address the possibility of silicate growth. We find that the mantle thickness decreases with increasing redshift, from about 120 to 20 layers for z varying from 0 to 8. Furthermore, the mantle composition is also a strong function of the cloud redshift, with the relative importance of CO, CO2, ammonia, methane and methanol highly varying with z. Finally, being Si-bearing species always a very minor component of the mantle, the formation of silicates in molecular clouds is practically impossible.

Towards the pre-stellar core L1544, the methanol (CH$_3$OH) emission forms an asymmetric ring around the core centre, where CH$_3$OH is mostly in solid form, with a clear peak 4000~au to the north-east of the dust continuum peak. As part of the NOEMA Large Project SOLIS (Seeds of Life in Space), the CH$_3$OH peak has been spatially resolved to study its kinematics and physical structure and to investigate the cause behind the local enhancement. We find that methanol emission is distributed in a ridge parallel to the main axis of the dense core. The centroid velocity increases by about 0.2~km~s$^{-1}$ and the velocity dispersion increases from subsonic to transonic towards the central zone of the core, where the velocity field also shows complex structure. This could be indication of gentle accretion of material onto the core or interaction of two filaments, producing a slow shock. We measure the rotational temperature and show that methanol is in local thermodynamic equilibrium (LTE) only close to the dust peak, where it is significantly depleted. The CH$_3$OH column density, $N_{tot}({\rm CH_3OH})$, profile has been derived with non-LTE radiative transfer modelling and compared with chemical models of a static core. The measured $N_{tot}({\rm CH_3OH})$ profile is consistent with model predictions, but the total column densities are one order of magnitude lower than those predicted by models, suggesting that the efficiency of reactive desorption or atomic hydrogen tunnelling adopted in the model may be overestimated; or that an evolutionary model is needed to better reproduce methanol abundance.

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.

We present here a systematic search for cyanopolyynes in the shock region L1157-B1 and its associated protostar L1157-mm in the framework of the Large Program "Astrochemical Surveys At IRAM" (ASAI), dedicated to chemical surveys of solar-type star forming regions with the IRAM 30m telescope. Observations of the millimeter windows between 72 and 272 GHz permitted the detection of HC$_3$N and its $^{13}$C isotopologues, and HC$_5$N (for the first time in a protostellar shock region). In the shock, analysis of the line profiles shows that the emission arises from the outflow cavities associated with L1157-B1 and L1157-B2. Molecular abundances and excitation conditions were obtained from analysis of the Spectral Line Energy Distributions under the assumption of Local Thermodynamical Equilibrium or using a radiative transfer code in the Large Velocity Gradient approximation. Towards L1157mm, the HC$_3$N emission arises from the cold envelope ($T_{rot}=10$ K) and a higher-excitation region ($T_{rot}$= $31$ K) of smaller extent around the protostar. We did not find any evidence of $^{13}$C or D fractionation enrichment towards L1157-B1. We obtain a relative abundance ratio HC$_3$N/HC$_5$N of 3.3 in the shocked gas. We find an increase by a factor of 30 of the HC$_3$N abundance between the envelope of L1157-mm and the shock region itself. Altogether, these results are consistent with a scenario in which the bulk of HC$_3$N was produced by means of gas phase reactions in the passage of the shock. This scenario is supported by the predictions of a parametric shock code coupled with the chemical model UCL_CHEM.

Glycolaldehyde (HOCH$_2$CHO) and ethylene glycol ((CH$_2$OH)$_2$) are two complex organic molecules detected in the hot cores and hot corinos of several star-forming regions. The ethylene glycol/glycolaldehyde abundance ratio seems to show an increase with the source luminosity. In the literature, several surface-chemistry formation mechanisms have been proposed for these two species. With the UCLCHEM chemical code, we explored the different scenarios and compared the predictions for a range of sources of different luminosities with the observations. None of the scenarios reproduce perfectly the trend. A better agreement is, however, found for a formation through recombination of two HCO radicals followed by successive hydrogenations. The reaction between HCO and CH$_2$OH could also contribute to the formation of glycolaldehyde in addition to the hydrogenation pathway. The predictions are improved when a trend of decreasing H$_2$ density within the core region with T $\geq$ 100 K as a function of luminosity, is included in the model. Destruction reactions of complex organic molecules in the gas phase would also need to be investigated, since they can affect the abundance ratios once the species have desorbed in the warm inner regions of the star-forming regions.

Phosphorus (P) is one of the essential elements for life due to its central role in biochemical processes. Recent searches have shown that P-bearing molecules (in particular PN and PO) are present in star-forming regions, although their formation routes remain poorly understood. In this Letter, we report observations of PN and PO towards seven molecular clouds located in the Galactic Center, which are characterized by different types of chemistry. PN is detected in five out of seven sources, whose chemistry is thought to be shock-dominated. The two sources with PN non-detections correspond to clouds exposed to intense UV/X-rays/cosmic-ray radiation. PO is detected only towards the cloud G+0.693$-$0.03, with a PO/PN abundance ratio of $\sim$1.5. We conclude that P-bearing molecules likely form in shocked gas as a result of dust grain sputtering, while are destroyed by intense UV/X-ray/cosmic ray radiation.

The high-sensitivity of the IRAM 30-m ASAI unbiased spectral survey in the mm-window allows us to detect NO emission towards both the Class I object SVS13-A and the protostellar outflow shock L1157-B1. We detect the hyperfine components of the $^2\Pi_{\rm 1/2}$ $J$ = 3/2 $\to$ 1/2 (at 151 GHz) and the $^2\Pi_{\rm 1/2}$ $J$ = 5/2 $\to$ 3/2 (250 GHz) spectral pattern. The two objects show different NO profiles: (i) SVS13-A emits through narrow (1.5 km s$^{-1}$) lines at the systemic velocity, while (ii) L1157-B1 shows broad ($\sim$ 5 km s$^{-1}$) blue-shifted emission. For SVS13-A the analysis leads to $T_{\rm ex}$ $\geq$ 4 K, $N(\rm NO)$ $\leq$ 3 $\times$ 10$^{15}$ cm$^{-2}$, and indicates the association of NO with the protostellar envelope. In L1157-B1, NO is tracing the extended outflow cavity: $T_{\rm ex}$ $\simeq$ 4--5 K, and $N(\rm NO)$ = 5.5$\pm$1.5 $\times$ 10$^{15}$ cm$^{-2}$. Using C$^{18}$O, $^{13}$C$^{18}$O, C$^{17}$O, and $^{13}$C$^{17}$O ASAI observations we derive an NO fractional abundance less than $\sim$ 10$^{-7}$ for the SVS13-A envelope, in agreement with previous measurements towards extended PDRs and prestellar objects. Conversely, a definite $X(NO)$ enhancement is measured towards L1157-B1, $\sim$ 6 $\times$ 10$^{-6}$, showing that the NO production increases in shocks. The public code UCLCHEM was used to interpret the NO observations, confirming that the abundance observed in SVS13-A can be attained in an envelope with a gas density of 10$^5$ cm$^{-3}$ and a kinetic temperature of 40 K. The NO abundance in L1157-B1 is reproduced with pre-shock densities of 10$^5$ cm$^{-3}$ subjected to a $\sim$ 45 km s$^{-1}$ shock.

Infrared-Dark Clouds (IRDCs) are cold, dense regions of high (optical and infrared) extinction, believed to be the birthplace of high-mass stars and stellar clusters. The physical mechanisms leading to the formation of these IRDCs are not completely understood and it is thus important to study their molecular gas kinematics and chemical content to search for any signature of the IRDCs formation process. Using the 30m-diameter antenna at the Instituto de Radioastronom\'ia Milim\'etrica, we have obtained emission maps of dense gas tracers (H$^{13}$CO$^{+}$ and HN$^{13}$C) and typical shock tracers (SiO and CH$_3$OH) toward three IRDCs, G028.37+00.07, G034.43+00.24 and G034.77-00.55 (clouds C, F and G, respectively). We have studied the molecular gas kinematics in these clouds and, consistent with previous works toward other IRDCs, the clouds show complex gas kinematics with several velocity-coherent sub-structures separated in velocity space by a few km s$^{-1}$. Correlated with these complex kinematic structures, widespread (parsec-scale) emission of SiO and CH$_3$OH is present in all the three clouds. For clouds C and F, known to be actively forming stars, widespread SiO and CH$_3$OH is likely associated with on-going star formation activity. However, for cloud G, which lacks either 8 $\mu$m or 24 $\mu$m sources and 4.5 $\mu$m H$_2$ shock-excited emission, the detected widespread SiO and CH$_3$OH emission may have originated in a large-scale shock interaction, although a scenario involving a population of low-mass stars driving molecular outflows cannot be fully ruled out.

Peptide bonds (N-C=O) play a key role in metabolic processes since they link amino acids into peptide chains or proteins. Recently, several molecules containing peptide-like bonds have been detected across multiple environments in the interstellar medium (ISM), growing the need to fully understand their chemistry and their role in forming larger pre-biotic molecules. We present a comprehensive study of the chemistry of three molecules containing peptide-like bonds: HNCO, NH$_2$CHO, and CH$_3$NCO. We also included other CHNO isomers (HCNO, HOCN), and C$_2$H$_3$NO isomers (CH$_3$OCN, CH$_3$CNO) to the study. We have used the \uclchem gas-grain chemical code and included in our chemical network all possible formation/destruction pathways of these peptide-like molecules recently investigated either by theoretical calculations or in laboratory experiments. Our predictions are compared to observations obtained toward the proto-star IRAS16293$-$2422 and the L1544 pre-stellar core. Our results show that some key reactions involving the CHNO and C$_2$H$_3$NO isomers need to be modified to match the observations. Consistently with recent laboratory findings, hydrogenation is unlikely to produce NH$_2$CHO on grain surfaces, while a combination of radical-radical surface reactions and gas-phase reactions is a better alternative. In addition, better results are obtained for NH$_2$CHO when a slightly higher activation energy of 25$\,$K is considered for the gas-phase reaction $\rm NH_2 + H_2CO \rightarrow NH_2CHO + H$. Finally, our modelling shows that the observed correlation between NH$_2$CHO and HNCO in star-forming regions may come from the fact that HNCO and NH$_2$CHO react to temperature in the same manner rather than from a direct chemical link between the two species.

Complex organic molecules have been observed for decades in the interstellar medium. Some of them might be considered as small bricks of the macromolecules at the base of terrestrial life. It is hence particularly important to understand organic chemistry in Solar-like star forming regions. In this article, we present a new observational project: SOLIS (Seeds Of Life In Space). This is a Large Project at the IRAM-NOEMA interferometer, and its scope is to image the emission of several crucial organic molecules in a sample of Solar-like star forming regions in different evolutionary stage and environments. Here, we report the first SOLIS results, obtained from analysing the spectra of different regions of the Class 0 source NGC1333-IRAS4A, the protocluster OMC-2 FIR4, and the shock site L1157-B1. The different regions were identified based on the images of formamide (NH2CHO) and cyanodiacetylene (HC5N) lines. We discuss the observed large diversity in the molecular and organic content, both on large (3000-10000 au) and relatively small (300-1000 au) scales. Finally, we derive upper limits to the methoxy fractional abundance in the three observed regions of the same order of magnitude of that measured in few cold prestellar objects, namely ~10^-12-10^-11 with respect to H2 molecules.

Ab initio grids of time dependent chemical models, varying in gas density, temperature, cosmic ray ionization rate, and radiation field, are used as input to RADEX calculations. Tables of abundances, column densities, theoretical line intensities, and line ratios for some of the most used dense gas tracers are provided. The degree of correlation as well as degeneracy inherent in molecular ratios is discussed. Comparisons of the theoretical intensities with example observations are also provided. We find that, within the parameters space explored, chemical abundances can be constrained by a well defined set of gas density-gas temperature-cosmic ray ionization rate for the species we investigate here. However, line intensities, as well as, more importantly, line ratios, from different chemical models can be very similar leading to a clear degeneracy. We also find that the gas subjected to a galactic cosmic ray ionization rate will not necessarily have reached steady state by 1 Myr. The species most affected by time dependency effects are HCN and CS, both high density tracers. We use our method to fit an example set of data from two galaxies. We find that (i) molecular line ratios can be easily matched even with erroneous individual line intensities; (ii) no set of species can be matched by a one-component ISM; (iii) a species may be a good tracer of an energetic process but only under specific density and temperature conditions. We show that by taking into consideration the chemistry behind each species and the individual line intensities, many degeneracies that arise by just using molecular line ratios can be avoided. Finally we show that using a species or a ratio as a tracer of an individual energetic process (e.g. cosmic rays, UV) ought to be done with caution.

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.

We have performed combined photoionization and photodissociation region (PDR) modelling of a Crab Nebula filament subjected to the synchrotron radiation from the central pulsar wind nebula, and to a high flux of charged particles; a greatly enhanced cosmic ray ionization rate over the standard interstellar value, $\zeta_0$, is required to account for the lack of detected [C I] emission in published Herschel SPIRE FTS observations of the Crab Nebula. The observed line surface brightness ratios of the OH$^+$ and ArH$^+$ transitions seen in the SPIRE FTS frequency range can only be explained with both a high cosmic ray ionization rate and a reduced ArH$^+$ dissociative recombination rate compared to that used by previous authors, although consistent with experimental upper limits. We find that the ArH$^+$/OH$^+$ line strengths and the observed H$_2$ vibration-rotation emission can be reproduced by model filaments with $n_{\rm{H}} = 2 \times 10^4$ cm$^{-3}$, $\zeta = 10^7 \zeta_0$ and visual extinctions within the range found for dusty globules in the Crab Nebula, although far-infrared emission from [O I] and [C II] is higher than the observational constraints. Models with $n_{\rm{H}} = 1900$ cm$^{-3}$ underpredict the H$_2$ surface brightness, but agree with the ArH$^+$ and OH$^+$ surface brightnesses and predict [O I] and [C II] line ratios consistent with observations. These models predict HeH$^+$ rotational emission above detection thresholds, but consideration of the formation timescale suggests that the abundance of this molecule in the Crab Nebula should be lower than the equilibrium values obtained in our analysis.

Context: Modern versions of the Miller-Urey experiment claim that formamide (NH$_2$CHO) could be the starting point for the formation of metabolic and genetic macromolecules. Intriguingly, formamide is indeed observed in regions forming Solar-type stars as well as in external galaxies. Aims: How NH$_2$CHO is formed has been a puzzle for decades: our goal is to contribute to the hotly debated question of whether formamide is mostly formed via gas-phase or grain surface chemistry. Methods: We used the NOEMA interferometer to image NH$_2$CHO towards the L1157-B1 blue-shifted shock, a well known interstellar laboratory, to study how the components of dust mantles and cores released into the gas phase triggers the formation of formamide. Results: We report the first spatially resolved image (size $\sim$ 9", $\sim$ 2300 AU) of formamide emission in a shocked region around a Sun-like protostar: the line profiles are blueshifted and have a FWHM $\simeq$ 5 km s$^{-1}$. A column density of $N_{\rm NH_2CHO}$ = 8 $\times$ 10$^{12}$ cm$^{-1}$, and an abundance (with respect to H-nuclei) of 4 $\times$ 10$^{-9}$ are derived. We show a spatial segregation of formamide with respect to other organic species. Our observations, coupled with a chemical modelling analysis, indicate that the formamide observed in L1157-B1 is formed by gas-phase chemical process, and not on grain surfaces as previously suggested. Conclusions: The SOLIS interferometric observations of formamide provide direct evidence that this potentially crucial brick of life is efficiently formed in the gas-phase around Sun-like protostars.

L1157-B1 is the brightest shocked region of the large-scale molecular outflow, considered the prototype of chemically rich outflows, being the ideal laboratory to study how shocks affect the molecular gas. Several deuterated molecules have been previously detected with the IRAM 30m, most of them formed on grain mantles and then released into the gas phase due to the shock. We aim to observationally investigate the role of the different chemical processes at work that lead to formation the of DCN and test the predictions of the chemical models for its formation. We performed high-angular resolution observations with NOEMA of the DCN(2-1) and H13CN(2-1) lines to compute the deuterated fraction, Dfrac(HCN). We detected emission of DCN(2-1) and H13CN(2-1) arising from L1157-B1 shock. Dfrac(HCN) is ~4x10$^{-3}$ and given the uncertainties, we did not find significant variations across the bow-shock. Contrary to HDCO, whose emission delineates the region of impact between the jet and the ambient material, DCN is more widespread and not limited to the impact region. This is consistent with the idea that gas-phase chemistry is playing a major role in the deuteration of HCN in the head of the bow-shock, where HDCO is undetected as it is a product of grain-surface chemistry. The spectra of DCN and H13CN match the spectral signature of the outflow cavity walls, suggesting that their emission result from shocked gas. The analysis of the time dependent gas-grain chemical model UCL-CHEM coupled with a C-type shock model shows that the observed Dfrac(HCN) is reached during the post-shock phase, matching the dynamical timescale of the shock. Our results indicate that the presence of DCN in L1157-B1 is a combination of gas-phase chemistry that produces the widespread DCN emission, dominating in the head of the bow-shock, and sputtering from grain mantles toward the jet impact region.

The interstellar delivery of carbon atoms locked into molecules might be one of the key ingredients for the emergence of life. Cyanopolyynes are carbon chains delimited at their two extremities by an atom of hydrogen and a cyano group, so that they might be excellent reservoirs of carbon. The simplest member, HC3N, is ubiquitous in the galactic interstellar medium and found also in external galaxies. Thus, understanding the growth of cyanopolyynes in regions forming stars similar to our Sun, and what affects it, is particularly relevant. In the framework of the IRAM/NOEMA Large Program SOLIS (Seeds Of Life In Space), we have obtained a map of two cyanopolyynes, HC3N and HC5N, in the protocluster OMC2-FIR4. Because our Sun is thought to be born in a rich cluster, OMC2-FIR4 is one of the closest and best known representatives of the environment in which the Sun may have been born. We find a HC3N/HC5N abundance ratio across the source in the range ~ 1 - 30, with the smallest values (< 10) in FIR5 and in the Eastern region of FIR4. The ratios < 10 can be reproduced by chemical models only if: (1) the cosmic-ray ionisation rate $\zeta$ is ~ $4 \times 10^{-14}$ s$^{-1}$; (2) the gaseous elemental ratio C/O is close to unity; (3) O and C are largely depleted. The large $\zeta$ is comparable to that measured in FIR4 by previous works and was interpreted as due to a flux of energetic (> 10 MeV) particles from embedded sources. We suggest that these sources could lie East of FIR4 and FIR5. A temperature gradient across FIR4, with T decreasing by about 10 K, could also explain the observed change in the HC3N/HC5N line ratio, without the need of a cosmic ray ionisation rate gradient. However, even in this case, a high constant cosmic-ray ionisation rate (of the order of $10^{-14}$ s$^{-1}$) is necessary to reproduce the observations.

We present a publicly available, open source version of the time dependent gas-grain chemical code UCLCHEM. UCLCHEM propagates the abundances of chemical species through a large network of chemical reactions in a variety of physical conditions. The model is described in detail along with its applications. As an example of possible uses, UCLCHEM is used to explore the effect of protostellar collapse on commonly observed molecules and to study the behaviour of molecules in C-type shocks. We find the collapse of a simple Bonnor-Ebert sphere successfully reproduces most of the behaviour of CO,CS and NH$_3$ from cores observed by \citet{Tafalla2004} but cannot predict the behaviour of N$_2$H$^+$. In the C-shock application, we find that molecules can be categorized so that they can be useful observational tracers of shocks and their physical properties. Whilst many molecules are enhanced in shocked gas, we identify two groups of molecules in particular. A small number of molecules are enhanced by the sputtering of the ices as the shock propagates and then remain high in abundance throughout the shock. A second, larger set are also enhanced by sputtering but then are destroyed as the gas temperature rises. Through these applications the general applicability of UCLCHEM is demonstrated.

The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in particular those of ESO (i.e., the VLT and ELT, ALMA and APEX), by producing a White Paper to identify and develop the: 1. needs to access ESO ground-based facilities to achieve the formulated Athena science objectives; 2. needs to access Athena to achieve the formulated science objectives of ESO facilities contemporary to Athena; 3. science areas where the synergetic use of Athena and ESO facilities in the late 2020s will result in scientific added value. Community input to the process happened primarily via a dedicated ESO - Athena Synergy Workshop that took place on Sept. 14 - 16, 2016 at ESO, Garching. This White Paper presents the results of the EAST's work, sorted by synergy area, and deals with the following topics: 1. the Hot Universe: Early groups and clusters and their evolution, Physics of the Intracluster medium, Missing baryons in cosmic filaments; 2. the Energetic Universe: Supermassive black hole (SMBH) history, SMBH accretion disks, Active Galactic Nuclei feedback - Molecular outflows, Ultra-fast outflows, Accretion Physics, Transient Science; 3. Observatory Science: Star Formation, Stars. It then discusses the optical-NIR-sub-mm perspective by providing details on VLT/MOONS, the E-ELT instruments, in particular the MOS, VISTA/4MOST, the ESO and ALMA archives, future ALMA and ESO developments, and finally the (likely) ESO - Athena astronomical scene in the 2020s. (abridged)

Nitrogen is one of the most abundant elements in the Universe and its 14N/15N isotopic ratio has the potential to provide information about the initial environment in which our Sun formed. Recent findings suggest that the Solar System may have formed in a massive cluster since the presence of short-lived radioisotopes in meteorites can only be explained by the influence of a supernova. The aim of this project is to determine the 14N/15N ratio towards a sample of cold, massive dense cores at the initial stages in their evolution. We have observed the J=1-0 transitions of HCN, H13CN, HC15N, HN13C and H15NC toward a sample of 22 cores in 4 Infrared-Dark Clouds (IRDCs). IRDCs are believed to be the precursors of high-mass stars and star clusters. Assuming LTE and a temperature of 15K, the column densities of HCN, H13CN, HC15N, HN13C and H15NC are calculated and their 14N/15N ratio is determined for each core. The 14N/15N ratio measured in our sample of IRDC cores range between ~70 and >763 in HCN and between ~161 and ~541 in HNC. They are consistent with the terrestrial atmosphere (TA) and protosolar nebula (PSN) values, and with the ratios measured in low-mass pre-stellar cores. However, the 14N/15N ratios measured in cores C1, C3, F1, F2 and G2 do not agree with the results from similar studies toward the same massive cores using nitrogen bearing molecules with nitrile functional group (-CN) and nitrogen hydrides (-NH) although the ratio spread covers a similar range. Amongst the 4 IRDCs we measured relatively low 14N/15N ratios towards IRDC G which are comparable to those measured in small cosmomaterials and protoplanetary disks. The low average gas density of this cloud suggests that the gas density, rather than the gas temperature, may be the dominant parameter influencing the initial nitrogen isotopic composition in young PSN.

We report on a systematic search for oxygen-bearing Complex Organic Molecules (COMs) in the Solar-like protostellar shock region L1157-B1, as part of the IRAM Large Program "Astrochemical Surveys At IRAM" (ASAI). Several COMs are unambiguously detected, some for the first time, such as ketene H$_2$CCO, dimethyl ether (CH$_3$OCH$_3$) and glycolaldehyde (HCOCH$_2$OH), and others firmly confirmed, such as formic acid (HCOOH) and ethanol (C$_2$H$_5$OH). Thanks to the high sensitivity of the observations and full coverage of the 1, 2 and 3mm wavelength bands, we detected numerous (10--125) lines from each of the detected species. Based on a simple rotational diagram analysis, we derive the excitation conditions and the column densities of the detected COMs. Combining our new results with those previously obtained towards other protostellar objects, we found a good correlation between ethanol, methanol and glycolaldehyde. We discuss the implications of these results on the possible formation routes of ethanol and glycolaldehyde.

We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae.

We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.

Methanol (CH$_3$OH) is found to be abundant and widespread towards the Central Molecular Zone, the inner few hundred parsecs of our Galaxy. Its origin is, however, not fully understood. It was proposed that the high cosmic ray ionisation rate in this region could lead to a more efficient non-thermal desorption of this species formed on grain surfaces, but it would also mean that this species is destroyed in a relatively short timescale. In a first step, we run chemical models with a high cosmic ray ionisation rate and find that this scenario can only reproduce the lowest abundances of methanol derived in this region ($\sim$10$^{-9}$-10$^{-8}$). In a second step, we investigate another scenario based on episodic explosions of grain mantles. We find a good agreement between the predicted abundances of methanol and the observations. We find that the dominant route for the formation of methanol is through hydrogenation of CO on the grains followed by the desorption due to the grain mantle explosion. The cyclic aspect of this model can explain the widespread presence of methanol without requiring any additional mechanism. We also model silicon monoxide (SiO), another species detected in several molecular clouds of the Galactic Centre. An agreement is found with observations for a high depletion of Si (Si/H $\sim$ 10$^{-8}$) with respect to the solar abundance.

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.

The [CII] 158$\mu$m emission line can arise in all phases of the ISM, therefore being able to disentangle the different contributions is an important yet unresolved problem when undertaking galaxy-wide, integrated [CII] observations. We present a new multi-phase 3D radiative transfer interface that couples Starburst99, a stellar spectrophotometric code, with the photoionisation and astrochemistry codes Mocassin and 3D-PDR. We model entire star forming regions, including the ionised, atomic and molecular phases of the ISM, and apply a Bayesian inference methodology to parametrise how the fraction of the [CII] emission originating from molecular regions, $f_{[CII],mol}$, varies as a function of typical integrated properties of galaxies in the local Universe. The main parameters responsible for the variations of $f_{[CII],mol}$ are specific star formation rate (sSFR), gas phase metallicity, HII region electron number density ($n_e$), and dust mass fraction. For example, $f_{[CII],mol}$ can increase from 60% to 80% when either $n_e$ increases from 10$^{1.5}$ to 10$^{2.5}$cm$^{-3}$, or SSFR decreases from $10^{-9.6}$ to $10^{-10.6}$ yr$^{-1}$. Our model predicts for the Milky Way that $f_{[CII],mol}$$=75.8\pm5.9$%, in agreement with the measured value of 75%. When applying the new prescription to a complete sample of galaxies from the Herschel Reference Survey (HRS), we find that anywhere from 60 to 80% of the total integrated [CII] emission arises from molecular regions.

Cyanopolyynes are chains of carbon atoms with an atom of hydrogen and a CN group on either side. They are detected almost everywhere in the ISM, as well as in comets. In the past, they have been used to constrain the age of some molecular clouds, since their abundance is predicted to be a strong function of time. We present an extensive study of the cyanopolyynes distribution in the solar-type protostar IRAS16293-2422 based on TIMASSS IRAM-30m observations. The goals are (i) to obtain a census of the cyanopolyynes in this source and of their isotopologues; (ii) to derive how their abundance varies across the protostar envelope; and (iii) to obtain constraints on the history of IRAS16293-2422. We detect several lines from HC3N and HC5N, and report the first detection of DC3N, in a solar-type protostar. We found that the HC3N abundance is roughly constant (~1.3x10^(-11)) in the outer cold envelope of IRAS16293-2422, and it increases by about a factor 100 in the inner region where Tdust>80K. The HC5N has an abundance similar to HC3N in the outer envelope and about a factor of ten lower in the inner region. The HC3N abundance derived in the inner region, and where the increase occurs, also provide strong constraints on the time taken for the dust to warm up to 80K, which has to be shorter than ~10^3-10^4yr. Finally, the cyanoacetylene deuteration is about 50\% in the outer envelope and <5$\% in the warm inner region. The relatively low deuteration in the warm region suggests that we are witnessing a fossil of the HC3N abundantly formed in the tenuous phase of the pre-collapse and then frozen into the grain mantles at a later phase. The accurate analysis of the cyanopolyynes in IRAS16293-2422 unveils an important part of its past story. It tells us that IRAS16293-2422 underwent a relatively fast (<10^5yr) collapse and a very fast (<10^3-10^4yr) warming up of the cold material to 80K.

The detection of complex organic molecules (COMs) toward cold sources such as pre-stellar cores (with T<10 K), has challenged our understanding of the formation processes of COMs in the interstellar medium. Recent modelling on COM chemistry at low temperatures has provided new insight into these processes predicting that COM formation depends strongly on parameters such as visual extinction and the level of CO freeze out. We report deep observations of COMs toward two positions in the L1544 pre-stellar core: the dense, highly-extinguished continuum peak with Av>=30 mag within the inner 2700 au; and a low-density shell with average Av~7.5-8 mag located at 4000 au from the core's center and bright in CH3OH. Our observations show that CH3O, CH3OCH3 and CH3CHO are more abundant (by factors ~2-10) toward the low-density shell than toward the continuum peak. Other COMs such as CH3OCHO, c-C3H2O, HCCCHO, CH2CHCN and HCCNC show slight enhancements (by factors <=3) but the associated uncertainties are large. This suggests that COMs are actively formed and already present in the low-density shells of pre-stellar cores. The modelling of the chemistry of O-bearing COMs in L1544 indicates that these species are enhanced in this shell because i) CO starts freezing out onto dust grains driving an active surface chemistry; ii) the visual extinction is sufficiently high to prevent the UV photo-dissociation of COMs by the external interstellar radiation field; and iii) the density is still moderate to prevent severe depletion of COMs onto grains.

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.

Dust growth via accretion of gas species has been proposed as the dominant process to increase the amount of dust in galaxies. We show here that this hypothesis encounters severe difficulties that make it unfit to explain the observed UV and IR properties of such systems, particularly at high redshifts. Dust growth in the diffuse ISM phases is hampered by (a) too slow accretion rates; (b) too high dust temperatures, and (c) the Coulomb barrier that effectively blocks accretion. In molecular clouds these problems are largely alleviated. Grains are cold (but not colder than the CMB temperature, (Tcmb = 20 K at redshift z=6). However, in dense environments accreted materials form icy water mantles, perhaps with impurities. Mantles are immediately photo-desorbed as grains return to the diffuse ISM at the end of the cloud lifetime, thus erasing any memory of the growth. We conclude that dust attenuating stellar light at high-$z$ must be ready-made stardust largely produced in supernova ejecta.

Sulfur-bearing molecules are highly reactive in the gas phase of the ISM. However, the form in which most of the sulfur is locked onto interstellar dust grains is unknown. By taking advantage of the short time-scales of shocks in young molecular outflows, one could track back the main form of sulfur in the ices. In this paper, six transitions of H$_2$S and its isotopologues in the L1157-B1 bowshock have been detected using data from the Herschel-CHESS survey and the IRAM-30m ASAI large program. These detections are used to calculate the properties of H$_2$S gas in L1157-B1 through use of a rotation diagram and to explore the possible carriers of sulfur on the grains. The isotopologue detections allow the first calculation of the H$_2$S deuteration fraction in an outflow from a low mass protostar. The fractional abundance of H$_2$S in the region is found to be 6.0$\times$10$^{-7}$ and the deuteration fraction is 2$\times$10$^{-2}$. In order to investigate the form of sulfur on the grains, a chemical model is run with four different networks, each with different branching ratios for the freeze out of sulfur bearing species into molecules such as OCS and H$_2$S. It is found that the model best fits the data when at least half of each sulfur bearing species hydrogenates when freezing. We therefore conclude that a significant fraction of sulfur in L1157-B1 is likely to be locked in H$_2$S on the grains.

As part of the Large Program ASAI (Astrochemical Surveys At IRAM), we have used the IRAM 30m telescope to lead a systematic search for the emission of rotational transitions of P-bearing species between 80 and 350 GHz towards L1157-B1, a shock position in the solar-type star forming region L1157. We report the detection of several transitions of PN and, for the first time, of prebiotic molecule PO. None of these species are detected towards the driving protostar of the outflow L1157-mm. Analysis of the line profiles shows that PN arises from the outflow cavity, where SiO, a strong shock tracer, is produced. Radiative transfer analysis yields an abundance of 2.5e-9 and 0.9e-9 for PO and PN, respectively. These results imply a strong depletion (approx 100) of Phosphorus in the quiescent cloud gas. Shock modelling shows that atomic N plays a major role in the chemistry of PO and PN. The relative abundance of PO and PN brings constraints both on the duration of the pre-shock phase, which has to be about 1 Myr, and on the shock parameters. The maximum temperature in the shock has to be larger than 4000K, which implies a shock velocity of 40 km/s.

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.

We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

Previous work by various authors has suggested that the detection by Herschel/HIFI of nitrogen hydrides along the low density lines of sight towards G10.6-0.4 (W31C) cannot be accounted for by gas-phase chemical models. In this paper we investigate the role of surface reactions on dust grains in diffuse regions, and we find that formation of the hydrides by surface reactions on dust grains with efficiency comparable to that for H$_2$ formation reconciles models with observations of nitrogen hydrides. However, similar surface reactions do not contribute significantly to the hydrides of O$^+$ ions detected by Herschel/HIFI present along many sight lines in the Galaxy. The O$^+$ hydrides can be accounted for by conventional gas-phase chemistry either in diffuse clouds of very low density with normal cosmic ray fluxes or in somewhat denser diffuse clouds with high cosmic ray fluxes. Hydride chemistry in dense dark clouds appears to be dominated by gas-phase ion-molecule reactions.

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

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.

It is now well established that chemistry in external galaxies is rich and complex. In this review I will explore whether one can use molecular emissions to determine their physical conditions. There are several considerations to bear in mind when using molecular emission, and in particular molecular ratios, to determine the densities, temperatures and energetics of a galaxy, which I will briefly summarise here. I will then present an example of a study that uses multiple chemical and radiative transfer analyses in order to tackle the too often neglected `degeneracies' implicit in the interpretation of molecular ratios and show that only via such analyses combined with multi-species and multi-lines high spatial resolution data one can truly make molecules into powerful diagnostics of the evolution and distribution of molecular gas.

We explore the potential of imaging vibrationally excited molecular emission at high angular resolution to better understand the morphology and physical structure of the dense gas in Arp~220 and to gain insight into the nature of the nuclear powering sources. Vibrationally excited emission of HCN is detected in both nuclei with a very high ratio relative to the total $L_{FIR}$, higher than in any other observed galaxy and well above what is observed in Galactic hot cores. HCN $v_2=1f$ is observed to be marginally resolved in $\sim60\times50$~pc regions inside the dusty $\sim100$~pc sized nuclear cores. Its emission is centered on our derived individual nuclear velocities based on HCO$^+$ emission ($V_{WN}=5342\pm4$ and $V_{EN}=5454\pm8$~\kms, for the western and eastern nucleus, respectively). With virial masses within $r\sim25-30$~pc based on the HCN~$v_2=1f$ line widths, we estimate gas surface densities (gas fraction $f_g=0.1$) of $3\pm0.3\times10^4~M_\odot~\rm pc^{-2}$ (WN) and $1.1\pm0.1\times10^4~M_\odot~\rm pc^{-2}$ (EN). The $4-3/3-2$ flux density ratio could be consistent with optically thick emission, which would further constrain the size of the emitting region to $>15$~pc (EN) and $>22$~pc (WN). The absorption systems that may hide up to $70\%$ of the HCN and HCO$^+$ emission are found at velocities of $-50$~\kms~(EN) and $6$, $-140$, and $-575$~\kms (WN) relative to velocities of the nuclei. Blueshifted absorptions are the evidence of outflowing motions from both nuclei. The bright vibrational emission implies the existence of a hot dust region radiatively pumping these transitions. We find evidence of a strong temperature gradient that would be responsible for both the HCN $v_2$ pumping and the absorbed profiles from the vibrational ground state as a result of both continuum and self-absorption by cooler foreground gas.

A fraction of the missing sulfur in dense clouds and circumstellar regions could be in the form of three species not yet de- tected in the interstellar medium: H2S2, HS2, and S2 according to experimental simulations performed under astrophysically relevant conditions. These S-S bonded molecules can be formed by the energetic processing of H2S-bearing ice mantles on dust grains, and subsequently desorb to the gas phase. The detection of these species could partially solve the sulfur depletion problem, and would help to improve our knowledge of the poorly known chemistry of sulfur in the interstellar medium. To this purpose we performed dedicated ground-based observations toward the low-mass warm core IRAS16293-2422. Observations in the submillimeter regime were obtained with the APEX 12 m telescope during 15 hours of observation, targeting a wide selection of the predicted rotational transitions of the three molecules. The 1{\sigma} noise rms values were extracted in the spectral regions where the targeted species should have been detected. These values were a factor of 2-7 lower than those reached by previous observations toward the same source, and allowed us to estimate a 1{\sigma} upper limit to their molecular abundances of 8.1E-9 , 1.1E-8 , and 2.9E-7 relative to H2, for H2S2, HS2, and S2, respectively. The upper limit abundances of the three molecules containing the S2 unit are up to two orders of magnitude lower than the H2S abundance in the source, and one order of magnitude lower than the expected abundances from the experimental simulations using ice analogs. Subsequent gas-phase chemistry after desorption could lower the abundances of the three species to undetectable levels in our observations.

We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

A gas-grain time dependent chemical code, UCL\_CHEM, has been used to investigate the possibility of using chemical tracers to differentiate between the possible formation mechanisms of brown dwarfs. In this work, we model the formation of a pre-brown dwarf core through turbulent fragmentation by following the depth-dependent chemistry in a molecular cloud through the step change in density associated with an isothermal shock and the subsequent freefall collapse once a bound core is produced. Trends in the fractional abundance of molecules commonly observed in star forming cores are then explored to find a diagnostic for identifying brown dwarf mass cores formed through turbulence. We find that the cores produced by our models would be bright in CO and NH$_3$ but not in HCO$^+$. This differentiates them from models using purely freefall collapse as such models produce cores that would have detectable transitions from all three molecules.

To date O2 has definitely been detected in only two sources, namely rho Oph A and Orion, reflecting the extremely low abundance of O2 in the interstellar medium. One of the sources in the HOP program is the +50 km/s Cloud in the Sgr A Complex in the centre of the Milky Way. The Herschel HIFI is used to search for the 487 and 774 GHz emission lines of O2. No O2 emission is detected towards the Sgr A +50 km/s Cloud, but a number of strong emission lines of methanol (CH3OH) and absorption lines of chloronium (H2Cl+) are observed. A 3 sigma upper limit for the fractional abundance ratio of (O2)/(H2) in the Sgr A +50 km/s Cloud is found to be X(O2) less than 5 x 10(-8). However, since we can find no other realistic molecular candidate than O2 itself, we very tentatively suggest that two weak absorption lines at 487.261 and 487.302 GHz may be caused by the 487 GHz line of O2 in two foreground spiral arm clouds. By considering that the absorption may only be apparent, the estimated upper limit to the O2 abundance of less than (10-20) x 10(-6) in these foreground clouds is very high. This abundance limit was determined also using Odin non-detection limits. If the absorption is due to a differential Herschel OFF-ON emission, the O2 fractional abundance may be of the order of (5-10) x 10(-6). With the assumption of pure absorption by foreground clouds, the unreasonably high abundance of (1.4-2.8) x 10(-4) was obtained. The rotation temperatures for CH3OH-A and CH3OH-E lines in the +50 km/s Cloud are found to be 64 and 79 K, respectively, and the fractional abundance of CH3OH is approximately 5 x 10(-7).

The circumnuclear disk (CND) of the Galactic Center is exposed to many energetic phenomena coming from the supermassive black hole Sgr A* and stellar activities. These energetic activities can affect the chemical composition in the CND by the interaction with UV-photons, cosmic-rays, X-rays, and shock waves. We aim to constrain the physical conditions present in the CND by chemical modeling of observed molecular species detected towards it. We analyzed a selected set of molecular line data taken toward a position in the southwest lobe of the CND with the IRAM 30m and APEX 12-meter telescopes and derived the column density of each molecule using a large velocity gradient (LVG) analysis. The determined chemical composition is compared with a time-dependent gas-grain chemical model based on the UCL\_CHEM code that includes the effects of shock waves with varying physical parameters. Molecules such as CO, HCN, HCO$^+$, HNC, CS, SO, SiO, NO, CN, H$_2$CO, HC$_3$N, N$_2$H$^+$ and H$_3$O$^+$ are detected and their column densities are obtained. Total hydrogen densities obtained from LVG analysis range between $2 \times 10^4$ and $1 \times 10^6\,$cm$^{-3}$ and most species indicate values around several $\times 10^5\,$cm$^{-3}$, which are lower than values corresponding to the Roche limit, which shows that the CND is tidally unstable. The chemical models show good agreement with the observations in cases where the density is $\sim10^4\,$cm$^{-3}$, the cosmic-ray ionization rate is high, $>10^{-15} \,$s$^{-1}$, or shocks with velocities $> 40\,$km s$^{-1}$ have occurred. Comparison of models and observations favors a scenario where the cosmic-ray ionization rate in the CND is high, but precise effects of other factors such as shocks, density structures, UV-photons and X-rays from the Sgr A* must be examined with higher spatial resolution data.

The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium either dynamically or chemically. In this paper we present TORUS-3DPDR, a new self-consistent code for treating the chemistry of three-dimensional photoionization and photodissociation regions. It is an integrated code coupling the two codes TORUS, a hydrodynamics and Monte Carlo radiation transport code, and 3D-PDR, a photodissociation regions code. The new code uses a Monte Carlo radiative transfer scheme to account for the propagation of the ionizing radiation including the diffusive component as well as a ray-tracing scheme based on the HEALPix package in order to account for the escape probability and column density calculations. Here, we present the numerical techniques we followed and we show the capabilities of the new code in modelling three-dimensional objects including single or multiple sources. We discuss the effects introduced by the diffusive component of the UV field in determining the thermal balance of PDRs as well as the effects introduced by a multiple sources treatment of the radiation field. We find that diffuse radiation can positively contribute to the formation of CO. With this new code, three-dimensional synthetic observations for the major cooling lines are possible, for making feasible a detailed comparison between hydrodynamical simulations and observations.

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-{\mu}m continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-{\sigma} surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average.

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($\nu_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $\nu_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($\nu$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

We obtained an ALMA Cycle 0 spectral scan of the dusty LIRG NGC 4418, spanning a total of 70.7 GHz in bands 3, 6, and 7. We use a combined local thermal equilibrium (LTE) and non-LTE (NLTE) fit of the spectrum in order to identify the molecular species and derive column densities and excitation temperatures. We derive molecular abundances and compare them with other Galactic and extragalactic sources by means of a principal component analysis. We detect 317 emission lines from a total of 45 molecular species, including 15 isotopic substitutions and six vibrationally excited variants. Our LTE/NLTE fit find kinetic temperatures from 20 to 350 K, and densities between 10$^5$ and 10$^7$ cm$^{-3}$. The spectrum is dominated by vibrationally excited HC$_3$N, HCN, and HNC, with vibrational temperatures from 300 to 450 K. We find high abundances of HC$_3$N, SiO, H$_2$S, and c-HCCCH and a low CH$_3$OH abundance. A principal component analysis shows that NGC 4418 and Arp 220 share very similar molecular abundances and excitation, which clearly set them apart from other Galactic and extragalactic environments. The similar molecular abundances observed towards NCG 4418 and Arp 220 are consistent with a hot gas-phase chemistry, with the relative abundances of SiO and CH$_3$OH being regulated by shocks and X-ray driven dissociation. The bright emission from vibrationally excited species confirms the presence of a compact IR source, with an effective diameter $<$5 pc and brightness temperatures $>$350 K. The molecular abundances and the vibrationally excited spectrum are consistent with a young AGN/starburst system. We suggest that NGC 4418 may be a template for a new kind of chemistry and excitation, typical of compact obscured nuclei (CON). Because of the narrow line widths and bright molecular emission, NGC 4418 is the ideal target for further studies of the chemistry in CONs.

Sulfur appears to be depleted by an order of magnitude or more from its elemental abundance in star-forming regions. In the last few years, numerous observations and experiments have been performed in order to to understand the reasons behind this depletion without providing a satisfactory explanation of the sulfur chemistry towards high-mass star-forming cores. Several sulfur-bearing molecules have been observed in these regions, and yet none are abundant enough to make up the gas-phase deficit. Where, then, does this hidden sulfur reside? This paper represents a step forward in our understanding of the interactions among the various S-bearing species. We have incorporated recent experimental and theoretical data into a chemical model of a hot molecular core in order to see whether they give any indication of the identity of the sulfur sink in these dense regions. Despite our model producing reasonable agreement with both solid-phase and gas-phase abundances of many sulfur-bearing species, we find that the sulfur residue detected in recent experiments takes up only ~6 per cent of the available sulfur in our simulations, rather than dominating the sulfur budget.

Aims. We observe the dense gas tracer CS in two nearby starburst galaxies to determine how the conditions of the dense gas varies across the circumnuclear regions in starburst galaxies. Methods. Using the IRAM-30m telescope, we mapped the distribution of the CS(2-1) and CS(3-2) lines in the circumnuclear regions of the nearby starburst galaxies NGC 3079 and NGC 6946. We also detected the formaldehyde (H2CO) and methanol (CH3OH) in both galaxies. We marginally detect the isotopologue C34S. Results. We calculate column densities under LTE conditions for CS and CH3OH. Using the detections accumulated here to guide our inputs, we link a time and depth dependent chemical model with a molecular line radiative transfer model; we reproduce the observations, showing how conditions where CS is present are likely to vary away from the galactic centres. Conclusions. Using the rotational diagram method for CH3OH, we obtain a lower limit temperature of 14 K. In addition to this, by comparing the chemical and radiative transfer models to observations, we determine the properties of the dense gas as traced by CS (and CH3OH). We also estimate the quantity of the dense gas. We find that, provided that there are a between 10^5 and 10^6 dense cores in our beam, for both target galaxies, emission of CS from warm (T = 100 - 400 K), dense (n(H2) = 10^5-6 cm-3) cores, possibly with a high cosmic ray ionisation rate (zeta = 100 zeta0) best describes conditions for our central pointing. In NGC 6946, conditions are generally cooler and/or less dense further from the centre, whereas in NGC 3079, conditions are more uniform. The inclusion of shocks allows for more efficient CS formation, leading to an order of magnitude less dense gas being required to replicate observations in some cases.

We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.

We report on the effects of cosmic rays (CRs) on the abundance of CO in $\rm H_2$ clouds under conditions typical for star-forming galaxies in the Universe. We discover that this most important molecule for tracing H$_2$ gas is very effectively destroyed in ISM environments with CR energy densities $\rm U_{CR}\sim(50-10^{3})\times U_{CR,Gal}$, a range expected in numerous star-forming systems throughout the Universe. This density-dependent effect operates volumetrically rather than only on molecular cloud surfaces (i.e. unlike FUV radiation that also destroys CO), and is facilitated by: a) the direct destruction of CO by CRs, and b) a reaction channel activated by CR-produced He$^{+}$. The effect we uncover is strong enough to render Milky-Way type Giant Molecular Clouds (GMCs) very CO-poor (and thus CO-untraceable), even in ISM environments with rather modestly enhanced average CR energy densities of $\rm U_{CR}\sim(10-50)\times\rm U_{CR,Gal}$. We conclude that the CR-induced destruction of CO in molecular clouds, unhindered by dust absorption, is perhaps the single most important factor controlling the CO-visibility of molecular gas in vigorously star-forming galaxies. We anticipate that a second order effect of this CO destruction mechanism will be to make the H$_2$ distribution in the gas-rich disks of such galaxies appear much clumpier in CO $J$=1--0, 2--1 line emission than it actually is. Finally we give an analytical approximation of the CO/H$_2$ abundance ratio as a function of gas density and CR energy density for use in galaxy-size or cosmological hydrodynamical simulations, and propose some key observational tests.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Molecular complexity builds up at each step of the Sun-like star formation process, starting from simple molecules and ending up in large polyatomic species. Complex organic molecules (COMs; such as methyl formate, HCOOCH$_3$, dymethyl ether, CH$_3$OCH$_3$, formamide, NH$_2$CHO, or glycoaldehyde, HCOCH$_2$OH) are formed in all the components of the star formation recipe (e.g. pre-stellar cores, hot-corinos, circumstellar disks, shocks induced by fast jets), due to ice grain mantle sublimation or sputtering as well as gas-phase reactions. Understanding in great detail the involved processes is likely the only way to predict the ultimate molecular complexity reached in the ISM, as the detection of large molecules is increasingly more difficult with the increase of the number of atoms constituting them. Thanks to the recent spectacular progress of astronomical observations, due to the Herschel (sub-mm and IR), IRAM and SMA (mm and sub-mm), and NRAO (cm) telescopes, an enormous activity is being developed in the field of Astrochemistry, extending from astronomical observatories to chemical laboratories. We are involved in several observational projects providing unbiased spectral surveys (in the 80-300 and 500-2000 GHz ranges) with unprecedented sensitivity of templates of dense cores and protostars. Forests of COM lines have been detected. In this chapter we will focus on the chemistry of both cold prestellar cores and hot shocked regions, (i) reviewing results and open questions provided by mm-FIR observations, and (ii) showing the need of carrying on the observations of COMs at lower frequencies, where SKA will operate. We will also emphasize the importance of analysing the spectra by the light of the experimental studies performed by our team, who is investigating the chemical effects induced by ionising radiation bombarding astrophysically relevant ices.

The formation of complex organic molecules (COMs) in protostellar environments is a hotly debated topic. In particular, the relative importance of the gas phase processes as compared to a direct formation of COMs on the dust grain surfaces is so far unknown. We report here the first high-resolution images of acetaldehyde (CH$_3$CHO) emission towards the chemically rich protostellar shock L1157-B1, obtained at 2 mm with the IRAM Plateau de Bure interferometer. Six blueshifted CH$_3$CHO lines with $E_{\rm u}$ = 26-35 K have been detected. The acetaldehyde spatial distribution follows the young ($\sim$ 2000 yr) outflow cavity produced by the impact of the jet with the ambient medium, indicating that this COM is closely associated with the region enriched by iced species evaporated from dust mantles and released into the gas phase. A high CH$_3$CHO relative abundance, 2-3 $\times$ 10$^{-8}$, is inferred, similarly to what found in hot-corinos. Astrochemical modelling indicates that gas phase reactions can produce the observed quantity of acetaldehyde only if a large fraction of carbon, of the order of 0.1%, is locked into iced hydrocarbons.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

We present a multiline CS survey towards the brightest bow-shock B1 in the prototypical chemically active protostellar outflow L1157. We made use of (sub-)mm data obtained in the framework of the Chemical HErschel Surveys of Star forming regions (CHESS) and Astrochemical Surveys at IRAM (ASAI) key science programs. We detected $^{12}$C$^{32}$S, $^{12}$C$^{34}$S, $^{13}$C$^{32}$S, and $^{12}$C$^{33}$S emissions, for a total of 18 transitions, with $E_{\rm u}$ up to $\sim$ 180 K. The unprecedented sensitivity of the survey allows us to carefully analyse the line profiles, revealing high-velocity emission, up to 20 km s$^{-1}$ with respect to the systemic. The profiles can be well fitted by a combination of two exponential laws that are remarkably similar to what previously found using CO. These components have been related to the cavity walls produced by the $\sim$ 2000 yr B1 shock and the older ($\sim$ 4000 yr) B2 shock, respectively. The combination of low- and high-excitation CS emission was used to properly sample the different physical components expected in a shocked region. Our CS observations show that this molecule is highlighting the dense, $n_{\rm H_2}$ = 1--5 $\times$ 10$^{5}$ cm$^{-3}$, cavity walls produced by the episodic outflow in L1157. In addition, the highest excitation (E$_u$ $\geq$ 130 K) CS lines provide us with the signature of denser (1--5 $\times$ 10$^{6}$ cm$^{-3}$) gas, associated with a molecular reformation zone of a dissociative J-type shock, which is expected to arise where the precessing jet impacting the molecular cavities. The CS fractional abundance increases up to $\sim$ 10$^{-7}$ in all the kinematical components. This value is consistent with what previously found for prototypical protostars and it is in agreement with the prediction of the abundances obtained via the chemical code Astrochem.

Understanding the physical conditions of dark molecular clouds and star forming regions is an inverse problem subject to complicated chemistry that varies non-linearly with time and the physical environment. In this paper we apply a Bayesian approach based on a Markov Chain Monte Carlo (MCMC) method for solving the non-linear inverse problems encountered in astrochemical modelling. We use observations for ice and gas species in dark molecular clouds and a time dependent, gas grain chemical model to infer the values of the physical and chemical parameters that characterize quiescent regions of molecular clouds. We show evidence that in high dimensional problems, MCMC algorithms provide a more efficient and complete solution than more classical strategies. The results of our MCMC method enable us to derive statistical estimates and uncertainties for the physical parameters of interest as a result of the Bayesian treatment.

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

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

New laboratory data of CH$_2$CHCN (vinyl cyanide) in its ground and vibrationally excited states at the microwave to THz domain allow searching for these excited state transitions in the Orion-KL line survey. Frequency-modulated spectrometers combined into a single broadband 50-1900 GHz spectrum provided measurements of CH$_2$CHCN covering a spectral range of 18-1893 GHz, whose assignments was confirmed by Stark modulation spectra in the 18-40 GHz region and by ab-initio anharmonic force field calculations. For analyzing the emission lines of CH$_2$CHCN species detected in Orion-KL we used the excitation and radiative transfer code (MADEX) at LTE conditions. The rotational transitions of the ground state of this molecule emerge from four cloud components of hot core nature which trace the physical and chemical conditions of high mass star forming regions in the Orion-KL Nebula. The total column density of CH$_2$CHCN in the ground state is (3.0$\pm$0.9)x10$^{15}$ cm$^{-2}$. We report on the first interstellar detection of transitions in the v10=1/(v11=1,v15=1) dyad in space, and in the v11=2 and v11=3 states in Orion-KL. The lowest energy vibrationally excited states of vinyl cyanide such as v11=1 (at 328.5 K), v15=1 (at 478.6 K), v11=2 (at 657.8 K), the v10=1/(v11=1,v15=1) dyad (at 806.4/809.9 K), and v11=3 (at 987.9 K) are populated under warm and dense conditions, so they probe the hottest parts of the Orion-KL source. Column density and rotational and vibrational temperatures for CH$_2$CHCN in their ground and excited states, as well as for the isotopologues, have been constrained by means of a sample of more than 1000 lines in this survey. Moreover, we present the detection of methyl isocyanide (CH$_3$NC) for the first time in Orion-KL and a tentative detection of vinyl isocyanide (CH$_2$CHNC) and give column density ratios between the cyanide and isocyanide isomers.

We present the results of a line identification analysis using data from the IRAM Plateau de Bure Inferferometer, focusing on six massive star-forming hot cores: G31.41+0.31, G29.96-0.02, G19.61-0.23, G10.62-0.38, G24.78+0.08A1 and G24.78+0.08A2. We identify several transitions of vibrationally excited methyl formate (HCOOCH$_3$) for the first time in these objects as well as transitions of other complex molecules, including ethyl cyanide (C$_2$H$_5$CN), and isocyanic acid (HNCO). We also postulate a detection of one transition of glycolaldehyde (CH$_2$(OH)CHO) in two new hot cores. We find G29.96-0.02, G19.61-0.23, G24.78+0.08A1 and 24.78+0.08A2 to be chemically very similar. G31.41+0.31, however, is chemically different: it manifests a larger chemical inventory and has significantly larger column densities. We suggest that it may represent a different evolutionary stage to the other hot cores in the sample, or it may surround a star with a higher mass. We derive column densities for methyl formate in G31.41+0.31, using the rotation diagram method, of $\times$10$^{17}$ cm$^{-2}$ and a T$_{rot}$ of $\sim$170 K. For G29.96-0.02, G24.78+0.08A1 and G24.78+0.08A2, glycolaldehyde, methyl formate and methyl cyanide all seem to trace the same material and peak at roughly the same position towards the dust emission peak. For G31.41+0.31, however, glycolaldehyde shows a different distribution to methyl formate and methyl cyanide and seems to trace the densest, most compact inner part of hot cores.

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 standard interstellar ratio of deuterium to hydrogen (D/H) atoms is $\sim 1.5 \times 10^{-5}$. However, the deuterium fractionation is in fact found to be enhanced, to different degrees, in cold, dark cores, hot cores around massive star forming regions, lukewarm cores, and warm cores ({\it hereafter}, hot corinos) around low-mass star forming regions. In this paper, we investigate the overall differences in the deuterium chemistry between hot cores and hot corinos. We have modelled the chemistry of dense gas around low-mass and massive star forming regions using a gas-grain chemical model. We investigate the influence of varying the core density, the depletion efficiency of gaseous species on to dust grains, the collapse mode and the final mass of the protostar on the chemical evolution of star forming regions. We find that the deuterium chemistry is, in general, most sensitive to variations of the depletion efficiency on to grain surfaces, in agreement with observations. In addition, the results showed that the chemistry is more sensitive to changes in the final density of the collapsing core in hot cores than in hot corinos. Finally, we find that ratios of deuterated sulphur bearing species in dense gas around hot cores and corinos may be good evolutionary indicators in a similar way as their non deuterated counterparts.

We study the sulphur chemistry evolution in the Orion KL along the gas and grain phases of the cloud. We investigate the processes that dominate the sulphur chemistry and to determine how physical and chemical parameters, such as the final star mass and the initial elemental abundances, influence the evolution of the hot core and of the surrounding outflows and shocked gas (the plateau). We independently modelled the chemistry evolution of both components using the time-dependent gas-grain model UCL_CHEM and considering two different phase calculations. Phase I starts with the collapsing cloud and the depletion of atoms and molecules onto grain surfaces. Phase II starts when a central protostar is formed and the evaporation from grains takes place. We show how the gas density, the gas depletion efficiency, the initial sulphur abundance, the shocked gas temperature and the different chemical paths on the grains leading to different reservoirs of sulphur on the mantles affect sulphur-bearing molecules at different evolutionary stages. We also compare the predicted column densities with those inferred from observations of the species SO, SO2, CS, OCS, H2S and H2CS. The models that reproduce the observations of the largest number of sulphur-bearing species are those with an initial sulphur abundance of 0.1 times the sulphur solar abundance and a density of at least n_H=5x10^6 cm^-3 in the shocked gas region. We conclude that most of the sulphur atoms were ionised during Phase I, consistent with an inhomogeneous and clumpy region where the UV interstellar radiation penetrates leading to sulphur ionisation. We also conclude that the main sulphur reservoir on the ice mantles was H2S. In addition, we deduce that a chemical transition currently takes place in the shocked gas, where SO and SO2 gas-phase formation reactions change from being dominated by O2 to being dominated by OH.

We present a model of the photodissociation regions of NGC 4038, which is part of the Antennae galaxies. We have considered one-dimensional slabs of uniform density all having a maximum $A_V=10\,{\rm mag}$, interacting with plane-parallel radiation. The density range in our simulations spans four orders of magnitude ($100\le n\le 10^6\,{\rm cm}^{-3}$) and the UV field strength spans more than three orders of magnitude ($10\le\chi\le10^{4.5}$ multiples of the ${\rm Draine}$ field), from which we generated a grid of about 1400 simulations. We compare our results with Herschel SPIRE-FTS, CSO and ISO-LWS observations of eight CO transition lines ($J=1-0$ to $8-7$) and the [CI] 609\mu m and [OI] 146\mu m fine structure lines. We find that the molecular and atomic emission lines trace different gas components of NGC 4038, thus single emission models are insufficient to reproduce the observed values. In general, low-$J$ CO transition lines correspond to either low density regions interacting with low UV field strengths, or high density regions interacting with high UV field strengths. Higher $J$ CO transition lines are less dependent on the UV field strength and are fitted by gas with density $n\sim10^{4.5}-10^{5.2}\,{\rm cm}^{-3}$. We find that the observed fine structure line ratio of [CI] 609\mu m/[OI] 146\mu m is reproduced by clouds subject to weaker UV fields compared to the CO lines. We make estimates of the $X_{\rm CO}$ factor which relates the CO emission with the column density of molecular hydrogen, and find that it is less than the canonical Milky Way value.

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.

Glycine (NH2CH2COOH) is the simplest amino acid relevant for life. Its detection in the interstellar medium is key to understand the formation mechanisms of pre-biotic molecules and their subsequent delivery onto planetary systems. Glycine has extensively been searched for toward hot molecular cores, although these studies did not yield any firm detection. In contrast to hot cores, low-mass star forming regions, and in particular their earliest stages represented by cold pre-stellar cores, may be better suited for the detection of glycine as well as more relevant for the study of pre-biotic chemistry in young Solar System analogs. We present 1D spherically symmetric radiative transfer calculations of the glycine emission expected to arise from the low-mass pre-stellar core L1544. Water vapour has recently been reported toward this core, indicating that a small fraction of the grain mantles in L1544 (~0.5%) has been injected into the gas phase. Assuming that glycine is photo-desorbed together with water in L1544, and considering a solid abundance of glycine on ices of ~1E-4 with respect to water, our calculations reveal that several glycine lines between 67 GHz and 80 GHz have peak intensities larger than 10 mK. These results show for the first time that glycine could reach detectable levels in cold objects such as L1544. This opens up the possibility to detect glycine, and other pre-biotic species, at the coldest and earliest stages in the formation of Solar-type systems with near-future instrumentation such as the Band 2 receivers of ALMA.

[Abridged] Ethylene oxide and its isomer acetaldehyde are important complex organic molecules because of their potential role in the formation of amino acids. Despite the fact that acetaldehyde is ubiquitous in the interstellar medium, ethylene oxide has not yet been detected in cold sources. We aim to understand the chemistry of the formation and loss of ethylene oxide in hot and cold interstellar objects (i) by including in a revised gas-grain network some recent experimental results on grain surfaces and (ii) by comparison with the chemical behaviour of its isomer, acetaldehyde. We test the code for the case of a hot core. The model allows us to predict the gaseous and solid ethylene oxide abundances during a cooling-down phase prior to star formation and during the subsequent warm-up phase. We can therefore predict at what temperatures ethylene oxide forms on grain surfaces and at what temperature it starts to desorb into the gas phase. The model reproduces the observed gaseous abundances of ethylene oxide and acetaldehyde towards high-mass star-forming regions. In addition, our results show that ethylene oxide may be present in outer and cooler regions of hot cores where its isomer has already been detected. Despite their different chemical structures, the chemistry of ethylene oxide is coupled to that of acetaldehyde, suggesting that acetaldehyde may be used as a tracer for ethylene oxide towards cold cores.

In order to study the temperature distribution of the extended gas within the Orion Kleinmann-Low nebula, we have mapped the emission by methyl cyanide (CH3CN) in its J=6_K-5_K, J=12_K-11_K, J=13_K-12_K, and J=14_K-13_K transitions at an average angular resolution of ~10 arcsec (22 arcsec for the 6_K-5_K lines), as part of a new 2D line survey of this region using the IRAM 30m telescope. These fully sampled maps show extended emission from warm gas to the northeast of IRc2 and the distinct kinematic signatures of the hot core and compact ridge source components. We have constructed population diagrams for the four sets of K-ladder emission lines at each position in the maps and have derived rotational excitation temperatures and total beam-averaged column densities from the fitted slopes. In addition, we have fitted LVG model spectra to the observations to determine best-fit physical parameters at each map position, yielding the distribution of kinetic temperatures across the region. The resulting temperature maps reveal a region of hot (T > 350 K) material surrounding the northeastern edge of the hot core, whereas the column density distribution is more uniform and peaks near the position of IRc2. We attribute this region of hot gas to shock heating caused by the impact of outflowing material from active star formation in the region, as indicated by the presence of broad CH3CN lines. This scenario is consistent with predictions from C-shock chemical models that suggest that gas-phase methyl cyanide survives in the post-shock gas and can be somewhat enhanced due to sputtering of grain mantles in the passing shock front.

We explore the utility of CI as an alternative high-fidelity gas mass tracer for Galactic molecular clouds. We evaluate the X$_{\rm CI}$-factor for the 609 $\mu$m carbon line, the analog of the CO X-factor, which is the ratio of the H$_2$ column density to the integrated $^{12}$CO(1-0) line intensity. We use 3D-PDR to post-process hydrodynamic simulations of turbulent, star-forming clouds. We compare the emission of CI and CO for model clouds irradiated by 1 and 10 times the average background and demonstrate that CI is a comparable or superior tracer of the molecular gas distribution for column densities up to $6 \times 10^{23}$ cm$^{-2}$. Our results hold for both reduced and full chemical networks. For our fiducial Galactic cloud we derive an average $X_{\rm CO}$ of $3.0\times 10^{20}$ cm$^{-2}$K$^{-1}$km$^{-1}$s and $X_{\rm CI}$ of $1.1\times 10^{21}$ cm$^{-2}$K$^{-1}$km$^{-1}$s.

We aim to quantify the molecular gas reservoir in a subset of 6 low-metallicity galaxies from the Herschel Dwarf Galaxy Survey with newly acquired CO data, and link this reservoir to the observed star formation activity. We present CO(1-0), CO(2-1), and CO(3-2) observations obtained at the ATNF Mopra 22-m, APEX, and IRAM 30-m telescopes, as well as [CII] 157um and [OI] 63um observations obtained with the Herschel/PACS spectrometer in the 6 galaxies: Haro11, Mrk1089, Mrk930, NGC4861, NGC625, and UM311. We derive molecular gas mass from several methods including the use of the CO-to-H2 conversion factor Xco (both Galactic and metallicity-scaled values) and of dust measurements. The molecular and atomic gas reservoirs are compared to the star formation activity. We also constrain the physical conditions of the molecular clouds using the non-LTE code RADEX and the spectral synthesis code Cloudy. We detect CO in 5 of the 6 galaxies, including first detections in Haro11 (Z~0.4 Zsun), Mrk930 (0.2 Zsun), and UM311 (0.5 Zsun), but CO remains undetected in NGC4861 (0.2 Zsun). The CO luminosities are low while [CII] is bright in these galaxies, resulting in [CII]/CO(1-0)>10000. Our dwarf galaxies are in relatively good agreement with the Schmidt-Kennicutt relation for total gas. They show short molecular depletion time scales, even when considering metallicity-scaled Xco factors. Those galaxies are dominated by their HI gas, except Haro11 which has high star formation efficiency and is dominated by ionized and molecular gas. We determine the mass of each ISM phase in Haro11 using Cloudy and estimate an equivalent Xco factor which is 10 times higher than the Galactic value. Overall, our results confirm the emerging picture that CO suffers from significant selective photodissociation in low-metallicity dwarf galaxies.

Within a few parsecs around the central Black Hole Sgr A*, chemistry in the dense molecular cloud material of the circumnuclear disk (CND) can be affected by many energetic phenomena such as high UV-flux from the massive central star cluster, X-rays from Sgr A*, shock waves, and an enhanced cosmic-ray flux. Recently, spectroscopic surveys with the IRAM 30 meter and the APEX 12 meter telescopes of substantial parts of the 80--500 GHz frequency range were made toward selected positions in and near the CND. These datasets contain lines from the molecules HCN, HCO$^+$, HNC, CS, SO, SiO, CN, H$_2$CO, HC$_3$N, N$_2$H$^+$, H$_3$O$^+$ and others. We conduct Large Velocity Gradient analyses to obtain column densities and total hydrogen densities, $n$, for each species in molecular clouds located in the southwest lobe of CND. The data for the above mentioned molecules indicate 10$^5\,$cm$^{-3} \lesssim n <10^6\,$cm$^{-3}$, which shows that the CND is tidally unstable. The derived chemical composition is compared with a chemical model calculated using the UCL\_CHEM code that includes gas and grain reactions, and the effects of shock waves. Models are run for varying shock velocities, cosmic-ray ionization rates, and number densities. The resulting chemical composition is fitted best to an extremely high value of cosmic-ray ionization rate $\zeta \sim 10^{-14}\,$s$^{-1}$, 3 orders of magnitude higher than the value in regular Galactic molecular clouds, if the pre-shock density is $n=10^5\,$cm$^{-3}$.

Molecular outflows powered by young protostars strongly affect the kinematics and chemistry of the natal molecular cloud through strong shocks resulting in substantial modifications of the abundance of several species. As part of the "Chemical Herschel Surveys of Star forming regions" guaranteed time key program, we aim at investigating the physical and chemical conditions of H20 in the brightest shock region B1 of the L1157 molecular outflow. We observed several ortho- and para-H2O transitions using HIFI and PACS instruments on board Herschel, providing a detailed picture of the kinematics and spatial distribution of the gas. We performed a LVG analysis to derive the physical conditions of H2O shocked material, and ultimately obtain its abundance. We detected 13 H2O lines probing a wide range of excitation conditions. PACS maps reveal that H2O traces weak and extended emission associated with the outflow identified also with HIFI in the o-H2O line at 556.9 GHz, and a compact (~10") bright, higher-excitation region. The LVG analysis of H2O lines in the bow-shock show the presence of two gas components with different excitation conditions: a warm (Tkin~200-300 K) and dense (n(H2)~(1-3)x10^6 cm-3) component with an assumed extent of 10" and a compact (~2"-5") and hot, tenuous (Tkin~900-1400 K, n(H2)~10^3-10^4 cm-3) gas component, which is needed to account for the line fluxes of high Eu transitions. The fractional abundance of the warm and hot H2O gas components is estimated to be (0.7-2)x10^{-6} and (1-3)x10^{-4}, respectively. Finally, we identified an additional component in absorption in the HIFI spectra of H2O lines connecting with the ground state level, probably arising from the photodesorption of icy mantles of a water-enriched layer at the edges of the cloud.

Glycolaldehyde, the simplest monosaccharide sugar, has recently been detected in low- and high-mass star-forming cores. Following on from our previous investigation into glycolaldehyde formation (Woods et al. 2012, ApJ, 750, 19), we now consider a further mechanism for the formation of glycolaldehyde that involves the dimerisation of the formyl radical, HCO. Quantum mechanical investigation of the HCO dimerisation process upon an ice surface is predicted to be barrierless and therefore fast. In an astrophysical context, we show that this mechanism can be very efficient in star-forming cores. It is limited by the availability of the formyl radical, but models suggest that only very small amounts of CO are required to be converted to HCO to meet the observational constraints.

We present a study of cyanoacetylene (HC3N) and cyanodiacetylene (HC5N) in Orion KL, through observations from two line surveys performed with the IRAM 30m telescope and the HIFI instrument on board the Herschel telescope. The frequency ranges covered are 80-280 GHz and 480-1906 GHz. We model the observed lines of HC3N, HC5N, their isotopologues (including DC3N), and vibrational modes, using a non-LTE radiative transfer code. To investigate the chemical origin of HC3N and DC3N in Orion KL, we use a time-dependent chemical model. We detect 40 lines of the ground state of HC3N and 68 lines of its 13C isotopologues. We also detect 297 lines of six vibrational modes of this molecule (nu_7, 2nu_7, 3nu_7, nu_6, nu_5, and nu_6+nu_7) and 35 rotational lines of the ground state of HC5N. We report the first tentative detection of DC3N in a giant molecular cloud with a DC3N/HC3N abundance ratio of 0.015. We provide column densities and isotopic and molecular abundances. We also perform a 2x2" map around Orion IRc2 and we present maps of HC3N lines and maps of lines of the HC3N vibrational modes nu_6 and nu_7. In addition, a comparison of our results for HC3N with those in other clouds allows us to derive correlations between the column density, the FWHM, the mass, and the luminosity of the clouds. The high column densities of HC3N obtained in the hot core, make this molecule an excellent tracer of hot and dense gas. In addition, the large frequency range covered reveals the need to consider a temperature and density gradient in the hot core in order to obtain better line fits. The high D/H ratio (comparable to that obtained in cold clouds) that we derive suggests a deuterium enrichment. Our chemical models indicate that the possible deuterated HC3N present in Orion KL is formed during the gas-phase. This fact provides new hints concerning the processes leading to deuteration.

We present high spatial resolution (750 AU at 250 pc) maps of the B1 shock in the blue lobe of the L1157 outflow in four lines: CS (3-2), CH3OH (3_K-2_K), HC3N (16-15) and p-H2CO (2_02-3_01). The combined analysis of the morphology and spectral profiles has shown that the highest velocity gas is confined in a few compact (~ 5 arcsec) bullets while the lowest velocity gas traces the wall of the gas cavity excavated by the shock expansion. A large velocity gradient model applied to the CS (3-2) and (2-1) lines provides an upper limit of 10^6 cm^-3 to the averaged gas density in B1 and a range of 5x10^3< n(H2)< 5x10^5 cm^-3 for the density of the high velocity bullets. The origin of the bullets is still uncertain: they could be the result of local instabilities produced by the interaction of the jet with the ambient medium or could be clump already present in the ambient medium that are excited and accelerated by the expanding outflow. The column densities of the observed species can be reproduced qualitatively by the presence in B1 of a C-type shock and only models where the gas reaches temperatures of at least 4000 K can reproduce the observed HC3N column density.

We present the first detection of N2H+ towards a low-mass protostellar outflow, namely the L1157-B1 shock, at about 0.1 pc from the protostellar cocoon. The detection was obtained with the IRAM 30-m antenna. We observed emission at 93 GHz due to the J = 1-0 hyperfine lines. The analysis of the emission coupled with the HIFI CHESS multiline CO observations leads to the conclusion that the observed N2H+(1-0) line originates from the dense (> 10^5 cm-3) gas associated with the large (20-25 arcsec) cavities opened by the protostellar wind. We find a N2H+ column density of few 10^12 cm-2 corresponding to an abundance of (2-8) 10^-9. The N2H+ abundance can be matched by a model of quiescent gas evolved for more than 10^4 yr, i.e. for more than the shock kinematical age (about 2000 yr). Modelling of C-shocks confirms that the abundance of N2H+ is not increased by the passage of the shock. In summary, N2H+ is a fossil record of the pre-shock gas, formed when the density of the gas was around 10^4 cm-3, and then further compressed and accelerated by the shock.

Complex organic molecules, such as propylene (CH3CHCH2), are detected in molecular clouds (such as TMC1) with high fractional abundances (~2.10^{-9}, relative to hydrogen) that cannot be explained by gas-phase chemical reactions under normal dark cloud conditions. To obtain such high abundances requires an efficient grain-driven chemistry to be operating, coupled with effective desorption of the complex organics back into the gas-phase. We propose that the mechanism that we have previously described; rapid high density gas-phase chemistry in the gas released following sudden, total, ice mantle sublimation - can explain the high abundances, without recourse to ill-defined surface chemical pathways. Our model provides a natural explanation for why it is that some sources harbour high abundances of propylene, whilst others do not; based on the age and level of dynamical activity within the source (which affects the ice composition) and the chemical composition of the ambient gas.

We present 74 MHz radio continuum observations of the Galactic center region. These measurements show nonthermal radio emission arising from molecular clouds that is unaffected by free-free absorption along the line of sight. We focus on one cloud, G0.13--0.13, representative of the population of molecular clouds that are spatially correlated with steep spectrum (alpha^{74MHz}_{327MHz}=1.3\pm0.3) nonthermal emission from the Galactic center region. This cloud lies adjacent to the nonthermal radio filaments of the Arc near l~0.2^0 and is a strong source of 74 MHz continuum, SiO (2-1) and FeI Kalpha 6.4 keV line emission. This three-way correlation provides the most compelling evidence yet that relativistic electrons, here traced by 74 MHz emission, are physically associated with the G0.13--0.13 molecular cloud and that low energy cosmic ray electrons are responsible for the FeI Kalpha line emission. The high cosmic ray ionization rate ~10-13 s-1 H-1 is responsible for heating the molecular gas to high temperatures and allows the disturbed gas to maintain a high velocity dispersion. LVG modeling of multi-transition SiO observations of this cloud implies H2 densities ~104-5 cm-3 and high temperatures. The lower limit to the temperature of G0.13-0.13 is ~100K, whereas the upper limit is as high as 1000K. Lastly, we used a time-dependent chemical model in which cosmic rays drive the chemistry of the gas to investigate for molecular line diagnostics of cosmic ray heating. When the cloud reaches chemical equilibrium, the abundance ratios of HCN/HNC and N2H+/HCO+ are consistent with measured values. In addition, significant abundance of SiO is predicted in the cosmic ray dominated region of the Galactic center. We discuss different possibilities to account for the origin of widespread SiO emission detected from Galactic center molecular clouds.

We use 3D-PDR, a three-dimensional astrochemistry code for modeling photodissociation regions (PDRs), to post-process hydrodynamic simulations of turbulent, star-forming clouds. We focus on the transition from atomic to molecular gas, with specific attention to the formation and distribution of H, C+, C, H2 and CO. First, we demonstrate that the details of the cloud chemistry and our conclusions are insensitive to the simulation spatial resolution, to the resolution at the cloud edge, and to the ray angular resolution. We then investigate the effect of geometry and simulation parameters on chemical abundances and find weak dependence on cloud morphology as dictated by gravity and turbulent Mach number. For a uniform external radiation field, we find similar distributions to those derived using a one-dimensional PDR code. However, we demonstrate that a three-dimensional treatment is necessary for a spatially varying external field, and we caution against using one-dimensional treatments for non-symmetric problems. We compare our results with the work of Glover et al. (2010), who self-consistently followed the time evolution of molecule formation in hydrodynamic simulations using a reduced chemical network. In general, we find good agreement with this in situ approach for C and CO abundances. However, the temperature and H2 abundances are discrepant in the boundary regions (Av < 5), which is due to the different number of rays used by the two approaches.

Recent laboratory experiments have demonstrated that, even though contribution from other reaction channels cannot be neglected, unsaturated hydrocarbons easily break their multiple C-C bonds to form CO after their interactions with atomic oxygen. Here we present an upgraded chemical modelling including a revision of the reactions between oxygen atoms and small unsaturated hydrocarbons for different astrochemical environments. A first conclusion is that towards hot cores/corinos atomic oxygen easily degrades unsaturated hydrocarbons directly to CO or to its precursor species (such as HCCO or HCO) and destroys the double or triple bond of alkenes and alkynes. Therefore, environments rich in atomic oxygen at a relatively high temperature are not expected to be rich in large unsaturated hydrocarbons or polycyclic aromatic hydrocarbons. On the contrary, in O-poor and C-rich objects, hydrocarbon growth can occur to a large extent. In addition, new radical species, namely ketyl and vinoxy radicals, generated from other reaction channels can influence the abundances of hydrocarbons towards hot cores. We, therefore, suggest they should be included in the available databases. Hydrocarbon column densities are calculated in the 10(13)-10(15) cm(-2) range, in good agreement with their observed values, despite the small number of data currently published in the literature.

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 present a model for the formation of large organic molecules in dark clouds. The molecules are produced in the high density gas-phase that exists immediately after ice mantles are explosively sublimated. The explosions are initiated by the catastrophic recombination of trapped atomic hydrogen. We propose that, in molecular clouds, the processes of freeze-out onto ice mantles, accumulation of radicals, explosion and then rapid (three-body) gas-phase chemistry occurs in a cyclic fashion. This can lead to a cumulative molecular enrichment of the interstellar medium. A model of the time-dependent chemistries, based on this hypothesis, shows that significant abundances of large molecular species can be formed, although the complexity of the species is limited by the short expansion timescale in the gas, immediately following mantle explosion. We find that this mechanism may be an important source of smaller organic species, such as methanol and formaldehyde, as well as precursors to bio-molecule formation. Most significantly, we predict the gas-phase presence of these larger molecular species in quiescent molecular clouds and not just dynamically active regions, such as hot cores. As such the mechanism that we propose complements alternative methods of large molecule formation, such as those that invoke solid-state chemistry within activated ice mantles.

We report the discovery of a widespread population of collisionally excited methanol J = 4_{-1} to 3$_0 E sources at 36.2 GHz from the inner 66'x18' (160x43 pc) of the Galactic center. This spectral feature was imaged with a spectral resolution of ~16.6 km/s taken from 41 channels of a VLA continuum survey of the Galactic center region. The revelation of 356 methanol sources, most of which are maser candidates, suggests a large abundance of methanol in the gas phase in the Galactic center region. There is also spatial and kinematic correlation between SiO (2--1) and CH3OH emission from four Galactic center clouds: the +50 and +20 km/s clouds and G0.13-0.13 and G0.25+0.01. The enhanced abundance of methanol is accounted for in terms of induced photodesorption by cosmic rays as they travel through a molecular core, collide, dissociate, ionize, and excite Lyman Werner transitions of H2. A time-dependent chemical model in which cosmic rays drive the chemistry of the gas predicts CH3OH abundance of 10^{-8} to 10^{-7} on a chemical time scale of 5x10^4 to 5x10^5 years. The average methanol abundance produced by the release of methanol from grain surfaces is consistent with the available data.

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

Ultra-high spectral resolution observations of time-varying interstellar absorption towards {\kappa} Vel are reported, using the Ultra-High Resolution Facility on the Anglo-Australian Telescope. Detections of interstellar Ca I, Ca II, K I, Na I and CH are obtained, whilst an upper limit on the column density is reported for C_2. The results show continued increases in column densities of K I and Ca I since observations ~ 4 yr earlier, as the transverse motion of the star carried it ~ 10 AU perpendicular to the line of sight. Line profile models are fitted to the spectra and two main narrow components (A & B) are identified for all species except CH. The column density N(K I) is found to have increased by 82 +10-9 % between 1994 and 2006, whilst N(Ca I) is found to have increased by 32 +- 5 % over the shorter period of 2002-2006. The line widths are used to constrain the kinetic temperature to T_k,A < 671 +18-17 K and T_k,B < 114 +15-14 K. Electron densities are determined from the Ca I / Ca II ratio, which in turn place lower limits on the total number density of n_A > 7 * 10^3 cm^-3 and n_B > 2 * 10^4 cm^-3. Calcium depletions are estimated from the Ca I / K I ratio. Comparison with the chemical models of Bell et al. (2005) confirms the high number density, with n = 5 * 10^4 cm^-3 for the best-fitting model. The first measurements of diffuse interstellar bands (DIBs) towards this star are made at two epochs, but only an upper limit of < 40 % is placed on their variation over ~ 9 years. The DIBs are unusually weak for the measured E(B-V) and appear to exhibit similar behaviour to that seen in Orion. The ratio of equivalent widths of the {\lambda}5780 to {\lambda}5797 DIBs is amongst the highest known, which may indicate that the carrier of {\lambda}5797 is more sensitive to UV radiation than to local density.

We aimed to study the molecular composition of the interstellar medium (ISM) surrounding an Active Galactic Nucleus (AGN), by making an inventory of molecular species and their abundances, as well as to establish a chemical differentiation between starburst galaxies and AGN. We used the IRAM-30 m telescope to observe the central 1.5-2 kpc region of NGC1068, covering the frequencies between 86.2 GHz and 115.6 GHz. Using Boltzmann diagrams, we calculated the column densities of the detected molecules. We used a chemical model to reproduce the abundances found in the AGN, to determine the origin of each detected species, and to test the influence of UV fields, cosmic rays, and shocks on the ISM. We identified 24 different molecular species and isotopologues, among which HC3N, SO, N2H+, CH3CN, NS, 13CN, and HN13C are detected for the first time in NGC1068. We obtained the upper limits to the isotopic ratios 12C/13C=49, 16O/18O=177 and 32S/34S=5. Our chemical models suggest that the chemistry in the nucleus of NGC1068 is strongly influenced by cosmic rays, although high values of both cosmic rays and far ultraviolet (FUV) radiation fields also explain well the observations. The gas in the nucleus of NGC1068 has a different chemical composition as compared to starburst galaxies. The distinct physical processes dominating galaxy nuclei (e.g. C-shocks, UV fields, X-rays, cosmic rays) leave clear imprints in the chemistry of the gas, which allow to characterise the nucleus activity by its molecular abundances.

Despite its potential reactivity due to ring strain, ethylene oxide (c-C2H4O) is a complex molecule that seems to be stable under the physical conditions of an interstellar dense core; indeed it has been detected towards several high-mass star forming regions with a column density of the order of 10e13cm-2 (Ikeda et al. 2001). To date, its observational abundances cannot be reproduced by chemical models and this may be due to the significant contribution played by its chemistry on grain surfaces. Recently, Ward and Price (2011) have performed experiments in order to investigate the surface formation of ethylene oxide starting with oxygen atoms and ethylene ice as reactants. We present a chemical model which includes the most recent experimental results from Ward and Price (2011) on the formation of c-C2H4O. We study the influence of the physical parameters of dense cores on the abundances of c-C2H4O. We verify that ethylene oxide can indeed be formed during the cold phase (when the ISM dense cores are formed), via addition of an oxygen atom across the C=C double bond of the ethylene molecule, and released by thermal desorption during the hot core phase. A qualitative comparison between our theoretical results and those from the observations shows that we are able to reproduce the abundances of ethylene oxide towards high-mass star-forming regions.

We surveyed the circumnuclear disk of the Seyfert galaxy NGC1068 between the frequencies 86.2 GHz and 115.6 GHz, and identified 17 different molecules. Using a time and depth dependent chemical model we reproduced the observational results, and show that the column densities of most of the species are better reproduced if the molecular gas is heavily pervaded by a high cosmic ray ionization rate of about 1000 times that of the Milky Way. We discuss how molecules in the NGC1068 nucleus may be influenced by this external radiation, as well as by UV radiation fields.

We present for the first time a detailed study of the properties of molecular gas in metal-rich environments such as early-type galaxies (ETGs). We have explored Photon-Dominated Region (PDR) chemistry for a wide range of physical conditions likely to be appropriate for these sources. We derive fractional abundances of the 20 most chemically reactive species as a function of the metallicity, as a function of the optical depth and for various volume number gas densities, Far-Ultra Violet (FUV) radiation fields and cosmic ray ionisation rates. We also investigate the response of the chemistry to the changes in $\alpha-$element enhancement as seen in ETGs. We find that the fractional abundances of CS, H$_{2}$S, H$_{2}$CS, H$_{2}$O, H$_{3}$O$^{+}$, HCO$^{+}$ and H$_{2}$CN seem invariant to an increase of metallicity whereas C$^{+}$, CO, C$_{2}$H, CN, HCN, HNC and OCS appear to be the species most sensitive to this change. The most sensitive species to the change in the fractional abundance of $\alpha-$elements are C$^{+}$, C, CN, HCN, HNC, SO, SO$_{2}$, H$_{2}$O and CS. Finally, we provide line brightness ratios for the most abundant species, especially in the range observable with ALMA. Discussion of favorable line ratios to use for the estimation of super-solar metallicities and $\alpha$-elements are also provided.

We present high angular resolution observations (0.5"x0.3") carried out with the Submillimeter Array (SMA) toward the AFGL2591 high-mass star forming region. Our SMA images reveal a clear chemical segregation within the AFGL2591 VLA 3 hot core, where different molecular species (Type I, II and III) appear distributed in three concentric shells. This is the first time that such a chemical segregation is ever reported at linear scales <3000 AU within a hot core. While Type I species (H2S and 13CS) peak at the AFGL2591 VLA 3 protostar, Type II molecules (HC3N, OCS, SO and SO2) show a double-peaked structure circumventing the continuum peak. Type III species, represented by CH3OH, form a ring-like structure surrounding the continuum emission. The excitation temperatures of SO2, HC3N and CH3OH (185+-11 K, 150+-20 K and 124+-12 K, respectively) show a temperature gradient within the AFGL2591 VLA 3 envelope, consistent with previous observations and modeling of the source. By combining the H2S, SO2 and CH3OH images, representative of the three concentric shells, we find that the global kinematics of the molecular gas follow Keplerian-like rotation around a 40 Mo-star. The chemical segregation observed toward AFGL2591 VLA 3 is explained by the combination of molecular UV photo-dissociation and a high-temperature (~1000 K) gas-phase chemistry within the low extinction innermost region in the AFGL2591 VLA 3 hot core.

The Gould Belt Legacy Survey will map star-forming regions within 500 pc, using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common-User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes HARP observations of the J = 3-2 transitions of 13CO and C18O towards Orion A. The 1500-resolution observations cover 5 pc of the Orion filament, including OMC1 (inc. BN-KL and Orion Bar), OMC 2/3 and OMC 4, and allow a comparative study of the molecular gas properties throughout the star-forming cloud. The filament shows a velocity gradient of ~1 km/s /pc between OMC 1, 2 and 3, and high velocity emission is detected in both isotopologues. The Orion Nebula and Bar have the largest masses and line widths, and dominate the mass and energetics of the high velocity material. Compact, spatially resolved emission from CH3CN, 13CH3OH, SO, HCOOCH3, C2H5OH, CH3CHO and CH3OCHO is detected towards the Orion Hot Core. The cloud is warm, with a median excitation temperature of ~24 K; the Orion Bar has the highest excitation temperature gas, at >80 K. The C18O excitation temperature correlates well with the dust temperature (to within 40%). The C18O emission is optically thin, and the 13CO emission is marginally optically thick; despite its high mass, OMC 1 shows the lowest opacities. A virial analysis indicates that Orion A is too massive for thermal or turbulent support, but is consistent with a model of a filamentary cloud that is threaded by helical magnetic fields. The variation of physical conditions across the cloud is reflected in the physical characteristics of the dust cores....continued

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

Glycolaldehyde is a simple monosaccharide sugar linked to prebiotic chemistry. Recently it was detected in a molecular core in the star-forming region G31.41+0.31 at a reasonably high abundance. We investigate the formation of glycolaldehyde at 10K to determine whether it can form efficiently under typical dense core conditions. Using an astrochemical model, we test five different reaction mechanisms that have been proposed in the astrophysical literature, finding that a gas-phase formation route is unlikely. Of the grain-surface formation routes, only two are efficient enough at very low temperatures to produce sufficient glycolaldehyde to match the observational estimates, with the mechanism culminating in CH3OH + HCO being favoured. However, when we consider the feasibility of these mechanisms from a reaction chemistry perspective, the second grain-surface route looks more promising, H3CO + HCO.

We aimed to study the chemistry of the circumnuclear molecular gas of NGC1068, and to compare it with those of the starburst galaxies M82 and NGC253. Using the IRAM-30m telescope, we observed the inner 2 kpc of NGC1068 between 86.2 GHz and 115.6 GHz. We identified 35 spectral features, corresponding to 24 different molecular species. Among them, HC3N, SO, N2H+, CH3CN, NS, 13CN, and HN13C are detected for the first time in NGC1068. Assuming local thermodynamic equilibrium (LTE), we calculated the column densities of the detected molecules, as well as the upper limits to the column densities of some undetected species. The comparison among the chemistries of NGC1068, M82, and NGC253, suggests that, apart from X-rays, shocks also determine the chemistry of NGC1068. We propose the column density ratio between CH3CCH and HC3N as a prime indicator of the imprints of starburst and AGN environments in the circumnuclear interstellar medium. This ratio is, at least, 64 times larger in M82 than in NGC1068, and, at least, 4 times larger in NGC253 than in NGC1068. Finally, we used the UCL_CHEM and UCL_PDR chemical codes to constrain the origin of the species, as well as to test the influence of UV radiation fields and cosmic rays on the observed abundances.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

We present the first detection of hydrogen chlorine in a protostellar shock, by observing the fundamental transition at 626 GHz with the Herschel HIFI spectrometer. We detected two of the three hyperfine lines, from which we derived a line opacity < 1. Using a non-LTE LVG code, we constrained the HCl column density, temperature and density of the emitting gas. The hypothesis that the emission originates in the molecular cloud is ruled out, as it would imply a too dense gas. Conversely, assuming that the emission originates in the 10"-15" size shocked gas previously observed at the IRAM PdB interferometer, we obtain: N(HCl)=0.7-2 x 10(13) cm-2, temperature > 15 K and density > 3 x 10(5) cm-3}. Combining with the Herschel HIFI CO(5-4) observations allows to further constrain the gas density and temperature, 10(5)-10(6) cm-3 and 120-250 K, as well as the HCl column density, 2 x 10(13) cm-2, and, finally, abundance: 3-6 x 10(-9). The estimated HCl abundance is consistent with that previously observed in low- and high- mass protostars. This puzzling result in the L1157-B1 shock, where species from volatile and refractory grains components are enhanced, suggests either that HCl is not the main reservoir of chlorine in the gas phase, against previous chemical models predictions, or that the elemental chlorine abundance is low in L1157-B1. Astrochemical modelling suggests that HCl is in fact formed in the gas phase, at low temperatures, prior to the occurance of the shock, and that the latter does not enhance its abundance.

We investigate the origin and nature of the profiles of water and ammonia observ ed toward the L1157-B1 clump as part of the HIFI CHESS survey (Ceccarelli et al. 2010) using a new code coupling a gas-grain chemical model with a parametric shock model. Fir st results from the unbiased survey (Lefloch et al. 2010, Codella et al. 2010) reveal different molecular components at different excitation conditions coexisting in the B1 bow shock structure, with NH$_3$, H$_2$CO and CH$_3$OH emitting only at relatively low outflow velocities whereas H$_2$O shows bright emission at high velocities. Our model suggests that these differences are purely chemical and can be explained by the presence of a C-type shock whose maximum temperature must be close to 4000 K along the B1 clump.

The molecular clouds Lupus 1, 3 and 4 were mapped with the Mopra telescope at 3 and 12 mm. Emission lines from high density molecular tracers were detected, i.e. NH$_3$ (1,1), NH$_3$ (2,2), N$_2$H$^+$ (1-0), HC$_3$N (3-2), HC$_3$N (10-9), CS (2-1), CH$_3$OH (2$_0-1_0$)A$^+$ and CH$_3$OH (2$_{-1}-1_{-1}$)E. Velocity gradients of more than 1 km s$^{-1}$ are present in Lupus 1 and 3 and multiple gas components are present in these clouds along some lines of sight. Lupus 1 is the cloud richest in high density cores, 8 cores were detected in it, 5 cores were detected in Lupus 3 and only 2 in Lupus 4. The intensity of the three species HC$_3$N, NH$_3$ and N$_2$H$^+$ changes significantly in the various cores: cores that are brighter in HC$_3$N are fainter or undetected in NH$_3$ and N$_2$H$^+$ and vice versa. We found that the column density ratios HC$_3$N/N$_2$H$^+$ and HC$_3$N/NH$_3$ change by one order of magnitude between the cores, indicating that also the chemical abundance of these species is different. The time dependent chemical code that we used to model our cores shows that the HC$_3$N/N$_2$H$^+$ and HC$_3$N/NH$_3$ ratios decrease with time therefore the observed column density of these species can be used as an indicator of the chemical evolution of dense cores. On this base we classified 5 out of 8 cores in Lupus 1 and 1 out of 5 cores in Lupus 3 as very young protostars or prestellar cores. Comparing the millimetre cores population with the population of the more evolved young stellar objects identified in the Spitzer surveys, we conclude that in Lupus 3 the bulk of the star formation activity has already passed and only a moderate number of stars are still forming. On the contrary, in Lupus 1 star formation is on-going and several dense cores are still in the pre--/proto--stellar phase. Lupus 4 is at an intermediate stage, with a smaller number of individual objects.

Localised regions of enhanced emission from HCO+, NH3 and other species near Herbig-Haro objects (HHOs) have been interpreted as arising in a photochemistry stimulated by the HHO radiation on high density quiescent clumps in molecular clouds. Static models of this process have been successful in accounting for the variety of molecular species arising ahead of the jet; however recent observations show that the enhanced molecular emission is widespread along the jet as well as ahead. Hence, a realistic model must take into account the movement of the radiation field past the clump. It was previously unclear as to whether the short interaction time between the clump and the HHO in a moving source model would allow molecules such as HCO+ to reach high enough levels, and to survive for long enough to be observed. In this work we model a moving radiation source that approaches and passes a clump. The chemical picture is qualitatively unchanged by the addition of the moving source, strengthening the idea that enhancements are due to evaporation of molecules from dust grains. In addition, in the case of several molecules, the enhanced emission regions are longer-lived. Some photochemically-induced species, including methanol, are expected to maintain high abundances for ~10,000 years.

Several complex organic molecules are routinely detected in high abundances towards hot cores and hot corinos. For many of them, their paths of formation in space are uncertain, as gas phase reactions alone seem to be insufficient. In this paper, we investigate a possible solid-phase route of formation for methyl formate (HCOOCH3). We use a chemical model updated with recent results from an experiment where simulated grain surfaces were irradiated with 200 keV protons at 16 K, to simulate the effects of cosmic ray irradiation on grain surfaces. We find that this model may be sufficient to reproduce the observed methyl formate in dark clouds, but not that found in hot cores and corinos.

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

We present the first observations of emission lines of CN(2-1), HCO$^{+}$(3-2) and C$_{2}$H(3-2) in the Perseus cluster. We observed at two positions: directly at the central galaxy, NGC 1275 and also at a position about 20$''$ to the east where associated filamentary structure has been shown to have strong CO emission. Clear detections in CN and HCO$^{+}$ transitions and a weak detection of the C$_{2}$H transition were made towards NGC 1275, while weak detections of CN and HCO$^{+}$ were made towards the eastern filamentary structure. Crude estimates of the column densities and fractional abundances (mostly upper limits) as functions of an unknown rotational temperature were made to both sources. These observational data were compared with the outputs of thermal/chemical models previously published by \citet{Baye10c} in an attempt to constrain the heating mechanisms in cluster gas. We find that models in which heating is dominated by cosmic rays can account for the molecular observations. This conclusion is consistent with that of \citet{Ferl09} in their study of gas traced by optical and infrared radiation. The cosmic ray heating rate in the regions probed by molecular emissions is required to be at least two orders of magnitude larger than that in the Milky Way.

Molecular line observations may serve as diagnostics of the degree to which the number density of cosmic ray protons, having energies of 10s to 100s of MeVs each, is enhanced in starburst galaxies and galaxies with active nuclei. Results, obtained with the UCL\_PDR code, for the fractional abundances of molecules as functions of the cosmic-ray induced ionisation rate, $\zeta$, are presented. The aim is not to model any particular external galaxies. Rather, it is to identify characteristics of the dependencies of molecular abundances on $\zeta$, in part to enable the development of suitable observational programmes for cosmic ray dominated regions (CRDRs) which will then stimulate detailed modelling. For a number density of hydrogen nuclei of of $10^4$ cm$^{-3}$, and high visual extinction, the fractional abundances of some species increase as $\zeta$ increases to $10^{-14}$ s$^{-1}$, but for much higher values of $\zeta$ the fractional abundances of all molecular species are significantly below their peak values. We show in particular that OH, H$_{2}$O, H$_{3}^{+}$, H$_{3}$O$^{+}$ and OH$^{+}$ attain large fractional abundances ($\geqslant 10^{-8}$) for $\zeta$ as large as $10^{-12}$ s$^{-1}$. HCO$^{+}$ is a poor tracer of CRDRs when $\zeta > 10^{-13}$ s$^{-1}$. Sulphur-bearing species may be useful tracers of CRDRs gas in which $\zeta \sim 10^{-16}$ s$^{-1}$. Ammonia has a large fractional abundance for $\zeta \leqslant 10^{-16}$ s$^{-1}$ and nitrogen appears in CN-bearing species at significant levels as $\zeta$ increases, even up to $\sim 10^{-14}$ s$^{-1}$. In this paper, we also discuss our model predictions, comparing them to recent detections in both galactic and extragalactic sources. We show that they agree well, to a first approximation, with the observational constraints.

We examine in detail the recent proposal that extreme Cosmic-Ray-Dominated-Regions (CRDRs) characterize the ISM of galaxies during events of high-density star formation, fundamentally altering its initial conditions (Papadopoulos 2010). Solving the coupled chemical and thermal state equations for dense UV-shielded gas reveals that the large cosmic ray energy densities in such systems (U_{CR} (few)x(10^3-10^4) U_{CR,Gal}) will indeed raise the minimum temperature of this phase (where the initial conditions of star formation are set) from ~10K (as in the Milky Way) to (50-100)K. Moreover in such extreme CRDRs the gas temperature remains fully decoupled from that of the dust, with T_{kin} >> T_{dust}, even at high densities (n(H_2)~10^5--10^6 cm^{-3}), quite unlike CRDRs in the Milky Way where T_k T_{dust} when n(H_2) >= 10^5 cm^{-3}. These dramatically different star formation initial conditions will: a) boost the Jeans mass of UV-shielded gas regions by factors of ~10--100 with respect to those in quiescent or less extreme star forming systems, and b) "erase" the so-called inflection point of the effective equation of state (EOS) of molecular gas. Both these effects occur across the entire density range of typical molecular clouds, and may represent {\it a new paradigm for all high-density star formation in the Universe}, with cosmic rays as the key driving mechanism, operating efficiently even in the high dust extinction environments of extreme starbursts...

We present a study of the Silicon-bearing species detected in a line confusion limited survey towards Orion KL performed with the IRAM 30-m telescope, which ranges from 80 to 280 GHz. Our aim is to derive physical and chemical conditions for each family taking into account all observed lines from all isotopologues of each species. We have modeled the lines of the detected molecules using a radiative transfer code, which permit us to choose between LVG and LTE approximations depending on the physical conditions of the gas. We have used appropriate collisional rates for the LVG calculations. For the v=1 state of SiO we have detected the J=2-1 line and, for the first time in this source, emission in the J=4-3 transition, both of them showing strong masering effect. For SiO v=0, we have detected 28SiO, 29SiO, and 30SiO; in addition, we have mapped the J = 5-4 SiO line. For SiS, we have detected the main species, 29SiS, and SiS v=1. Unlikely other species detected in Orion KL (IRc2), the emission peak of SiS appears at a velocity of ~15.5 km s-1; a study of the 5-4 SiO line around IRc2 shows this feature as an extended component that probably arises from the interaction of the outflow with the ambient cloud. We derive a SiO/SiS column density ratio of ~13 in the plateau component. Besides, we provide upper limits to the column density of several non-detected Silicon-bearing species. The results of our chemical models show that while it is possible to reproduce SiO in the gas phase (as well as on the grains), SiS is a product of surface reactions, most likely involving direct reactions of S with Si.

We suggest that three-body chemistry may occur in warm high density gas evaporating in transient co\textendash desorption events on interstellar ices. Using a highly idealised computational model we explore the chemical conversion from simple species of the ice to more complex species containing several heavy atoms, as a function of density and of adopted three body rate coefficients. We predict that there is a wide range of densities and rate coefficients in which a significant chemical conversion may occur. We discuss the implications of this idea for the astrochemistry of hot cores.

We present a theoretical study of CS line profiles in archetypal hot cores. We provide estimates of line fluxes from the CS(1-0) to the CS(15-14) transitions and present the temporal variation of these fluxes. We find that \textit{i)} the CS(1-0) transition is a better tracer of the Envelope of the hot core whereas the higher-J CS lines trace the ultra-compact core; \textit{ii)} the peak temperature of the CS transitions is a good indicator of the temperature inside the hot core; \textit{iii)} in the Envelope, the older the hot core the stronger the self-absorption of CS; \textit{iv)} the fractional abundance of CS is highest in the innermost parts of the ultra-compact core, confirming the CS molecule as one of the best tracers of very dense gas.

We present a theoretical study of the deuterated species detectability in various types of extragalactic star-forming regions based on our predictions of chemical abundances. This work is motivated by the past and current attempts at observing deuterated species in external galaxies such as NGC~253, IC~342 and the LMC. Here, we investigate the influence of the density, the temperature, the FUV radiation field, the cosmic ray ionisation, and the metallicity on the fractional abundances and D/H abundance ratios of about 20 deuterated species. Without modelling any particular source, we determined how the deuterium chemistry behaves in different physical environments such as starburst, cosmic-rays enhanced environments, low metallicity and high redshift galaxies. In general, our predicted column densities seem in good agreement with those derived from the current limited dataset of observations in external galaxies. We provide, for the first time, a list of key deuterated species whose abundances are high enough to be possibly detectable by the Atacama Large Millimeter Array (ALMA) and Herschel, as a function of galactic nuclear activity and redshift.

Many theoretical and laboratory studies predict H2 to be formed in highly excited ro-vibrational states. The consequent relaxation of excited levels via a cascade of infrared transitions might be observable in emission from suitable interstellar regions. In this work, we model H2 formation pumping in standard dense clouds, taking into account the H/H2 transition zone, through an accurate description of chemistry and radiative transfer. The model includes recent laboratory data on H2 formation, as well as the effects of the interstellar UV field, predicting the populations of gas-phase H2 molecules and their IR emission spectra. Calculations suggest that some vibrationally excited states of H2 might be detectable towards lines of sight where significant destruction of H2 occurs, such as X-ray sources, and provide a possible explanation as to why observational attempts resulted in no detections reported to date.

We aim at investigating with high angular resolution the NH3/N2H+ ratio toward the high-mass star-forming region AFGL 5142 in order to study whether this ratio behaves similarly to the low-mass case, for which the ratio decreases from starless cores to cores associated with YSOs. CARMA was used to observe the 3.2 mm continuum and N2H+(1-0) emission. We used NH3(1,1) and (2,2), HCO+(1-0) and H13CO+(1-0) data from the literature and we performed a time-dependent chemical modeling of the region. The 3.2 mm continuum emission reveals a dust condensation of ~23 Msun associated with the massive YSOs, deeply embedded in the strongest NH3 core (hereafter central core). The N2H+ emission reveals two main cores, the western and eastern core, located to the west and to the east of the mm condensation, and surrounded by a more extended and complex structure of ~0.5 pc. Toward the central core the N2H+ emission drops significantly, indicating a clear chemical differentiation in the region. We found low values of the NH3/N2H+ ratio ~50-100 toward the western/eastern cores, and high values up to 1000 in the central core. The chemical model indicates that density, and in particular temperature, are key parameters in determining the NH3/N2H+ ratio. The high density and temperature reached in the central core allow molecules like CO to evaporate from grain mantles. The CO desorption causes a significant destruction of N2H+, favoring the formation of HCO+. This result is supported by our observations, which show that N2H+ and HCO+ are anticorrelated in the central core. The observed values of the NH3/N2H+ ratio in the central core can be reproduced by our model for times t~4.5-5.3x10^5 yr (central) and t~10^4-3x10^6 yr (western/eastern). The NH3/N2H+ ratio in AFGL 5142 does not follow the same trend as in regions of low-mass star formation mainly due to the high temperature reached in hot cores.

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.

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.

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

Optical emission is detected from filaments around the central galaxies of clusters of galaxies. These filaments have lengths of tens of kiloparsecs. The emission is possibly due to heating caused by the dissipation of mechanical energy and by cosmic ray induced ionisation. CO millimeter and submillimeter line emissions as well as H$_{2}$ infrared emission originating in such filaments surrounding NGC~1275, the central galaxy of the Perseus cluster, have been detected. Our aim is to identify those molecular species, other than CO, that may emit detectable millimeter and submillimeter line features arising in these filaments, and to determine which of those species will produce emissions that might serve as diagnostics of the dissipation and cosmic ray induced ionisation. The time-dependent UCL photon-dominated region modelling code was used in the construction of steady-state models of molecular filamentary emission regions at appropriate pressures, for a range of dissipation and cosmic ray induced ionisation rates and incident radiation fields.HCO$^+$ and C$_2$H emissions will potentially provide information about the cosmic ray induced ionisation rates in the filaments. HCN and, in particular, CN are species with millimeter and submillimeter lines that remain abundant in the warmest regions containing molecules. Detections of the galaxy cluster filaments in HCO$^{+}$, C$_{2}$H, and CN emissions and further detections of them in HCN emissions would provide significant constraints on the dissipation and cosmic ray induced ionisation rates.

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.

The deuterium fractionation, Dfrac, has been proposed as an evolutionary indicator in pre-protostellar and protostellar cores of low-mass star-forming regions. We investigate Dfrac, with high angular resolution, in the cluster environment surrounding the UCHII region IRAS 20293+3952. We performed high angular resolution observations with the IRAM Plateau de Bure Interferometer (PdBI) of the ortho-NH2D 1_{11}-1_{01} line at 85.926 GHz and compared them with previously reported VLA NH3 data. We detected strong NH2D emission toward the pre-protostellar cores identified in NH3 and dust emission, all located in the vicinity of the UCHII region IRAS 20293+3952. We found high values of Dfrac~0.1-0.8 in all the pre-protostellar cores and low values, Dfrac<0.1, associated with young stellar objects. The high values of Dfrac in pre-protostellar cores could be indicative of evolution, although outflow interactions and UV radiation could also play a role.

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.

The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region

Warm cores (or hot corinos) around low-mass protostellar objects show a rich chemistry with strong spatial variations. This chemistry is generally attributed to the sublimation of icy mantles on dust grains initiated by the warming effect of the stellar radiation. We have used a model of the chemistry in warm cores in which the sublimation process is based on extensive laboratory data; these data indicate that sublimation from mixed ices occurs in several well-defined temperature bands. We have determined the position of these bands for the slow warming by a solar-mass star. The resulting chemistry is dominated by the sublimation process and by subsequent gas-phase reactions; strong spatial and temporal variations in certain molecular species are found to occur, and our results are, in general, consistent with observational results for the well-studied source IRAS 16293-2422. The model used is similar to one that describes the chemistry of hot cores. We infer that the chemistry of both hot cores and warm cores may be described by the same model (suitably adjusted for different physical parameters).

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.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

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.

Shock modelling predicts an electron density enhancement within the magnetic precursor of C-shocks. Previous observations of SiO, H13CO+, HN13C and H13CN toward the young L1448-mm outflow showed an over-excitation of the ion fluid that was attributed to an electron density enhancement in the precursor. We re-visit this interpretation and test if it still holds when we consider different source morphologies and kinetic temperatures for the observed molecules, and also give some insight on the spatial extent of the electron density enhancement around L1448-mm. We estimate the opacities of H13CO+ and HN13C by observing the J=3\to2 lines of rarer isotopologues to confirm that the emission is optically thin. To model the excitation of the molecules, we use the large velocity gradient (LVG) approximation with updated collisional coefficients to i) re- analyse the observations toward the positions where the over-excitation of H13CO+ has previously been observed [i.e. toward L1448- mm at offsets (0,0) and (0,-10)], and ii) to investigate if the electron density enhancement is still required for the cases of extended and compact emission, and for kinetic temperatures of up to 400 K. We also report several lines of SiO, HN13C and H13CO+ toward new positions around this outflow, to investigate the spatial extent of the over-excitation of the ions in L1448-mm. From the isotopologue observations, we find that the emission of H13CO+ and HN13C from the precursor is optically thin if this emission is extended. Using the new collisional coefficients, an electron density enhancement is still needed to explain the excitation of H13CO+ for extended emission and for gas temperatures of\le 400 K toward L1448-mm (0,-10), and possibly also toward L1448-mm (0,0). For compact emission the data cannot be fitted. We do not find any evidence for the over-excitation of the ion fluid toward the newly observed positions around L1448-mm. The observed line emission of SiO, H13CO+ and HN13C toward L1448-mm (0,0) and (0,-10) is consistent with an electron density enhancement in the precursor component, if this emission is spatially extended. This is also true for the case of high gas temperatures (\le400 K) toward the (0,-10) offset. The electron density enhancement seems to be restricted to the southern, redshifted lobe of the L1448-mm outflow. Interferometric images of the line emission of these molecules are needed to confirm the spatial extent of the over-excitation of the ions and thus, of the electron density enhancement in the magnetic precursor of L1448-mm.

We present chemical and radiative transfer models for the many far-IR ortho- and para-H2O lines that were observed from the Orion-KL region in high resolution Fabry-Perot (FP) mode by the Long Wavelength Spectrometer (LWS) on board the Infrared Space Observatory (ISO). The chemistry of the region was first studied by simulating the conditions in the different known components of Orion-KL: chemical models for a hot core, a plateau and a ridge were coupled with an accelerated Lambda$-iteration (ALI) radiative transfer model to predict H2O line fluxes and profiles. Our models include the first 45 energy levels of ortho- and para-H2O. We find that lines arising from energy levels below 560 K were best reproduced by a gas of density 3x10^5 cm^-3 at a temperature of 70-90 K, expanding at a velocity of 30 km s^-1 and with a H2O/H2 abundance ratio of the order of 2 - 3 x 10^-5, similar to the abundance derived by Cernicharo et al. (2006). However, the model that best reproduces the fluxes and profiles of H2O lines arising from energy levels above 560 K has a significantly higher H2O/H2 abundance, 1 - 5 x 10^-4, originating from gas of similar density, in the Plateau region, that has been heated to 300 K, relaxing to 90-100 K. We conclude that the observed water lines do not originate from high temperature shocks.

As part of a JCMT Legacy Survey of star formation in the Gould Belt, we present early science results for Taurus. CO J=3-2 maps have been secured along the north-west ridge and bowl, collectively known as L 1495, along with deep 13CO and C18O J=3-2 maps in two sub-regions. With these data we search for molecular outflows, and use the distribution of flows, HH objects and shocked H2 line emission features, together with the population of young stars, protostellar cores and starless condensations to map star formation across this extensive region. In total 21 outflows are identified. It is clear that the bowl is more evolved than the ridge, harbouring a greater population of T Tauri stars and a more diffuse, more turbulent ambient medium. By comparison, the ridge contains a much younger, less widely distributed population of protostars which, in turn, is associated with a greater number of molecular outflows. We estimate the ratio of the numbers of prestellar to protostellar cores in L 1495 to be ~ 1.3-2.3, and of gravitationally unbound starless cores to (gravitationally bound) prestellar cores to be ~ 1. If we take previous estimates of the protostellar lifetime of ~ 5 x 10^5 yrs, this indicates a prestellar lifetime of 9(+/-3) x 10^5 yrs. From the number of outflows we also crudely estimate the star formation efficiency in L 1495, finding it to be compatible with a canonical value of 10-15 %. We note that molecular outflow-driving sources have redder near-IR colours than their HH jet-driving counterparts. We also find that the smaller, denser cores are associated with the more massive outflows, as one might expect if mass build-up in the flow increases with the collapse and contraction of the protostellar envelope.

Far-infrared [C ii] 158 micron, [O i] 145 micron and [O i] 63 micron fine structure emission line fluxes were measured from archival Infrared Space Observatory Long Wavelength Spectrometer spectra of 46 extragalactic sources, with 28 sources providing detections in all three lines. For 12 of the sources, the contribution to the [C ii] 158 micron line flux from H ii regions could be estimated from their detected [N ii] 122 micron line fluxes. The measured [C ii]/[O i] and [O i] 63/145 micron line flux ratios were compared with those from a grid of PDR models previously computed using the UCL PDR code. Persistent offsets between the observed and modelled line ratios could be partly attributed to the effects of [O i] 63 micron self-absorption. Using the SMMOL code, we calculated model [O i] line profiles and found that the strength of the [O i] 63 micron line was reduced by 20-80%, depending on the PDR parameters. We conclude that high PDR densities and radiation field strengths, coupled with the effects of [O i] 63 micron self-absorption, are likely to provide the best match to the observed line flux ratios.

We report the first detection of PDR molecular tracers, namely HOC+, and CO+, and confirm the detection of the also PDR tracer HCO towards the starburst galaxy NGC 253, claimed to be mainly dominated by shock heating and in an earlier stage of evolution than M 82, the prototypical extragalactic PDR. Our CO+ detection suffers from significant blending to a group of transitions of 13CH3OH, tentatively detected for the first time in the extragalactic interstellar medium. These species are efficiently formed in the highly UV irradiated outer layers of molecular clouds, as observed in the late stage nuclear starburst in M 82. The molecular abundance ratios we derive for these molecules are very similar to those found in M 82. This strongly supports the idea that these molecules are tracing the PDR component associated with the starburst in the nuclear region of NGC 253. A comparison with the predictions of chemical models for PDRs shows that the observed molecular ratios are tracing the outer layers of UV illuminated clouds up to two magnitudes of visual extinction. Chemical models, which include grain formation and photodissociation of HNCO, support the scenario of a photo-dominated chemistry as an explanation to the abundances of the observed species. From this comparison we conclude that the molecular clouds in NGC 253 are more massive and with larger column densities than those in M 82, as expected from the evolutionary stage of the starbursts in both galaxies.

Some well-studied Herbig Haro objects have associated with them one or more cold, dense, and quiescent clumps of gas. We propose that such clumps near an HH object can be used as a general measure of clumpiness in the molecular cloud that contains that HH object. Our aim is to make a survey of clumps around a sample of HH objects, and to use the results to make an estimate of the clumpiness in molecular clouds. All known cold, dense, and quiescent clumps near HH objects are anomalously strong HCO+ emitters. Our method is, therefore, to search for strong HCO+ emission as an indicator of a clump near to an HH object. The searches were made using JCMT and SEST in the HCO+ 3-2 and also H13CO+ 1-0 lines, with some additional searches for methanol and sulphur monoxide lines. The sources selected were a sample of 22 HH objects in which no previous HCO+ emission had been detected. We find that half of the HH objects have clumps detected in the HCO+ 3-2 line and that all searches in H13CO$+ 1-0 lines show evidence of clumpiness. All condensations have narrow linewidths and are evidently unaffected dynamically by the HH jet shock. We conclude that the molecular clouds in which these HH objects are found must be highly heterogeneous on scales of less than 0.1 pc. An approximate calculation based on these results suggests that the area filling factor of clumps affected by HH objects is on the order of 10%. These clumps have gas number densities larger than 3e4 cm-2.

We present a coherent and homogeneous multi-line study of the CS molecule in nearby (D$<$10Mpc) galaxies. We include, from the literature, all the available observations from the $J=1-0$ to the $J=7-6$ transitions towards NGC 253, NGC 1068, IC 342, Henize~2-10, M~82, the Antennae Galaxies and M~83. We have, for the first time, detected the CS(7-6) line in NGC 253, M~82 (both in the North-East and South-West molecular lobes), NGC 4038, M~83 and tentatively in NGC 1068, IC 342 and Henize~2-10. We use the CS molecule as a tracer of the densest gas component of the ISM in extragalactic star-forming regions, following previous theoretical and observational studies by Bayet et al. (2008a,b and 2009). In this first paper out of a series, we analyze the CS data sample under both Local Thermodynamical Equilibrium (LTE) and non-LTE (Large Velocity Gradient-LVG) approximations. We show that except for M~83 and Overlap (a shifted gas-rich position from the nucleus NGC 4039 in the Antennae Galaxies), the observations in NGC 253, IC 342, M~82-NE, M~82-SW and NGC 4038 are not well reproduced by a single set of gas component properties and that, at least, two gas components are required. For each gas component, we provide estimates of the corresponding kinetic temperature, total CS column density and gas density.

The Gould Belt Legacy Survey will survey nearby star-forming regions (within 500 pc), using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common- User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes the initial data obtained using HARP to observe 12CO, 13CO and C18O J = 3 - 2 towards two regions in Orion B, NGC 2024 and NGC 2071. We describe the physical characteristics of the two clouds, calculating temperatures and opacities utilizing all three isotopologues. We find good agreement between temperatures calculated from CO and from dust emission in the dense, energetic regions. We determine the mass and energetics of the clouds, and of the high-velocity material seen in 12CO emission, and compare the relative energetics of the high- and low-velocity material in the two clouds. We present a CLUMPFIND analysis of the 13CO condensations. The slope of the condensation mass functions, at the high-mass ends, is similar to the slope of the initial mass function.

Banerji et al. (2009) suggested that top-heavy stellar Initial Mass Functions (IMFs) in galaxies may arise when the interstellar physical conditions inhibit low-mass star formation, and they determined the physical conditions under which this suppression may or may not occur. In this work, we explore the sensitivity of the chemistry of interstellar gas under a wide range of conditions. We use these results to predict the relative velocity-integrated antenna temperatures of the CO rotational spectrum for several models of high redshift active galaxies which may produce both top-heavy and unbiased IMFs. We find that while active galaxies with solar metallicity (and top-heavy IMFs) produce higher antenna temperatures than those with sub-solar metallicity (and unbiased IMFs) the actual rotational distribution is similar. The high-J to peak CO ratio however may be used to roughly infer the metallicity of a galaxy provided we know whether it is active or quiescent. The metallicity strongly influences the shape of the IMF. High order CO transitions are also found to provide a good diagnostic for high far-UV intensity and low metallicity counterparts of Milky Way type systems both of which show some evidence for having top-heavy IMFs. We also compute the relative abundances of molecules known to be effective tracers of high density gas in these galaxy models. We find that the molecules CO and CS may be used to distinguish between solar and sub-solar metallicity in active galaxies at high redshift whereas HCN, HNC and CN are found to be relatively insensitive to the IMF shape at the large visual magnitudes typically associated with extragalactic sources.

We present a deep and wide field-of-view (4'x 7') image of the planetary nebula (PN) NGC 7293 (the Helix Nebula) in the 2.12 micron H2 v=1-0 S(1) line. The excellent seeing (0.4'') at the Subaru Telescope, allows the details of cometary knots to be examined. The knots are found at distances of 2.2'-6.4' from the central star (CS). At the inner edge and in the inner ring (up to 4.5' fromthe CS), the knot often show a `tadpole' shape, an elliptical head with a bright crescent inside and a long tail opposite to the CS. In detail, there are variations in the tadpole shapes, such as narrowing tails, widening tails, meandering tails, or multi-peaks within a tail. In the outer ring (4.5'-6.4' from the CS), the shapes are more fractured, and the tails do not collimate into a single direction. The transition in knot morphology from the inner edge to the outer ring is clearly seen. The number density of knots governs the H2 surface brightness in the inner ring: H2 exists only within the knots. Possible mechanisms which contribute to the shaping of the knots are discussed, including photo-ionization and streaming motions. A plausible interpretation of our images is that inner knots are being overrun by a faster wind, but that this has not (yet) reached the outer knots. Based on H2 formation and destruction rates, H2 gas can survive in knots from formation during the late asymptotic giant branch (AGB) phase throughout the PN phase. These observations provide new constraints on the formation and evolution of knots, and on the physics of molecular gas embedded within ionized gas.

Photon-dominated regions (PDRs) are powerful molecular line emitters in external galaxies. They are expected in galaxies with high rates of massive star formation due to either starburst (SB) events or starburst coupled with active galactic nuclei (AGN) events. We have explored the PDR chemistry for a range of physical conditions representing a variety of galaxy types. Our main result is a demonstration of the sensitivity of the chemistry to changes in the physical conditions. We adopt crude estimates of relevant physical parameters for several galaxy types and use our models to predict suitable molecular tracers of those conditions. The set of recommended molecular tracers differs from that which we recommended for use in galaxies with embedded massive stars. Thus, molecular observations can in principle be used to distinguish between excitation by starburst and by SB+AGN in distant galaxies. Our recommendations are intended to be useful in preparing Herschel and ALMA proposals to identify sources of excitation in galaxies.

We present a Berkeley-Illinois-Maryland Association (BIMA) Array molecular survey of the star forming core ahead of HH 80N, the optically obscured northern counterpart of the Herbig-Haro objects HH 80/81. Continuum emission at 1.4 mm and 8 micron is detected at the center of the core, which confirms the presence of an embedded very young stellar object in the core. All detected molecular species arise in a ring-like structure, which is most clearly traced by CS(2-1) emission. This molecular ring suggests that strong molecular depletion occurs in the inner part of the core (at a radius of ~0.1 pc and densities higher than ~5 x 10^4 cm^{-3}). Despite of the overall morphology and kinematic similarity between the different species, there is significant molecular differentiation along the ring-like structure. The analysis of the chemistry along the core shows that part of this differentiation may be caused by the UV irradiation of the nearby HH 80N object, that illuminates the part of the core facing HH 80N, which results in an abundance enhancement of some of the detected species.

We present the first results of the expected variations of the molecular line emission arising from material recently affected by C-shocks (shock precursors). Our parametric model of the structure of C-shocks has been coupled with a radiative transfer code to calculate the molecular excitation and line profiles of shock tracers such as SiO, and of ion and neutral molecules such as H13CO+ and HN13C, as the shock propagates through the unperturbed medium. Our results show that the SiO emission arising from the early stage of the magnetic precursor typically has very narrow line profiles slightly shifted in velocity with respect to the ambient cloud. This narrow emission is generated in the region where the bulk of the ion fluid has already slipped to larger velocities in the precursor as observed toward the young L1448-mm outflow. This strongly suggests that the detection of narrow SiO emission and of an ion enhancement in young shocks, is produced by the magnetic precursor of C-shocks. In addition, our model shows that the different velocity components observed toward this outflow can be explained by the coexistence of different shocks at different evolutionary stages, within the same beam of the single-dish observations.

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

Glycolaldehyde is the simplest of the monosaccharide sugars and is directly linked to the origin of life. We report on the detection of glycolaldehyde (CH2OHCHO) towards the hot molecular core G31.41+0.31 through IRAM PdBI observations at 1.4, 2.1, and 2.9 mm. The CH2OHCHO emission comes from the hottest (> 300 K) and densest (>2E8 cm^-3) region closest (< 10^4 AU) to the (proto)stars. The comparison of data with gas-grain chemical models of hot cores suggests for G31.41+0.31 an age of a few 10^5 yr. We also show that only small amounts of CO need to be processed on grains in order for existing hot core gas-grain chemical models to reproduce the observed column densities of glycolaldehyde, making surface reactions the most feasible route to its formation.

We investigate the physical and chemical conditions necessary for low-mass star formation in extragalactic environments by calculating various characteristic timescales associated with star formation for a range of initial conditions. The balance of these timescales indicates whether low-mass star formation is enhanced or inhibited under certain physical conditions. In this study, we consider timescales for free-fall, cooling, freeze-out, desorption, chemistry and ambipolar diffusion and their variations with changes in the gas density, metallicity, cosmic ray ionisation rate and FUV radiation field strength. We find that extragalactic systems with high FUV radiation field strengths and high cosmic ray fluxes considered at a range of metallicities, are likely to have enhanced low-mass star formation unless the magnetic pressure is sufficient to halt collapse. Our results indicate that this is only likely to be the case for high-redshift galaxies approaching solar metallicities. Unless this is true for all high-redshift sources, this study finds little evidence for a high-mass biased IMF at high redshifts.

We present the 3.5m SPICA space telescope, a proposed Japanese-led JAXA-ESA mission scheduled for launch around 2017. The actively cooled (<5 K), single aperture telescope and monolithic mirror will operate from ~3.5 to ~210 um and will provide superb sensitivity in the mid- and far-IR spectral domain (better than JWST at lambda > 18 um). SPICA is one of the few space missions selected to go to the next stage of ESA's Cosmic Vision 2015-2025 selection process. In this White Paper we present the main specifications of the three instruments currently baselined for SPICA: a mid-infrared (MIR) coronagraph (~3.5 to ~27 um) with photometric and spectral capabilities (R~200), a MIR wide-field camera and high resolution spectrometer (R~30,000), and a far-infrared (FIR ~30 to ~210 um) imaging spectrometer - SAFARI - led by a European consortium. We discuss their capabilities in the context of MIR direct observations of exo-planets (EPs) and multiband photometry/high resolution spectroscopy observations of transiting exo-planets. We conclude that SPICA will be able to characterize the atmospheres of transiting exo-planets down to the super-Earth size previously detected by ground- or space-based observatories. It will also directly detect and characterize Jupiter/Neptune-size planets orbiting at larger separation from their parent star (>5-10 AU), by performing quantitative atmospheric spectroscopy and studying proto-planetary and debris disks. In addition, SPICA will be a scientific and technological precursor for future, more ambitious, IR space missions for exo-planet direct detection as it will, for example, quantify the prevalence exo-zodiacal clouds in planetary systems and test coronographic techniques, cryogenic systems and lightweight, high quality telescopes. (abridged)

We present the first detection of CS in the Antennae galaxies towards the NGC 4038 nucleus, as well as the first detections of two high-J (5-4 and 7-6) CS lines in the center of M 82. The CS(7-6) line in M 82 shows a profile that is surprisingly different to those of other low-J CS transitions we observed. This implies the presence of a separate, denser and warmer molecular gas component. The derived physical properties and the likely location of the CS(7-6) emission suggests an association with the supershell in the centre of M 82.

We present high resolution (R = 90,000) mid-infrared spectra of M dwarfs. The mid infrared region of the spectra of cool low mass stars contain pure rotational water vapour transitions that may provide us with a new methodology in the determination of the effective temperatures for low mass stars. We identify and assign water transitions in these spectra and determine how sensitive each pure rotational water transition is to small (25 K) changes in effective temperature. We find that, of the 36 confirmed and assigned pure rotational water transitions, at least 10 should be sensitive enough to be used as temperature indicators.

The evolution of star-forming core analogues undergoing inside-out collapse is studied with a multi-point chemodynamical model which self-consistently computes the abundance distribution of chemical species in the core. For several collapse periods the output chemistry of infall tracer species such as HCO+, CS, and N2H+, is then coupled to an accelerated Lambda-iteration radiative transfer code, which predicts the emerging molecular line profiles using two different input gas/dust temperature distributions. We investigate the sensitivity of the predicted spectral line profiles and line asymmetry ratios to the core temperature distribution, the time-dependent model chemistry, as well as to ad hoc abundance distributions. The line asymmetry is found to be strongly dependent on the adopted chemical abundance distribution. In general, models with a warm central region show higher values of blue asymmetry in optically thick HCO+ and CS lines than models with a starless core temperature profile. We find that in the formal context of Shu-type inside-out infall, and in the absence of rotation or outflows, the relative blue asymmetry of certain HCO+ and CS transitions is a function of time and, subject to the foregoing caveats, can act as a collapse chronometer. The sensitivity of simulated HCO+ line profiles to linear radial variations, subsonic or supersonic, of the internal turbulence field is investigated in the separate case of static cores.

As part of the first high resolution far-IR spectral survey of the Orion KL region (Lerate et al. 2006), we observed 20 CO emission lines with Jup=16 to Jup=39 (upper levels from approx 752 K to 4294 K above the ground state). Observations were taken using the Long Wavelength Spectrometer (LWS) on board the Infrared Space Observatory (ISO), in its high resolution Fabry-Perot (FP) mode (approx 33 km s$^{-1}$). We present here an analysis of the final calibrated CO data, performed with a more sophisticated modelling technique than hitherto, including a detailed analysis of the chemistry, and discuss similarities and differences with previous results. The inclusion of chemical modelling implies that atomic and molecular abundances are time-predicted by the chemistry. This provides one of the main differences with previous studies in which chemical abundances needed to be assumed as initial condition. The chemistry of the region is studied by simulating the conditions of the different known components of the KL region: chemical models for a hot core, a plateau and a ridge are coupled with an accelerated Lambda-iteration (ALI)radiative transfer model to predict line fluxes and profiles. We conclude that the CO transitions with 18<Jup<25 mainly arise from a hot core of diameter 0.02 pc and a density of 10$^{7}$ cm$^{-3}$ rather from the plateau as previous studies had indicated.

Hot core molecules should be detectable in external active galaxies out to high redshift. We present here a detailed study of the chemistry of star-forming regions under physical conditions that differ significantly from those likely to be appropriate in the Milky Way Galaxy. We examine, in particular, the trends in molecular abundances as a function of time with respect to changes in the relevant physical parameters. These parameters include metallicity, dust:gas mass ratio, the H$_{2}$ formation rate, relative initial elemental abundances, the cosmic ray ionization rate, and the temperature of hot cores. These trends indicate how different tracers provide information on the physical conditions and on evolutionary age. We identify hot core tracers for several observed galaxies that are considered to represent spirals, active galaxies, low-metallicity galaxies, and high-redshift galaxies. Even in low-metallicity examples, many potential molecular tracers should be present at levels high enough to allow unresolved detection of active galaxies at high redshift containing large numbers of hot cores.

Knots are commonly found in nearby planetary nebulae (PNe) and star forming regions. Within PNe, knots are often found to be associated with the brightest parts of the nebulae and understanding the physics involved in knots may reveal the processes dominating in PNe. As one of the closest PNe, the Helix Nebula (NGC 7293) is an ideal target to study such small-scale (~300 AU) structures. We have obtained infrared integral spectroscopy of a comet-shaped knot in the Helix Nebula using SINFONI on the Very Large Telescope at high spatial resolution (50-125 mas). With spatially resolved 2 micron spectra, we find that the H2 rotational temperature within the cometary knots is uniform. The rotational-vibrational temperature of the cometary knot (situated in the innermost region of the nebula, 2.5 arcmin away from the central star), is 1800 K, higher than the temperature seen in the outer regions (5-6 arcmin from the central star) of the nebula (900 K), showing that the excitation temperature varies across the nebula. The obtained intensities are reasonably well fitted with 27 km s-1 C-type shock model. This ambient gas velocity is slightly higher than the observed [HeII] wind velocity of 13 km s-1. The gas excitation can also be reproduced with a PDR (photo dominant region) model, but this requires an order of magnitude higher UV radiation. Both models have limitations, highlighting the need for models that treats both hydrodynamical physics and the PDR.

We have carried out a systematic search for the molecular ion CO+ in a sample of 8 protoplanetary and planetary nebulae in order to determine the origin of the unexpectedly strong HCO+ emission previously detected in these sources. An understanding of the HCO+ chemistry may provide direct clues to the physical and chemical evolution of planetary nebulae. We find that the integrated intensity of the CO+ line may be correlated with that of HCO+, suggesting that the reaction of CO+ with molecular hydrogen may be an important formation route for HCO+ in these planetary nebulae.

The desorption of molecular species from ice mantles back into the gas phase in molecular clouds results from a variety of very poorly understood processes. We have investigated three mechanisms; desorption resulting from H_2 formation on grains, direct cosmic ray heating and cosmic ray induced photodesorption. Whilst qualitative differences exist between these processes (essentially deriving from the assumptions concerning the species-selectivity of the desorption and the assumed threshold adsorption energies, E_t) all three processes are found to be potentially very significant in dark cloud conditions. It is therefore important that all three mechanisms should be considered in studies of molecular clouds in which freeze-out and desorption are believed to be important. Employing a chemical model of a typical static molecular core and using likely estimates for the quantum yields of the three processes we find that desorption by H_2 formation probably dominates over the other two mechanisms. However, the physics of the desorption processes and the nature of the dust grains and ice mantles are very poorly constrained. We therefore conclude that the best approach is to set empirical constraints on the desorption, based on observed molecular depletions - rather than try to establish the desorption efficiencies from purely theoretical considerations. Applying this method to one such object (L1689B) yields upper limits to the desorption efficiencies that are consistent with our understanding of these mechanisms.

We present high spatial resolution maps, obtained with the Plateau de Bure Interferometer, of the blue lobe of the L1157 outflow. We observed four lines at 3 mm, namely CH3OH (2_K-1_K), HC3N (11-10), HCN (1-0) and OCS (7-6). Moreover, the bright B1 clump has also been observed at better spatial resolution in CS (2-1), CH3OH (2_1-1_1)A-, and 34SO (3_2-2_1). These high spatial resolution observations show a very rich structure in all the tracers, revealing a clumpy structure of the gas superimposed to an extended emission. In fact, the three clumps detected by previous IRAM-30m single dish observations have been resolved into several sub-clumps and new clumps have been detected in the outflow. The clumps are associated with the two cavities created by two shock episodes driven by the precessing jet. In particular, the clumps nearest the protostar are located at the walls of the younger cavity with a clear arch-shape form while the farthest clumps have slightly different observational characteristics indicating that they are associated to the older shock episode. The emission of the observed species peaks in different part of the lobe: the east clumps are brighter in HC3N (11-10), HCN (1-0) and CS (2-1) while the west clumps are brighter in CH3OH(2_K-1_K), OCS (7-6) and 34SO (3_2-2_1). This peak displacement in the line emission suggests a variation of the physical conditions and/or the chemical composition along the lobe of the outflow at small scale, likely related to the shock activity and the precession of the outflow. In particular, we observe the decoupling of the silicon monoxide and methanol emission, common shock tracers, in the B1 clump located at the apex of the bow shock produced by the second shock episode.

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

We present SEST observations of the molecular environment ahead of the southern Herbig-Haro object 188 (HH188), associated with the low-mass protostar Re4 IRS. We have also used the SuperCosmos Halpha survey to search for Halpha emission associated with the Re4 IRS - HH188 region. The aim of the present work is to study the properties of the molecular gas and to better characterize this southern star forming region. We mapped the HCO+ 3-2 and H13CO+ 1-0 emission around the YSO and took spectra of the CH3OH 2(0)-1(0)A+ and 2(-1)-1(-1)E and SO 6(5)-5(4) towards the central source. Column densities are derived and different scenarios are considered to explain the origin of the molecular emission. HCO+ arises from a relatively compact region around the YSO; however, its peak emission is displaced to the south following the outflow direction. Our chemical analysis indicates that a plausible scenario is that most of the emission arises from the cold, illuminated dense gas ahead of the HH188 object. We have also found that HH188, a high excitation object, seems to be part of a parsec scale and highly collimated HH system. Re4 IRS is probably a binary protostellar system, in the late Class 0 or Class I phase. One of the protostars, invisible in the near-IR, seems to power the HH188 system.

The sub-millimetre wavelength regime is perhaps the most poorly explored over large areas of the sky, despite the considerable effort that has been expended in making deep maps over small regions. As a consequence the properties of the sub-millimetre sky as a whole, and of rare bright objects in particular, remains largely unknown. Here we describe a forthcoming survey (the SCUBA-2 ``All-Sky'' Survey, or SASSy) designed to address this issue by making a large-area map of approximately one-fifth of the sky visible from the JCMT (4800 square degrees) down to a 1 sigma noise level of 30 mJy/beam. This map forms the pilot for a much larger survey, which will potentially map the remaining sky visible from the JCMT, with the region also visible to ALMA as a priority. SASSy has been awarded 500 hours for the 4800 square degree pilot phase and will commence after the commissioning of SCUBA-2, expected in early 2008.

We present theoretically-established values of the CO-to-H2 and C-to-H2 conversion factors that may be used to estimate the gas masses of external galaxies. We consider four distinct galaxy types, represented by M51, NGC 6946, M82 and SMC N27. The physical parameters that best represent the conditions within the molecular clouds in each of the galaxy types are estimated using a chi^2 analysis of several observed atomic fine structure and CO rotational lines. This analysis is explored over a wide range of density, radiation field, extinction, and other relevant parameters. Using these estimated physical conditions in methods that we have previously established, CO-to-H2 conversion factors are then computed for CO transitions up to J=9-8. For the conventional CO(1-0) transition, the computed conversion factor varies significantly below and above the canonical value for the Milky Way in the four galaxy types considered. Since atomic carbon emission is now frequently used as a probe of external galaxies, we also present, for the first time, the C-to-H2 conversion factor for this emission in the four galaxy types considered.

We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.

We report experimental oscillator strengths for 88 Ti I transitions covering the wavelength range 465 to 3892 nm, 67 of which had no previous experimental values. Radiative lifetimes for thirteen energy levels, including the low energy levels 3d2 (3F) 4s4p (3P) z5Dj, have been measured using time resolved laser induced fluorescence. Intensity calibrated Ti I spectra have been measured using Fourier transform spectroscopy to determine branching fractions for the decay channels of these levels. The branching fractions are combined with the radiative lifetimes to yield absolute transition probabilities and oscillator strengths. Our measurements include 50 transitions in the previously unobserved infrared region lambda > 1000 nm, a region of particular interest to the analysis of cool stars and brown dwarfs.

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

Magnetic precursors of C-shocks accelerate, compress and heat molecular ions, modifying the kinematics and the physical conditions of the ion fluid with respect to the neutral one. Electron densities are also expected to be significantly enhanced in shock precursors. In this Letter, we present observations of strongly polar ion and neutral molecules such as SiO, H13CO+, HN13C and H13CN, which reveal the electron density enhancements associated with the precursor of the young L1448-mm outflow. While in the ambient gas the excitation of the ions and neutrals is explained by collisional excitation by H2 with a single density of 10E5cmE-3, H13CO+ shows an over excitation in the shock precursor component that requires H2 densities of a factor of >10 larger than those derived from the neutral species. This over excitation in H13CO+ can be explained if we consider an additional excitation by collisions with electrons and an electron density enhancement in the precursor stage by a factor of 500, i.e. a fractional ionization of 5x10E-5. These results show that multiline observations can be used to study the evolution of the ion and electron fluids at the first stages of the C-shock interaction.

We present the detection and modeling of more than 70 far-IR pure rotational lines of water vapor, including the 18O and 17O isotopologues, towards Orion KL. Observations were performed with the Long Wavelength Spectrometer Fabry-Perot (LWS/FP; R~6800-9700) on board the Infrared Space Observatory (ISO) between ~43 and ~197 um. The water line profiles evolve from P-Cygni type profiles (even for the H2O18 lines) to pure emission at wavelengths above ~100 um. We find that most of the water emission/absorption arises from an extended flow of gas expanding at 25+-5 kms^-1. Non-local radiative transfer models show that much of the water excitation and line profile formation is driven by the dust continuum emission. The derived beam averaged water abundance is 2-3x10^-5. The inferred gas temperature Tk=80-100 K suggests that: (i) water could have been formed in the "plateau" by gas phase neutral-neutral reactions with activation barriers if the gas was previously heated (e.g. by shocks) to >500 K and/or (ii) H2O formation in the outflow is dominated by in-situ evaporation of grain water-ice mantles and/or (iii) H2O was formed in the innermost and warmer regions (e.g. the hot core) and was swept up in ~1000 yr, the dynamical timescale of the outflow.

Clouds containing molecular dark matter in quantities relevant for star formation may exist in minihaloes of the type of cold dark matter included in many cosmological simulations or in the regions of some galaxies extending far beyond their currently known boundaries. We have systematically explored parameter space to identify conditions under which plane-parallel clouds contain sufficient column densities of molecular dark matter that they could be reservoirs for future star formation. Such clouds would be undetected or at least appear by current observational criteria to be uninteresting from the perspective of star formation. We use a time-dependent PDR code to produce theoretical models of the chemistry and emission arising in clouds for our chosen region of parameter space. We then select a subset of model clouds with levels of emission that are low enough to be undetectable or at least overlooked by current surveys. The existence of significant column densities of cold molecular dark matter requires that the background radiation field be several or more orders of magnitude weaker than that in the solar neighbourhood. Lower turbulent velocities and cosmic ray induced ionization rates than typically associated with molecular material within a kpc of the Sun are also required for the molecular matter to be dark. We find that there is a large region within the parameter space that results in clouds that might contain a significant mass of molecular gas whilst remaining effectively undetectable or at least not particularly noticeable in surveys. We note briefly conditions under which molecular dark matter may contain a dynamically interesting mass.

A reliable estimate of the molecular gas content in galaxies plays a crucial role in determining their dynamical and star-forming properties. However, H2, the dominant molecular species, is difficult to observe directly, particularly in the regions where most molecular gas is thought to reside. Its mass is therefore commonly inferred by assuming a direct proportionality with the integrated intensity of the CO(J=1-0) emission line, using a CO-to-H2 conversion factor, X. Although a canonical value for X is used extensively in such estimates, there is increasing evidence, both theoretical and observational, that the conversion factor may vary by over an order of magnitude under conditions different to those of the local neighbourhood. In an effort to understand the influence of changing environmental conditions on the conversion factor, we derive theoretical estimates of X for a wide range of physical parameters using a photon-dominated region (PDR) time-dependent chemical model, benchmarking key results against those of an independent PDR code to ensure reliability. Based on these results, the sensitivity of the X factor to change in each physical parameter is interpreted in terms of the chemistry and physical processes within the cloud. In addition to confirming previous observationally derived trends, we find that the time-dependence of the chemistry, often neglected in such models, has a considerable influence on the value of the conversion factor.

Recent observations reveal galaxies in the early Universe (2<z<6.4) with large reservoirs of molecular gas and extreme star formation rates. For a very large range of sources, a tight relationship exists between star formation rate and the luminosity of the HCN J=1-0 spectral line, but sources at redshifts of z~2 and beyond do not follow this trend. The deficit in HCN is conventionally explained by an excess of infrared (IR) radiation due to active galactic nuclei (AGN). We show in this letter not only that the presence of AGN cannot account for the excess of IR over molecular luminosity, but also that the observed abundance of HCN is in fact consistent with a population of stars forming from near-primordial gas.

We have carried out a high spectral resolution line survey towards the Orion Kleinmann-Low (KL) cluster from 44-188 um. The observations were taken with the Long Wavelength Spectrometer (LWS) in Fabry-Perot mode, on board the Infrared Space Observatory (ISO). A total of 152 lines are clearly detected and a further 34 features are present as possible detections. The spectrum is dominated by the molecular species H2O, OH and CO, along with [OI] and [CII] lines from PDR or shocked gas and [OIII], [NIII] lines from the foreground M42 HII region. Several isotopic species, as well as NH3, are also detected. HDO and H3O+ are tentatively detected for the first time in the far-infrared range towards Orion-KL. A basic analysis of the line observations is carried out, by comparing with previous measurements and published models and deriving rotational temperatures and column densities in the case of the molecular species. The complexity of the region requires more sophisticated models for the interpretation of all the line observations.

We present new observations of the CepA-East region of massive star formation and describe an extended and dynamically distinct feature not previously recognised. This feature is present in emission from H2CS, OCS, CH3OH, and HDO at -5.5 km/s, but is not traced by conventional tracers of star forming regions H2S, SO2, SO, CS. The feature is extended up to at least 0.1 pc. We show that the feature is neither a hot core nor a shocked outflow. However, the chemistry of the feature is consistent with predictions of a model of an eroding interface between a fast wind and a dense core; mixing between the two media occurs in the interface on a timescale of 10-50 years. If these observations are confirmed by detailed maps and by detections in species also predicted to be abundant (e.g. HCO+, H2CO, and NH3) this feature would be the first detection of such an interface in regions of massive star formation. An important implication of the model is that a significant reservoir of sulfur in grain mantles is required to be in the form of OCS.

In order to investigate the origin and the structure of the low velocity, chemically rich clumps observed along the lobes of low- and intermediate-mass outflows, we construct a detailed model of the S1 clump along the CB3 outflow. We use a time-dependent chemical model coupled with a radiative transfer model to reproduce the observed line profile for a direct comparison with previous observations of this clump. We find that the simultaneous fitting of multiple species and transitions is a powerful tool in constraining the physical parameters of the gas. Different scenarios for the clump formation have been investigated. The models that better reproduce all the observed lines are those where the clump is formed, at least partially, before the advent of the outflow; with the advent of the outflow the clump undergoes a short period of non-dissociative shock and the consequent release of the icy mantle together with the high temperature chemistry leads to the observed chemical enrichment. Our results also suggest the presence of substructure within the clump: a more extended component traced by CS, SO and the lower energy transitions (3_k-2_k and 2_k-1_k) of CH_3OH, and a more compact component traced by H_2CO, SO_2 and the higher energy transitions (5_k-4_k) of CH_3OH.

As part of the first far-IR line survey towards Orion KL, we present the detection of seven new rotationally excited OH Lambda-doublets (at 48, 65, 71, 79, 98 and 115 um). Observations were performed with the Long Wavelength Spectrometer (LWS) Fabry-Perots on board the Infrared Space Observatory (ISO). In total, more than 20 resolved OH rotational lines, with upper energy levels up to 620 K, have been detected at an angular and velocity resolutions of 80$'' and 33 km s^-1 respectively. OH line profiles show a complex behavior evolving from pure absorption, P-Cygni type to pure emission. We also present a large scale 6' declination raster in the OH ^2\Pi_3/2 J=5/2^+-3/2^- and ^2\Pi_3/2 J=7/2^-5/2^+ lines (at 119.441 and 84.597 um) revealing the decrease of excitation outside the core of the cloud. From the observed profiles, mean intrinsic line widths and velocity offsets between emission and absorption line peaks we conclude that most of the excited OH arises from Orion outflow(s), i.e. the ``plateau'' component. We determine an averaged OH abundance relative to H_2 of X(OH)=(0.5-1.0)x10^-6, a kinetic temperature of 100 K and a density of n(H_2)=5x10^5 cm^-3. Even with these conditions, the OH excitation is heavily coupled with the strong dust continuum emission from the inner hot core regions and from the expanding flow itself.

Context. A molecular survey recently performed ahead of HH~2 supports the idea that the observed molecular enhancement is due to UV radiation from the HH object. Aims. The aim of the present work is to determine whether all HH objects with enhanced HCO$^+$ emission ahead of them also exhibit the same enhanced chemistry as HH~2. We thus observed several molecular lines at several positions ahead of five Herbig-Haro objects where enhanced HCO$^+$ emission was previously observed. Methods. We mapped the five Herbig-Haro objects using the IRAM-30 m. For each position we searched for more than one molecular species, and where possible for more than one transition per species. We then estimated the averaged beam column densities for all species observed and also performed LVG analyses to constrain the physical properties of the gas. Results. The chemically richest quiescent gas is found ahead of the HH~7-11 complex, in particular at the HH~7-11 A position. In some regions we also detected a high velocity gas component. We find that the gas densities are always higher than those typical of a molecular cloud while the derived temperatures are always quite low, ranging from 10 to 25 K. The emission of most species seems to be enhanced with respect to that of a typical dense clump, probably due to the exposure to a high UV radiation from the HH objects. Chemical differentiation among the positions is also observed. We attempt a very simple chemical analysis to explain such differentiation.

We present models for the chemistry in gas moving towards the ionization front of an HII region. When it is far from the ionization front, the gas is highly depleted of elements more massive than helium. However, as it approaches the ionization front, ices are destroyed and species formed on the grain surfaces are injected into the gas phase. Photodissociation removes gas phase molecular species as the gas flows towards the ionization front. We identify models for which the OH column densities are comparable to those measured in observations undertaken to study the magnetic fields in star forming regions and give results for the column densities of other species that should be abundant if the observed OH arises through a combination of the liberation of H2O from surfaces and photodissociation. They include CH3OH, H2CO, and H2S. Observations of these other species may help establish the nature of the OH spatial distribution in the clouds, which is important for the interpretation of the magnetic field results.

The effect of a newly born star cluster inside a giant molecular cloud (GMC) is to produce a hot bubble and a thin, dense shell of interstellar gas and dust swept up by the H II expansion, strong stellar winds, and repeated supernova explosions. Lying at the inner side of the shell is the photodissociation region (PDR), the origin of much of the far-infrared/sub-millimeter/millimeter (FIR/sub-mm/mm) radiation from the interstellar medium (ISM). We present a model for the expanding shell at different stages of its expansion which predict mm/sub-mm and far-IR emission line intensities from a series of key molecular and atomic constituents in the shell. The kinematic properties of the swept-up shell predicted by our model are in very good agreement with the measurements of the supershell detected in the nearby starburst galaxy M 82. We compare the modeling results with the ratio-ratio plots of the FIR/sub-mm/mm line emission in the central 1.0 kpc region to investigate the mechanism of star forming activity in M 82. Our model has yielded appropriate gas densities, temperatures, and structure scales compared to those measured in M 82, and the total H2 content is compatible with the observations. This implies that the neutral ISM of the central star-forming region is a product of fragments of the evolving shells.

We present the results of a multiline survey at mm-wavelengths of the Cepheus A star forming region. Four main flows have been identified: three pointing in the SW, NE, and SE directions and accelerating high density CS clumps. The fourth outflow, revealed by high-sensitivity HDO observations, is pointing towards South and is associated with conditions particularly favourable to a chemical enrichment. The analysis of the line profiles shows that the SiO molecule dominates at the highest velocities and at the highest excitation conditions, confirming its close association with shocks. H2S, SO2, and SO preferentially trace more quiescent regions than SiO, and in particular a lack of bright H2S emission at the highest velocities is found. OCS and H2CS emit at quite high velocities, where the abundances of three shock tracers like SiO, CH3OH, and HDO are higher. These results may indicate that H2S is not the only major sulphur carrier in the grain mantles, and that OCS and H2CS may probably play an important role on the grains; or that alternatively they rapidly form once the mantle is evaporated after the passage of a shock. Finally, the outflow peak emission has been compared with recent time-dependent sulphur chemistry models: the results indicate that, if associated with accurate measurements of the physical conditions, the CH3OH/H2CS column density ratio can be used as an effective chemical clock to date the age of shocked gas.

The very high rates of second generation star formation detected and inferred in high redshift objects should be accompanied by intense millimetre-wave emission from hot core molecules. We calculate the molecular abundances likely to arise in hot cores associated with massive star formation at high redshift, using several independent models of metallicity in the early Universe. If the number of hot cores exceeds that in the Milky Way Galaxy by a factor of at least one thousand, then a wide range of molecules in high redshift hot cores should have detectable emission. It should be possible to distinguish between independent models for the production of metals and hence hot core molecules should be useful probes of star formation at high redshift.

We present the results of an optical spectroscopic study of the massive Wolf-Rayet binary WR137. These data cover the dust-formation maximum in 1997. Combining all available measurements of radial velocities, we derive, for the first time, a spectroscopic orbit with period 13.05 +/- 0.18 years. The resulting masses, adopting i=67degrees, are M(O)= 20 +/- 2 Mo and M(WR)= 4.4+/- 1.5 Mo. These appear, respectively, around normal and on the low side for the given spectral types. Analysis of the intense multi-site spectroscopic monitoring in 1999 shows that the CIII5696 and CIV5802/12 lines have the highest intrinsic variability levels. The periodogram analysis yields a small-amplitude modulation in the absorption troughs of the CIV5802/12 and HeI5876 lines with a period of 0.83 days, which could be related either to pulsations or large-scale rotating structures as seen in the WN4 star EZ Canis Majoris (WR6).Wavelet analysis of the strong emission lines of CIII5696 and CIV5802/12 enabled us to isolate and follow for several hours small structures associated with density enhancements within the wind of the Wolf-Rayet star. Cross-correlating the variability patterns seen in different lines, we find a weak but significant correlation between the varability in emission lines with different ionization potential, i.e. in lines formed at different distances from the WR stellar core. Adopting a beta wind-velocity law, from the motion of individual subpeaks we find beta around 5, which is significantly larger than the canonical value of 1 usually found in O-star winds.

We have carried out an extensive observational study (from BIMA data) and made a preliminary theoretical investigation of the molecular gas around HH2. The molecular maps show a very complex morphological, kinematical and chemical structure. The overall main conclusion of this work confirms the findings of Paper I and II, by demonstrating that in addition to the strong photochemical effects caused by penetration of the UV photons from HH2 into molecular cloud, a range of complex radiative and dynamical interactions occur. Thus, despite the apparent `quiescent' nature of the molecular cloud ahead of HH2, the kinematical properties observed within the field of view suggest that it is possibly being driven out by powerful winds from the VLA 1 protostar.

We compare high resolution infrared observations of the CO 3-1 bands in the 2.297-2.310 micron region of M dwarfs and one L dwarf with theoretical expectations. We find a good match between the observational and synthetic spectra throughout the 2000-3500K temperature regime investigated. Nonetheless, for the 2500-3500 K temperature range the temperatures that we derive from synthetic spectral fits are higher than expected from more empirical methods by several hundred K. In order to reconcile our findings with the empirical temperature scale it is necessary to invoke warming of the model atmosphere used to construct the synthetic spectra. We consider that the most likely reason for the back-warming is missing high temperature opacity due to water vapour. We compare the water vapour opacity of the Partridge & Schwenke (1997) line list used for the model atmosphere with the output from a preliminary calculation by Barber & Tennyson (2004). While the Partridge & Schwenke line list is a reasonable spectroscopic match for the new line list at 2000 K, by 4000 K it is missing around 25% of the water vapour opacity. We thus consider that the offset between empirical and synthetic temperature scales is explained by the lack of hot water vapour used for computation of the synthetic spectra. For our coolest objects with temperatures below 2500 K we find best fits when using synthetic spectra which include dust emission. Our spectra also allow us to constrain the rotational velocities of our sources, and these velocities are consistent with the broad trend of rotational velocities increasing from M to L.

Transient microstructure in the diffuse interstellar medium has been observed towards galactic and extragalactic sources for decades, usually in lines of atoms and ions, and, more recently, in molecular lines. Evidently, there is a molecular component to the transient microstructure. In this paper we explore the chemistry that may arise in such microstructure. We use a PDR code to model the conditions of relatively high density, low temperature, very low visual extinction and very short elapsed time that are appropriate for these objects. We find that there is a well defined region of parameter space where detectable abundances of molecular species might be found. The best matching models are those where the interstellar microstructure is young (<100 yr), small (~100 AU), and dense (>10^4 cm^-3).

Recent observations suggest that the behaviour of tracer species such as N_2H+ and CS is significantly different in regions of high and low mass star formation. In the latter, N_2H+ is a good tracer of mass, while CS is not. Observations show the reverse to be true in high-mass star formation regions. We use a computational chemical model to show that the abundances of these and other species may be significantly altered by a period of accelerated collapse in high mass star forming regions. We suggest these results provide a potential explanation of the observations, and make predictions for the behaviour of other species.

Hot cores and their precursors contain an integrated record of the physics of the collapse process in the chemistry of the ices deposited during that collapse. In this paper, we present results from a new model of the chemistry near high mass stars in which the desorption of each species in the ice mixture is described as indicated by new experimental results obtained under conditions similar to those hot cores. Our models show that provided there is a monotonic increase in the temperature of the gas and dust surrounding the protostar, the changes in the chemical evolution of each species due to differential desorption are important. The species H$_2$S, SO, SO$_2$, OCS, H$_2$CS, CS, NS, CH$_3$OH, HCOOCH$_3$, CH$_2$CO, C$_2$H$_5$OH show a strong time dependence that may be a useful signature of time evolution in the warm-up phase as the star moves on to the Main Sequence. This preliminary study demonstrates the consequences of incorporating reliable TPD data into chemical models.

We investigate the origin of the small, chemically rich molecular clumps observed along the main axis of chemically rich outflows such as CB3 and L1157. We develop a chemical model where we explore the chemical evolution of these clumps, assuming they are partially pre-existing to the outflow, or alternatively newly formed by the impact of the outflow on the surrounding medium. The effects of the impact of the outflow are reproduced by density and temperature changes in the clump. We find that the observed abundances of CH3OH, SO and SO2 are best reproduced by assuming a scenario where the dense molecular gas observed is probably pre-existing in the interstellar medium before the formation of their exciting (proto)stars and that the clumpiness and the rich chemistry of the clumps are a consequence of a pre-existing density enhancement and of its interaction with the outflow.

We provide a set of atomic lines which are suitable for the description of ultracool dwarf spectra from 10000 to 25000 \AA. This atomic linelist was made using both synthetic spectra calculations and existing atlases of infrared spectra of Arcturus and Sunspot umbra. We present plots, which show the comparison of synthetic spectra and observed Arcturus and Sunspot umbral spectra for all atomic lines likely to be observable in high resolution infrared spectra.

We study the chemistry in the slow wind during the transition from AGB star to pre-planetary nebula (PPN) and planetary nebula (PN). We show that there is a very rich chemistry of degradation products created by photoprocessing, driven by the gradually hardening radiation field of the central star. Most of these products are, however, photodissociated during the PPN phase if the wind is smooth. By contrast, if the wind is clumpy, possibly because of clumpiness in the AGB atmosphere, then many of these degradation products survive into the PN phase. Thus, chemistry may be used to infer the existence of clumpiness in the AGB phase. We identify potential molecular tracers, and we note that, in the case of clumpiness, large molecules may survive the transport from the stellar atmosphere to the interstellar medium. We compare between our model results with observations of three objects at differing evolutionary stages: CRL618, NGC 7027 and the Helix nebula (NGC 7293).

There is now strong evidence suggesting that the 12CO J = 1-0 transition, widely used to trace H2 gas, significantly underestimates its mass in metal-poor regions. In spiral disks such regions are found in large galactocentric distances where we show that any unaccounted H2 gas phase is likely to be diffuse (~5-20 cm^-3) and warmer (T(kin) ~ 50-100 K) than the cool (T(kin) ~ 15-20 K) CO-luminous one. Moreover we find that a high value of the H2 formation rate on grains, suggested by recent observational work, can compensate for the reduction of the available grain surface in the metal-poor part of typical galactic disks and thus enhance this CO-poor H2 component which may be contributing significantly to the mass and pressure of spiral disks beyond their optical radius.

We present a CSO and BIMA molecular line survey of the dense, quiescent molecular environment ahead of HH2. The molecular gas is cold, 13 K, and moderately dense, 3x10^5 cm^-3. A total of 14 species has been detected (including different isotopes and deuterated species). The relative abundances of the clump are compared with other dense molecular environments, including quiescent dark clouds, and active low and high mass star forming regions. This comparison confirms the peculiar chemical composition of the quiescent gas irradiated by the HH objects. Thus, from this comparison, we found that the HCO+, CH3OH and H2CO are strongly enhanced. SO and SO2 are weakly enhanced, whereas HCN and CS are underabundant. The CN abundance is within the range of value found in starless dark clouds, but it is low with respect to high mass star forming regions. Finally, the chemical composition of HH2 confirms the qualitative results of the Viti & Williams (1999) complex chemical model that follows the chemical behavior of a molecular clump irradiated by a HH object.

M star spectra, at wavelengths beyond 1.35 $\mu$m, are dominated by water vapour yet terrestrial water vapour makes it notoriously difficult to make accurate measurement from ground-based observations. We have used the short wavelength spectrometer on the Infrared Space Observatory at four wavelength settings to cover the 2.5--3.0 $\mu$m region for a range of M stars. The observations show a good match with previous ground-based observations and with synthetic spectra based on the Partridge & Schwenke (1997) line list though not with the SCAN (Jorgensen et al. 2001) line list. We used a least-squared minimisation technique to systematically find best fit parameters for the sample of stars. The temperatures that we find indicate a relatively hot temperature scale for M dwarfs. We consider that this could be a consequence of problems with the Partridge & Schwenke linelist which leads to synthetic spectra predicting water bands which are too strong for a given temperature. Such problems need to be solved in the next generation of water vapour line lists which will extend the calculation of water vapour to higher energy levels with the good convergence necessary for reliable modelling of hot water vapour. Then water bands can assume their natural role as the primary tool for the spectroscopic analysis of M stars.

HH80N is one of the Herbig-Haro objects that have associated quiescent dense clumps. We report CO and CS BIMA observations that reveal star formation within the HH80N dense clump. The CO emission reveals clearly a bipolar molecular outflow centered on the dense clump. The CS emission traces a ring-like structure of radius ~0.24 pc. The CS kinematics shows that the ring is collapsing with an infall speed of ~0.6 km/s. The required mass to produce the collapse is in agreement with previous ammonia observations of the 20 solar mass core, which is embedded within the CS structure. However, we cannot discard that the ring structure is expanding driven by protostellar winds, if the CS abundance if unusually high and the CO momentum rate is much higher than that measured, due to inclination and optical depth effects. The properties of the molecular outflow and of the dense core suggest that it harbors a Class 0 object. There are also signatures of interaction of the HH 80/81/80N outflow with the dense gas. In particular it is possible that the HH 80/81/80N outflow has triggered or at least speed up the star formation in this region.