Large millimeter interferometers are revealing a growing number of rotating outflows, which are suggested to trace magneto-centrifugal disk winds (MHD DWs). However, their impact on disk accretion is not yet well quantified. Here we identify systematic biases in retrieving the true launch zone, magnetic lever arm, and angular momentum flux of an MHD DW from apparent rotation signatures. Synthetic position-velocity cuts are constructed from self-similar MHD DWs over a broad range of parameters, and three different methods are applied for estimating the specific angular momentum. We find that the launch radius inferred using the well-known relation from Anderson et al. (2006) can markedly differ from the true outermost launch radius $r_{out}$ of the DW. The "double-peak separation" and "flow width" methods provide only a strict lower limit to $r_{out}$. This bias is independent of angular resolution and can reach a factor ten. In contrast, the "rotation curve" method gives a good estimate of $r_{out}$ when the flow is well resolved, and an upper limit otherwise. The magnetic lever arm is always underestimated. Only comparison with synthetic predictions can take into account properly all observational effects. As an application, we present a comparison with ALMA observations of HH212 at resolutions from 250 au to 16 au, which represents the most stringent observational test of MHD DW to date. This comparison confirms our predicted biases for the double-peak separation method, and the large $r_{out}\sim40~$au and small magnetic lever arm first suggested by Tabone et al. (2017). We also derive the first accurate analytical expression for the fraction of disk angular momentum extracted by an MHD disk wind of given radial extent, magnetic lever arm, and mass flux. Application to HH212 confirms that MHD DWs are serious candidates for the steady angular momentum extraction process in young disks.

Context. Low-mass protostars drive powerful molecular outflows that can be observed with mm and sub-mm telescopes. Various sulfuretted species are known to be bright in shocks and could be used to infer the physical and chemical conditions throughout the observed outflows. Aims. The evolution of sulfur chemistry is studied along the outflows driven by the NGC1333-IRAS4A protobinary system located in the Perseus cloud to constrain the physical and chemical processes at work in shocks. Methods. We observed various transitions from OCS, CS, SO, and SO$_2$ towards NGC1333-IRAS4A in the 1.3, 2, and 3mm bands using the IRAM NOEMA array and we interpreted the observations through the use of the Paris-Durham shock model. Results. The targeted species clearly show different spatial emission along the two outflows driven by IRAS4A. OCS is brighter on small and large scales along the south outflow driven by IRAS4A1, whereas SO$_2$ is detected rather along the outflow driven by IRAS4A2 that is extended along the north east - south west (NE-SW) direction. Column density ratio maps estimated from a rotational diagram analysis allowed us to confirm a clear gradient of the OCS/SO$_2$ column density ratio between the IRAS4A1 and IRAS4A2 outflows. SO is detected at extremely high radial velocity up to 25 km/s relative to the source velocity, clearly allowing us to distinguish the two outflows on small scales. Conclusions. The observed chemical differentiation between the two outflows of the IRAS4A system could be explained by a different chemical history. The outflow driven by IRAS4A1 is likely younger and more enriched in species initially formed in interstellar ices, such as OCS, and recently sputtered into the shock gas. In contrast, the longer and likely older outflow triggered by IRAS4A2 is more enriched in species that have a gas phase origin, such as SO$_2$.

The accretion history of protostars remains widely mysterious even though it represents one of the best ways to understand the protostellar collapse that leads to the formation of stars. Molecular outflows are here used to characterize the protostellar accretion phase in W43-MM1. The W43-MM1 protocluster host a sufficient number of protostars to statistically investigate molecular outflows in a single, homogeneous region. We used the CO(2-1) and SiO(5-4) line datacubes, taken as part of an ALMA mosaic with a 2000 AU resolution, to search for protostellar outflows, evaluate the influence that the environment has on these outflows' characteristics and put constraints on outflow variability in W43-MM1. We discovered a rich cluster of 46 outflow lobes, driven by 27 protostars with masses of 1-100 Msun. The complex environment inside which these outflow lobes develop has a definite influence on their length, limiting the validity of using outflows' dynamical timescales as a proxy of the ejection timescale in clouds with high dynamics and varying conditions. We performed a detailed study of Position-Velocity (PV) diagrams of outflows that revealed clear events of episodic ejection. The time variability of W43-MM1 outflows is a general trend and is more generally observed than in nearby, low- to intermediate-mass star-forming regions. The typical timescale found between two ejecta, about 500 yr, is consistent with that found in nearby protostars. If ejection episodicity reflects variability in the accretion process, either protostellar accretion is more variable or episodicity is easier to detect in high-mass star-forming regions than in nearby clouds. The timescale found between accretion events could be resulting from disk instabilities, associated with bursts of inflowing gas arising from the dynamical environment of high-mass star-forming cores.

The investigation of star forming regions have enormously benefited from the recent advent of the ALMA interferometer. More specifically, the unprecedented combination of high-sensitivity and high-angular resolution provided by ALMA allows one to shed light on the jet/disk systems associated with a Sun-like mass protostar. Also astrochemistry enjoyed the possibility to analyze complex spectra obtained using large bandwidths: several interstellar Complex Organic Molecules (iCOMs; C-bearing species with at least 6 atoms) have been imaged around protostars. This in turn boosted the study of the astrochemistry at work during the earliest phases of star formation paving the way to the chemical complexity in planetary systems where Life could emerge. There is mounting evidence that the observations of iCOMs can be used as unique tool to shed light, on Solar System scales (< 50 au), on the molecular content of protostellar disk. The increase of iCOMs abundances occur only under very selective physical conditions, such as those associated low-velocity shocks found where the infalling envelope is impacting the rotating accretion disk. The imaging of these regions with simpler molecules such as CO or CS is indeed paradoxically hampered by their high abundances and consequently high line opacities which do not allow the observers to disentangle all the emitting components at these small scales. In this respect, we review the state-of-the art of the ALMA analysis about the standard Sun-like star forming region in Orion named HH 212. We show (i) how all the physical components involved in the formation of a Sun-like star can be revealed only by observing different molecular tracers, and (ii) how the observation of iCOMs emission, observed to infer the chemical composition of star forming regions, can be used also as unique tracer to image protostellar disks on Solar System scales.

Supernovae from core-collapse of massive stars drive shocks into the molecular clouds from which the stars formed. Such shocks affect future star formation from the molecular clouds, and the fast-moving, dense gas with compressed magnetic fields is associated with enhanced cosmic rays. This paper presents new theoretical modeling, using the Paris-Durham shock model, and new observations, using the Stratospheric Observatory for Infrared Astronomy (SOFIA), of the H$_2$ S(5) pure rotational line from molecular shocks in the supernova remnant IC443. We generate MHD models for non-steady-state shocks driven by the pressure of the IC443 blast wave into gas of densities $10^3$ to $10^5$ cm$^{-3}$. We present the first detailed derivation of the shape of the velocity profile for emission from H$_2$ lines behind such shocks, taking into account the shock age, preshock density, and magnetic field. For preshock densities $10^3$-$10^5$ cm$^{-3}$, the the predicted shifts of line centers, and the line widths, of the H$_2$ lines range from 20-2, and 30-4 km/s, respectively. The a priori models are compared to the observed line profiles, showing that clumps C and G can be explained by shocks into gas with density 10$^3$ to $2\times 10^4$ cm$^{-3}$ and strong magnetic fields. For clump B2 (a fainter region near clump B), the H$_2$ spectrum requires a J-type shock into moderate density (~100 cm$^{-3}$) with the gas accelerated to 100 km/s from its pre-shock location. Clump B1 requires both a magnetic-dominated C-type shock (like for clumps C and G) and a J-type shock (like for clump B1) to explain the highest observed velocities. The J-type shocks that produce high-velocity molecules may be locations where the magnetic field is nearly parallel to the shock velocity, which makes it impossible for a C-type shock (with ions and neutrals separated) to form.

NGC 5128 (Centaurus A) is one of the best targets to study AGN-feedback in the local Universe. Optical filaments located at 16 kpc from the galaxy along the radio jet direction show recent star formation, likely triggered by the interaction of the jet with an HI shell. A large reservoir of molecular gas has been discovered outside the HI. In this reservoir, lies the Horseshoe complex: a filamentary structure seen in CO with ALMA and in Halpha with MUSE. The ionised gas is mostly excited by shocks, with only a minor contribution of star formation. We used the Atacama Pathfinder EXperiment (APEX) to observe the 12CO(3-2) and 12CO(4-3) transitions, as well as dense gas tracers in the Horseshoe complex. 12CO(3-2) and 12CO(4-3) are detected for the first time in the northern filaments of Centaurus A, with integrated intensity line ratios R32~0.2 and R43~0.1, compared to the 12CO(1-0) emission. We also derived a line ratio R21~0.6, based on the previous 12CO(2-1) observations of Salom\'e et al. (2016). We used the non-LTE radiative transfer code RADEX and determined that the molecular gas in this region has a temperature of 55-70 K and densities between 2-6x10^2 cm^-3. Such densities are also in agreement with results from the Paris-Durham shock code that predicts a post-shock density of a few 100 cm^-3. However, we need more observations of emission lines at a better angular resolution in order to place tighter constraints on our radiative models, whether they are used as a stand-alone tool (LVG codes) or combined with a shock model.

