J-type Shocks Research Articles

Laboratory investigation of shock-induced dissociation of buckminsterfullerene and astrophysical insights

Fullerene C60 is one of the most iconic forms of carbon found in the interstellar medium (ISM). The interstellar chemistry of carbon-rich components, including fullerenes, is driven by a variety of energetic processes including UV and X-ray irradiation, cosmic-ray (CR) bombardment, electron impact, and shock waves. These violent events strongly alter the particle phase and lead to the release of new molecular species in the gas phase. Only a few experimental studies on the shock processing of cosmic analogs have been conducted so far. We explored in the laboratory the destruction of buckminsterfullerene C60 using a pressure-driven shock tube coupled with optical diagnostics. Our efforts were first devoted to probing in situ the shock-induced processing of C60 at high temperatures (≤ 4500 K) by optical emission spectroscopy. The analysis of the spectra points to the massive production of C2 units. A broad underlying continuum was observed as well and was attributed to the collective visible emission of carbon clusters, generated similarly in large amounts. This proposed assignment was performed with the help of calculated emission spectra of various carbon clusters. The competition between dissociation and radiative relaxation, determined by statistical analysis, alludes to a predominance of clusters with less than 40 carbon atoms. Our laboratory experiments, supported by molecular dynamics simulations performed in the canonical ensemble, suggest that C60 is very stable, and that high-energy input is required to process it under interstellar low-density conditions and to produce C2 units and an abundance of intermediate-sized carbon clusters. These results provide some insights into the life cycle of carbon in space. Our findings hint that only J-type shocks with velocities above ~100 km s−1 or C-type shocks with velocities above 9 km s−1 can lead to the destruction of fullerenes. Observational tracers of this process remain elusive, however. Our work confirms the potential of shock tubes for laboratory astrophysics.

Shock excitation of H2 in the James Webb Space Telescope era

Context. Molecular hydrogen, H2, is the most abundant molecule in the Universe. Thanks to its widely spaced energy levels, it predominantly lights up in warm gas, T ≳ 102 K, such as shocked regions externally irradiated or not by interstellar UV photons, and it is one of the prime targets of James Webb Space Telescope (JWST) observations. These may include shocks from protostellar outflows, supernova remnants impinging on molecular clouds, all the way up to starburst galaxies and active galactic nuclei. Aims. Sophisticated shock models are able to simulate H2 emission from such shocked regions. We aim to explore H2 excitation using shock models, and to test over which parameter space distinct signatures are produced in H2 emission. Methods. We here present simulated H2 emission using the Paris-Durham shock code over an extensive grid of ~14 000 plane-parallel stationary shock models, a large subset of which are exposed to a semi-isotropic external UV radiation field. The grid samples six input parameters: the preshock density, shock velocity, transverse magnetic field strength, UV radiation field strength, the cosmic-ray-ionization rate, and the abundance of polycyclic aromatic hydrocarbons, PAHs. Physical quantities resulting from our self-consistent calculations, such as temperature, density, and width, have been extracted along with H2 integrated line intensities. These simulations and results are publicly available on the Interstellar Medium Services platform. Results. The strength of the transverse magnetic field, as quantified by the magnetic scaling factor, b, plays a key role in the excitation of H2. At low values of b (≲0.3, J-type shocks), H2 excitation is dominated by vibrationally excited lines; whereas, at higher values (b ≳ 1, C-type shocks), rotational lines dominate the spectrum for shocks with an external radiation field comparable to (or lower than) the solar neighborhood. Shocks with b ≥ 1 can potentially be spatially resolved with JWST for nearby objects. H2 is typically the dominant coolant at lower densities (≲104 cm−3); at higher densities, other molecules such as CO, OH, and H2O take over at velocities ≲20 km s−1 and atoms, for example, H, O, and S, dominate at higher velocities. Together, the velocity and density set the input kinetic energy flux. When this increases, the excitation and integrated intensity of H2 increases similarly. An external UV field mainly serves to increase the excitation, particularly for shocks where the input radiation energy is comparable to the input kinetic energy flux. These results provide an overview of the energetic reprocessing of input kinetic energy flux and the resulting H2 line emission.

Modeling accretion shocks at the disk–envelope interface

