Free Access
Issue
A&A
Volume 602, June 2017
Article Number A8
Number of page(s) 10
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201730454
Published online 24 May 2017

© ESO, 2017

1. Introduction

Protostellar outflows constitute the most prominent signposts of star formation. In the earliest phases, when the young protostar builds up most of its final mass by accreting material from its parental cloud, the ejection of magneto-centrifugal jets is believed to help the protostellar system to release angular momentum from its inner disk system. Slower wide-angle winds surrounding the fast jets stir up the ambient cloud environment and create molecular outflow cavities (e.g. Arce et al. 2007; Frank et al. 2014). These protostellar feedback processes play a major role in protostellar mass accretion, and in the interaction of the protostar with its parental environment. Their understanding is still subject to many open questions, for instance about the nature and universality of jets, their temporal activity, and their properties in terms of energy and momentum transfer. It is therefore of paramount importance to study the properties of protostellar jets in various protostellar sources, as has been done recently towards L1157 (e.g. Podio et al. 2016; Tafalla et al. 2015) and BHR71 (e.g. Benedettini et al. 2017), for example.

Jets have been indirectly traced by bow-shock emission (or working surfaces, e.g. Gueth et al. 1998), outflow cavities (or swept-up material, Gueth et al. 1996; Bachiller et al. 2001), and compact dissociative (e.g. Kristensen et al. 2013; Nisini et al. 2015) or high-velocity molecular (e.g. Tafalla et al. 2015) shocks. However, a direct probe of jets in outflows requires observations of typical tracers that are spatially or spectroscopically resolved. One such typical jet tracer is the fine-structure line of atomic oxygen at 63 μm, which is one of the main coolants in jets. The German Receiver for Astronomy at TeraHertz frequencies (GREAT) in the Stratospheric Observatory for Infrared Astronomy (SOFIA) now allows spectroscopically resolved observations of [OI]63μm and thereby offers this transition as a valuable tool for studying protostellar jet systems (for a detailed study on the potential of this tracer in the context of star formation we refer to Leurini et al. 2015).

thumbnail Fig. 1

Left: a finder’s chart for the positions of interest in the Cep E protostellar outflow: Spitzer-IRAC band-two (4.5 μm) image, retrieved from the Spitzer archive. Grey circles mark the diffraction-limited SOFIA beams at the frequencies of [OI]63μm (smaller circles) and CO (16–15) or OH at 1838 GHz (larger circles) lines. The three observed positions are indicated (Cep E-mm, red star, BI and BII, white rectangles), as well as others in the red lobe (RI and RII). The knots denomination is from Gómez-Ruiz et al. (2012). Right: zoom on the southern lobe of the outflow showing the spatial components identified by L15 in CO (2–1) (see Sect. 3): the jet (cyan contours), the terminal bow-shock (yellow contours distant from the protostar) and the cavity walls (coloured background). The three observed positions are indicated with a white star, and the SOFIA beams at the frequencies of [OI]63μm (smaller circles) and CO (16–15) or OH at 1838 GHz (larger circles) lines are shown in BII.

In order to fully exploit the potential of [OI]63μm as a diagnostic tool, its interpretation relies on the use of detailed shock modelling. The Paris-Durham shock model (Flower & Pineau des Forêts 2015) describes slow (νs ≲ 50 km s-1), C-, J-, and time-dependent magnetohydrodynamic (MHD) shocks, propagating in dark molecular environments. It internally computes the [OI]63μm emissivities in optically thin local thermodynamical equilibrium (LTE) conditions. In faster shocks, the high temperature reached in the shock front triggers UV emission from the gas. This UV emission can propagate to the pre-shock medium and cause a radiative precursor that leads to photodissociation and ionization in that gas (Hollenbach & McKee 1989). For such faster shocks, with velocities up to 150 km s-1, these authors provide [OI]63μm emissivities for a range of pre-shock densities. A combination of these different models can be used to constrain the physical and dynamical characteristics of jet-driven shocks, as well as the conditions in the pre-shock environment by comparing model predictions with the observed intensities of [OI]63μm and other available tracers.

In the present study, we applied this strategy to understand the nature of shocks in the Cep E protostellar outflow. This source is an intermediate-mass Class 0 protostar (L = 100 L, Lefloch et al. 1996; Chini et al. 2001) in the Cepheus OB3 association at a distance of 730 pc (Sargent 1977). Its prominent and well-studied outflow is driven by a jet of a size of 1.7′′ × 21′′, as identified in several CO transitions (Gómez-Ruiz et al. 2012) up to CO (16–15) and interpreted by means of time-dependent MHD shock models by (Lefloch et al. 2015; hereafter L15). In particular, the CO emission of shocked gas along the jet was found to be consistent with very young (220–740 yr) shocks with velocities of 20–30 km s-1, propagating into a pre-shock medium of density nH = (0.51.0) × 105 cm-3. However, this study already demonstrated the need for additional observations in order to achieve a more precise determination of the jet shock parameters. In this work, we aim at further constraining the nature of this protostellar jet using [OI]63μm observations, together with an extended grid of MHD Paris-Durham shock models. This article is structured as follows: observations are introduced in Sect. 2 and direct results are presented in Sect. 3. In Sect. 4, we infer abundances from the kinematical components identified in our observations. We compare our observations to sophisticated models of interstellar shocks in Sect. 5 to understand the nature of such shocks, while Sect. 6 contains our concluding remarks.

2. Observations

The observations of the Cep E protostellar outflow were conducted with the GREAT1 receiver (Heyminck et al. 2012) during three SOFIA flights on 15, 16, and 18 December 2015 (legs 12, 8, and 8, respectively), as part of the Cycle 3 Community science program. Three positions were observed: Cep E mm, BI, and BII (see Fig. 1 and Table 1).