OMC-2/3 is one of the nearest embedded cluster-forming region that includes intermediate-mass protostars at early stages of evolution. A previous CO (3--2) mapping survey towards this region revealed outflow activity related to sources at different evolutionary phases. The present work aims to study the warm gas in the high-velocity emission from several outflows found in CO (3--2) emission by previous observations, determine their physical conditions, and make comparison with previous results in low-mass star-forming regions. We use the CHAMP+ heterodyne array on the APEX telescope to map the CO (6--5) and CO (7--6) emission in the OMC-2 FIR 6 and OMC-3 MMS 1-6 regions, and to observe $^{13}$CO (6--5) at selected positions. We analyze these data together with previous CO (3--2) observations. In addition, we mapped the SiO (5--4) emission in OMC-2 FIR 6. The CO (6--5) emission was detected in most of the outflow lobes in the mapped regions, while the CO (7--6) was found mostly from the OMC-3 outflows. In the OMC-3 MMS 5 outflow, a previously undetected extremely high velocity gas was found in CO (6--5). This extremely high velocity emission arises from the regions close to the central object MMS 5. Radiative transfer models revealed that the high-velocity gas from MMS 5 outflow consists of gas with $n_{\rm H_2}=$10$^4$--10$^{5}$ cm$^{-3}$ and $T>$ 200 K, similar to what is observed in young Class 0 low-mass protostars. For the other outflows, values of $n_{\rm H_2}>$10$^{4}$ cm$^{-3}$ were found. The physical conditions and kinematic properties of the young intermediate-mass outflows presented here are similar to those found in outflows from Class 0 low-mass objects. Due to their excitation requirements, mid$-J$ CO lines are good tracers of extremely high velocity gas in young outflows likely related to jets.

High-mass analogues of low-mass prestellar cores are searched for to constrain the models of high-mass star formation. Several high-mass cores, at various evolutionary stages, have been recently identified towards the massive star-forming region W43-MM1 and amongst them a high-mass prestellar core candidate. We aim to characterise the chemistry in this high-mass prestellar core candidate, referred to as W43-MM1 core #6, and its environment. Using ALMA high-spatial resolution data of W43-MM1, we have studied the molecular content of core #6 and a neighbouring high-mass protostellar core, referred to as #3, which is similar in size and mass to core #6. We first subtracted the continuum emission using a method based on the density distribution of the intensities on each pixel. Then, from the distribution of detected molecules, we identified the molecules centred on the prestellar core candidate (core #6) and those associated to shocks related to outflows and filament formation. Then we constrained the column densities and temperatures of the molecules detected towards the two cores. While core #3 appears to contain a hot core with a temperature of about 190 K, core #6 seems to have a lower temperature in the range from 20 K to 90 K from a rotational diagram analysis. We have considered different source sizes for core #6 and the comparison of the abundances of the detected molecules towards the core with various interstellar sources shows that it is compatible with a core of size 1000 au with $T = 20-90$ K or a core of size 500 au with $T \sim 80$ K. Core #6 of W43-MM1 remains one of the best high-mass prestellar core candidates even if we cannot exclude that it is at the very beginning of the protostellar phase of high-mass star formation.

With an aim of probing the physical conditions and excitation mechanisms of warm molecular gas in individual star-forming regions, we performed Herschel SPIRE FTS observations of 30 Doradus in the LMC. In our FTS observations, important FIR cooling lines in the ISM, including CO J=4-3 to 13-12, [CI] 370 micron, and [NII] 205 micron, were clearly detected. In combination with ground-based CO data, we then constructed CO spectral line energy distributions (SLEDs) on 10 pc scales over a 60 pc x 60 pc area and found that the shape of the observed CO SLEDs considerably changes across 30 Doradus, e.g., the peak transition varies from J=6-5 to 10-9, while the slope characterized by the high-to-intermediate J ratio ranges from 0.4 to 1.8. To examine the source(s) of these variations in CO transitions, we analyzed the CO observations, along with [CII] 158 micron, [CI] 370 micron, [OI] 145 micron, H2 0-0 S(3), and FIR luminosity data, using state-of-the-art models of PDRs and shocks. Our detailed modeling showed that the observed CO emission likely originates from highly-compressed (thermal pressure ~ 1e7-1e9 K cm-3) clumps on 0.7-2 pc scales, which could be produced by either UV photons (UV radiation field ~ 1e3-1e5 Mathis fields) or low-velocity C-type shocks (pre-shock medium density ~ 1e4-1e6 cm-3 and shock velocity ~ 5-10 km s-1). Considering the stellar content in 30 Doradus, however, we tentatively excluded the stellar origin of CO excitation and concluded that low-velocity shocks driven by kpc scale processes (e.g., interaction between the Milky Way and the Magellanic Clouds) are likely the dominant source of heating for CO. The shocked CO-bright medium was then found to be warm (temperature ~ 100-500 K) and surrounded by a UV-regulated low pressure component (a few (1e4-1e5) K cm-3) that is bright in [CII] 158 micron, [CI] 370 micron, [OI] 145 micron, and FIR dust continuum emission.

Spectrally-resolved observations of three pure rotational lines of H$_2$, conducted with the EXES instrument on SOFIA toward the classic bow shock HH7, reveal systematic velocity shifts between the S(5) line of ortho-H$_2$ and the two para-H$_2$ lines [S(4) and S(6)] lying immediately above and below it on the rotational ladder. These shifts, reported here for the first time, imply that we are witnessing the conversion of para-H$_2$ to ortho-H$_2$ within a shock wave driven by an outflow from a young stellar object. The observations are in good agreement with the predictions of models for non-dissociative, C-type molecular shocks. They provide a clear demonstration of the chemical changes wrought by interstellar shock waves, in this case the conversion of para-H$_2$ to ortho-H$_2$ in reactive collisions with atomic hydrogen, and provide among the most compelling evidence yet obtained for C-type shocks in which the flow velocity changes continuously.

Aims. The goal of the paper is to present a detailed study of the propagation of low velocity (5 to 25 km s-1) stationary molecular shocks in environments illuminated by an external ultraviolet (UV) radiation field. In particular, we intend to show how the structure, dynamics, energetics, and chemical properties of shocks are modified by UV photons and to estimate how efficiently shocks can produce line emission. Methods. We implemented several key physico-chemical processes in the Paris-Durham shock code to improve the treatment of the radiative transfer and its impact on dust and gas particles. We propose a new integration algorithm to find the steady-state solutions of magnetohydrodynamics equations in a range of parameters in which the fluid evolves from a supersonic to a subsonic regime. We explored the resulting code over a wide range of physical conditions, which encompass diffuse interstellar clouds and hot and dense photon-dominated regions (PDR). Results. We find that C-type shock conditions cease to exist as soon as G0 > 0.2 (nH/cm-3)^1/2. Such conditions trigger the emergence of another category of stationary solutions, called C*-type and CJ-type shocks, in which the shocked gas is momentarily subsonic along its trajectory. These solutions are shown to be unique for a given set of physical conditions and correspond to dissipative structures in which the gas is heated up to temperatures comprised between those found in C-type and adiabatic J-type shocks. High temperatures combined with the ambient UV field favour the production or excitation of a few molecular species to the detriment of others, hence leading to specific spectroscopic tracers such as rovibrational lines of H2 and rotational lines of CH+. Unexpectedly, the rotational lines of CH+ may carry as much as several percent of the shock kinetic energy.

We present a characterization of the binary protostar system that is forming within a dense core in the isolated dark cloud BHR71. The pair of protostars, IRS1 and IRS2, are both in the Class 0 phase, determined from observations that resolve the sources from 1 um out to 250 um and from 1.3 mm to 1.3cm. The resolved observations enable the luminosities of IRS1 and IRS2 to be independently measured (14.7 and 1.7L_sun, respectively), in addition to the bolometric temperatures 68~K, and 38~K, respectively. The surrounding core was mapped in NH3 (1,1) with the Parkes radio telescope, and followed with higher-resolution observations from ATCA in NH3 (1,1) and 1.3cm continuum. The protostars were then further characterized with ALMA observations in the 1.3~mm continuum along with N2D+ (J=3-2), 12CO, 13CO, and C18O (J=2-1) molecular lines. The Parkes observations find evidence for a velocity gradient across the core surrounding the two protostars, while ATCA reveals more complex velocity structure toward the protostars within the large-scale gradient. The ALMA observations then reveal that the two protostars are at the same velocity in C18O, and N2H+ exhibits a similar velocity structure as NH3. However, the C18O kinematics reveal that the rotation on scales $<$1000~AU around IRS1 and IRS2 are in opposite directions. Taken with the lack of a systematic velocity difference between the pair, it is unlikely that their formation resulted from rotational fragmentation. We instead conclude that the binary system most likely formed via turbulent fragmentation of the core.

Aims. To constrain the physical processes that lead to the birth of high-mass stars it is mandatory to study the very first stages of their formation. We search for high-mass analogs of low-mass prestellar cores in W43-MM1. Methods. We conducted a 1.3 mm ALMA mosaic of the complete W43-MM1 cloud, which has revealed numerous cores with ~ 2000 au FWHM sizes. We investigated the nature of cores located at the tip of the main filament, where the clustering is minimum. We used the continuum emission to measure the core masses and the $^{13}$CS(5-4) line emission to estimate their turbulence level. We also investigated the prestellar or protostellar nature of these cores by searching for outflow signatures traced by CO(2-1) and SiO(5-4) line emission, and for molecular complexity typical of embedded hot cores. Results. Two high-mass cores of ~ 1300 au diameter and ~ $60~M_\odot$ mass are observed to be turbulent but gravitationally bound. One drives outflows and is associated with a hot core. The other core, W43-MM1\#6, does not yet reveal any star formation activity and thus is an excellent high-mass prestellar core candidate.