Context.As material from an infalling protostellar envelope hits the forming disk, an accretion shock may develop which could (partially) alter the envelope material entering the disk. Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) indicate that emission originating from warm SO and SO2might be good tracers of such accretion shocks.Aims.The goal of this work is to test under what shock conditions the abundances of gas-phase SO and SO2increase in an accretion shock at the disk–envelope interface.Methods.Detailed shock models including gas dynamics were computed using the Paris-Durham shock code for nonmagnetizedJ-type accretion shocks in typical inner envelope conditions. The effect of the preshock density, shock velocity, and strength of the ultraviolet (UV) radiation field on the abundance of warm SO and SO2is explored. Compared with outflows, these shocks involve higher densities (~107cm−3), lower shock velocities (~few km s−1), and large dust grains (~0.2 μm) and thus probe a different parameter space.Results.Warm gas-phase chemistry is efficient in forming SO under mostJ-type shock conditions considered. In lower-velocity (~3 km s−1) shocks, the abundance of SO is increased through subsequent reactions starting from thermally desorbed CH4toward H2CO and finally SO. In higher velocity (≳4 km s−1) shocks, both SO and SO2are formed through reactions of OH and atomic S. The strength of the UV radiation field is crucial for SO and in particular SO2formation through the photodissociation of H2O. Thermal desorption of SO and SO2ice is only relevant in high-velocity (≳5 km s−1) shocks at high densities (≳107cm−3). Both the composition in the gas phase, in particular the abundances of atomic S and O, and in ices such as H2S, CH4, SO, and SO2play a key role in the abundances of SO and SO2that are reached in the shock.Conclusions.Warm emission from SO and SO2is a possible tracer of accretion shocks at the disk–envelope interface as long as a local UV field is present. Observations with ALMA at high-angular resolution could provide further constraints given that other key species for the gas-phase formation of SO and SO2, such as H2S and H2CO, are also covered. Moreover, theJames WebbSpace Telescope will give access to other possible slow, dense shock tracers such as H2, H2O, and [S I] 25 μm.

Probing the hidden atomic gas in Class I jets with SOFIA

Context. We present SOFIA/FIFI-LS observations of five prototypical, low-mass Class I outflows (HH111, SVS13, HH26, HH34, HH30) in the far-infrared [O I]63μm and [O I]145μm transitions. Aims. Spectroscopic [O I]63μm,145μm maps enable us to study the spatial extent of warm, low-excitation atomic gas within outflows driven by Class I protostars. These [O I] maps may potentially allow us to measure the mass-loss rates (Ṁjet) of this warm component of the atomic jet. Methods. A fundamental tracer of warm (i.e. T ~ 500–1500 K), low-excitation atomic gas is the [O I]63μm emission line, which is predicted to be the main coolant of dense dissociative J-type shocks caused by decelerated wind or jet shocks associated with protostellar outflows. Under these conditions, the [O I]63μm line can be directly connected to the instantaneous mass ejection rate. Thus, by utilising spectroscopic [O I]63μm maps, we wish to determine the atomic mass flux rate Ṁjet ejected from our target outflows. Results. Strong [O I]63μm emission is detected at the driving sources HH111IRS, HH34IRS, SVS13, as well as at the bow shock region, HH7. The detection of the [O I]63μm line at HH26A and HH8/HH10 can be attributed to jet deflection regions. The far-infrared counterpart of the optical jet is detected in [O I]63μm only for HH111, but not for HH34. We interpret the [O I]63μm emission at HH111IRS, HH34IRS, and SVS13 to be coming primarily from a decelerated wind shock, whereas multiple internal shocks within the HH111 jet may cause most of the [O I]63μm emission seen there. At HH30, no [O I]63μm,145μm was detected. The [O I]145μm line detection is at noise level almost everywhere in our obtained maps. The observed outflow rates of our Class I sample are to the order of Ṁjet ~ 10−6M⊙ yr−1, if proper shock conditions prevail. Independent calculations connecting the [O I]63μm line luminosity and observable jet parameters with the mass -loss rate are consistent with the applied shock model and lead to similar mass-loss rates. We discuss applicability and caveats of both methods. Conclusions. High-quality spectroscopic [O I]63μm maps of protostellar outflows at the jet driving source potentially allow a clear determination of the mass ejection rate.

Tracing shock type with chemical diagnostics

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

Supernova Shocks in Molecular Clouds: Velocity Distribution of Molecular Hydrogen