The observation of the [OI]63μm line is possible in the H channel of GREAT using a hot electron bolometer heterodyne mixer (Büchel et al. 2015), with a 2.5 GHz bandwidth receiver and a spectral resolution of 0.99 km s-1. For the tuning to the [OI]63μm line, a quantum-cascade laser was used as local oscillator (Richter et al. 2015). At the Cep E-mm position, the zero-level width observed for example by Lefloch et al. (2011) in SiO or H2O lines is wider than 200 km s-1. This is more than the bandwidth available at the 4744.778 GHz frequency of the [OI]63μm line (~160 km s-1). For this reason, we covered the emission of this line in two spectral windows. For the BI and BII positions, we prepared our tuning setup based on the CO line profiles showed by L15 at the BII position, which were about 150 km s-1 wide (at baseline level). A single USB tuning was sufficient to cover the [−150; 10] km s-1 interval at the frequency of the [OI]63μm line.

For OH and CO observations of all three positions, the lower frequency, L2 channel was connected to 4 GHz wide digital back-ends described in Risacher et al. (2016), providing respective spectral resolutions of 1.00 and 0.99 km s-1 (given in Table 2). This channel was tuned in LSB to the frequency of the OH triplet between Π1/22\hbox{$^2\Pi_{1/2}$}J = 3/2 and J = 1/2 at 1837.817 GHz. This setup allows us to pick up the CO (16–15) in the upper sideband at 1841.346 GHz (see Table 2 for spectroscopic parameters).

The resulting rms uncertainties are given in Table 1. The observations were performed in double beam-switching mode with an amplitude of 40′′ (80′′ throw) at the position angle of 45° and a phase time of 0.5 s. The nominal focus position was updated regularly against temperature drifts of the telescope structure. The pointing was established with the optical guide cameras to an accuracy of better than ~1′′. The beam widths and efficiencies are indicated in Table 2. The data were calibrated with the KOSMA/GREAT calibrator (Guan et al. 2012), which removes residual telluric lines, and were further processed (mostly by removing linear baseline) with the CLASS software2.

3. Results

3.1. Line profiles

Table 1

Observed positions and sensitivity (rms of baseline noise) at 1 km s-1 spectral resolution for each observed line.

The spectra of the CO (16–15), [OI]63μm, and OH 1838 GHz lines obtained at the three positions are shown in Fig. 2. L15 mapped the Cep E protostellar outflow in CO (2–1) at high angular resolution with the Plateau de Bure Interferometer (PdBI) and performed pointed CO observations with numerous single-dish telescopes at the Cep E-BII position. This allowed them to spectrally and spatially identify three components in the outflow, which we refer to throughout this article:

  • a narrow jet emitting between 140 and 110, centred on 125 km s-1, with a typical line width of 15 km s-1, well separated from the bipolar outflow wing emission;

  • the outflow cavity, composed of the gas contained in the cavity walls, that is, of the cavity walls plus the entrained gas that is not part of the jet that is emitted from the ambient cloud velocity (10.9 km s-1, Lefloch et al. 1996) up to 100 km s-1;

  • the terminal bow-shock labelled ‘HH377’, which is emitted between 90 and 50 km s-1 and appears as a bump between the outflow cavity and the jet in the line profile, and whose CO (16–15) intensity is similar to the entrained gas at 20 km s-1.

At the protostar position, all lines are undetected except for a faint 3σ emission of the CO (16–15) line around the systemic velocity of Cep E (at 10.9 km s-1). Hints of emission can be seen in the blue- and red-shifted wing, between 90 and 60, and 40 and 60 km s-1, respectively.

At the Cep E-BI and BII positions, only OH is undetected. At the Cep E-BI position, only the spectral component that was found by L15 to be associated with the walls of the outflow cavity is detected in CO and O i. In CO, this component exhibits a wing that extends up to 75 km s-1, whereas the corresponding O i emission is found to be much narrower.

Table 2

Characteristics and observational parameters of the OH transitions between the Π1/22\hbox{$^2\Pi_{1/2}$}J = 3/2 and J = 1/2 states, and of the [OI]63μm and CO (16–15) lines.

thumbnail Fig. 2

CO (16–15) (grey histograms), [OI]63μm (red line), and OH at 1838 GHz (blue line) observations obtained in Cep E-mm (top panel), Cep E-BI (middle panel), and Cep E-BII (lower panel). The vertical dashed lines represent the central value for the velocity components identified by L15 in CO transitions: 125 (cyan), 68 (pink), and 10.9 km s-1 (black; systemic velocity of Cep E). The spectral resolution is 1 km s-1 for the Cep E-BI and BII spectra, and 4 km s-1 for Cep E-mm.

Finally, at the Cep E-BII position, the three components identified by L15, that is, the jet, the terminal bow-shock, and the outflow cavity, appear in the CO (16–15) and in [OI]63μm line profiles. At the BI and BII positions, a slight shift in velocity of the order of 5 km s-1 might exist between the peak of the CO (16–15) emission and that of [OI]63μm. This shift could be a consequence of the two lines tracing a different material, or it could be an effect of self-absorption of [OI]63μm near the systemic velocity. In particular, the bright bump in the entrained gas seen at 20 km s-1 in CO (16–15) appears between 30 and 20 km s-1 in [OI]63μm. Additionally, this bright bump is weaker than the HH377 component in [OI]63μm, in contrast to the CO (16–15) case. This could be an effect of different filling factors resulting from the different telescope beams at 4744.778 GHz and 1841.3455 GHz (se Fig. 1). As there is no [OI]63μm detection at the protostar position, and as the BI and BII observations of CO (16–15) are not independent (see Fig. 1), we focused on the spectra obtained at the BII position.

3.2. Line ratios