The HH 212 protostellar system, in Orion B, has been mapped thanks to ALMA-Band 7 Cycle 1 and Cycle 4 observations of dueterated water (HDO) and acetaldehyde (CH$_3$CHO) emission with an angular resolution down to $\sim$0.15 arcsec (60 au). Many emission lines due to 14 CH$_3$CHO and 1 HDO transitions at high excitation ($E_{\rm u}$ between 163 K and 335 K) have been imaged in the inner $\sim$ 70 au region. The local thermal equilibrium analysis of the CH$_3$CHO emission leads to a temperature of 78$\pm$14 K and a column density of 7.6$\pm$3.2 $\times$ 10$^{15}$ cm$^{-2}$, which, when $N_{\rm H_2}$ of 10$^{24}$ cm$^{-2}$ is assumed, leads to an abundance of $X_{\rm CH_3CHO}$ $\simeq$ 8 $\times$ 10$^{-9}$. The large velocity gradient analysis of the HDO emission also places severe constraints on the volume density, n$_{\rm H_2}$ $\geq$ 10$^8$ cm$^{-3}$. The line profiles are 5--7 km s$^{-1}$ wide, and CH$_3$CHO and HDO both show a $\pm$ 2 km s$^{-1}$ velocity gradient over a size of $\sim$ 70 au (blue-shifted emission towards the north-west and red-shifted emission towards the south-east) along the disc equatorial plane, in agreement with what was found so far using other molecular tracers. The kinematics of CH$_3$CHO and HDO are consistent with the occurrence of a centrifugal barrier, that is, the infalling envelope-rotating disc ring, which is chemically enriched through low-velocity accretion shocks. The emission radius is $\sim$ 60 au, in good agreement with what was found before for another interstellar complex organic molecule such as NH$_2$CHO. We support a vertical structure for the centrifugal barrier, suggesting the occurrence of two outflowing, expanding, and rotating rings above and below (of about 40-45 au) the optically thick equatorial disc plane. It is tempting to speculate that these rings could probe the basis of a wind launched from this region.

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

Understanding the processes that determine the stellar Initial Mass Function (IMF) is a critical unsolved problem, with profound implications for many areas of astrophysics. In molecular clouds, stars are formed in cores, gas condensations which are sufficiently dense that gravitational collapse converts a large fraction of their mass into a star or small clutch of stars. In nearby star-formation regions, the core mass function (CMF) is strikingly similar to the IMF, suggesting that the shape of the IMF may simply be inherited from the CMF. Here we present 1.3 mm observations, obtained with ALMA, the world's largest interferometer, of the active star-formation region W43-MM1, which may be more representative of the Galactic-disk regions where most stars form. The unprecedented resolution of these observations reveals, for the first time, a statistically robust CMF at high masses, with a slope that is markedly shallower than the IMF. This seriously challenges our understanding of the origin of the IMF.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

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

When a fast moving star or a protostellar jet hits an interstellar cloud, the surrounding gas gets heated and illuminated: a bow shock is born which delineates the wake of the impact. In such a process, the new molecules that are formed and excited in the gas phase become accessible to observations. In this article, we revisit models of H2 emission in these bow shocks. We approximate the bow shock by a statistical distribution of planar shocks computed with a magnetized shock model. We improve on previous works by considering arbitrary bow shapes, a finite irradiation field, and by including the age effect of non-stationary C-type shocks on the excitation diagram and line profiles of H2. We also examine the dependence of the line profiles on the shock velocity and on the viewing angle: we suggest that spectrally resolved observations may greatly help to probe the dynamics inside the bow shock. For reasonable bow shapes, our analysis shows that low velocity shocks largely contribute to H2 excitation diagram. This can result in an observational bias towards low velocities when planar shocks are used to interpret H2 emission from an unresolved bow. We also report a large magnetization bias when the velocity of the planar model is set independently. Our 3D models reproduce excitation diagrams in BHR71 and Orion bow shocks better than previous 1D models. Our 3D model is also able to reproduce the shape and width of the broad H2 1-0S(1) line profile in an Orion bow shock.

We wish to constrain the possible contribution of a magnetohydrodynamic disk wind (DW) to the HH212 molecular jet. We mapped the flow base with ALMA Cycle 4 at 0.13"$\sim$ 60 au resolution and compared these observations with synthetic DW predictions. We identified, in SO/SO$_2$, a rotating flow that is wider and slower than the axial SiO jet. The broad outflow cavity seen in C$^{34}$S is not carved by a fast wide-angle wind but by this slower agent. Rotation signatures may be fitted by a DW of a moderate lever arm launched out to $\sim$ 40 au with SiO tracing dust-free streamlines from 0.05-0.3 au. Such a DW could limit the core-to-star efficiency to $\leq$50%.

Context: Methanol is thought to be mainly formed during the prestellar phase and its deuterated form keeps memory of the conditions at that epoch. Thanks to the unique combination of high angular resolution and sensitivity provided by ALMA, we wish to measure methanol deuteration in the planet formation region around a Class 0 protostar and to understand its origin. Aims: We mapped both the $^{13}$CH$_3$OH and CH$_2$DOH distribution in the inner regions ($\sim$100 au) of the HH212 system in Orion B. To this end, we used ALMA Cycle 1 and Cycle 4 observations in Band 7 with angular resolution down to $\sim$0.15$"$. Results: We detected 6 lines of $^{13}$CH$_3$OH and 13 lines of CH$_2$DOH with upper level energies up to 438 K in temperature units. We derived a rotational temperature of (171 $\pm$ 52) K and column densities of 7$\times$10$^{16}$ cm$^{-2}$ ($^{13}$CH$_3$OH) and 1$\times$10$^{17}$ cm$^{-2}$ (CH$_2$DOH), respectively. Consequently, the D/H ratio is (2.4 $\pm$ 0.4)$\times$10$^{-2}$, a value lower by an order of magnitude with respect to what was previously measured using single dish telescopes toward protostars located in Perseus. Our findings are consistent with the higher dust temperatures in Orion B with respect to that derived for the Perseus cloud. The emission is tracing a rotating structure extending up to 45 au from the jet axis and elongated by 90 au along the jet axis. So far, the origin of the observed emission appears to be related with the accretion disk. Only higher spatial resolution measurements however, will be able to disentangle between different possible scenarios: disk wind, disk atmosphere, or accretion shocks.

(Abridged) Through spectrally unresolved observations of high-J CO transitions, Herschel-PACS has revealed large reservoirs of warm (300 K) and hot (700 K) molecular gas around low-mass protostars. We aim to shed light on the excitation and origin of the CO ladder observed toward protostars, and on the water abundance in different physical components using spectrally resolved Herschel-HIFI data. Observations are presented of the highly excited CO line J=16-15 with Herschel-HIFI toward 24 low-mass protostellar objects. The spectrally resolved profiles show two distinct velocity components: a broad component with an average FWHM of 20 km/s, and a narrower component with a FWHM of 5 km/s that is often offset from the source velocity. The average rotational temperature over the entire profile, as measured from comparison between CO J=16-15 and 10-9 emission, is ~300 K. A radiative-transfer analysis shows that the average H2O/CO column-density ratio is ~0.02, suggesting a total H2O abundance of ~2x10^-6. Two distinct velocity profiles observed in the HIFI line profiles suggest that the CO ladder observed with PACS consists of two excitation components. The warm component (300 K) is associated with the broad HIFI component, and the hot component (700 K) is associated with the offset HIFI component. The former originates in either outflow cavity shocks or the disk wind, and the latter in irradiated shocks. The ubiquity of the warm and hot CO components suggests that fundamental mechanisms govern the excitation of these components; we hypothesize that the warm component arises when H2 stops being the dominant coolant. In this scenario, the hot component arises in cooling molecular H2-poor gas just prior to the onset of H2 formation. High spectral resolution observations of highly excited CO transitions uniquely shed light on the origin of warm and hot gas in low-mass protostellar objects.

The shock L1157-B1 driven by the low-mass protostar L1157-mm is an unique environment to investigate the chemical enrichment due to molecules released from dust grains. IRAM-30m and Plateau de Bure Interferometer observations allow a census of Si-bearing molecules in L1157-B1. We detect SiO and its isotopologues and, for the first time in a shock, SiS. The strong gradient of the [SiO/SiS] abundance ratio across the shock (from >=180 to ~25) points to a different chemical origin of the two species. SiO peaks where the jet impacts the cavity walls ([SiO/H2] ~ 1e-6), indicating that SiO is directly released from grains or rapidly formed from released Si in the strong shock occurring at this location. In contrast, SiS is only detected at the head of the cavity opened by previous ejection events ([SiS/H2] ~ 2e-8). This suggests that SiS is not directly released from the grain cores but instead should be formed through slow gas-phase processes using part of the released silicon. This finding shows that Si-bearing molecules can be useful to distinguish regions where grains or gas-phase chemistry dominates.

In molecular outflows from forming low-mass protostars, most oxygen is expected to be locked up in water. However, Herschel observations have shown that typically an order of magnitude or more of the oxygen is still unaccounted for. To test if the oxygen is instead in atomic form, SOFIA-GREAT observed the R1 position of the bright molecular outflow from NGC1333-IRAS4A. The [OI] 63 um line is detected and spectrally resolved. From an intensity peak at +15 km/s, the intensity decreases until +50 km/s. The profile is similar to that of high-velocity (HV) H2O and CO 16-15, the latter observed simultaneously with [OI]. A radiative transfer analysis suggests that ~15% of the oxygen is in atomic form toward this shock position. The CO abundance is inferred to be ~10^-4 by a similar analysis, suggesting that this is the dominant oxygen carrier in the HV component. These results demonstrate that a large portion of the observed [OI] emission is part of the outflow. Further observations are required to verify whether this is a general trend.

Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]$_{\rm 63 \mu m}$ emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. We present observations of the OI $^3$P$_1 \rightarrow$ $^3$P$_2$, OH between $^2\Pi_{1/2}$ $J = 3/2$ and $J = 1/2$ at 1837.8 GHz, and CO (16-15) lines with SOFIA-GREAT at three positions in the Cep E outflow: mm (the driving protostar), BI (in the southern lobe), and BII (the terminal position in the southern lobe). The CO line is detected at all three positions. The OI line is detected in BI and BII, whereas the OH line is not detected. In BII, we identify three kinematical components in OI and CO, already detected in CO: the jet, the HH377 terminal bow-shock, and the outflow cavity. The OI column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is where the abundance ratio of OI to CO is the lowest (about 0.2), whereas the jet component is atomic (ratio $\sim$2.7). In the jet, we compare the OI observations with shock models that successfully fit the integrated intensity of 10 CO lines: these models do not fit the OI data. The high intensity of OI emission points towards the propagation of additional dissociative or alternative FUV-irradiated shocks, where the illumination comes from the shock itself. From the sample of low-to-high mass protostellar outflows where similar observations have been performed, the effects of illumination seem to increase with the mass of the protostar.