Abstract 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 at high spectral resolution, using the Stratospheric Observatory for Infrared Astronomy, of the H2 S(5) pure rotational line from molecular shocks in the supernova remnant IC 443. We generate MHD models for nonsteady-state shocks driven by the pressure of the IC 443 blast wave into gas of densities 103–105 cm−3. We present the first detailed derivation of the shape of the velocity profile for emission from H2 lines behind such shocks, taking into account the shock age, preshock density, and magnetic field. For preshock densities 103–105 cm−3, the H2 emission arises from layers that extend 0.01–0.0003 pc behind the shock, respectively. The predicted shifts of line centers, and the line widths, of the H2 lines range from 20–2 and 30–4 km s−1, 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 103 to 2 × 104 cm−3 and strong magnetic fields. Two positions in clump B were observed. For clump B2 (a fainter region near clump B), the H2 spectrum requires a J-type shock into moderate-density gas (∼102 cm−3) with the gas accelerated to 100 km s−1 from its preshock 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 B2) 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.

Models of irradiated molecular shocks

Context. The recent discovery of excited molecules in starburst galaxies observed with ALMA and the Herschel space telescope has highlighted the necessity to understand the relative contributions of radiative and mechanical energies in the formation of molecular lines and explore the conundrum of turbulent gas bred in the wake of galactic outflows. Aims. The goal of the paper is to present a detailed study of the propagation of low velocity (5–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. 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. Conclusions. Ultraviolet photons are found to strongly modify the way the mechanical energy of interstellar shocks is processed and radiated away. In spite of what intuition dictates, a strong external UV radiation field boosts the efficiency of low velocity interstellar shocks in the production of several molecular lines which become evident tracers of turbulent dissipation.

The VLA Nascent Disk and Multiplicity Survey of Perseus Protostars (VANDAM). IV. Free–Free Emission from Protostars: Links to Infrared Properties, Outflow Tracers, and Protostellar Disk Masses

Abstract Emission from protostars at centimeter radio wavelengths has been shown to trace the free–free emission arising from ionizing shocks as a result of jets and outflows driven by protostars. Therefore, measuring properties of protostars at radio frequencies can provide valuable insights into the nature of their outflows and jets. We present a C-band (4.1 and 6.4 cm) survey of all known protostars (Class 0 and Class I) in Perseus as part of the VLA Nascent Disk and Multiplicity (VANDAM) Survey. We examine the known correlations between radio flux density and protostellar parameters, such as bolometric luminosity and outflow force, for our sample. We also investigate the relationship between radio flux density and far-infrared line luminosities from Herschel. We show that free–free emission most likely originates from J-type shocks; however, the large scatter indicates that those two types of emission probe different time and spatial scales. Using C-band fluxes, we removed an estimation of free–free contamination from the corresponding Ka-band (9 mm) flux densities that primarily probe dust emission from embedded disks. We find that the compact (<1″) dust emission is lower for Class I sources (median dust mass 96 M ⊕) relative to Class 0 (248 M ⊕), but several times higher than in Class II (5–15 M ⊕). If this compact dust emission is tracing primarily the embedded disk, as is likely for many sources, this result provides evidence of decreasing disk masses with protostellar evolution, with sufficient mass for forming giant planet cores primarily at early times.

MHOs toward HMOs: A Search for Molecular Hydrogen Emission-Line Objects toward High-mass Outflows

Abstract We present the results of a narrow-band near-infrared imaging survey for Molecular Hydrogen emission-line Objects (MHOs) toward 26 regions containing high-mass protostellar candidates and massive molecular outflows. We have detected a total of 236 MHOs, 156 of which are new detections, in 22 out of the 26 regions. We use H2 2.12 μm/H2 2.25 μm flux ratios, together with morphology, to separate the signatures of fluorescence associated with photo-dissociation regions (PDRs) from shocks associated with outflows in order to identify the MHOs. PDRs have typical low flux ratios of ∼1.5–3, while the vast majority of MHOs display flux ratios typical of C-type shocks (∼6–20). A few MHOs exhibit flux ratios consistent with expected values for J-type shocks (∼3–4), but these are located in regions that may be contaminated with fluorescent emission. Some previously reported MHOs have low flux ratios, and are likely parts of PDRs rather than shocks indicative of outflows. We identify a total of 36 outflows across the 22 target regions where MHOs were detected. In over half these regions, MHO arrangements and fluorescent structures trace features present in CO outflow maps, suggesting that the CO emission traces a combination of dynamical effects, which may include gas entrained in expanding PDRs as well as bipolar outflows. Where possible, we link MHO complexes to distinct outflows and identify candidate driving sources.

The shocked gas of the BHR71 outflow observed byHerschel: indirect evidence for an atomic jet

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 ORION FINGERS: NEAR-IR SPECTRAL IMAGING OF AN EXPLOSIVE OUTFLOW

ABSTRACT We present near-IR (1.1–2.4 μm) position–position–velocity cubes of the 500 year old Orion BN/KL explosive outflow with spatial resolution 1″ and spectral resolution 86 km s−1. We construct integrated intensity maps free of continuum sources of 15 H2 and [Fe ii] lines while preserving kinematic information of individual outflow features. Included in the detected H2 lines are the 1-0 S(1) and 1-0 Q(3) transitions, allowing extinction measurements across the outflow. Additionally, we present dereddened flux ratios for over two dozen outflow features to allow for the characterization of the true excitation conditions of the BN/KL outflow. All of the ratios show the dominance of the shock excitation of the H2 emission, although some features exhibit signs of fluorescent excitation from stellar radiation or J-type shocks. We also detect tracers of the PDR/ionization front north of the Trapezium stars in [O i] and [Fe ii] and analyze other observed outflows not associated with the BN/KL outflow.

The density structure of the L1157 molecular outflow★†

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

DISCOVERY OF A PRE-EXISTING MOLECULAR FILAMENT ASSOCIATED WITH SUPERNOVA REMNANT G127.1+0.5

We performed millimeter observations in CO lines toward the supernova remnant (SNR) G127.1+0.5. We found a molecular filament at 4--13 km s$^{-1}$ consisting of two distinct parts: a straight part coming out of the remnant region and a curved part in the remnant region. The curved part is coincides well with the bright SNR shell detected in 1420 MHz radio continuum and mid-infrared observations in the northeastern region. In addition, redshifted line wing broadening is found only in the curved part of the molecular filament, which indicates a physical interaction. These provide strong evidences, for the first time, to confirm the association between an SNR and a pre-existing long molecular filament. Multi-band observations in the northeastern remnant shell could be explained by the interaction between the remnant shock and the dense molecular filament. RADEX radiative transfer modeling of the quiet and shocked components yield physical conditions consistent with the passage of a non-dissociative J-type shock. We argue that the curved part of the filament is fully engulfed by the remnant's forward shock. A spatial correlation between aggregated young stellar objects (YSOs) and the adjacent molecular filament close to the SNR is also found, which could be related to the progenitor's activity.

Dust, ice and gas in time (DIGIT):HerschelandSpitzerspectro-imaging of SMM3 and SMM4 in Serpens

We report on spectro-imaging observations employing Spitzer IRS and Herschel PACS, aiming to constrain the physical conditions around SMM3 and SMM4 in Serpens. The combined power of both instruments provides an almost complete wavelength coverage between 5 and 200 micron at an angular resolution of 10". We detect line emission from all major molecular (H2, CO, H2O and OH) and many atomic ([OI], [CII], [FeII], [SiII] and [SI]) coolants. Line emission tends to peak at distances of 10 - 20" from the protostellar sources, at positions of known outflow shocks. The only exception is [CII] which likely traces a PDR excited from the neighboring source SMM6. Excitation analysis indicates that H2 and CO originate from gas at two distinct rotational temperatures of 300 K and 1000 K, while H2O and OH emission corresponds to rotational temperatures of 100 - 200 K. The morphological and physical association between CO and H2 suggests a common excitation mechanism which allows direct comparisons between the two molecules. The CO/H2 abundance ratio varies from 10^-5 in the warm gas up to 10^-4 in the hotter regions. The occurrence of J-shocks is suggested by the strong atomic/ionic (except for [CII]) emission as well as a number of line ratio diagnostics. Both C- and J-shocks can account for the observed molecular emission, however J-shocks are strongly advocated by the atomic emission and provide simpler and more homogeneous solutions for CO and H2. C-shocks describe better the emission from H2O and OH. The variations in the CO/H2 abundance ratio for gas at different temperatures can be interpreted by their reformation rates in dissociative J-type shocks, or the simultaneous influence of both C and J shocks.

Sputtering in oblique C-type shocks

We present the first results for the sputtering of grain mantles and cores obtained with self-consistent multifluid hydromagnetic models of C-type shocks propagating through dusty media. The threshold shock speed for mantle sputtering is about 10 km/s and is independent of density. The mantles are completely vapourised in shocks with speeds of 20-25 km/s. At such shock speeds core sputtering commences and gas-phase SiO forms. Core destruction is not total in any C-type shock because grains are not completely destroyed in shocks with speeds near the minimum speeds at which J-type shocks appear. Due to the density-dependence of the critical shock speed for this transition, higher gas-phase SiO fractional abundances are produced behind shocks propagating in lower density gas. For shock speeds near the threshold speeds for both core and mantle sputtering, sputtering is much greater for shock velocities at smaller angles relative to the upstream magnetic field. At higher shock speeds, the angular variation is still present but less pronounced.

The CHESS survey of the L1157-B1 shock: the dissociative jet shock as revealed byHerschel–PACS

Outflows generated by protostars heavily affect the kinematics and chemistry of the hosting molecular cloud through strong shocks that enhance the abundance of some molecules. L1157 is the prototype of chemically active outflows, and a strong shock, called B1, is taking place in its blue lobe between the precessing jet and the hosting cloud. We present the Herschel-PACS 55--210 micron spectra of the L1157-B1 shock, showing emission lines from CO, H2O, OH, and [OI]. The spatial resolution of the PACS spectrometer allows us to map the warm gas traced by far-infrared (FIR) lines with unprecedented detail. The rotational diagram of the high-Jup CO lines indicates high-excitation conditions (Tex ~ 210 +/- 10 K). We used a radiative transfer code to model the hot CO gas emission observed with PACS and in the CO (13-12) and (10-9) lines measured by Herschel-HIFI. We derive 200<Tkin<800 K and n>10^5 cm-3. The CO emission comes from a region of about 7 arcsec located at the rear of the bow shock where the [OI] and OH emission also originate. Comparison with shock models shows that the bright [OI] and OH emissions trace a dissociative J-type shock, which is also supported by a previous detection of [FeII] at the same position. The inferred mass-flux is consistent with the "reverse" shock where the jet is impacting on the L1157-B1 bow shock. The same shock may contribute significantly to the high-Jup CO emission.

TheHerschelHIFI water line survey in the low-mass proto-stellar outflow L1448

As part of the WISH (Water In Star-forming regions with Herschel) key project, we report on the observations of several ortho- and para-H2O lines performed with the HIFI instrument towards two bright shock spots (R4 and B2) along the outflow driven by the L1448 low-mass proto-stellar system, located in the Perseus cloud. These data are used to identify the physical conditions giving rise to the H2O emission and infer any dependence with velocity. These observations provide evidence that the observed water lines probe a warm (T_kin~400-600 K) and very dense (n 10^6 - 10^7 cm^-3) gas, not traced by other molecules, such as low-J CO and SiO, but rather traced by mid-IR H2 emission. In particular, H2O shows strong differences with SiO in the excitation conditions and in the line profiles in the two observed shocked positions, pointing to chemical variations across the various velocity regimes and chemical evolution in the different shock spots. Physical and kinematical differences can be seen at the two shocked positions. At the R4 position, two velocity components with different excitation can be distinguished, with the component at higher velocity (R4-HV) being less extended and less dense than the low velocity component (R4-LV). H2O column densities of about 2 10^13 and 4 10^14 cm^-2 have been derived for the R4-LV and the R4-HV components, respectively. The conditions inferred for the B2 position are similar to those of the R4-HV component, with H2O column density in the range 10^14 - 5 10^14 cm^-2, corresponding to H2O/H2 abundances in the range 0.5 - 1 10^-5. The observed line ratios and the derived physical conditions seem to be more consistent with excitation in a low velocity J-type shock with large compression rather than in a stationary C-shock, although none of these stationary models seems able to reproduce all the characteristics of the observed emission.

High-resolution radio study of SNR IC 443 at low radio frequencies

We investigated in detail the morphology at low radio frequencies of the supernova remnant IC443 and accurately established its radio continuum spectral properties. We used the VLA in multiple configurations to produce high resolution radio images of IC443 at 74 and 330 MHz. The changes with position in the radio spectral index were correlated with data in near infrared from 2MASS, in gamma-rays from VERITAS, and with the molecular 12^CO line emission. The new image at 74 MHz has HPBW=35", rms=30 mJy/beam and at 330 MHz HPBW= 17" and rms=1.7 mJy/beam. The integrated flux densities for the whole SNR are S_74MHz=470+/-51 Jy and S_330MHz=248+/-15 Jy. For the pulsar wind nebula associated with the compact source CXOUJ061705.3+222127, we calculated S_330MHz=0.23+/-0.05 Jy, S_1420MHz=0.20+/-0.04 Jy, and alpha~0.0. Substantial variations are observed in spectral index between 74 and 330 MHz across IC443. The flattest spectral components (-0.25< alpha<-0.05) coincide with the brightest parts of the SNR along the eastern border, with an impressive agreement with ionic lines as observed in the 2MASS J and H bands. This result strongly suggests that the passage of a fast dissociating J-type shock across the interacting molecular cloud dissociated the molecules and ionized the gas. We therefore conclude that thermal absorption at 74 MHz is responsible for the localized spectral index flattening observed along the eastern border of IC443. The diffuse interior of IC443 has a spectrum steeper than found anywhere in the SNR (-0.85<alpha<-0.6), while the southern ridge again has a flatter spectrum (alpha -0.4), consequence in this case of the strong shock/molecular cloud interaction. At the available statistics the VERITAS gamma-ray emission strikingly matches the CO distribution, but no clear evidence is found for a morphological correlation between the TeV distribution and radio emission.

Herschelobservations of the Herbig-Haro objects HH 52-54

We are aiming at the observational estimation of the relative contribution to the cooling by CO and H2O, as this provides decisive information for the understanding of the oxygen chemistry behind interstellar shock waves. Methods. The high sensitivity of HIFI, in combination with its high spectral resolution capability, allows us to trace the H2O outflow wings at unprecedented signal-to-noise. From the observation of spectrally resolved H2O and CO lines in the HH52-54 system, both from space and from ground, we arrive at the spatial and velocity distribution of the molecular outflow gas. Solving the statistical equilibrium and non-LTE radiative transfer equations provides us with estimates of the physical parameters of this gas, including the cooling rate ratios of the species. The radiative transfer is based on an ALI code, where we use the fact that variable shock strengths, distributed along the front, are naturally implied by a curved surface. Based on observations of CO and H2O spectral lines, we conclude that the emission is confined to the HH54 region. The quantitative analysis of our observations favours a ratio of the CO-to-H2O-cooling-rate >> 1. From the best-fit model to the CO emission, we arrive at an H2O abundance close to 1e-5. The line profiles exhibit two components, one of which is triangular and another, which is a superposed, additional feature. This additional feature likely originates from a region smaller than the beam where the ortho-water abundance is smaller than in the quiescent gas. Comparison with recent shock models indicate that a planar shock can not easily explain the observed line strengths and triangular line profiles.We conclude that the geometry can play an important role. Although abundances support a scenario where J-type shocks are present, higher cooling rate ratios than predicted by these type of shocks are derived.

First hyperfine resolved far-infrared OH spectrum from a star-forming region

OH is an important molecule in the H2O chemistry and the cooling budget of star-forming regions. The goal of the Herschel key program `Water in Star-forming regions with Herschel' (WISH) is to study H2O and related species during protostellar evolution. Our aim in this letter is to assess the origin of the OH emission from star-forming regions and constrain the properties of the emitting gas. High-resolution observations of the OH 2Pi1/2 J = 3/2-1/2 triplet at 1837.8 GHz (163.1 micron) towards the high-mass star-forming region W3 IRS 5 with the Heterodyne Instrument for the Far-Infrared (HIFI) on Herschel reveal the first hyperfine velocity-resolved OH far-infrared spectrum of a star-forming region. The line profile of the OH emission shows two components: a narrow component (FWHM approx. 4-5 km/s) with partially resolved hyperfine structure resides on top of a broad (FWHM approx. 30 km/s) component. The narrow emission agrees well with results from radiative transfer calculations of a spherical envelope model for W3 IRS 5 with a constant OH abundance of approx. 8e-9. Comparison with H2O yields OH/H2O abundance ratios of around 1e-3 for T > 100 K and around unity for T < 100K, consistent with the current picture of the dense cloud chemistry with freeze-out and photodesorption. The broad component is attributed to outflow emission. An abundance ratio of OH/H2O > 0.028 in the outflow is derived from comparison with results of water line modeling. This ratio can be explained by a fast J-type shock or a slower UV-irradiated C-type shock.