We extracted integrated intensities from the [OI]63μm and CO (16–15) spectra at the BII position. We based our analysis on the three kinematical components identified by L15 (see previous section): the jet at high velocities, the terminal bow-shock at intermediate velocities, and the outflow cavity at lower velocities. Our method was the same for the [OI]63μm and CO (16–15) lines, given the similarity of their profiles. The method that we used to extract the integrated intensity from each kinematical component is fully consistent with that of L15:

  • for the jet component, we integrated the emission between 140 and 100 km s-1;

  • for the HH377 component, we considered the velocity interval between 100 and 25 km s-1, in which we removed the outflow cavity component by fitting a first-order baseline, and then integrated the emission of the residual between 90 and 40 km s-1;

  • for the outflow cavity, we integrated the emission between 100 and 14 km s-1, from which we subtracted the integrated intensity associated with the HH377 terminal shock.

The final values that we extracted from this analysis can be found in Table 3. We used conservative error bars of ±10%, as the 3σ value probably does not reflect the uncertainty intrinsic to the extraction we performed.

The table indicates similar integrated intensities for the jet and the terminal bow-shock components in [OI]63μm and CO (16–15). These integrated intensities are significantly lower than the contribution from the outflow cavity. The line ratios found for the three components are similar.

4. Abundances

4.1. Filling factors

We then estimated filling factor corrections to the integrated intensities. Based on their CO (2–1) interferometric observations, L15 inferred the size of the emission region for the jet, bow-shock, and outflow cavity components: 1.7′′ × 21′′, 4.5′′, and 22′′ × 10′′. The authors were then able to infer filling factors for the jet, shock, and cavity wall components of their SOFIA/GREAT observations of CO (16–15) based on the assumption that the emission region was the same. As our observations were pointed in exactly the same BII position, we used the filling factor values derived from their observations for CO (16–15): 0.09 for the jet, 0.06 for the shock, and 0.15 for the cavity walls (see Table 4).

For the [OI]63μm line emission, our method relies on the combined use of the Herschel-PACS footprint (see Appendix A and Fig. A.1) obtained at this wavelength, of the Spitzer-IRAC band-two (4.5 μm) image (see Fig. 1), and of the CO (2–1) interferometric map of L15 (also see Fig. 1):

  • the outflow cavity most likely dominates the emission in thePACS footprint. In this map, the emission region has a typicaldeconvolved size of 10′′. Given that the map is not Nyquist-sampled, and given the significantly different beam sizes for CO (16–15) and [OI]63μm lines, we adopted a simple filling factor of 1 for this component in [OI]63μm;

  • the size of HH377 is measured in CO (2–1) (~4.5′′, L15) and in Spitzer-IRAC band-two emission (deconvolved size of about 10′′). As HH377 emission is visible in the low-excitation CO (2–1) and the high-excitation CO (16–15) line, we consider that its typical size in [OI]63μm is between 4.5 and 10′′, yielding a filling factor of 0.55 ± 0.15;

  • finally, we inferred the filling factor for the jet emission from the interferometric CO (2–1) map because this component is similar in CO (2–1) and [OI]63μm: 0.25 ± 0.05.

The filling factors are indicated in Table 4 together with the corrected integrated intensities and fluxes (in erg cm-2 s-1 sr-1) obtained using the [c3/(2 kν3)] conversion factor. The integrated intensity ratios strongly depend on the filling factor correction, as very different corrections affect both lines. This new set of results indicates that the [OI]63μm/CO (16–15) integrated intensity ratio is significantly lower in the HH377 bow-shock than in the jet and arguably in the outflow cavity, as well. We note that valuable information would be gained from a fully sampled map of this region with SOFIA-GREAT. This would allow us to study the extension of the emission spectral component per spectral component, although the spatial resolution of SOFIA might not suffice to reach constraining conclusions.

Table 3

Raw integrated intensities extracted from the [OI]63μm and CO (16–15) spectra, in units of K km s-1, and corresponding line ratio, kinematic component by kinematic component (jet, bow-shock, outflow cavity).

Table 4

Filling-factor-corrected, integrated intensities extracted from the [OI]63μm and CO (16–15) spectra, in units of K km s-1 and erg cm-2 s-1 sr-1, and corresponding line ratio, kinematic component by kinematic component (jet, bow-shock, outflow cavity).

4.2. Abundances

We finally computed the O i column density. We assumed optically thin emission and LTE conditions. Both assumptions are discussed at the end of this section. Respectively denoting with h, k, c, ν, Eu, Aul, and gu Planck’s constant, Boltzmann’s constant, the speed of light, the line frequency, upper level energy, Einstein’s coefficient, and the upper level statistical weight of the observed transition (see Table 2 for these last parameters), we hence used the following formula: NOI(cm-2)=8πkν2hc3105AulQTexexp(Eu/kTex)guTMBdν(Kkms-1).\begin{equation} N_{\rm OI}{\rm(cm^{-2})}\! =\! \frac{8\pi k \nu^2}{hc^3} \frac{10^5}{A_{\rm ul}} \frac{Q_{\rm Tex} \exp(E_{\rm u}/kT_{\rm ex})} {g_{\rm u}} \!\int\! T_{\rm MB} {\rm d} \varv\, \rm{(K~km~s^{-1})}. \label{equation1} \end{equation}(1)In their analysis of CO lines excitation, L15 identified a low-excitation and a high-excitation component in the jet and the outflow cavity, and a high-excitation component was associated with HH377. The main properties inferred by L15 for each of these components are listed in Table 5. In order to calculate a column density and column density ratio for O i, it is necessary to understand the component of origin for the [OI]63μm emission. Unfortunately, the three kinematic components (jet, HH377, and outflow cavity) are all seen in low- and high-excitation lines (CO 2–1 and 16–15, respectively), and the upper level of our observed O i transition has spectroscopic parameters (Eup = 227.7 K, Aij = 8.9 × 10-5 s-1) similar to those of the CO Jup = 9 transition (Eup = 248.9 K, Aij = 7.3 × 10-5 s-1), where the transition from low- to high-excitation occurs in the CO rotational diagram of the jet and the outflow cavity. We consequently calculated the column density of O i assuming that its emission originates from both the low- and high-excitation components for the jet and the outflow cavity. Our results are given in Table 5.