Protoplanetary disks undergo substantial mass-loss by photoevaporation, a mechanism which is crucial to their dynamical evolution. However, the processes regulating the gas energetics have not been well constrained by observations so far. We aim at studying the processes involved in disk photoevaporation when it is driven by far-UV photons. We present a unique Herschel survey and new ALMA observations of four externally-illuminated photoevaporating disks (a.k.a. proplyds). For the analysis of these data, we developed a 1D model of the photodissociation region (PDR) of a proplyd, based on the Meudon PDR code and computed the far infrared line emission. We successfully reproduce most of the observations and derive key physical parameters, i.e. densities at the disk surface of about $10^{6}$ cm$^{-3}$ and local gas temperatures of about 1000 K. Our modelling suggests that all studied disks are found in a transitional regime resulting from the interplay between several heating and cooling processes that we identify. These differ from those dominating in classical PDRs, i.e. grain photo-electric effect and cooling by [OI] and [CII] FIR lines. This energetic regime is associated to an equilibrium dynamical point of the photoevaporation flow: the mass-loss rate is self-regulated to set the envelope column density at a value that maintains the temperature at the disk surface around 1000 K. From our best-fit models, we estimate mass-loss rates - of the order of $10^{-7}$ $\mathrm{M}_\odot$/yr - that are in agreement with earlier spectroscopic observation of ionised gas tracers. This holds only if we assume an evaporation flow launched from the disk surface at sound speed (supercritical regime). We have identified the energetic regime regulating FUV-photoevaporation in proplyds. This regime could be implemented into models of the dynamical evolution of protoplanetary disks.

In the BHR71 region, two low-mass protostars drive two distinguishable outflows. They constitute an ideal laboratory to investigate the effects of shock chemistry and the mechanisms that led to their formation. We aim to define the morphology of the warm gas component of the BHR 71 outflow and at modelling its shocked component. We present the first far infrared Herschel images of the BHR71 outflow in the CO(14-13), H$_2$O (2$_{21}$-1$_{10}$), H$_2$O (2$_{12}$-1$_{01}$) and [OI] 145 $\mu$m, lines, revealing the presence of several knots of warm, shocked gas associated with fast outflowing gas. In two of these knots we performed a detailed study of the physical conditions by comparing a large set of transitions from several molecules to a grid of shock models. Herschel lines ratios in the outflow knots are quite similar, showing that the excitation conditions of the fast moving gas do not change significantly within the first $\sim$ 0.068 pc of the outflow, apart at the extremity of the southern blue-shifted lobe that is expanding outside the molecular cloud. Rotational diagram, spectral line profile and LVG analysis of the CO lines in knot A show the presence of two gas components: one extended, cold ($T\sim$80 K) and dense ($n$(H$_2$) = 3$\times$10$^5$-4$\times$10$^6$ cm$^{-3}$) and another compact (18 arcsec), warm ($T$ = 1700-2200 K) with slightly lower density ($n$(H$_2$) = (2-6)$\times$10$^4$ cm$^{-3}$). In the two brightest knots (where we performed shock modelling) we found that H$_2$ and CO are well fitted with non-stationary (young) shocks. These models, however, significantly underestimate the observed fluxes of [OI] and OH lines, but are not too far off those of H$_2$O, calling for an additional, possibly dissociative, J-type shock component. Our modelling indirectly suggests that an additional shock component exists, possibly a remnant of the primary jet

The relative populations in rotational transitions of CO can be useful for inferring gas conditions and excitation mechanisms at work in the interstellar medium. We present CO emission lines from rotational transitions observed with Herschel/HIFI in the star-forming cores Orion S, Orion KL, Sgr B2(M), and W49N. Integrated line fluxes from these observations are combined with those from Herschel/PACS observations of the same sources to construct CO spectral line energy distributions (SLEDs) from $5\leq J_u\leq48$. These CO SLEDs are compared to those reported in other galaxies, with the intention of empirically determining which mechanisms dominate excitation in such systems. We find that CO SLEDs in Galactic star-forming cores cannot be used to reproduce those observed in other galaxies, although the discrepancies arise primarily as a result of beam filling factors. The much larger regions sampled by the Herschel beams at distances of several Mpc contain significant amounts of cooler gas which dominate the extragalactic CO SLEDs, in contrast to observations of Galactic star-forming regions which are focused specifically on cores containing primarily hot molecular gas.

The mechanisms leading to the formation of disks around young stellar objects (YSOs) and to the launching of the associated jets are crucial to the understanding of the earliest stages of star and planet formation. HH 212 is a privileged laboratory to study a pristine jet-disk system. Therefore we investigate the innermost region ($<100$ AU) around the HH 212-MM1 protostar through ALMA band\,7 observations of methanol. The 8 GHz bandwidth spectrum towards the peak of the continuum emission of the HH 212 system reveals at least 19 transitions of methanol. Several of these lines (among which several vibrationally excited lines in the v$_{\rm t}=1,2$ states) have upper energies above 500 K. They originate from a compact ($<135$ AU in diameter), hot ($\sim 295$ K) region elongated along the direction of the SiO jet. We performed a fit in the $uv$ plane of various velocity channels of the strongest high-excitation lines. The blue- and red-shifted velocity centroids are shifted roughly symmetrically on either side of the jet axis, indicating that the line-of-sight velocity beyond 0.7 km s$^{-1}$ from systemic is dominated by rotational motions. The velocity increases moving away from the protostar further indicating that the emission of methanol is not associated with a Keplerian disk or rotating-infalling cavity, and it is more likely associated with outflowing gas. We speculate that CH$_3$OH traces a disk wind gas accelerated at the base. The launching region would be at a radius of a few astronomical units from the YSO.

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 Herschel SPIRE Fourier Transform Spectrometer (FTS) observations of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In our observations, a number of far-infrared cooling lines including CO(4-3) to CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected. With an aim of investigating the physical conditions and excitation processes of molecular gas, we first construct CO spectral line energy distributions (SLEDs) on 10 pc scales by combining the FTS CO transitions with ground-based low-J CO data and analyze the observed CO SLEDs using non-LTE radiative transfer models. We find that the CO-traced molecular gas in N159W is warm (kinetic temperature of 153-754 K) and moderately dense (H2 number density of (1.1-4.5)e3 cm-3). To assess the impact of the energetic processes in the interstellar medium on the physical conditions of the CO-emitting gas, we then compare the observed CO line intensities with the models of photodissociation regions (PDRs) and shocks. We first constrain the properties of PDRs by modelling Herschel observations of [OI] 145, [CII] 158, and [CI] 370 micron fine-structure lines and find that the constrained PDR components emit very weak CO emission. X-rays and cosmic-rays are also found to provide a negligible contribution to the CO emission, essentially ruling out ionizing sources (ultraviolet photons, X-rays, and cosmic-rays) as the dominant heating source for CO in N159W. On the other hand, mechanical heating by low-velocity C-type shocks with ~10 km/s appears sufficient enough to reproduce the observed warm CO.

Aims: Using the unprecedented combination of high resolution and sensitivity offered by ALMA, we aim to investigate whether and how hot corinos, circumstellar disks, and ejected gas are related in young solar-mass protostars. Methods: We observed CH$_3$CHO and deuterated water (HDO) high-excitation ($E_{\rm u}$ up to 335 K) lines towards the Sun-like protostar HH212--MM1. Results: For the first time, we have obtained images of CH$_3$CHO and HDO emission in the inner $\simeq$ 100 AU of HH212. The multifrequency line analysis allows us to contrain the density ($\geq$ 10$^{7}$ cm$^{-3}$), temperature ($\simeq$ 100 K), and CH$_3$CHO abundance ($\simeq$ 0.2--2 $\times$ 10$^{-9}$) of the emitting region. The HDO profile is asymmetric at low velocities ($\leq$ 2 km s$^{-1}$ from $V_{\rm sys}$). If the HDO line is optically thick, this points to an extremely small ($\sim$ 20--40 AU) and dense ($\ge$ 10$^{9}$ cm$^{-3}$) emitting region. Conclusions: We report the first detection of a hot corino in Orion. The HDO asymmetric profile indicates a contribution of outflowing gas from the compact central region, possibly associated with a dense disk wind.

Shocks are ubiquitous in the interstellar medium of galaxies, where they contribute to the energetic balance and to the cycle of matter, and where they are thought to be the primary sites for cosmic rays acceleration. Most of the time: in jets and outflows, supernova remnants, or colliding flows, they are linked with star formation. The study of shocks is hence a powerful tool to probe the evolution of the interstellar medium and to better understand star formation. To these aims, the most precise observations must be compared with the most precise models of shocks. The SOFIA/GREAT instrument represents a powerful observational tool to support our progresses, as it allows to observe numerous shock tracers in the far-infrared range.