Table 5

O i abundance calculations.

These calculations show that the O i column density is higher in the outflow cavity than in the jet, itself presenting a higher O i column density than the terminal shock of the Cep E outflow. In terms of O i quantity over the entire outflow, this effect is accentuated by the larger extension of the outflow cavity with respect to the jet, let alone to the HH377 region. The abundance ratio of O i to CO is also lowest in the terminal shock region (about 0.2, see Table 5). Finally, our strongest conclusion from Table 5 is that the jet component is atomic (N(O i)/N(CO) ~ 2.7) regardless of the assumption that we made on the origin of the [OI]63μm emission.

The LTE assumption is only validated in the jet, bow-shock, and outflow cavity under the high-excitation assumption by the local density values constrained by L15. These are indeed greater than the critical density of the [OI]63μm line, about 5 × 105 cm-3 (see Table 5). The LTE approximation is not so well validated under the low-excitation assumption, where the local densities are lower. This means that the O i column densities obtained in this way are most likely lower limits to the column densities and not definitive values. This does not crucially modify the conclusions of the above paragraph.

The optically thin assumption could also be questioned. We adopted a line width of 15 km s-1 for the jet component, consistent with our observations and those of L15. Using RADEX online (van der Tak et al. 2007), we found that the [OI]63μm line only becomes thick if N(O i) exceeds 3.1 × 1018 and 6.3 × 1018 cm-2 in the low- and high-excitation conditions, respectively, that is, a factor 10 to 100 times higher than the O i column densities we found.

The adequate remedy to improve our column density measurement would be to map the region also in the 3P03P1 transition of O i (at 145 μm). This transition is indeed likely to be optically thin, and we could hence apply an LVG treatment to these two lines. We would thus overcome the limitations intrinsic to the LTE approach and be able to study the origin of the [OI]63μm emission (from the low- or high-excitation CO component). SOFIA/upGREAT will allow this observation in a near future. The [OI]63μm/[OI]145μm line ratio predicted by Nisini et al. (2015) would make this detection difficult. They are in the 25–50 range (depending on the kinetic temperature and main collision partner, see their Fig. 10) for the values of the local density of the high-excitation component inferred by L15. On the other hand, the ratios we found with RADEX online are significantly lower (of the order of 1.5–3), and the 3P03P1 transition of O i at 145 μm has been successfully detected by Herschel-PACS.

5. Nature of shocks in Cepheus E

5.1. Shock modelling of the jet component

In L15, we successfully reproduced the CO emission in the high-excitation component of the jet by means of shock models computed with the Paris-Durham code. We compared the observations to a grid of such models and found several solutions depending on filling factor assumptions of this high-excitation component. We recall the solutions here and compare them to our [OI]63μm observations in an attempt to lift the degeneracy of the solution by fitting a maximum number of observational data.

Without a filling factor correction in addition to the correction that is based on the interferometric observation of CO (2–1) emission, six models were found to provide a satisfying fit to the CO observations of the jet. All of them were young non-stationary (Lesaffre et al. 2004a,b) shock models with the following parameters: pre-shock density nH between 5 × 104 or 105 cm-3, magnetic field parameter b = 1 (B(μG) = b [nH(cm-3)] 1/2, where B is the transverse magnetic field strength), shock velocity νs between 20 and 30 km s-1, and age between 225 and 435 yr. When we applied an additional filling factor of 0.25 to the CO data (effectively multiplying the integrated intensities by a factor 4) to account for the possibility that the high-excitation component is less extended than the bulk of the CO emission, six older models fit the data, with nH = 105 cm-3, b = 1, νs = 25–30 km s-1, and age between 595 and 740 yr. These models and their associated [OI]63μm flux are given in Table 6.

Table 6

Predictions for the [OI]63μm line emission from the models that were found to be a satisfying fit to the high-excitation CO emission from the jet by L15.

The first conclusion is that the predicted value is very dependent on the shock parameters, as our results are spread over two orders of magnitude, between 2.9 × 10-6 and 2.8 × 10-4 erg cm-2 s-1 sr-1. The most important parameter to explain this spread is the age of the shock, which determines the strength of the J-type contribution, and that of the [OI]63μm line emission in our non-stationary models. The observed value in the jet component is 2.3 ± 0.5 × 10-3 erg cm-2 s-1 sr-1. If an additional filling factor of 0.25 is applied, similarly to what was considered for CO by L15, these values should be multiplied by 4. In any case, we find that none of our shock models approaches the observed value, as they predict at least a factor 10 less O i emission than observed. The first explanation is that the bulk of the [OI]63μm line emission comes from the (not shock-modelled) low-excitation component instead of the high-excitation component.

However, if the [OI]63μm line emission comes from the high-excitation component, a number of explanations could account of the discrepancy between observations and shock models. A first explanation is that our model is inaccurate regarding the oxygen chemistry. The main oxygen carriers are CO, H2O, OH, O2, and O i. The models that best fit the CO emission predict that the local fractional abundance of both CO and H2O are about 10-4 in the post-shock region and that the local fractional abundance of O2 is greater than that of O i. The high value of the local fractional abundance of O2 is a well-documented problem in our simulation and in chemical models in general (see e.g. Viti et al. 2001; and Hollenbach et al. 2009, for a discussion of this specific aspect). Meaningful constraints on the oxygen chemistry could then be obtained through the observation of H2O and O2. Unfortunately, no velocity-resolved observations of water have been performed in Cep E by Herschel, and no instrument can allow them at the moment. O2 is best observed from space observatories (the NJ = 33−12 transition at 487.249 GHz is close to an atmospheric absorption feature), but it has not been observed in Cep E, and its symmetry makes its detection difficult from the outset.

Another possibility is that the shock that we observed is illuminated by the protostar that drives the outflow. In such irradiated shocks (e.g. Lesaffre et al. 2013; Melnick & Kaufman 2015), both H2O and CO can be photo-dissociated, increasing the gas-phase abundance of O i. This effect is not taken into account in the Paris-Durham code. However, the possibility that shocks at position BII could be affected by the distant protostar is relatively low: the BII and mm position are 8.4 × 10-2 pc away. The Cep E-mm protostar has a luminosity of about 100 L, and a T-Tauri star typically emits between 10-1.5 and 10-3.5 of its total luminosity in the UV range between 3200 and 5200 Å (Gullbring et al. 1998). Conservatively assuming that Cep E-mm is radiating 10% of its total luminosity in this UV range and taking dilution into account, we found that the radiation field in BII is about 4.8 G0 in units of Draine’s (Draine 1978) or Mathis’s (Mathis et al. 1983) interstellar radiation field (ISRF). This value has then to be corrected for the absorption by dust grains. We assumed an average local density in the outflow cavity of 10-4 cm-3, which translates into a column density of NH ≃ 2.6 × 1021 cm-2 and into a visual extinction of Av ≃ 1.4 mag (Lefloch et al. 1996 provided local estimates for this value: 61 mag around the driving protostar, and 3.2–3.4 mag in the H2 knots along the flow). Using the exponential attenuation approximation I = I0 × eτ, where τ = NH × σ, and σ is the effective attenuation cross section from Draine & Bertoldi (1996), we found that the radiation field really felt in BII is about 2.6 × 10-3G0 in Draine’s or Mathis’s ISRF units.

Another explanation is that several shock structures are present in our beam. In this case, the missing O i emission could be accounted for by the presence of another shock structure. This additional shock layer must present different physical characteristics, as simply combining shock layers similar to those we used would simultaneously increase the [OI]63μm emission insufficiently and increase the CO emission too much, producing another type of discrepancy between observations and models. One such shock structure could be a stationary dissociative J-type shock. In particular, a J-type shock with nH = 104 cm-3, magnetic field parameter b = 0.1, and shock velocity νs ≥ 30 km s-1 (Flower & Pineau des Forêts 2015) could be combined with our current best-fit solution, providing a satisfying fit for the O i and CO emission.

The additional shock structure could also be a shock with a radiative precursor, or in other words, a shock structure whose high temperatures generate FUV-illumination that affects both pre- and post-shock regions (Hollenbach & McKee 1989). Figure 7 of Hollenbach & McKee 1989 shows that these types of shocks (especially with a pre-shock density nH = 104 cm-3) could reproduce the observed levels of the [OI]63μm line emission. However, these shocks also produce CO emission (see Fig. 6 of Hollenbach & McKee 1989): combining this solution with our models would perhaps prove successful at fitting the [OI]63μm emission, but would generate a discrepancy in the CO comparisons. On the other hand, it is possible that the CO and [OI]63μm line emission could both arise from the same shock structure with a radiative precursor.

All these explanations (including the possibility that the [OI]63μm emission comes from the low-excitation component) deserve to be explored in more detail by comparing a maximum of observables with the models. In this respect, the important role of SOFIA/upGREAT can only be stressed, as the observation of the O i 3P03P1 line at 145 μm and the 2P3/22P1/2 C ii line at 158 μm will allow placing more meaningful constraints on our comparisons between observations and models.

5.2. Comparison with other star-forming regions

We recall that the BII spectrum has been obtained at a position that is distant from the protostar (r = 0.084 pc), which is located at d = 730 pc from the Earth and has a luminosity L = 100 L. To compare our observations with other regions where [OI]63μm and CO (16–15) have been observed in similar conditions, we plot the two spectra and their channel-by-channel ratio in the upper panel of Fig. 3. In this panel, the ratio is shown as green points for every channel where the two lines have been detected above the 3σ level. We extracted the same spectra and ratio from two other star-forming regions. The first is the R1 shock spot of the NGC1333-IRAS4A bipolar outflow, which is driven by a low-mass protostar (for which r = 0.011 pc, d = 235 pc, L = 9.1L; Kristensen et al. 2017). In this region, the [OI]63μm filling factor is of the order of unity, as shown by the half-maximum contour of Fig. 2 of Nisini et al. (2015). The second region is the central position of G5.89–0.39, a complex high-mass star-forming region with several outflows and spectacular [OI]63μm emission (for which r = 0 pc, d = 1.28 kpc, L = 1.3 × 105L; Leurini et al. 2015). In this region, the [OI]63μm filling factor is also of the order of unity given the size of the emiting region. The comparison between the three datasets is shown in Fig. 3.

thumbnail Fig. 3

[OI]63μm (red line) and CO (16–15) (grey line) spectra, and {Tmb [OI63μm]}/{Tmb [CO (16–15)]} ratio (green points) obtained at nominal spatial resolutions, channel by channel, in NGC1333-IRAS4A RI (upper panel, from Kristensen et al. 2017), Cep E-BII (middle panel), the central position of the G5.89–0.39 massive star-forming region (from Leurini et al. 2015). The line ratio is only shown for channels where both lines have been detected above the 3σ threshold. Horizontal black lines at 0.5, 1, 1.5, and 2 values are shown in relation with this ratio as a guide for the eye in both panels.