The fine-structure line of [OI] at 63micron is an important diagnostic tool in different fields of astrophysics. However, our knowledge of this line relies on observations with low spectral resolution, and the real contribution of each component (PDR, jet) in complex environment of star-forming regions (SFRs) is poorly understood. We investigate the contribution of jet and PDR emission, and of absorption to the [OI]63micron line towards the ultra-compact H{\sc ii} region G5.89--0.39 and study its far-IR line luminosity in different velocity regimes through [OI], [CII], CO, OH, and H2O. We mapped G5.89--0.39 in [OI] and in CO(16--15) with the GREAT receiver onboard SOFIA. We observed the central position of the source in the OH^2\Pi_{3/2}, J=5/2\toJ=3/2 and ^2\Pi_{1/2}, J=3/2\toJ=1/2 lines. These data were complemented with APEX CO(6-5) and CO(7-6) and HIFI maps and single-pointing observations in [CII], H2O, and HF. The [OI] spectra in G5.89--0.39 are severely contaminated by absorptions from the envelope and from different clouds along the line of sight. Emission is detected only at HV, clearly associated with the compact north-south outflows traced by extremely HV low-J CO. The mass-loss rate and energetics of derived from [OI] agree well with estimates from CO, suggesting that the molecular outflows in G5.89--0.39 are driven by the jet system seen in [OI]. The far-IR line luminosity of G5.89--0.39 is dominated by [OI] at HV; the second coolant in this velocity regime is CO, while [CII], OH and H2O are minor contributors to the cooling in the outflow. Our study shows the importance of spectroscopically resolved data of [OI]63micron for using this line as diagnostic of SFRs. While this was not possible until now, the GREAT receiver onboard SOFIA has recently opened the possibility of detailed studies of this line to investigate its potential for probing different environments.

The interstellar medium (ISM) is constantly evolving due to unremitting injection of energy in various forms. Energetic radiation transfers energy to the ISM: from the UV photons, emitted by the massive stars, to X- and $\gamma$-ray ones. Cosmic rays are another source of energy. Finally, mechanical energy is injected through shocks or turbulence. Shocks are ubiquitous in the interstellar medium of galaxies. They are associated to star formation (through jets and bipolar outflows), life (via stellar winds), and death (in AGB stellar winds or supernovae explosion). The dynamical processes leading to the formation of molecular clouds also generate shocks where flows of interstellar matter collide. Because of their ubiquity, the study of interstellar shocks is also a useful probe to the other mechanisms of energy injection in the ISM. This study must be conducted in order to understand the evolution of the ISM as a whole, and to address various questions: what is the peculiar chemistry associated to shocks, and what is their contribution to the cycle of matter in galaxies ? What is the energetic impact of shocks on their surroundings on various scales, and hence what is the feedback of stars on the galaxies ? What are the scenarios of star formation, whether this star formation leads to the propagation of shocks, or whether it is triggered by shock propagation ? What is the role of shocks in the acceleration of cosmic rays ? Can they shed light on their composition and diffusion processes ? In order to progress on these questions, it is paramount to interpret the most precise observations with the most precise models of shocks. From the observational point of view, the James Webb Space Telescope represents a powerful tool to better address the above questions, as it will allow to observe numerous shock tracers in the infrared range at an unprecedented spatial and spectral resolution.

OH is a key molecule in H2O chemistry, a valuable tool for probing physical conditions, and an important contributor to the cooling of shock regions. OH participates in the re-distribution of energy from the protostar towards the surrounding ISM. Our aim is to assess the origin of the OH emission from the Cepheus A massive star-forming region and to constrain the physical conditions prevailing in the emitting gas. We thus want to probe the processes at work during the formation of massive stars. We present spectrally resolved observations of OH towards the outflows of Cepheus A with the GREAT spectrometer onboard the SOFIA telescope. Three triplets were observed at 1834.7 GHz, 1837.8 GHz, and 2514.3 GHz (163.4, 163.1, and 119.2 microns), at angular resolutions of 16.3", 16.3", and 11.9", respectively. We present the CO (16-15) spectrum at the same position. We compared the integrated intensities in the redshifted wings to shock models. The two triplets near 163 microns are detected in emission with blending hyperfine structure unresolved. Their profiles and that of CO can be fitted by a combination of 2 or 3 Gaussians. The observed 119.2 microns triplet is seen in absorption, since its blending hyperfine structure is unresolved, but with three line-of-sight components and a blueshifted emission wing consistent with that of the other lines. The OH line wings are similar to those of CO, suggesting that they emanate from the same shocked structure. Under this common origin assumption, the observations fall within the model predictions and within the range of use of our model only if we consider that four shock structures are caught in our beam. Our comparisons suggest that the observations might be consistently fitted by a J-type model with nH > 1e5 cm-3, v > 20 km/s, and with a filling factor of ~1. Such a high density is generally found in shocks associated to high-mass protostars.

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

After the death of a runaway massive star, its supernova shock wave interacts with the bow shocks produced by its defunct progenitor, and may lose energy, momentum, and its spherical symmetry before expanding into the local interstellar medium (ISM). We investigate whether the initial mass and space velocity of these progenitors can be associated with asymmetric supernova remnants. We run hydrodynamical models of supernovae exploding in the pre-shaped medium of moving Galactic core-collapse progenitors. We find that bow shocks that accumulate more than about 1.5 Mo generate asymmetric remnants. The shock wave first collides with these bow shocks 160-750 yr after the supernova, and the collision lasts until 830-4900 yr. The shock wave is then located 1.35-5 pc from the center of the explosion, and it expands freely into the ISM, whereas in the opposite direction it is channelled into the region of undisturbed wind material. This applies to an initially 20 Mo progenitor moving with velocity 20 km/s and to our initially 40 Mo progenitor. These remnants generate mixing of ISM gas, stellar wind and supernova ejecta that is particularly important upstream from the center of the explosion. Their lightcurves are dominated by emission from optically-thin cooling and by X-ray emission of the shocked ISM gas. We find that these remnants are likely to be observed in the [OIII] lambda 5007 spectral line emission or in the soft energy-band of X-rays. Finally, we discuss our results in the context of observed Galactic supernova remnants such as 3C391 and the Cygnus Loop.

To investigate the disk formation and jet launch in protostars is crucial to comprehend the earliest stages of star and planet formation. We aim to constrain the properties of the molecular jet and the disk of the HH 212 protostellar system at unprecedented angular scales through ALMA observations of sulfur-bearing molecules, SO 9(8)-8(7), SO 10(11)-10(10), SO2 8(2,6)-7(1,7). SO 9(8)-8(7) and SO2 8(2,6)-7(1,7) show broad velocity profiles. At systemic velocity they probe the circumstellar gas and the cavity walls. Going from low to high blue-/red-shifted velocities the emission traces the wide-angle outflow and the fast (~100-200 km/s) and collimated (~90 AU) molecular jet revealing the inner knots with timescales <50 years. The jet transports a mass loss rate >0.2-2e-6 Msun/yr, implying high ejection efficiency (>0.03-0.3). The SO and SO2 abundances in the jet are ~1e-7-1e-6. SO 10(11)-10(10) emission is compact and shows small-scale velocity gradients indicating that it originates partly from the rotating disk previously seen in HCO+ and C17O, and partly from the base of the jet. The disk mass is >0.002-0.013 Msun, and the SO abundance in the disk is ~1e-8-1e-7. SO and SO2 are effective tracers of the molecular jet in the inner few hundreds AU from the protostar. Their abundances indicate that 1% - 40% of sulfur is in SO and SO2 due to shocks in the jet/outflow and/or to ambipolar diffusion at the wind base. The SO abundance in the disk is 3-4 orders of magnitude larger than in evolved protoplanetary disks. This may be due to an SO enhancement in the accretion shock at the envelope-disk interface or in spiral shocks if the disk is partly gravitationally unstable.

Supernova remnants (SNRs) are considered as being the sources of galactic cosmic rays. In order to understand the origin, acceleration, and composition of these cosmic rays, detailed knowledge of the physical conditions in the local interstellar medium is needed. The shock interaction of SNRs with molecular clouds that gives rise to strong molecular emission in the far-IR and sub-mm wavelength regimes can be used as a highly valuable tracer of these conditions. The application of MHD shock models in the interpretation of the resulting line emission can yield information on the energetic and chemical impact of supernova remnants. We have mapped two regions in the supernova remnant W44 with the APEX telescope in ${}^{12}$CO (3-2), (4-3), (6-5), (7-6) and ${}^{13}$CO (3-2). The extraction of integrated intensities on five different positions, corresponding to local maxima of CO emission, allows to compare these intensities to the outputs of a grid of models, which combine an MHD shock code with a radiative transfer module based on the 'large velocity gradient' approximation. We find that the observed CO line emission is compatible with non-stationary shocks and a pre-shock density of $10^4$ cm${}^{-3}$. Our models furthermore allow to constrain shock ages, velocities, the pre-shock magnetic field strength components perpendicular to the line-of-sight, and the full ladder of CO transitions. Finally, our analysis can be used to estimate the contribution of such SNRs to, e.g. the galactic energy balance and the momentum-injection into the surrounding interstellar medium.

The formation of massive stars play a crucial role in galaxies from numerous points of view. The protostar generates a strong ultraviolet radiation field that ionizes its surroundings, and it drives powerful shock waves in the neighbouring medium in the form of jets and bipolar outflows, whose structure can be partially organized by local, strong magnetic field. Such an ejection activity locally modifies the interstellar chemistry, contributing to the cycle of matter. It also significantly participates to the energetic balance of galaxies. In the latter stages of massive star formation, the protostar is surrounded by an ultra-compact HII region, and irradiates its bipolar outflows, where an intrinsically strong magnetic field structure is associated to the generally high densities. In the HII region, or in the bipolar outflows, the question of in situ cosmic rays acceleration can then be raised by the simultaneous presence of strong magnetic fields, significant ionization of the matter, and mechanical energy available in large quantities. In this contribution, we will only summarise the results of our study of potential in situ cosmic rays acceleration the irradiated shocks in the W28 A2 massive star forming region, based on observations from the APEX, IRAM 30m, and Herschel telescopes, and presented in Gusdorf et al., recently submitted to A&A.