The line temperature ratio is very different in the three cases. In NGC1333-IRAS4A, the [OI]63μm/CO (16–15) temperature ratio does not exceed the value of 0.6. In this case, only four points are shown that correspond to the four velocity channels where [OI]63μm emission is detected above 3σ. In Cep E-BII, the great majority of the ratio values is between 0.5 and 1.5, with lower values found in the outflow cavity component near the systemic velocity. Finally, in G5.89–0.39, the ratio is lower than 2 only at velocities close to the systemic velocity (between −5 and 40 km s-1 for a νlsr ~ 10 km s-1). Given that the [OI]63μm emission is self-absorbed at these velocities (Leurini et al. 2015), these points do not reflect the local shock conditions. Moreover, it is possible that self-absorption affects the emission for velocities between 20 to 50 km s-1. It follows that the ratio is systematically in the 5–20 range. Although not corrected for filling factor effects, this is significantly higher than what is seen in Cep E and NGC1333-IRAS4A.

Depending on filling factor corrections, this [OI]63μm/CO (16–15) temperature ratio could be reflective of the nature of shocks in the different sources:

  • in the low-mass bipolar outflow NGC1333-IRAS4A theobserved line width seems to suggest that the shock propagates ata moderate velocity and is hence likely molecular. Additionally,the region is offset from the (low-luminosity) protostar, thus notirradiated by its far-UV-illumination. In such conditions,non-dissociative C-, J- or CJ-type shock models could be used to interpret the observations;

  • in Cep E-BII, we have seen that those models are probably not sufficient to fit the observations in the jet, where additional dissociative shocks or alternative solutions with radiative precursors must be invoked;

  • finally, in G5.89–0.39, where the observations were centred on the very luminous protostar, the very high intensity of [OI]63μm emission could be a signature of the effect of the far-UV field of the central protostar on the various outflows that propagate in the region, all the more so because the filling-factor-corrected ratio is in the 1–5 range.

This global picture of the [OI]63μm/CO (16–15) interpretation needs to be tested through thorough comparisons between additional observations and various models. In particular, the observation of the 2P3/22P1/2 C ii line at 158 μm with SOFIA/upGREAT, combined with ground-based observations of CO and C i, will allow us to quantify dissociation effects in the various targeted regions.

6. Concluding remarks

  • 1.

    We have used the SOFIA/GREAT instrument to observe theemission of the O i 3P13P2, OH between Π1/22\hbox{$^2\Pi_{1/2}$}J = 3/2 and J = 1/2 at 1837.8 GHz and CO (16–15) lines towards three positions in the Cepheus E outflow (namely Cep E-mm, Cep E BI, and Cep E BII). The CO (16–15) line was detected at all three positions. The [OI]63μm line was detected in Cep E BI and BII, whereas the OH line was not detected.

  • 2.

    We focused our study on the most significant detection in Cep E BII. Three kinematical components were detected in the [OI]63μm spectrum, which are remarkably similar to those seen in CO (16–15). These are themselves similar to those seen by L15 in other CO transitions, and related to various spatial components: the jet, the HH377 terminal bow-shock, and the cavity walls all contribute to the line emission. This allowed us to estimate “most realistic” filling factor values for the [OI]63μm emission in each of these components based on the observations presented in L15: 0.25 ± 0.05, 0.55 ± 0.15, and 1.

  • 3.

    We used the results of L15 and analysed the [OI]63μm emission to infer O i column densities and N(O i)/N(CO) ratios in each of these components. The two quantities depend on the filling factor correction and on the characteristics of the region of origin for the [OI]63μm emission. The O i column density is higher in the outflow cavity than in the jet, which itself presents a higher O i column density than the terminal shock. In terms of O i quantity over the entire outflow, this effect is accentuated by the larger extension of the outflow cavity with respect to the jet, let alone to the HH377 region. The terminal shock is also the region where the abundance ratio of O i to CO is the lowest (about 0.2). The jet component is atomic (N(O i)/N(CO) ~ 2.7) regardless of the assumption that we made about the origin of the [OI]63μm emission.

  • 4.

    We compared the observed [OI]63μm flux in the jet component with shock models from the Paris-Durham code that proved successful to fit the high-excitation emission of the jet component of 10 CO lines. All of these models are at least a factor 10 below the observed [OI]63μm flux. This discrepancy might indicate that the O i emission comes from the low-excitation component. This discrepancy could also be a sign that the current chemical modelling of oxygen is not accurate enough, or that the shocks are irradiated by the protostar. Two more compelling explanations exist for this discrepancy: several shock structures (including dissociative J-type shocks) could be present in the beam of our observations, or the CO and O i emission could arise from a J-type shock with a radiative precursor.

  • 5.

    We compared our observations in the Cep E-BII region with similar observations of other star-forming regions: one shock position in a bipolar outflow driven by a low-mass protostar, NGC1333-IRAS4A, and one position centred on a massive protostar at the centre of several outflow structures in G5.89–0.39. The trend observed in the channel-by-channel [OI]63μm/CO (16–15) temperature ratio (which increases with the mass of the central protostar) could be a signature of various types of shocks propagating in these regions.

  • 6.

    Most of the limitations that we have encountered in the present study will find their remedy in the next round of SOFIA/upGREAT observations that could allow us to map the entire outflow in O i lines (the 3P13P2 line at 63 μm and the 3P03P1 line at 145 μm). This will allow us to improve the accuracy of all our measurements. We would then better be able to constrain the size of emitting regions, the excitation conditions for O i, and also conclude more precisely about energetic impacts by extending our conclusions to the full outflow. Finally, the observation of the 2P3/22P1/2 C ii line at 158 μm with SOFIA/upGREAT is also an interesting prospect, as it will allow us to more tightly constrain the types of shocks that can fit the data.


1