During the formation of a star, material is ejected along powerful jets that impact the ambient material. This outflow regulates star formation by e.g. inducing turbulence and heating the surrounding gas. Understanding the associated shocks is therefore essential to the study of star formation. We present comparisons of shock models with CO, H2, and SiO observations in a 'pure' shock position in the BHR71 bipolar outflow. These comparisons provide an insight into the shock and pre-shock characteristics, and allow us to understand the energetic and chemical feedback of star formation on Galactic scales. New CO (Jup = 16, 11, 7, 6, 4, 3) observations from the shocked regions with the SOFIA and APEX telescopes are presented and combined with earlier H2 and SiO data (from the Spitzer and APEX telescopes). The integrated intensities are compared to a grid of models that were obtained from a magneto-hydrodynamical shock code which calculates the dynamical and chemical structure of these regions combined with a radiative transfer module based on the 'large velocity gradient' approximation. The CO emission leads us to update the conclusions of our previous shock analysis: pre-shock densities of 1e4 cm-3 and shock velocities around 20-25 km s-1 are still constrained, but older ages are inferred ( 4000 years). We evaluate the contribution of shocks to the excitation of CO around forming stars. The SiO observations are compatible with a scenario where less than 4% of the pre-shock SiO belongs to the grain mantles. We infer outflow parameters: a mass of 1.8x1e-2 Msun was measured in our beam, in which a momentum of 0.4 Msun km s-1 is dissipated, for an energy of 4.2x1e43erg. We analyse the energetics of the outflow species by species. Comparing our results with previous studies highlights their dependence on the method: H2 observations only are not sufficient to evaluate the mass of outflows.

Ionized carbon is the main gas-phase reservoir of carbon in the neutral diffuse interstellar medium and its 158 micron fine structure transition [CII] is the most important cooling line of the diffuse interstellar medium (ISM). We combine [CII] absorption and emission spectroscopy to gain an improved understanding of physical conditions in the different phases of the ISM. We present high resolution [CII] spectra obtained with the Herschel/HIFI instrument towards bright dust continuum sources regions in the Galactic plane, probing simultaneously the diffuse gas along the line of sight and the background high-mass star forming regions. These data are complemented by observations of the 492 and 809 GHz fine structure lines of atomic carbon and by medium spectral resolution spectral maps of the fine structure lines of atomic oxygen at 63 and 145 microns with Herschel/PACS. We show that the presence of foreground absorption may completely cancel the emission from the background source in medium spectral resolution data and that high spectral resolution spectra are needed to interpret the [CII] and [OI] emission and the [CII]/FIR ratio. This phenomenon may explain part of the [CII]/FIR deficit seen in external luminous infrared galaxies. The C+ and C excitation in the diffuse gas is consistent with a median pressure of 5900 Kcm-3 for a mean TK ~100 K. The knowledge of the gas density allows us to determine the filling factor of the absorbing gas along the selected lines of sight: the median value is 2.4 %, in good agreement with the CNM properties. The mean excitation temperature is used to derive the average cooling due to C+ in the Galactic plane : 9.5 x 10^{-26} erg/s/H. Along the observed lines of sight, the gas phase carbon abundance does not exhibit a strong gradient as a function of Galacto-centric radius and has a weighted average of C/H = 1.5 +/- 0.4 x 10^{-4}.

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

Molecular outflows are a direct consequence of accretion, and therefore they represent one of the best tracers of accretion processes in the still poorly understood early phases of high-mass star formation. Previous studies suggested that the SiO abundance decreases with the evolution of a massive young stellar object probably because of a decay of jet activity, as witnessed in low-mass star-forming regions. We investigate the SiO excitation conditions and its abundance in outflows from a sample of massive young stellar objects through observations of the SiO(8-7) and CO(4-3) lines with the APEX telescope. Through a non-LTE analysis, we find that the excitation conditions of SiO increase with the velocity of the emitting gas. We also compute the SiO abundance through the SiO and CO integrated intensities at high velocities. For the sources in our sample we find no significant variation of the SiO abundance with evolution for a bolometric luminosity-to-mass ratio of between 4 and 50 $L_\odot/M_\odot$. We also find a weak increase of the SiO(8-7) luminosity with the bolometric luminosity-to-mass ratio. We speculate that this might be explained with an increase of density in the gas traced by SiO. We find that the densities constrained by the SiO observations require the use of shock models that include grain-grain processing. For the first time, such models are compared and found to be compatible with SiO observations. A pre-shock density of $10^5\, $cm$^{-3}$ is globally inferred from these comparisons. Shocks with a velocity higher than 25 km s$^{-1}$ are invoked for the objects in our sample where the SiO is observed with a corresponding velocity dispersion. Our comparison of shock models with observations suggests that sputtering of silicon-bearing material (corresponding to less than 10% of the total silicon abundance) from the grain mantles is occurring.

Aims. The interaction of supernova remnants (SNRs) with molecular clouds gives rise to strong molecular emission in the far-IR and sub-mm wavelength regimes. The application of MHD shock models in the interpretation of this line emission can yield valuable information on the energetic and chemical impact of supernova remnants. Methods. New mapping observations with the APEX telescope in CO (3-2), (4-3), (6-5), (7-6) and 13CO (3-2) towards two regions in the supernova remnant W44 are presented. Integrated intensities are extracted on five different positions, corresponding to local maxima of CO emission. The integrated intensities are compared to the outputs of a grid of models, which combine an MHD shock code with a radiative transfer module based on the large velocity gradient approximation. Results. All extracted spectra show ambient and line-of-sight components as well as blue- and red-shifted wings indicating the presence of shocked gas. Basing the shock model fits only on the highest-lying transitions that unambiguously trace the shock-heated gas, we find that the observed CO line emission is compatible with non-stationary shocks and a pre-shock density of 10^4 cm-3. The ages of the modelled shocks scatter between values of \sim1000 and \sim3000 years. The shock velocities in W44F are found to lie between 20 and 25 km/s, while in W44E fast shocks (30-35 km/s) as well as slower shocks (\sim20 km/s) are compatible with the observed spectral line energy diagrams. The pre-shock magnetic field strength components perpendicular to the line of sight in both regions have values between 100 and 200 \muG. Our best-fitting models allow us to predict the full ladder of CO transitions, the shocked gas mass in one beam as well as the momentum- and energy injection.

We used PdBI observations of SiO (2-1) to investigate the morphology and profile of the SiO emission within several massive dense clumps (MDCs) in Cygnus-X. We find that most molecular outflows are detected in both SiO and CO, although there are some cases of CO outflows with no SiO counterpart. We find a significant amount of narrow line SiO emission that appears to be unrelated to outflows. The fraction of the total SiO luminosity that is not associated with outflows is highly variable in the different MDCs (from 10% to 90%); this might be a problem when extrapolating outflow properties from SiO luminosities without resolving individual outflows. The extent of the narrow SiO emission varies from rather compact (~ 0.03 pc) to widespread (~0.2 pc), and its kinematics often differs from those found by other high-density tracers such as H13CO+. We find that the least centrally concentrated clumps with the least massive protostellar cores have the most widespread narrow SiO emission. In line with previous evidence of SiO emission associated with low-velocity shocks, we propose an evolutionary picture to explain the existence and distribution of narrow SiO line profiles. In this scenario, the least centrally condensed MDCs are at an early stage where the SiO emission traces shocks from the large-scale collapse of material onto the MDC (e.g. CygX-N40). As the MDC collapses, the SiO emission becomes more confined to the close surroundings of cores, tracing the post-shock material from the infalling MDC against the dense cores (e.g. CygX-N3, N12, and N48). At later stages, when single massive protostars are formed, the SiO luminosity is largely dominated by powerful outflows, and the weaker narrow component shows perhaps the last remnants of the initial collapse (e.g. CygX-N53 and N63).

Context: Because it is viewed simply edge-on, the HH212 protostellar system is an ideal laboratory for studying the interplay of infall, outflow, and rotation in the earliest stages of low-mass star formation. Aims: We wish to exploit the unmatched combination of high angular resolution, high sensitivity, high-imaging fidelity, and spectral coverage provided by ALMA to shed light on the complex kinematics of the innermost central regions of HH212. Methods: We mapped the inner 10" (4500 AU) of the HH212 system at about 0.5 arcsec resolution in several molecular tracers and in the 850 $\mu$m dust continuum using the ALMA interferometer in band 7 in the extended configuration of the Early Science Cycle 0 operations. Results: Within a single ALMA spectral set-up, we simultaneously identify all the crucial ingredients known to be involved in the star formation recipe namely: (i) the fast, collimated bipolar SiO jet driven by the protostar, (ii) the large-scale swept-up CO outflow, (iii) the flattened rotating and infalling envelope, with bipolar cavities carved by the outflow (in C$^{17}$O(3--2)), and (iv) a rotating wide-angle flow that fills the cavities and surrounding the axial jet (in C$^{34}$S(7--6)). In addition, the compact high-velocity C$^{17}$O emission ($\pm$ 1.9--3.5 km s$^{-1}$ from systemic) shows a velocity gradient along the equatorial plane consistent with a rotating disk of about 0.2 arcsec = 90 AU around a $\simeq 0.3 \pm 0.1 M_{\rm \odot}$ source. The rotating disk is possibly Keplerian. Conclusions: HH212 is the third Class 0 protostar with possible signatures of a Keplerian disk of radius $\geq 30 AU$. The warped geometry in our CS data suggests that this large keplerian disk might result from misaligned magnetic and rotation axes during the collapse phase. The wide-angle CS flow suggests that disk winds may be present in this source.