GREAT is a development by the MPI für Radioastronomie and the KOSMA/Universität zu Köln, in cooperation with the MPI für Sonnensystemforschung and the DLR Institut für Planetenforschung.

Acknowledgments

We thank an anonymous referee for meaningful comments that helped to improve the clarity and quality of this study. We thank the SOFIA operations and the GREAT instrument teams, whose support has been essential for the GREAT accomplishments, and the DSI telescope engineering team. Based [in part] on observations made with the NASA/DLR Stratospheric Observatory for Infrared Astronomy. SOFIA Science Mission Operations are conducted jointly by the Universities Space Research Association, Inc., under NASA contract NAS2-97001, and the Deutsches SOFIA Institut, under DLR contract 50 OK 0901. B. Lefloch acknowledges support from a grant from LabeX Osug2020 (Investissements d’avenir – ANR10LABX56). This work was supported by the French program “Physique et Chimie du Milieu Interstellaire” (PCMI) funded by the Centre National de la Recherche Scientifique (CNRS) and Centre National d’Études Spatiales (CNES).

References

  1. Arce, H. G., Shepherd, D., Gueth, F., et al. 2007, Protostars and Planets V, 245 [Google Scholar]
  2. Bachiller, R., Pérez Gutiérrez, M., Kumar, M. S. N., & Tafalla, M. 2001, A&A, 372, 899 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Benedettini, M., Gusdorf, A., Nisini, B., et al. 2017, A&A, 598, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Büchel, D., Pütz, P., Jacobs, K., et al. 2015, IEEE Trans. Terahertz Science and Technology, 5, 207 [NASA ADS] [CrossRef] [Google Scholar]
  5. Chini, R., Ward-Thompson, D., Kirk, J. M., et al. 2001, A&A, 369, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Draine, B. T. 1978, ApJS, 36, 595 [NASA ADS] [CrossRef] [Google Scholar]
  7. Draine, B. T., & Bertoldi, F. 1996, ApJ, 468, 269 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Flower, D. R., & Pineau des Forêts, G. 2015, A&A, 578, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Frank, A., Ray, T. P., Cabrit, S., et al. 2014, Protostars and Planets VI, 451 [Google Scholar]
  10. Gómez-Ruiz, A. I., Gusdorf, A., Leurini, S., et al. 2012, A&A, 542, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Guan, X., Stutzki, J., Graf, U. U., et al. 2012, A&A, 542, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Gueth, F., Guilloteau, S., & Bachiller, R. 1996, A&A, 307, 891 [NASA ADS] [Google Scholar]
  13. Gueth, F., Guilloteau, S., & Bachiller, R. 1998, A&A, 333, 287 [NASA ADS] [Google Scholar]
  14. Gullbring, E., Hartmann, L., Briceño, C., & Calvet, N. 1998, ApJ, 492, 323 [NASA ADS] [CrossRef] [Google Scholar]
  15. Heyminck, S., Graf, U. U., Güsten, R., et al. 2012, A&A, 542, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Hollenbach, D., & McKee, C. F. 1989, ApJ, 342, 306 [NASA ADS] [CrossRef] [Google Scholar]
  17. Hollenbach, D., Kaufman, M. J., Bergin, E. A., & Melnick, G. J. 2009, ApJ, 690, 1497 [CrossRef] [Google Scholar]
  18. Kristensen, L. E., van Dishoeck, E. F., Benz, A. O., et al. 2013, A&A, 557, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Kristensen, L. E., Gusdorf, A., Mottram, J. C., et al. 2017, A&A, 601, L4 [CrossRef] [EDP Sciences] [Google Scholar]
  20. Lefloch, B., Eisloeffel, J., & Lazareff, B. 1996, A&A, 313, L17 [NASA ADS] [Google Scholar]
  21. Lefloch, B., Cernicharo, J., Pacheco, S., & Ceccarelli, C. 2011, A&A, 527, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Lefloch, B., Gusdorf, A., Codella, C., et al. 2015, A&A, 581, A4 [CrossRef] [EDP Sciences] [Google Scholar]
  23. Lesaffre, P., Chièze, J.-P., Cabrit, S., & Pineau des Forêts, G. 2004a, A&A, 427, 147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Lesaffre, P., Chièze, J.-P., Cabrit, S., & Pineau des Forêts, G. 2004b, A&A, 427, 157 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Lesaffre, P., Pineau des Forêts, G., Godard, B., et al. 2013, A&A, 550, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Leurini, S., Wyrowski, F., Wiesemeyer, H., et al. 2015, A&A, 584, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, 128, 212 [NASA ADS] [Google Scholar]
  28. Melnick, G. J., & Kaufman, M. J. 2015, ApJ, 806, 227 [NASA ADS] [CrossRef] [Google Scholar]
  29. Nisini, B., Santangelo, G., Giannini, T., et al. 2015, ApJ, 801, 121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Ott, S. 2010, in Astronomical Data Analysis Software and Systems XIX, eds. Y. Mizumoto, K.-I. Morita, & M. Ohishi, ASP Conf. Ser., 434, 139 [Google Scholar]
  31. Pickett, H. M., Poynter, R. L., Cohen, E. A., et al. 1998, J. Quant. Spectr. Rad. Transf., 60, 883 [Google Scholar]
  32. Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
  33. Podio, L., Codella, C., Gueth, F., et al. 2016, A&A, 593, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Richter, H., Wienold, M., Schrottke, L., et al. 2015, IEEE Trans. Terahertz Science and Technology, 5, 539 [Google Scholar]
  36. Risacher, C., Güsten, R., Stutzki, J., et al. 2016, IEEE Trans. Terahertz Science and Technology, 6, 199 [Google Scholar]
  37. Sargent, A. I. 1977, ApJ, 218, 736 [Google Scholar]
  38. Tafalla, M., Bachiller, R., Lefloch, B., et al. 2015, A&A, 573, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. van der Tak, F. F. S., Black, J. H., Schöier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, A&A, 468, 627 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Viti, S., Roueff, E., Hartquist, T. W., Pineau des Forêts, G., & Williams, D. A. 2001, A&A, 370, 557 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Herschel-PACS footprint of the [OI]63μm emission in the southern lobe of the Cep E outflow