Context: Star formation efficiency (SFE) theories are currently based on statistical distributions of turbulent cloud structures and a simple model of star formation from cores. They remain poorly tested, especially at the highest densities. Aims: We investigate the effects of gas density on the SFE through measurements of the core formation efficiency (CFE). With a total mass of $\sim2\times10^4$ M$_\odot$, the W43-MM1 ridge is one of the most convincing candidate precursor of starburst clusters and thus one of the best place to investigate star formation. Methods: We used high-angular resolution maps obtained at 3 mm and 1 mm within W43-MM1 with the IRAM Plateau de Bure Interferometer to reveal a cluster of 11 massive dense cores (MDCs), and, one of the most massive protostellar cores known. An Herschel column density image provided the mass distribution of the cloud gas. We then measured the 'instantaneous' CFE and estimated the SFE and the star formation rate (SFR) within subregions of the W43-MM1 ridge. Results: The high SFE found in the ridge ($\sim$6% enclosed in $\sim$8 pc$^3$) confirms its ability to form a starburst cluster. There is however a clear lack of dense cores in the northern part of the ridge, which may be currently assembling. The CFE and the SFE are observed to increase with volume gas density while the SFR steeply decreases with the virial parameter, $\alpha_{vir}$. Statistical models of the SFR may well describe the outskirts of the W43-MM1 ridge but struggle to reproduce its inner part, which corresponds to measurements at low $\alpha_{vir}$. It may be that ridges do not follow the log-normal density distribution, Larson relations, and stationary conditions forced in the statistical SFR models.

Supernovae constitute a critical source of energy input to the interstellar medium (ISM). In this short review, we focus on their latest phase of evolution, the supernova remnants (SNRs). We present observations of three old SNRs that have reached the phase where they interact with the ambient ISM: W28, IC443, and 3C391. We show that such objects make up clean laboratories to constrain the physical and chemical processes at work in molecular shock environments. Our studies subsequently allow us to quantify the impact of SNRs on their environment in terms of mass, momentum, and energy dissipation. In turn, their contribution to the energy balance of galaxies can be assessed. Their potential to trigger a further generation of star formation can also be investigated. Finally, our studies provide strong support for the interpretation of gamma-ray emission in SNRs, a crucial step to answer questions related to cosmic rays population and acceleration.

We investigate the physical and chemical processes at work during the formation of a massive protostar based on the observation of water in an outflow from a very young object previously detected in H2 and SiO in the IRAS 17233-3606 region. We estimated the abundance of water to understand its chemistry, and to constrain the mass of the emitting outflow. We present new observations of shocked water obtained with the HIFI receiver onboard Herschel. We detected water at high velocities in a range similar to SiO. We self-consistently fitted these observations along with previous SiO data through a state-of-the-art, one-dimensional, stationary C-shock model. We found that a single model can explain the SiO and H2O emission in the red and blue wings of the spectra. Remarkably, one common area, similar to that found for H2 emission, fits both the SiO and H2O emission regions. This shock model subsequently allowed us to assess the shocked water column density, N(H2O)=1.2x10^{18} cm^{-2}, mass, M(H2O)=12.5 M_earth, and its maximum fractional abundance with respect to the total density, x(H2O)=1.4x10^{-4}. The corresponding water abundance in fractional column density units ranges between 2.5x10^{-5} and 1.2x10^{-5}, in agreement with recent results obtained in outflows from low- and high-mass young stellar objects.

We have investigated the water abundance in shock-heated molecular gas, making use of $Herschel$ measurements of far-infrared CO and H$_2$O line emissions in combination with $Spitzer$ measurements of mid-IR H$_2$ rotational emissions. We present far-infrared line spectra obtained with $Herschel$'s PACS instrument in range spectroscopy mode towards two positions in the protostellar outflow NGC 2071 and one position each in the supernova remnants W28 and 3C391. These spectra provide unequivocal detections, at one or more positions, of 12 rotational lines of water, 14 rotational lines of CO, 8 rotational lines of OH (4 lambda doublets), and 7 fine-structure transitions of atoms or atomic ions. We first used a simultaneous fit to the CO line fluxes, along with H$_2$ rotational line fluxes measured previously by $Spitzer$, to constrain the temperature and density distribution within the emitting gas; and we then investigated the water abundances implied by the observed H$_2$O line fluxes. The water line fluxes are in acceptable agreement with standard theoretical models for nondissociative shocks that predict the complete vaporization of grain mantles in shocks of velocity $v \sim 25$ km/s, behind which the characteristic gas temperature is $\sim 1300$ K and the H$_2$O/CO ratio is 1.2

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.

Context: We investigate non-Zeeman circular polarization and linear polarization levels of up to 1% of $^{12}$CO spectral line emission detected in a shocked molecular clump around the supernova remnant (SNR) IC 443, with the goal of understanding the magnetic field structure in this source. Aims: We examine our polarization results to confirm that the circular polarization signal in CO lines is caused by a conversion of linear to circular polarization, consistent with anisotropic resonant scattering. In this process background linearly polarized CO emission interacts with similar foreground molecules aligned with the ambient magnetic field and scatters at a transition frequency. The difference in phase shift between the orthogonally polarized components of this scattered emission can cause a transformation of linear to circular polarization. Methods: We compared linear polarization maps from dust continuum, obtained with PolKa at APEX, and $^{12}$CO ($J=2\rightarrow1$) and ($J=1\rightarrow0$) from the IRAM 30-m telescope and found no consistency between the two sets of polarization maps. We then reinserted the measured circular polarization signal in the CO lines across the source to the corresponding linear polarization signal to test whether before this linear to circular polarization conversion the linear polarization vectors of the CO maps were aligned with those of the dust. Results: After the flux correction for the two transitions of the CO spectral lines, the new polarization vectors for both CO transitions aligned with the dust polarization vectors, establishing that the non-Zeeman CO circular polarization is due to a linear to circular polarization conversion.

The formation of high-mass stars is tightly linked to that of their parental clouds. We here focus on the high-density parts of W43, a molecular cloud undergoing an efficient event of formation. The cloud structure is studied with a column density image derived from Herschel continuum maps obtained at 70, 160, 250, 350, and 500 micron. We identify two high-column density filamentary clouds, quoted as the W43-MM1 and W43-MM2 ridges, which both account for 1.5x10^4 Msun gas mass above 10^23 cm-2 and within areas of 5 and 14pc^2, respectively. We used the N_2H^+ 1--0 line to confirm that the W43-MM1 and W43-MM2 ridges are structures coherent in velocity and gravitationally bound, despite their large velocity dispersion and ~5 kms line widths. The most intriguing result of the W43 large program is the bright wide-spread SiO 2--1 emission: 1--11 K kms$ stretching an area of ~28 pc^2. Concentrated toward the W43-MM1 and W43-MM2 ridges and their immediate surroundings, it leads to a total luminosity of L_SiO 2-1 ~4 10^4 K kms kpc^2pc^2. We measured a steep relation between the luminosity and velocity extent of the SiO~2--1 lines and propose to use it to distinguish the low-velocity shocks observed here from the more classical high-velocity ones associated with outflows of high-mass young stellar objects. We used state-of-the-art shock models to demonstrate that low-velocity (<10 kms^-1) shocks with a small amount (10%) of Si atoms initially in gas phase or in grain mantles can explain the observed SiO column density in W43. The spatial and velocity overlaps between the ridges high-density gas (n_H2>10^4-10^5 cm^-3) and the shocked SiO gas suggests that ridges could be forming via colliding flows driven by gravity and accompanied by low-velocity shocks. This mechanism may be the initial conditions for the formation of young massive clusters in these ridges.

Studies of molecular outflows in high-mass young stellar objects reveal important information about the formation process of massive stars. We therefore selected the close-by IRAS 17233-3606 massive star-forming region to perform SiO observations with the SMA interferometer in the (5-4) line and with the APEX single-dish telescope in the (5-4) and (8-7) transitions. In this paper, we present a study of one of the outflows in the region, OF1, which shows several properties similar to jets driven by low-mass protostars, such as HH211 and HH212. It is compact and collimated, and associated with extremely high velocity CO emission, and SiO emission at high velocities. We used a state-of-the-art shock model to constrain the pre-shock density and shock velocity of OF1. The model also allowed us to self-consistently estimate the mass of the OF1 outflow. The shock parameters inferred by the SiO modelling are comparable with those found for low-mass protostars, only with higher pre-shock density values, yielding an outflow mass in agreement with those obtained for molecular outflows driven by early B-type young stellar objects. Our study shows that it is possible to model the SiO emission in high-mass star-forming regions in the same way as for shocks from low-mass young stellar objects.

Previous SiO maps of the innermost regions of HH212 set strong constraints on the structure and origin of this jet. They rule out a fast wide-angle wind, and tentatively favor a magneto-centrifugal disk wind launched out to 0.6 AU. We aim to assess the SiO content at the base of the HH212 jet to set an independent constraint on the location of the jet launch zone with respect to the dust sublimation radius. We present the first sub-arcsecond (0"44x0"96) CO map of the HH212 jet base, obtained with the IRAM Plateau de Bure Interferometer. Combining this with previous SiO(5-4) data, we infer the CO(2-1) opacity and mass-flux in the high-velocity jet and arrive at a much tighter lower limit to the SiO abundance than possible from the (optically thick) SiO emission alone. Gas-phase SiO at high velocity contains at least 10% of the elemental silicon if the jet is dusty, and at least 40% if the jet is dust-free, if CO and SiO have similar excitation temperatures. Such a high SiO content is challenging for current chemical models of both dust-free winds and dusty interstellar shocks. Updated chemical models (equatorial dust-free winds, highly magnetized dusty shocks) and observations of higher J CO lines are required to elucidate the dust content and launch radius of the HH212 high-velocity jet.