In this appendix, we show a footprint of the [OI]63μm emission obtained in the southern lobe of the Cep E outflow (Fig. A.1) with the PACS receiver (Poglitsch et al. 2010) onboard the Herschel space observatory (Pilbratt et al. 2010). The image corresponds to the ObsId 1342262543 and was obtained in line spectroscopy mode. It was created using HIPE version 14.2.1 (Ott 2010). We solely use the image for the purpose of discussing our filling factor assumptions. The more thorough presentation and deeper analysis of all PACS data (including water line emission) will be the subject of a forthcoming publication.

thumbnail Fig. A.1

Spitzer-IRAC band-one (3.6 μm) image of the southern lobe of the Cep E outflow (in colours), retrieved from the Spitzer archive, overlaid with the CO (2–1) cyan contours of the jet (PdBI data, from L15), and with the [OI]63μm emission (yellow contours) seen by Herschel-PACS. The HPBW contour is plotted in purple and roughly corresponds to a 14′′ circle shown as a white dashed contour. The red squares indicate the position of the spaxels, and the black star indicates the position of HH377, the terminal bow-shock.

The half-power contour shown in purple in Fig. A.1 has a typical diameter of 14′′ (also see the dashed white contour in the figure). This indicates a typical size of 10′′ for the emitting region, after rudimentary deconvolution from the under-sampled PACS footprint, assuming each spaxel has a Gaussian size of 9\hbox{$\farcs$}4. This value is indicative of the emission size that is most likely dominated at this position by the outflow cavity.

All Tables

Table 1

Observed positions and sensitivity (rms of baseline noise) at 1 km s-1 spectral resolution for each observed line.

Table 2

Characteristics and observational parameters of the OH transitions between the Π1/22\hbox{$^2\Pi_{1/2}$}J = 3/2 and J = 1/2 states, and of the [OI]63μm and CO (16–15) lines.

Table 3

Raw integrated intensities extracted from the [OI]63μm and CO (16–15) spectra, in units of K km s-1, and corresponding line ratio, kinematic component by kinematic component (jet, bow-shock, outflow cavity).

Table 4

Filling-factor-corrected, integrated intensities extracted from the [OI]63μm and CO (16–15) spectra, in units of K km s-1 and erg cm-2 s-1 sr-1, and corresponding line ratio, kinematic component by kinematic component (jet, bow-shock, outflow cavity).

Table 5

O i abundance calculations.

Table 6

Predictions for the [OI]63μm line emission from the models that were found to be a satisfying fit to the high-excitation CO emission from the jet by L15.

All Figures

thumbnail Fig. 1

Left: a finder’s chart for the positions of interest in the Cep E protostellar outflow: Spitzer-IRAC band-two (4.5 μm) image, retrieved from the Spitzer archive. Grey circles mark the diffraction-limited SOFIA beams at the frequencies of [OI]63μm (smaller circles) and CO (16–15) or OH at 1838 GHz (larger circles) lines. The three observed positions are indicated (Cep E-mm, red star, BI and BII, white rectangles), as well as others in the red lobe (RI and RII). The knots denomination is from Gómez-Ruiz et al. (2012). Right: zoom on the southern lobe of the outflow showing the spatial components identified by L15 in CO (2–1) (see Sect. 3): the jet (cyan contours), the terminal bow-shock (yellow contours distant from the protostar) and the cavity walls (coloured background). The three observed positions are indicated with a white star, and the SOFIA beams at the frequencies of [OI]63μm (smaller circles) and CO (16–15) or OH at 1838 GHz (larger circles) lines are shown in BII.

In the text
thumbnail Fig. 2

CO (16–15) (grey histograms), [OI]63μm (red line), and OH at 1838 GHz (blue line) observations obtained in Cep E-mm (top panel), Cep E-BI (middle panel), and Cep E-BII (lower panel). The vertical dashed lines represent the central value for the velocity components identified by L15 in CO transitions: 125 (cyan), 68 (pink), and 10.9 km s-1 (black; systemic velocity of Cep E). The spectral resolution is 1 km s-1 for the Cep E-BI and BII spectra, and 4 km s-1 for Cep E-mm.

In the text
thumbnail Fig. 3

[OI]63μm (red line) and CO (16–15) (grey line) spectra, and {Tmb [OI63μm]}/{Tmb [CO (16–15)]} ratio (green points) obtained at nominal spatial resolutions, channel by channel, in NGC1333-IRAS4A RI (upper panel, from Kristensen et al. 2017), Cep E-BII (middle panel), the central position of the G5.89–0.39 massive star-forming region (from Leurini et al. 2015). The line ratio is only shown for channels where both lines have been detected above the 3σ threshold. Horizontal black lines at 0.5, 1, 1.5, and 2 values are shown in relation with this ratio as a guide for the eye in both panels.

In the text
thumbnail Fig. A.1

Spitzer-IRAC band-one (3.6 μm) image of the southern lobe of the Cep E outflow (in colours), retrieved from the Spitzer archive, overlaid with the CO (2–1) cyan contours of the jet (PdBI data, from L15), and with the [OI]63μm emission (yellow contours) seen by Herschel-PACS. The HPBW contour is plotted in purple and roughly corresponds to a 14′′ circle shown as a white dashed contour. The red squares indicate the position of the spaxels, and the black star indicates the position of HH377, the terminal bow-shock.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.