We report on single-dish radio CO observations towards the inter-galactic medium (IGM) of the Stephan's Quintet (SQ) group of galaxies. Extremely bright mid-IR H2 rotational line emission from warm molecular gas has been detected by Spitzer in the kpc-scale shock created by a galaxy collision. We detect in the IGM CO(1-0), (2-1) and (3-2) line emission with complex profiles, spanning a velocity range of 1000 km/s. The spectra exhibit the pre-shock recession velocities of the two colliding gas systems (5700 and 6700 km/s), but also intermediate velocities. This shows that much of the molecular gas has formed out of diffuse gas accelerated by the galaxy-tidal arm collision. A total H2 mass of 5x10^9 Msun is detected in the shock. The molecular gas carries a large fraction of the gas kinetic energy involved in the collision, meaning that this energy has not been thermalized yet. The turbulent kinetic energy of the H2 gas is at least a factor of 5 greater than the thermal energy of the hot plasma heated by the collision. The ratio between the warm H2 mass derived from Spitzer IRS spectroscopy and the H2 mass derived from CO fluxes is ~0.3 in the IGM of SQ, which is 10-100 times higher than in star-forming galaxies. In the shocked region, the ratio of the PAH-to-CO surface luminosities, commonly used to measure the star formation efficiency of the H2 gas, is lower (up to a factor 75) than the observed values in star-forming galaxies. We suggest that turbulence fed by the galaxy-tidal arm collision maintains a high heating rate within the H2 gas. This interpretation implies that the velocity dispersion on the scale of giant molecular clouds in SQ is one order of magnitude larger than the Galactic value. The high amplitude of turbulence may explain why this gas is not forming stars efficiently. [abridged version]

Context. Observing supernova remnants (SNRs) and modelling the shocks they are associated with is the best way to quantify the energy SNRs re-distribute back into the Interstellar Medium (ISM). Aims. We present comparisons of shock models with CO observations in the F knot of the W28 supernova remnant. These comparisons constitute a valuable tool to constrain both the shock characteristics and pre-shock conditions. Methods. New CO observations from the shocked regions with the APEX and SOFIA telescopes are presented and combined. The integrated intensities are compared to the outputs of a grid of models, which were combined from an MHD shock code that calculates the dynamical and chemical structure of these regions, and a radiative transfer module based on the 'large velocity gradient' (LVG) approximation. Results. We base our modelling method on the higher J CO transitions, which unambiguously trace the passage of a shock wave. We provide fits for the blue- and red-lobe components of the observed shocks. We find that only stationary, C-type shock models can reproduce the observed levels of CO emission. Our best models are found for a pre-shock density of 104 cm-3, with the magnetic field strength varying between 45 and 100 {\mu}G, and a higher shock velocity for the so-called blue shock (\sim25 km s-1) than for the red one (\sim20 km s-1). Our models also satisfactorily account for the pure rotational H2 emission that is observed with Spitzer.

We present and analyze two spectrally resolved high-J CO lines towards the molecular outflow Cep E, driven by an intermediate-mass class 0 protostar. Using the GREAT receiver on board SOFIA, we observed the CO (12--11) and (13--12) transitions (E_u ~ 430 and 500 K, respectively) towards one position in the blue lobe of this outflow, that had been known to display high-velocity molecular emission. We detect the outflow emission in both transitions, up to extremely high velocities (~ 100 km/s with respect to the systemic velocity). We divide the line profiles into three velocity ranges that each have interesting spectral features: standard, intermediate, and extremely high-velocity. One distinct bullet is detected in each of the last two. A large velocity gradient analysis for these three velocity ranges provides constraints on the kinetic temperature and volume density of the emitting gas, >~ 100 K and > ~ 10^4 cm^-3, respectively. These results are in agreement with previous ISO observations and are comparable with results obtained by Herschel for similar objects. We conclude that high-J CO lines are a good tracer of molecular bullets in protostellar outflows. Our analysis suggests that different physical conditions are at work in the intermediate velocity range compared with the standard and extremely high-velocity gas at the observed position.

Carbon monoxide is an excellent tracer of the physical conditions of gas in molecular outflows from young stars. To understand the outflow mechanism we need to investigate the origin of the molecular emission and the structure and interaction of the outflowing molecular gas. Deriving the physical parameters of the gas will help us to trace and understand the various gas components in the flow. We observed CO(12-11) line emission at various positions along the L1157 bipolar outflow with GREAT aboard SOFIA. Comparing these new data with CO(2-1), we find basically constant line ratios along the outflow and even at the position of the source. These line ratios lead us to estimates of 10^5 to 10^6 cm^-3 for the gas density and 60 to 100 K for the gas temperature of the outflowing gas. The constrained density and temperature values indicate that we are mostly tracing a low-velocity gas component everywhere along the outflow, which is intermediate between the already known cold gas component, which gets entrained into the flow, and the hot gas, which gets shocked in the outflow.

During the formation of a star, material is ejected along powerful jets that impact the ambient material. This outflow phenomenon plays an important role in the regulation of star formation. Understanding the associated shocks and their energetic effects is therefore essential to the study of star formation. We present comparisons of shock models with observations of H$_2$ and SiO emission in the bipolar outflow BHR71, and predict water emission, under the basic assumption that the emission regions of the considered species coincide, at the resolution of currently available observations. New SiO observations from APEX are presented, and combined with \textit{Spitzer} and ground-based observations of H$_2$ to constrain shock models. The shock regions associated with targeted positions in the molecular outflow are studied by means of a 1D code that generates models of the propagation of stationary shock waves, and approximations to non-stationary ones. The SiO emission in the inner part of the outflow is concentrated near the apex of the corresponding bow-shock that is also seen in the pure rotational transitions of H$_2$. Simultaneous modelling is possible for H$_2$ and SiO and leads to constraints on the silicon pre-shock distribution on the grain mantles and/or cores. The best-fitting models are found to be of the non-stationary type, but the degeneracy of the solutions is still large. Pre-shock densities of 10$^4$ and 10$^5$ cm$^{-3}$ are investigated, and the associated best-model candidates have rather low velocity (respectively, 20-30 and 10-15 km s$^{-1}$) and are not older than 1000 years. We provide emission predictions for water, focusing on the brightest transitions, to be observed with the PACS and HIFI instruments of the \textit{Herschel} Telescope.

Spitzer-IRS maps of H_2 pure rotational lines from S(0) to S(7) in 3 outflows from Class 0 sources - L1448, BHR71, and NGC2071 - are presented. These lines are used, in conjunction with available rovibrational, near-IR H_2 lines, to probe the physical conditions of the warm gas between hundreds and thousands of Kelvin. We have constructed maps of the molecular hydrogen column density, ortho-to-para ratio and volume density, together with the index beta of the power law describing the distribution of gas temperature. In all three outflows, the present ortho-to-para ratio significantly deviates from the high temperature equilibrium of 3, being on average between 2.0 and 2.3. These low values, that reflect the young age of these flows, are found also in regions of relatively high temperature (~ 1000 K), likely indicating that shocks are occurring in a time shorter than that needed for a complete para to ortho conversion. Density maps indicate upper limits close to LTE conditions, i.e. between 10^6-10^7 cm^-3; moreover we demonstrate, based on the detections of HD emission spots (R(3) and R(4) lines), that a density stratification does exist, with the low density components (10^4-10^5 cm^-3) associated with the coldest gas. The beta index is found in all flows to be above 3.8: this value is consistent with predictions of multiple C-type bow shocks with a range of velocities, some of which insufficient to achieve the temperature at which H_2 partially dissociates. The contribution of H_2 to the total cooling is quantitatively similar to that of other abundant molecules emitting in the far-infrared, such as water and CO; moreover, the luminosity radiated away is comparable with the estimated kinetic energy of the swept-out outflow: this supports the scenario in which the shock from which H_2 is emitted is also capable of accelerating the molecular outflow driven at the shock working surface.

We study the production of SiO in the gas phase of molecular outflows, through the sputtering of Si--bearing material in refractory grain cores, which are taken to be olivine; we calculate also the rotational line spectrum of the SiO. The sputtering is driven by neutral particle impact on charged grains, in steady--state C-type shock waves, at the speed of ambipolar diffusion. The emission of the SiO molecule is calculated by means of an LVG code. A grid of models has been generated. We compare our results with those of an earlier study (Schilke et al. 1997). Improvements in the treatment of the coupling between the charged grains and the neutral fluid lead to narrower shock waves and lower fractions of Si being released into the gas phase. More realistic assumptions concerning the initial fractional abundance of O2 lead to SiO formation being delayed, so that it occurs in the cool, dense postshock flow. Good agreement is obtained with recent observations of SiO line intensities in the L1157 and L1448 molecular outflows. The inferred temperature, opacity, and SiO column density in the emission region differ significantly from those estimated by means of LVG `slab' models. The fractional abundance of SiO is deduced. Observed line profiles are wider than predicted and imply multiple, unresolved shock regions within the beam.

The bipolar HH 212 outflow has been mapped in SiO using the extended configuration of the Plateau de Bure Interferometer (PdBI), revealing a highly collimated SiO jet closely associated with the H2 jet component. We study at unprecedented resolution (0.34" across the jet axis) the properties of the innermost SiO ``microjet'' within 1000 AU of this young Class 0 source, to compare it with atomic microjets from more evolved sources and to constrain its origin. The SiO channel maps are used to investigate the microjet collimation and velocity structure. A large velocity gradient analysis is applied to SiO (2-1), (5-4) and (8-7) data from the PdBI and the Submillimeter Array to constrain the SiO opacity and abundance. The HH212 Class 0 microjet shows striking similarities in collimation and energetic budget with atomic microjets from T Tauri sources. Furthermore, the SiO lines appear optically thick, unlike what is generally assumed. We infer T(kin) ~ 50-500 K and an SiO/H2 abundance greater than 4 10(-8)-6 10(-5) for n(H2) = 10(7)-10(5) cm(-3), i.e. 0.05-90% of the elemental silicon. This similar jet width, regardless of the presence of a dense envelope, definitely rules out jet collimation by external pressure, and favors a common MHD self-collimation (and possibly acceleration) process at all stages of star formation. We propose that the more abundant SiO in Class 0 jets could mainly result from rapid (less than 25 yrs) molecular synthesis at high jet densities.