Spitzer spectral line mapping of the HH211 outflow
O. Dionatos, B. Nisini, S. Cabrit, L. Kristensen, G. Pineau des Forets
AAstronomy & Astrophysics manuscript no. HH211˙astro˙ph c (cid:13)
ESO 2018November 5, 2018
Spitzer spectral line mapping of the HH211 outflow
O. Dionatos , , B. Nisini , S. Cabrit , L. Kristensen , and G. Pineau des Forˆets , Natural History Museum, University of Copenhagen, Øster Voldgade 5-7 1350 Copenhagen, Denmarke-mail: [email protected] INAF - Osservatorio Astronomico di Roma, Via di Frascati, 33 00040 Monte Porzio Catone (RM), Italye-mail: [email protected] LERMA, Observatoire de Paris, UMR 8112 du CNRS, 61 Avenue de Observatoire, 75014 Paris, Francee-mail: [email protected] Leiden Observatory, Leiden University, Niels Bohrweg 2, 2300 CA Leiden, The Netherlandse-mail: [email protected] Institut Astrophysique Spatiale (IAS), UMR 8617, CNRS, Universite Paris-Sud 11, Batiment 121, 91405 Orsay Cedex, Francee-mail: [email protected]
Preprint online version: November 5, 2018
ABSTRACT
Context.
Jets from the youngest protostars are often detected only at mm wavelengths, through line emission of CO and SiO. However,it is not yet clear if such jets are mostly molecular or atomic, nor if they trace ejected gas or an entrained layer around an embeddedatomic jet.
Aims.
We investigate the warm gas content in the HH211 protostellar outflow to assess the jet mass-flux in the form of H and probefor the existence of an embedded atomic jet. Methods.
We employ archival Spitzer slit-scan observations of the HH211 outflow over 5.2 – 37 µ m obtained with the low resolutionIRS modules. Detected molecular and atomic lines are interpreted by means of emission line diagnostics and an existing grid ofmolecular shock models. The physical properties of the warm gas are compared against other molecular jet tracers and to the resultsof a similar study towards the L1448-C outflow. Results.
We have detected and mapped the v = lines as well as fine-structure lines of S, Fe + , and Si + . The H isdetected down to 5 (cid:48)(cid:48) from the source and is characterized by a ”cool” T ∼ ∼ ±
300 K component, with anextinction A V ∼ towards the jet agrees with that estimated from CO assuming fully molecular gas. Thewarm component is well fitted by C–type shocks with a low beam filling factor ∼ .The fine-structure line emission arises from dense gas with ionization fraction ∼ . − × − , suggestive of dissociative shocks.Line ratios to sulfur indicate that iron and silicon are depleted compared to solar abundances by a factor ∼ Conclusions.
Spitzer spectral mapping observations reveal for the first time a cool H component towards the CO jet of HH211consistent with the CO material being fully molecular and warm at (cid:39)
300 K. The maps also reveal for the first time the existence ofan embedded atomic jet in the HH211 outflow that can be traced down to the central source position. Its significant iron and silicondepletion excludes an origin from within the dust sublimation zone around the protostar. The momentum-flux seems insu ffi cient toentrain the CO jet, although current uncertainties on jet speed and shock conditions are too large for a definite conclusion. Key words.
STARS: FORMATION, ISM: JETS AND OUTFLOWS, ISM: INDIVIDUAL OBJECTS: HH211-mm, INFRARED:ISM: LINES AND BANDS
1. Introduction
The process of mass accretion leading to the formation of lowmass protostars is always associated with ejection of materialin the form of well collimated high-velocity jets and / or widerslower bipolar outflows (Reipurth & Bachiller 1997; Arce et al.2007). Accretion and ejection phenomena are believed to be in-timately connected through the existence of a magnetized accre-tion disk around the forming protostar, but there is no consensusso far about the exact ejection mechanism (Pudritz et al. 2007;Shang et al. 2007; Ferreira et al. 2006). The study of the physicalproperties of jets (e.g. temperature, density, abundances, massflux) is therefore crucial to obtain indirect constraints on themass accretion and ejection processes, especially in the case ofprotostars in the earliest evolutionary stages - the so called Class 0 protostars - where the launching zone is heavily embedded ina dense cocoon of dust and gas.Due to this high extinction, jets from the youngest Class 0sources cannot be seen in the optical and are mostly traced innear-IR H and in mm CO and SiO lines (Guilloteau et al. 1992;McCaughrean et al. 1994; Davis & Smith 1996). However, it isstill unclear whether the molecules trace a sheath of entrainedambient material around an unseen underlying atomic jet, or ifthey trace the primary jet material itself. In the latter case, thechemical composition of the jet (molecular vs. atomic fraction,dust depletion) holds crucial clues to the ejection zone, as thestrong UV flux generated by the accretion shock should destroydust grains and H molecules in the innermost regions of thedisk. A wind from such inner zones would be mostly atomic even a r X i v : . [ a s t r o - ph . GA ] J un O. Dionatos et al.: Spitzer spectral line mapping of the HH211 outflow if CO and SiO are abundant — unless the mass flux is unusuallylarge (Ruden et al. 1990).Mid-infrared spectroscopic observations are essential toprogress on these questions. A Spitzer study of the innermostjet regions of the Class 0 L 1448 outflow (Dionatos et al. 2009)revealed an underlying deeply embedded atomic / ionic compo-nent, not seen in the near-IR range, as well as mid-IR H emis-sion tracing the warm molecular content. Here we present a moreextensive Spitzer mapping study of the HH211 outflow, a partic-ularly interesting target that appears as a ”text-book” example ofjet-driven flow.HH211 is located in the IC 348 complex in Perseus at adistance estimated to be between 250 pc (Enoch et al. 2006),adopted in the current paper, and ∼
320 pc (Herbig 1998; Ladaet al. 2006). The driving source HH211-mm is detected in mmcontinuum emission (e.g. Gueth & Guilloteau 1999; Lee et al.2007), and has been classified as a low-mass and low-luminosityClass 0 young stellar object (Froebrich 2005).The outflow was discovered by McCaughrean et al. (1994)in near-IR H emission tracing 2000 K hot, shocked gas in twobright symmetric bowshocks separated by 0.13pc with a faintchain of knots in between. A few compact knots of atomic emis-sion at optical and near-IR wavelengths (i.e. H α , [SII]6730Åand [FeII]1.64 µ m) are also seen near the brightest H peaks(O’Connell et al. 2005; Caratti o Garatti et al. 2006; Walawenderet al. 2005, 2006).The absence of any further shock structures on larger scalesmakes it the youngest outflow known so far, with a kinematicalage of only 1000 yr × ( V / − (Eisl¨o ff el et al. 2003a;Gueth & Guilloteau 1999).Low velocity CO J = = bow-shocksand its flanks connecting back to the driving source, while thehigh-velocity components of the same transitions reveal an in-ner well-collimated CO jet extending out to ± (cid:48)(cid:48) from thesource. The CO jet is also traced in SiO through interferomet-ric observations of the J = J levels implies that the SiOjet is warmer and denser than the outflow cavity traced by low-velocity CO emission.All observations reveal a highly bipolar structure, with blueand red-shifted lobes pointing southeast and northwest from thedriving source. The two lobes are well separated and such mor-phology combined with the low observed radial velocities ( ∼ − , Gueth & Guilloteau 1999) implies that the inclination ofthe outflow is less than 10 o from the plane of the sky (Gueth &Guilloteau 1999; Chandler & Richer 2001; Lee et al. 2007). Thissimple geometry is favorable to modelling. Indeed, the CO cav-ity shape is well fitted by a dynamical model of jet-driven bow-shock propagating into a medium of decreasing density (Gueth& Guilloteau 1999). This interpretation is supported by fittingthe near-IR line fluxes and morphologies at the apex of the red-shifted lobe with a series of bowshocks (O’Connell et al. 2005).A Spitzer map obtained at the tip the blue lobe indicates a higherexcitation 40 km / s dissociative shock on that side, with a sub-stantial UV flux (Tappe et al. 2008).In this paper, we present archival Spitzer-IRS slit-scan ob-servations covering most of the HH211 flow with low spectralresolution (Section 3), which we analyse to address the physicalconditions and the molecular and dust content. In Section 4, thedetected rotational H transitions are employed as probes of the
10 20 3000.0500.511.500.5 ambientoutflowon-source F l u x ( J y ) λ (μm) Fig. 1.
Continuum subtracted spectrum, integrated over a regionof 612 square arcseconds covering the HH211 flow ( middle ),shown in comparison with a spectrum of the o ff -outflow regionof the same area ( top ). H emission from the 0–0 S(0)– S(7)lines is detected in the outflow, together with atomic emissionfrom the ground state transitions of [FeII],[SiII] and [SI]. TheH S(0) and S(1) lines are detected also in the o ff -outflow spec-trum, indicating the presence of a di ff use line emission compo-nent not related to the HH211 flow. Di ff use emission from PAHfeatures, such as the 11.3 µ m feature, is also evidenced. ( Bottom: )On source spectrum extracted from an area equal to the LL pixelscale (110.25 square arcseconds), after removing the zodiacallight and ambient contribution presenting the intrinsic contin-uum emission from the outflow source.physical conditions in the jet close to the driving source and theshocked gas further out, and compared with an existing grid of Cand J-type shock models. Atomic lines present in our spectra aretreated separately as probes for the existence of a deeply embed-ded atomic jet, and used to measure dust depletion. Mass-fluxesare inferred for the H and atomic components, and compared tothose measured in CO and SiO, as well as to a similar study ofthe inner jet of L1448. Our conclusions are presented in Section5.
2. Observations and data reduction
Observations were obtained with the
Spitzer satellite (Werneret al. 2004) and retrieved from the Spitzer Public Data Archiveusing the Leopard software. They were performed as part of the”Shock dissipation in Nearby Star Forming Regions” programconducted by J. Bally (P.I.). In these, the low resolution modules(R ∼ Spitzer Infrared Spectrograph (IRS, Houcket al. 2004) were used in slit-scan mode to cover an area of57 (cid:48)(cid:48) × (cid:48)(cid:48) and 157 (cid:48)(cid:48) × (cid:48)(cid:48) for the Short Low (SL) and LongLow (LL) modules respectively, centered on the Class 0 sourceHH211-mm ( α J :03h43m56.5s, δ J : + . Dionatos et al.: Spitzer spectral line mapping of the HH211 outflow 3 H S(2) H S(4) H S(6)H S(7)H S(5)H S(3) + + ++++
B1 R1 R2B2B3
Fig. 2. H v = Spitzer
SL module, reconstructed with apixel scale of 3.5 (cid:48)(cid:48) . Contours are from 10 − W cm − sr − with an increment of 8 10 − W cm − sr − for the para lines (even J ),and from 2 10 − W cm − sr − with an increment of 2.5 10 − W cm − sr − for the ortho transitions (odd J ). Crosses indicate thedriving source position at α J :03h43m56.5s, δ J : + S(5) map, points B1-B3 and R1, R2 indicate peaksof emission, where A V and ortho to para ratio are estimated and excitation analysis is carried out.of both IRS low resolution modules gives a complete wavelengthcoverage of 5.2 – 37.0 µ m.Initial data processing was performed at the Spitzer ScienceCenter using Version 15 of the processing pipeline. Spectraldata-cubes were built using the CUBISM software package(Smith et al. 2007) and bad / rogue pixels were masked throughvisual inspection.Figure 1 presents the extracted continuum-subtracted spec-trum, integrated over a region encompassing the HH211 flow.The full series of H pure rotational lines (S(0)-S(7)) were de-tected, along with atomic and ionic lines from the fundamentaltransitions of [SI], [FeII], and [SiII] at 25 µ m, 26 µ m, and 35 µ mrespectively. For comparison, Fig. 1 shows also a spectrum of aregion encompassing the same area but in a direction perpendic-ular to the outflow. The H S(0) and S(1) lines are detected alsoin this o ff -outflow spectrum, indicating the presence of a di ff useline emission component not related to the HH211 flow. Di ff useemission from PAH features, such as the 11.3 µ m feature, is alsoevidenced. In Figure 1 we also plot the on source spectrum foran area equal to the LL pixel size, after subtracting the contribu-tion from an o ff -source, free of line-emission position. This onlyremoves zodiacal light and reveals the intrinsic mid-IR contin-uum from the central source, which is typical of low-luminosityClass 0 sources.Subsequent analysis consisted in the construction of individ-ual line emission maps, using a home-built pipeline; in this, foreach spatial pixel of the data-cube, the brightness of each spec-tral line of interest was calculated by Gaussian fitting after sub-tracting a local second order polynomial baseline. The resultingline intensity maps have a square pixel of side equal to the slitwidth of the IRS module, namely 3.5 (cid:48)(cid:48) and 10.5 (cid:48)(cid:48) for the SL andLL modules respectively, while the di ff raction limit of the tele-scope is 2.4 (cid:48)(cid:48) and 6.0 (cid:48)(cid:48) at λ = µ m and 25 µ m respectively. The astrometric accuracy of the maps was found to be good withinthe limits imposed by the pixel size of each IRS module.
3. Spectral line maps emission Figure 2 presents the emission line maps of the S(2) - S(7) ro-tational H lines observed with the SL - IRS modules with a3.5 (cid:48)(cid:48) sampling. Contours shape a characteristic bipolar outflowpattern, where the H emission is detected down to a projectedangular distance ∼ (cid:48)(cid:48) from the driving source. Further down-wind, peaks of emission which can be attributed to shocked gasare observed. In the H S(5) map of Fig. 2 we label these peaksas B1-B3 in the southeast, blue-shifted lobe and R1-R2 in thenorthwest, red-shifted lobe, respectively.Figure 3 shows the H S(5) emission map overlaid on thehigh-velocity CO J = → J = → J = → v = lines is quite goodexcept towards the broadest parts of the outflow cavity, wheremidIR emission appears to trail behind the near-IR one.In the southeastern - blue lobe, the peak of mid-IR H emis-sion closest to the source (B1) coincides with the 2.12 µ m H knot G in the nomenclature of McCaughrean et al. (1994), theBII knot of the CO jet (Gueth & Guilloteau 1999) and knotsBK4-6 in the SiO jet (Lee et al. 2007) (see lower panel of Fig.3). Therefore it appears to trace mostly the jet beam. Further out-wards, Peak B2 is centered in a hole behind the bowshock rimdelineated in the near-IR image, while Peak B3 appears to lie onthe bow rim (note however that the spatial extent of our SL mapdoes not fully cover the tip of this region). O. Dionatos et al.: Spitzer spectral line mapping of the HH211 outflow
B1B1 R1R1B2B2
B1B1B2B2B3B3 R1R1 R2R2
B1B1 R1R1B2B2
Fig. 3. (Ipper panel) H v = J = (Lower panels) Inner part of the H v = J = J = v = µ m image.In the blue-shifted jet lobe, point B1 is coincident with peak BIIin CO J = emis-sion close to the source forms an extended curving ”finger” trac-ing both the jet and the cavity wall north of it, as seen in the near-IR (see Fig. 3) and peak R1 coincides with the near-IR knot F ofMcCaughrean et al. (1994). However, peak R2 is again centeredon the cavity behind the bright bow rim traced in the near-IR.The spatial o ff set between B2 and R2 with respect to the near-IRcounterparts suggests that mid-IR emission in these broad re-gions may be tracing the outer bow-shock wings where gas isexpected to be in lower excitation conditions.The top row of Figure 4 presents the emission maps of theS(0) - S(1) rotational H lines observed with the LL - IRS mod-ules with a 10.5 (cid:48)(cid:48) sampling. Despite the lower angular resolutioncompared to the SL module, the H S (1) line clearly delineatesthe same outflow pattern as in Figure 3. As discussed in Sect.2, faint extended di ff use emission is also present, and becomesmostly evident at the north end of the map; di ff use emission iseven more apparent in the H S(0) line map of Fig. 4.This di ff use brightness could be caused by the presence of anextended photo-dissociation region (PDR) created at the cloudsurface by the illumination of FUV radiation from nearby stars.The existence of a di ff use PDR is strongly supported by the de-tection of di ff use PAH emission, which signals the presence ofUV radiation in the field. The H di ff use emission has a bright-ness of the order of 10 − erg / s / cm for both the S(0) and S(1)lines: according to the Kaufman et al. (1999) model, such an in-tensity is compatible with a PDR with a density of 10 cm − orlarger and a FUV field of the order of 10 times the average in-terstellar field, measured in units of G (i.e. in units of 1.6 × − erg cm − s − ). The most likely candidate for producing the dif- fuse PDR is the ø Per B0.5 star, located at a projected distanceof about 15 (cid:48) north of the HH211 system. The FUV field of thisstar (assuming an e ff ective temperature of 25 000 K, Hern´andezet al. (2005)), diluted for the distance, can account for a G ∼ emission also originates from other out-flows located north and south of HH211-mm, seen in both opti-cal and near-IR narrow band imaging of the region (Walawenderet al. 2006; Eisl¨o ff el et al. 2003b). However, no close spatial cor-respondence is evident at our spatial resolution. The bottom row of Figure 4 presents the emission line maps of[SiII]35 µ m, [SI] 25 µ m, and [FeII]26 µ m obtained from the LL -IRS modules with a 10.5 (cid:48)(cid:48) sampling. The emission from atomicand ionic lines is substantially weaker in comparison to the H lines. However, unlike H lines which are not detected closeto the driving source, all three fine structure transition contoursare continuous, evidencing emission very close to the protostarHH211-mm. Morphologically, the [SI] and [FeII] maps peak atthe blue lobe bow-shock . Emission from [SiII] instead showsstrong peaks on both the blue and red bow-shocks.The enhanced atomic - ionic emission on the blue-lobe bow-shock evidences a highly energetic interaction of the underlyingjet with the ambient material. In the same region, a wealth ofother molecular and ionic lines such as OH, H O, HD, [NeII]were mapped with Spitzer by Tappe et al. (2008), suggesting astrong UV flux typical of a 40 km s − dissociative shock. Thepresence of such high velocity shocks in this region is also con-sistent with the detection of optical H α and [SII] emission spots(Walawender et al. 2006, 2005). In the following, we will showthat rotational H lines over the HH 211 outflow appear to traceslower shocks, either internal to the jet or in bowshock wings.
4. Excitation conditions and mass fluxes emission The H pure rotational transitions are easily thermalized in pro-tostellar outflow environments owing to their low critical den-sities ( ∼ -10 cm − ) and can be used to probe temperaturesin the range between ∼
300 and 1500 K. The most direct wayto estimate the temperature of the emitting gas is by means ofan excitation diagram (EXD); this involves plotting the quanti-ties ln(N u , J / g u ) against E u , J , where N u , J and E u , J are the columndensity and the energy of the upper level, and g u the statisticalweight including the spin degeneracy (2S + . Considering LTEconditions, the values of these quantities for the observed H transitions should fall on a straight line, from the inverse slopeof which the excitation temperature can be derived. Deviationsfrom this linear correlation can be used to constrain other physi-cal properties, such as the H ortho-to-para ratio (hereafter, OTP)and the optical extinction A V , as described below. In the [FeII] 26 µ m map, contours are interrupted just before theactual blue-shifted bowshock peak; this is an artifact resulting from arogue pixel in the LL-IRS module that had to be masked out since itwas introducing unreliable brightness values forming a bright verticalstripe across the map.. Dionatos et al.: Spitzer spectral line mapping of the HH211 outflow 5 + ++ + R1+R2LL2LL1B2+B3 + H S(1) H S(0)[SiII] [SI] [FeII]
Fig. 4.
Intensity maps of the HH211 flow in the H S(1) and S(0) lines ( top row ) and the atomic [SiII]35 µ m, [SI] 25 µ m, and[FeII]26 µ m lines ( bottom row ) obtained with the Spitzer
IRS - LL module (pixel scale ∼ (cid:48)(cid:48) ). Crosses mark the position ofthe driving source HH211-mm. The dashed circle in the H maps masks the area around a bright star (IC348 IR) where line inten-sity measurements are unreliable. In the [FeII] 26 µ m map, contours are interrupted just before the eastern bowshock peak due toa rogue pixel. Contour levels in the H S(1) map start from 10 − W cm − sr − with 10 − increments, and from 7 10 − W cm − sr − with 5 × − increments for the S(0) map. Contours of the ionic line maps start at 10 − W cm − sr − with 10 − W cm − sr − steps for [SiII] and [FeII], and at 10 − W cm − sr − with 2 × − W cm − sr − steps for [SI].Deviations of the OTP from its LTE value are reflected asvertical displacements between the ortho and para transitionsin the EXD, forming a ”saw-tooth” pattern between the twoH species. The observed OTP was estimated by examining thealignment of the S(5) data point with the neighboring S(4) andS(6) transitions, following the method outlined in Wilgenbuset al. (2000); in all cases the spatial OTP variations are small,and values lie very close, within the statistical error limits, to thehigh-temperature LTE value of 3.The H S(3) transition at 9.7 µ m is sensitive to the amount ofdust along the line of sight, as it is located within a wide-band sil-icate absorption feature at the same wavelength; consequently avisual extinction value can be estimated by examining the align-ment of the S(3) point in comparison to the S(2) and S(4), havingpreviously corrected for any deviations of the OTP ratio from theequilibrium value. Visual extinction is then estimated assumingan A . / A V ratio equal to 0.087 (Rieke & Lebofsky 1985). A V is found to exhibit no regular pattern along the outflow and torange between 7 and 9 mag, being consistent with previous es-timates of 7–15 mag by McCaughrean et al. (1994); Caratti oGaratti et al. (2006) .Considering an average value of A V equal to 8 mag acrossthe map, we have dereddened all the detected H lines falling inthe SL range (ie. S(2) to S(7)) using the extinction law of Rieke& Lebofsky (1985). Consequently, for all the points of the SLdata cube where at least 3 H lines are detected above 3-sigma,we have calculated the excitation temperature from the slope ofa least square fit on the dereddened data points in the EXD. An O’Connell et al. (2005) reported a higher extinction of 25 mag in theredshifted lobe but considered it rather tentative given their observingconditions example of these operations is shown in the left panel of Fig. 5in the case of the R1 peak.The results over the whole HH211 region are summarized inFig. 6 in the form of an excitation temperature map. Commonmorphological characteristics between the excitation tempera-ture map and the H S(2)-S(7) emissions are evident; peaks inexcitation temperature tend to coincide with peaks of emission.At those peaks, temperature reaches values up to 1200 K, whilein the more di ff use outflow regions the temperature is as low as700 K. In Table 1 we report the derived excitation temperatureswith this method for the peaks of the S(5) emission along theblue and the red lobes of the outflow.In addition to the ”warm” H component at (cid:39) + R2 region in the LL 10.5 (cid:48)(cid:48) pixel, in-cluding the S(0) and S(1) datapoints. For consistency, the fluxesof S(2) to S(7) lines from the SL module were re-extracted atthe coarser LL module sampling. Due to the extended di ff useemission, the S(0) line flux was obtained after subtracting theambient contribution estimated in adjacent o ff -outflow regions.While the resampled S(3) - S(7) emission keeps the same”warm” temperature ∼ ∼
300 K not probed by the higher J transitions.Another interesting result is that column densities obtained fromthe S(2)-S(7) lines in the LL pixel are smaller than in the SLpixel of 3.5 (cid:48)(cid:48) , indicating that the warm component does not fillthe LL pixel. O. Dionatos et al.: Spitzer spectral line mapping of the HH211 outflow
Fig. 5.
Excitation diagrams of H in various pixel sizes. ( Left ) H columns in peak R1 from transitions in the SL module in 3.5 (cid:48)(cid:48) pixels (open circles). Ortho to para ratio and extinction are estimated from the alignment of the S(5) and S(3) data points with theirneighboring S(2), S(4) and S(4), S(6) respectively. Dereddened columns (filled red circles) are linearly fitted (solid line) to derivethe excitation temperature. Derived values are listed on the upper right corner of the diagram. ( Center, right ) Excitation diagrams ofpeaks R1 + R2 and LL1 in 10.5 (cid:48)(cid:48) pixels, combining the H transitions observed in both the LL and SL modules. The SL line mapswere re-extracted at the LL map pixel scale. A second cooler component is apparent in S(0) – S(2) with T ∼
300 K; the higher H transitions trace a warmer component with temperature T ∼ B1 R2R1B2B3
Fig. 6.
Excitation temperature map of the ”warm” H componentdeduced from the SL module map (3.5 (cid:48)(cid:48) sampling). Temperatureis derived by a linear fit to the dereddened columns from thev = transitions are inblack.A similarly resampled EXD towards the LL1 peak also re-veals a cool component at 300 K, in addition to the warm com-ponent traced by the J > + B3, as part of the corresponding LL pixel falls beyondthe SL map coverage. This double temperature structure di ff ersfrom a similar Spitzer analysis of the L1448 outflow where allv = The total H column density of the ”cool” and ”warm” compo-nents can be measured from the intersection of the linear fit to the data points in the EXD with the N / g axis, multiplied by thepartition function Z ( T ) = (cid:80) g i exp ( − E i / kT ). The later is esti-mated for the given temperatures, summing up for the first 35energy levels of H . Derived column densities for the two tem-perature components are reported in Tables 1 and 2 and rangefrom ∼ for the warm component (averaged over the SLpixel size) to 8 × cm − for the cool component (averagedover the LL pixel size).As already mentioned in the previous section, the inner peakof H emission at B1 / LL1 coincides with peaks of the underlyingCO and SiO jet; it is therefore instructive to compare H columndensities with those deduced from CO, taking into account thedi ff ering beam sizes, to investigate possible physical connectionsbetween the di ff erent tracers.Lee et al. (2007) estimate an H column towards the innerCO jet of 4 − × cm − in their 1 (cid:48)(cid:48) beam, assuming an ex-citation temperature of 100 K and a standard CO / H abundanceratio of 8 . × − (ie. a fully molecular gas). If the H mid-IRemission arises from the same narrow jet, and is relatively uni-form along the jet axis, our derived column densities have to becorrected for the di ff erent beam sizes across the jet by a factor = pixel size / (cid:48)(cid:48) . This increases the H columns to ∼ . × cm − for the ”warm” component, and to ∼ × cm − forthe ”cool” component. The latter value is in excellent agree-ment with the CO-derived one. This coincidence suggests thatthe ”cool” H gas detected by Spitzer towards LL1 could be trac-ing the same material as the CO jet, which would then be indeed mostly molecular . Note also that the temperature of 300 K ofthe cool H matches that inferred from SiO line ratios in the jet(Hirano et al. 2006).In contrast, the warm H component column density cor-rected for beam dilution across the jet remains 10 times smallerthan that of the CO jet. This may indicate that the warm H emis-sion close to the source arises only from a lower-density outerenvelope around the CO jet. Alternatively, the warm H com-ponent could trace thin shock-heated zones within the CO jet,thus filling only a fraction of the 3.5 (cid:48)(cid:48) SL pixel length along thejet and su ff ering further beam dilution. We investigate the shockhypothesis in more detail in the next section ( § components, assuming a uni- . Dionatos et al.: Spitzer spectral line mapping of the HH211 outflow 7 Table 1.
Physical properties of the warm H component from H v = (cid:48)(cid:48) ) scale Position O ff sets H ) ˙M a ( (cid:48)(cid:48) ) (10 − W cm − sr − ) (K) ( 10 cm − ) (10 − M (cid:12) yr − )R1 [-20.1, 13.3] 3.3 1100 ±
70 1.15 1.1R2 [-29.9, 17.0] 2.3 1050 ±
75 1.3 1.2B1 [12.2, -4.7] 1.3 860 ±
40 1.1 1.1B2 [26.9, -10.3] 2.0 950 ±
70 1.5 1.4B3 [35.1, -10.6] 1.9 980 ±
75 1.4 1.4 a ˙M is proportional to ( V jet /
100 km s − ) × (3.5 (cid:48)(cid:48) / l t ) and is not corrected for postshock compression. Table 2.
Physical properties of the cool H component from H v = (cid:48)(cid:48) ) scale Position O ff sets H ) ˙M a ( (cid:48)(cid:48) ) (10 − W cm − sr − ) (K) ( 10 cm − ) (10 − M (cid:12) yr − )R1 + R2 [-19.3, 15.1] 1.4 330 ±
50 7.5 2.0LL1 [6.3, 2.0] 1.3 300 ±
45 7.8 2.8 a ˙M is proportional to ( V jet /
100 km s − ) × (10.5 (cid:48)(cid:48) / l t ) and is not corrected for postshock compression. form laminar flow across the corresponding pixel, using the re-lationship from Dionatos et al. (2009):˙ M = µ m H × (cid:104) N ( H ) A (cid:105) × ( V t / l t ) (1)where µ = m H isthe proton mass, N ( H ) the column density averaged over thearea A of the pixel, l t is the projected emitting length along theflow direction (assumed in Tables 1,2 to be equal to the pixellength), and V t the tangential flow speed.The value of V t is quite uncertain. Taking into account radialvelocities in the range 5 - 20 km s − measured for the near-IR H knots in the study of Salas et al. (2003), and assuming an inclina-tion angle between 5 o and 10 o (see § − , introducing an absolute uncertainty ofa factor 3 either way in the inferred ˙ M values. In the following,we adopt V t =
100 km s − for ease of comparison with earlierwork, bearing in mind that ˙ M values scale proportionally to theassumed velocity. In column 6 of Tables 1 and 2 we report theresulting ˙ M for the peaks of emission along the outflow, for bothtemperature components.For the cool H gas component, the laminar mass-flux is ∼ − × − M (cid:12) yr − . For comparison, the one-sided jetmass-loss rate based on CO(3-2) emission obtained by Lee et al.(2007) assuming the same velocity of 100 km s − and a com-pression factor of 3 is 3 . − × − M (cid:12) yr − . Adjusting theirestimations for Helium, and uncorrecting for compression, thelaminar mass flux of the CO jet rises to 15 − × − M (cid:12) yr − ;this is in excellent agreement with the cool H laminar mass fluxestimated here, again consistent with both tracing the same fullymolecular gas.The warm H laminar mass-flux is 15-30 times lower at ∼ − M (cid:12) yr − × ( V t / − ). Similar mass flux estimationswere obtained from the pure rotational H lines along the outflowof L1448 (Dionatos et al. 2009, 0 . − . × − M (cid:12) yr − ), wheresimilar H columns were measured . Note however that if thewarm H does not trace a laminar flow as assumed, but arises ina single shock within the pixel, then the relevant emitting length l t along the flow will be generally narrower than the SL pixel Note that a typographical error was introduced in Table 2 ofDionatos et al. (2009); the values of N ( H ) are ∼ × cm − at CSand ∼ . − . × cm − for the other positions. size, and the true mass flux will be higher than listed in Table 1.Mass-fluxes in the warm component in the shock hypothesis arepresented in the next section. component As pointed out in § J ≥ ) H emission coin-cide with peaks of temperature (Fig. 6) and possibly correspondto regions of shocked gas. In order to constrain the shock condi-tions in these regions, we employ the existing shock model griddescribed in the work of Kristensen et al. (2007) which is basedon the MHD-VODE multi-fluid steady shock code of Flower &Pineau des Forˆets (2003). The grid includes both continuous (C)and jump (J) type shocks and predicts the H lines brightness forvarious values of pre-shock density ( n iniH ), shock velocity ( V s ),initial OTP ratio ( otp ini ) and transverse magnetic field density( ϕ ). The ranges of parameters investigated here are as follows: – n iniH = to 10 cm − by alternating factors of 5 and 2. – V s =
10 - 50 km s − with a step of 1 km s − – otp ini = – ϕ ( µ G ) = b × (cid:112) n H ( cm − ) where we fix b = b = cosmic ray ionization rate of 5 × − s − , yielding an ionizationfraction ∼ × − ( n iniH / cm s − ) − . . The H level populationsand line brightnesses are integrated along the postshock coolingzone down to a temperature of 50 K, which is reached in a cool-ing time ∼
100 (10 cm − / n H ) yrs in C-shocks (much shorter inJ shocks). For more information on both the model and the griddata, the interested reader is referred to the articles of Flower &Pineau des Forˆets (2003) and Kristensen et al. (2007).In the following, we have employed a χ fitting method tooptimally reproduce the observed H emission selecting the bestmatching C or J-type steady shock models from the grid. The χ - fitting is performed with and without considering the S(2) tran- O. Dionatos et al.: Spitzer spectral line mapping of the HH211 outflow sym. n H v S b ff χ χ
38 0.031 15.7 6.7 510
23 0.03 14.5 6.5 10
19 0.03 12.5 5.7 510
13 0.023 6.6 3.2 10
12 0.02 6.1 2.7 510
10 0.016 5.4 2.5 10
10 0.012 5.2 2.3
Fig. 7.
Steady 1D C-shock model fits to the excitation diagramof the S(2)–S(7) lines at the peak of emission R1 in a 3.5 (cid:48)(cid:48) pixel. Observed points (filled circles with 10% error bars) aredisplayed along with best fitting χ models for pre-shock den-sities between 10 and 10 cm − . The top panel lists the corre-sponding shock velocities and beam filling factors, along withthe reduced χ fit values to the S(2) - S(7) (J >
1) and to the S(3)- S(7) (J >
2) H transitions. A strong degeneracy among modelsis found (see text).sition, since the latter may be dominated by the cool component(see Fig.5). As a general trend we found that steady C-shock models are ableto reproduce the observed S(2)-S(7) brightnesses only if we as-sume a small beam filling factor. For this reason we have per-formed our χ fits on relative line brightnesses, scaled to the H S(5) line; for the best fitting models we have then calculateda beam filling factor as the average ratio of the observed overmodel line brightness.In Figure 7 we plot the EXDs of the best fitting χ models foreach preshock density, scaled by the corresponding beam fillingfactor, along with the observed data (filled circles) for the peakof emission R1; note that beam filling factors decrease for in-creasing values of density. There is a strong degeneracy amongthe models, with a very small spread in column densities for n iniH in the range between 5 10 and 10 cm − and velocities in therange 10-13 km s − , however models for lower values of n H ≤ and velocities up to 30 km s − cannot be excluded; each oneof these models can account fairly well for the observed emis-sions. Figure 8 presents EXD of the best matching C–type shocks(red squares, always scaled by the corresponding beam fillingfactor) along with the observed data points (filled circles), forthe peaks of emission as observed along the outflow. Input pa-rameters for models are reported in Table 3. All C-shock modelspredict moderate shock velocities (10-15 km s − ) and an initialortho to para ratio at the high temperature LTE value of 3.Steady C-shock models reproduce the data only if we con-sider low beam filling factors b f f ∼ A shock = A pixel × b f f , corresponding to a typical sizeof √ A shock ∼ (cid:48)(cid:48) - 0.7 (cid:48)(cid:48) . These dimensions are comparable tothe reported jet width < (cid:48)(cid:48) in the interferometric CO and SiOmaps of Lee et al. (2007). High post-shock densities in the samerange as the best fitting C-shock models (5 10 - 10 cm − ) havebeen retrieved in the jet by SiO emission studies (Nisini et al.2002b; Hirano et al. 2006; Lee et al. 2007), supporting our mod-eling results. The shock models predict thin H post-shock cool-ing zones of 0.015 (cid:48)(cid:48) - 0.25 (cid:48)(cid:48) , much smaller than the SL pixel sizeof 3.5 (cid:48)(cid:48) , explaining the small warm column densities derived inthe previous section. The cooling zone is also smaller than theshock surface size inferred from the filling factor, as required forthe 1-D planar shock approximation to be locally valid.If the warm H emission indeed traces low-velocity internalshocks within a fast jet beam, we may infer the jet mass-fluxfrom the pre-shock density and the shock area as:˙ M j , shock = µ m H × b f f A pixel × n initH × V jet . (2)Considering a pre-shock density of 5 × cm − with a beamfilling factor b f f of 0.023, and a jet speed of 100 km s − , weobtain for the warm H ˙ M j , shock = × − M (cid:12) yr − . This is 7times higher than the value obtained in Section 4.1.2 assuminga laminar flow filling the SL pixel. The mass-flux entering theshock is only a factor 3-4 below the laminar mass-flux for thecool H (and CO). Such a discrepancy is expected in a time-variable flow undergoing internal shocks (Hartigan et al. 1994).Therefore, the ”warm” H component seems compatible withlocalized heating in slow internal C-shocks within the molecularjet beam, and may not trace a separate lower density flow. Notethat pre-shock densities higher than 5 × cm − seem to beexcluded, as the pre-shock mass-flux rate from Equ. 2 wouldbecome comparable to the accretion rate onto the HH211-mmdriving source (8 × − M (cid:12) yr − Lee et al. 2007).However, SiO is observed only within ∼ (cid:48)(cid:48) from the source,therefore pre-shock densities of 5 × cm − might apply only tothe B1 region. In addition, in our maps the H emission remainsunresolved only close to the driving source. Further out at thebow-shocks, the emission appears slightly resolved at the 3.5 (cid:48)(cid:48) scale. A small b f f < .
04 would therefore imply that the bow-shock surface is highly patchy. This is a reasonable possibility,since HST observations in several outflows have shown that whatappears as an individual bow-shock at low resolution is in factresolved in many sub-structures of sub-arcsec scale (e.g. HH2 inOrion, Bally et al. (2002)). In HH211, SMA observations, per-formed with a resolution of ∼ (cid:48)(cid:48) (Lee et al. 2007), show indeedthat the CO emission is highly clumped, both in the high velocityjet and along the low velocity cavity ahead of the bow-shocks.An alternative explanation for the small beam filling factor inthe outer HH211 bow-shocks might be that the pre-shock densitydrops to ≤ cm − , so that the cooling length becomes compa-rable to the 3.5 (cid:48)(cid:48) pixel size; the hypothesis underlining our modelcomparison, that the whole cooling zone fits inside the pixel,would then become invalid. Indeed, O’Connell et al. (2005) de-rived a low pre-shock density in the R1 and R2 regions by fitting . Dionatos et al.: Spitzer spectral line mapping of the HH211 outflow 9 Fig. 8.
Excitation diagrams for all emission peaks mapped with the SL module; observed points (filled circles with error bars), aredisplayed along with best χ fitting C (red squares) and J (blue triangles) shock models. C- shock models can account fairly wellfor the observed emissions considering a small beam filling factor, whereas J-shock models cannot reproduce well the slopes of theobserved data points, except marginally in R1 and R2.a 3D C-type bow-shocks to the near-IR ro-vibrational H lines.However, no predictions for the mid-IR H emissions were pro-vided. Extension to 3D bow-shock geometry of the model ap-plied here (Gustafsson et al. 2010) will be needed to performmore reliable modeling of the outer bowshocks of HH211 thanpossible in the present analysis.A third possibility might be that the preshock gas ahead ofthe leading bowshocks is more ionized than we have assumed,due to the presence of a UV flux from nearby dissociative shockssuch as that observed by Tappe et al. (2008) in the blue-shiftedapex. This would increase the ion-coupling in the C-shock, re-ducing its width and column density, and yielding a larger fillingfactor. Indeed, O’Connell et al. (2005) find that a high pre-shockionization fraction of 10 − and a high atomic hydrogen fractionof 20% provide a good fit to the near-IR lines in the redshiftedbow. However, given the already high degeneracy of our modelfits, and the lack of information on shock geometry and thick-ness in low angular resolution Spitzer maps, introducing the UVfield as an additional free parameter is not feasible in practice.Additional constraints from other line tracers or from higher res-olution mid-IR maps (eg. with JWST) will be required to furtherexplore this alternative. In the case of J-shocks, absolute brightness can be fairly wellmatched by models and beam filling factors < ∼ − × cm − are lower than for the bestfitting C-shock ones. As a general trend though, the best fittingJ-shock models, presented with blue triangles in Fig. 8 and listedin Table 4, cannot reproduce at the same time the slope of the ob-served data points, except for the red-shifted bowshock (R1 andR2).An analogous behavior has been noted in the work ofGiannini et al. (2006), where models produced by the same codeof Flower & Pineau des Forˆets (2003) are employed to studythe conditions of the HH54 outflow. Considering only the ro-tational H emission in reported excitation diagrams, Gianniniet al. (2006) find that steady C-type shocks overpredict the ob-served v = + J shock) optimallyreproduces the observed H emission. The same conclusion wasreached also in the case of L1157 (Gusdorf et al. 2008a). Indeed,the overall e ff ect as reflected in excitation diagrams is to lowerthe column densities of low excitation energy levels (truncatedC-precursor) and increase the column densities of the rovibra-tional energy levels of H , excited in the J-type front (Floweret al. 2003). In our case, such a model might potentially fit ourmid-IR data in the outer bowshocks without the need to con-sider a beam filling factor <
1. A more detailed comparison withunsteady C + J-shock models would be needed to test this pos-sibility, however such models were not included in our adoptedgrid as the age of the shock would have, again, added an extra
Table 3.
Parameters of shock models best fitting the S(3) to S(7) lines
Peak C-shock (b =
1) J-shock (b =
0, b f f = n iniH (cm − ) V s (km s − ) opr ini b f f n iniH (cm − ) V s (km s − ) opr ini R1 5 10 - 10
10 - 13 3 0.01 - 0.02 5 10
14 3R2 5 10 - 10
10 - 11 3 0.02 - 0.04 5 10
14 3B1 5 10 - 10
10 3 0.02 - 0.03 10
10 3B2 5 10 - 10
10 - 11 3 0.01 - 0.03 10
10 3B3 5 10 - 10
10 - 11 3 0.02 - 0.03 10
11 3 free parameter to already degenerate fits. The shock ages wouldhave to be much shorter than the outflow dynamical age of about1000 yrs in order to reduce by a factor 50 the warm H column(typically a shock age of 1 /
50 of the H cooling time, ie. ∼ cm − / n H ) yrs). As pointed out in section §
3, all atomic / ionic emission fallswithin the LL IRS module and consequently these maps are ofinferior spatial resolution compared to the SL ones. As atomicemission is found to be enhanced at the bow shocks but remainssignificant also close to the driving source, we have focused ouranalysis all along the jet axis including the intermediate pointsLL1 and LL2 (see Fig. 4) between the bow shocks, where the[FeII] 26 µ m line is detected. For the analysis of the ionic emis-sion we have followed the techniques employed by Dionatoset al. (2009), and the reader is referred to this article for a de-tailed description. As a first step, the observed emissions of [SiII] and [FeII] are em-ployed in order to constrain the electron density. In Figure 9 theratio of [SiII]34.8 µ m over [FeII]26.0 µ m as a function of the elec-tron density is presented, for temperatures of 500K, 1000K and4000K (solid and dashed lines). A statistical equilibrium modelwas employed that uses the first 16 levels for [FeII] (Nisini et al.2002a), and a two level system was considered for [SiII], withradiative and collisional rates with electrons taken from Dufton& Kingston (1991); a solar abundance ratio of Si + / Fe + was con-sidered, taking into account the similar ionization potential ofthe two species as well as the similar e ffi ciency of sputtering ofSi and Fe atoms from olivine grains in C-shocks, according toFigure 2 of Gusdorf et al. (2008b). In the diagram we plot theobserved ratio limits at the bow shocks for clarity, but the re-sults for the other points along the jet are similar and are listedin Table 4. Lower limits for the electron density vary between 40and 150 cm − for temperatures of 4000K and 500K respectively,while upper limits converge for the considered temperatures to ∼
400 cm − . Temperature limits lower than a few thousand K arereported in the work of Tappe et al. (2008) for the blue-shiftedbow shock considering a non-LTE analysis of the three [FeII]lines (18 µ m, 26 µ m and 35 µ m) that are detected with the highresolution IRS modules. Based on this, we restrict the electrondensity lower limit to 100 cm − by assuming a maximum tem-perature of 3000 K. Fig. 9.
Diagnostic diagram of the [SiII]35 µ m / [FeII26] µ m ratioversus electron density, for temperatures of 500K, 1000K and4000K (solid and dashed lines) and assuming a solar Si + / Fe + abundance ratio. The hatched area represent the observed ratiofor the north and south bow-shock regions, from which electrondensities between ∼
100 and 400 cm − are found . We investigate the depletion onto grains following the analysisof Dionatos et al. (2009). We calculate the gas phase abundanceof refractory elements such as iron and silicon by comparingtheir emissions against the emission of a non-refractory species,in the present case sulphur. Theoretical emissivities for sulphurwere calculated considering a five level statistical equilibriumcode for an upper temperature limit of 3000 K assumed fromthe work of (Tappe et al. 2008) and a lower limit of 700 K asimposed from the H analysis. For the excitation of sulphur, wehave also calculated emissivities due to collisions with atomichydrogen, taking into account that in environments with low ion-ization fraction and high total density such excitation may proveto be significant, despite the fact that collisional de-excitationrates for atomic hydrogen are lower than the electronic ones(Hollenbach & McKee 1989, γ H /γ e ∼ .
02) for the [SI]25 µ mline. For these calculations we have considered a medium withn(H) = cm − . Excitation of [FeII] and [SiII] by H-collisionswas not taken into account because neutral hydrogen collisional . Dionatos et al.: Spitzer spectral line mapping of the HH211 outflow 11 Table 4.
Physical properties of the atomic / ionic line component Position O ff sets [SI] [FeII] / [SI] [SiII] / [FeII] ˙M([SI]) a , b n b , cH ( (cid:48)(cid:48) ) (10 − Watt cm − sr − ) (10 − M (cid:12) yr − ) (10 cm − )R1 + R2 [-19.3, 15.1] 5.34 ± ± ± ± ± ± ± ± ± + B3 [31.9, -11.0] 16.1 ± ± ± a ˙M is proportional to ( V jet /
100 km s − ) × (10.5 (cid:48)(cid:48) / l t ) and is not corrected for postshock compression. b values calculated for collisions with electrons at T = =
700 K (upper value); including collisions with hydrogen forn(H) = cm − would decrease the minimum and maximum values by a factor of ∼ c average proton density assuming an emitting volume of 10.5 (cid:48)(cid:48) × (cid:48)(cid:48) × (cid:48)(cid:48) Fig. 10.
Diagnostic diagram of [FeII]26 µ m / [SI]25 µ m ratio ver-sus temperature assuming solar relative abundances. Electrondensities of 100 and 400 cm − are considered (full and dot-ted lines, respectively). The e ff ect of collisions with atomic hy-drogen on the [SI] excitation is also shown (long-dashed line, n ( H ) = cm − ). The observed ratios are always below thepredicted values for solar abundances, which is interpreted asiron being depleted onto dust-grains in the considered regions.rates are much less important for ions (i.e. ∼ -10 weakerthan electronic collisional rates), and are not available for all thelevels considered in our model.Results of these calculations are illustrated in Fig. 10 wherewe plot the ratio of [FeII] 26 µ m over [SI] 25 µ m against gas tem-perature for electron densities of 100 and 400 cm − , assumingthat iron is fully ionised (which is the case in most dissociativeshocks) and has a solar abundance. The observed ratios for thepoints R1 + R2 and B1 + B2 are way below the expected values,giving rise to the hypothesis that iron is heavily depleted ontodust grains in the considered regions, more so if the emitting gasis cooler. The case including atomic hydrogen collisions withn(H) = cm − is separately plotted in Fig. 10. The iron de-pletion is less severe in this case but is still substantial, about afactor 10. Table 5.
Gas phase abundance of Fe + and Si + from line ratios tosulfur Position: R1 + R2 LL1 LL2 B2 + B3[Fe + gas ] / [Fe (cid:12) ] a − − − − [Fe + gas ] / [Fe (cid:12) ] b − − − − [Si + gas ] / [Si (cid:12) ] a − − − − [Si + gas ] / [Si (cid:12) ] b − − − − a assuming only collisions with electrons ( n e = −
400 cm − ) and T e =
700 K (lower value) and 3000 K (upper value) b including excitation of sulfur by atomic hydrogen with n ( H ) = cm − In Table 5 we give the inferred gas phase depletion for tem-peratures of 700 and 3000 K, both without and with atomic hy-drogen excitation of [SI]. The gas phase abundance of Fe + isabout 2–20% of the solar value for T = + , but since we assumeda solar Fe + / Si + ratio in the first place to derive our n e values, thisis just a self-consistency check. Interestingly, the SiO abundancein the HH211 molecular jet, estimated as ∼ − by (Nisini et al.2002b; Chandler & Richer 2001) corresponds to 3% of the solarabundance of silicon, therefore the fraction of silicon releasedin the gas phase seems comparable in the atomic and molecularcomponents of HH211. The total number of H nuclei in the emitting volume V can bederived from the line luminosity according to the relation: n H V = L ( line ) (cid:18) h ν A i f i (cid:20) XH (cid:21)(cid:19) − , (3)where A i and f i are the spontaneous radiative decay rate and frac-tional population of the upper level, and [X / H] is the gas phaseabundance of the atom / ion under consideration with respect to Hnuclei. Line luminosities were integrated over the LL pixel sizearea, and the fraction of atoms at the upper level f i was calcu-lated for n e between 100 and 400 cm − and T = Values of n H at individual positions are listed in column 7 ofTable 4 assuming an emitting volume in the LL pixel of 10.5 (cid:48)(cid:48) × (cid:48)(cid:48) × (cid:48)(cid:48) (uniform narrow jet). The inferred n H depends stronglyon the adopted excitation conditions: it is about 0 . − . × cm − for T = n e =
100 cm − , and 8 times higher forT =
700 K, n e =
400 cm − . The ionisation fraction is then 1 . − × − at 3000 K, or twice smaller at 700 K. Such an ionizationlevel seems to favor the upper temperature, and therefore thelower density range. If the fraction of H atoms is important with n ( H ) (cid:39) cm − , the n H values would be further reduced bya factor 2, down to a few 10 cm − . In the case again, that theemission comes from a smaller volume, eg. a small 1 (cid:48)(cid:48) knot asobserved in the near-IR [FeII]1.64 µ m line, the density would riseup to 10 cm − .The mass flux in the atomic / ionic component assuming alaminar flow can also be derived applying the relationship givenin Nisini et al. (2005) as applied to the mid-IR lines in Dionatoset al. (2009):˙ M = µ m H × ( n H V ) × ( V t / l t ) , (4)where µ = m H is the protonmass, n H V is the total number of protons in the emitting regionas given by Eq. 3, l t is the projected emitting length along theflow and V t the tangential flow velocity.Derived mass flux values at the various emission peaks arelisted in column 6 of Table 4 for a tangential speed of 100 kms − and l t = . (cid:48)(cid:48) (LL pixel size). Like n H V , they again dependstrongly on the adopted excitation conditions. From the ioniza-tion fraction considerations made above, we favor the high tem-perature, smaller ˙ M values of ∼ . − × − M (cid:12) yr − . Theseare 20–50 times smaller than the cool H jet mass-flux, assumingthe same flow speed of 100 km s − .However, we stress again that ˙ M estimates assuming a lam-inar jet flow may be in error if the [SI] emission arises fromshocks, which is likely given the ionization fraction inferredabove. This is di ffi cult to quantify without appropriate shockmodels.The atomic / ionic emission appears to require higher excita-tion shocks than those producing the mid-IR H lines. The lowvelocity 10–15 km s − shock models that best fit the warm H emission are unable to reproduce also the observed [FeII] and[SiII] line intensities unless the shocks are of J-type and theseatoms are essentially undepleted. As discussed in Section 4.3.2there is overwhelming evidence of depletion of these refractoryspecies, so [FeII] and [SiII] probably arise from faster shocksthan those dominating the mid-IR H lines. It is indeed likelythat the mid-IR atomic lines originates from the same dissocia-tive shocks at v s ≥
30 km s − that give rise to the optical and NIRemission of H α , [SII] and [FeII] (O’Connell et al. 2005; Carattio Garatti et al. 2006; Walawender et al. 2006, 2005).Unfortunately, our shock models do not yet include ioniza-tion and dissociation by the shock UV flux, so we cannot explorethe relevant excitation range. However, we note that the J–shockmodels of Hollenbach & McKee (1989), which include UV flux,suggest that the [SI]25 µ m brightness is relatively independent ofshock speed over the range 30–100 km s − and is roughly pro-portional to n iniH . Our measured flux inside the LL 10.5 (cid:48)(cid:48) pixelwould suggest n iniH × b f f (cid:39) − − . From Eq. 2 wewould infer a preshock mass-flux of ∼ . − × − M (cid:12) yr − ,still smaller by a factor 4–40 than the cool H jet mass-flux.These results would suggest that the atomic component inthe HH211 jet does not have enough momentum flux to entrainthe CO / SiO / H jet, if their velocities are comparable and the shock compression in the cool H / CO component does not ex-ceed a factor 3. The molecular jet would then have to trace mate-rial ejected from the accretion disk, while the atomic componentwould trace a separate ejection, e.g. from hotter more internalregions of the accretion disk. On the other hand, a higher com-pression factor in the cool H / CO emitting zone, or a faster ionicjet would su ffi ce to remove the discrepancy, therefore a definiteconclusion on whether the molecular jet is ejected or entrainedcannot be reached from the present data alone.
5. Conclusions
We have carried out Spitzer spectral mapping observationstowards the jet driven by the Class 0 source HH211-mm.Molecular lines (pure rotational H ) as well as fundamentalatomic and ionic lines ([SI], [SiII], [FeII]) were detected, andtheir maps follow the characteristic bipolar outflow pattern astraced by near-IR H and CO lines. H emission becomes im-portant only 5 (cid:48)(cid:48) away from the driving source while atomic andionic lines are detected very close to the driving source. In theinner part of the blue-shifted lobe, the H emission is spatiallycoincident with the high velocity jet observed in CO and SiO.Analysis of the observed H lines reveal two temperaturecomponents: ”cool” gas at T ∼ ∼ towards the CO jet is compatible with hydro-gen being mostly molecular. The warmer component traced bythe S(2) to S(7) lines is well fitted by C-shock models of highdensity (cid:39) × and a small shock cross section of 0.5 (cid:48)(cid:48) , com-patible with the density and width of the CO jet quoted by Leeet al. (2007). The warm H emission could then trace thin layersof warm post-shock gas within the CO jet.Similarly, high-density C-shock models can also account forthe brightness of the H mid-IR emission further downwind be-yond the CO and SiO jet, but the small shock surface filling fac-tor of 0.01-0.04 is not easy to reconcile with the spatial exten-sion visible in Spitzer maps, unless the bow-shock wings arevery clumpy. Lower density shocks of n H ≤ cm − where thecooling zone spreads over several pixels and / or higher preshockionization may need to be considered in these regions, as previ-ously invoked to model near-IR H excitation in the red-shiftedbow-shock of HH211 (O’Connell et al. 2005).The detected fine structure lines mapped very close to thedriving source signify the presence of an embedded atomic jet.Line ratio diagnostics indicate a gas-phase depletion of iron andsilicon of at least a factor 10, and lower excitation conditionsthan in optically visible jets. An atomic jet of similar proper-ties has been also detected by Spitzer in the outflow of L1448-C(Dionatos et al. 2009). As suggested in that case, the detectedatomic gas in the HH211-mm outflow may represent the equiv-alent for the Forbidden Emission Line (FEL) region observed inmore evolved ClassI / II sources. The gas depletion of iron and sil-icon indicates that dust grains have survived in the atomic flow,ruling out an origin from within the dust sublimation zone closeto the protostar. The excitation conditions require faster shocksthan those producing the H mid-IR lines.Estimations of the molecular and atomic mass-flux rateshave been performed using both a laminar flow assumption andshock models. The cool H mass-flux is comparable to that in-ferred from CO observations, with a value uncorrected for com-pression of ∼ × − M (cid:12) yr − ( V /
100 km s − ); a similar value isfound for the warm H if it arises from dense C-shocks as sug-gested by best fit models. The atomic component mass-flux is . Dionatos et al.: Spitzer spectral line mapping of the HH211 outflow 13 uncertain by a factor 8 due to its uncertain temperature, but con-siderations of ionization fraction as well as published dissocia-tive shock models suggest a mass-flux smaller than in the molec-ular jet. However, the momentum fluxes could become compa-rable if the atomic jet is faster than the molecular jet, or if thecool H and CO su ff ered shock compression. Given the uncer-tainties, it is still unclear whether the molecular CO jet is tracingambient gas entrained by the atomic jet, or if it has to be ejected,for example in a molecular MHD disk wind as recently modeledby Panoglou et al. (2010, A&A, submitted). Acknowledgements.
This work is based on archival data obtained with theSpitzer Space Telescope, which is operated by the Jet Propulsion Laboratory,California Institute of Technology under a contract with NASA. It was supportedin part by the European Community’s Marie Curie Actions - Human Resourceand Mobility within the JETSET (Jet Simulations, Experiments and Theory) net-work under contract MRTN-CT-2004 05592. Financial contribution from con-tract ASI I / / / References
Arce, H. G., Shepherd, D., Gueth, F., et al. 2007, in Protostars and Planets V, ed.B. Reipurth, D. Jewitt, & K. Keil, 245–260Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Astronomical Society ofthe Pacific Conference Series, Vol. 336, Cosmic Abundances as Records ofStellar Evolution and Nucleosynthesis, ed. T. G. Barnes, III & F. N. Bash,25– + Bally, J., Heathcote, S., Reipurth, B., et al. 2002, AJ, 123, 2627Caratti o Garatti, A., Giannini, T., Nisini, B., & Lorenzetti, D. 2006, A&A, 449,1077Chandler, C. J. & Richer, J. S. 2001, ApJ, 555, 139Davis, C. J. & Smith, M. D. 1996, A&A, 309, 929Dionatos, O., Nisini, B., Garcia Lopez, R., et al. 2009, ApJ, 692, 1Dufton, P. L. & Kingston, A. E. 1991, MNRAS, 248, 827Eisl¨o ff el, J., Froebrich, D., Stanke, T., & McCaughrean, M. J. 2003a, ApJ, 595,259Eisl¨o ff el, J., Froebrich, D., Stanke, T., & McCaughrean, M. J. 2003b, ApJ, 595,259Enoch, M. L., Young, K. E., Glenn, J., et al. 2006, ApJ, 638, 293Ferreira, J., Dougados, C., & Cabrit, S. 2006, A&A, 453, 785Flower, D. R., Le Bourlot, J., Pineau des Forˆets, G., & Cabrit, S. 2003, MNRAS,341, 70Flower, D. R. & Pineau des Forˆets, G. 2003, MNRAS, 343, 390Froebrich, D. 2005, ApJS, 156, 169Giannini, T., McCoey, C., Nisini, B., et al. 2006, A&A, 459, 821Gibb, A. G., Richer, J. S., Chandler, C. J., & Davis, C. J. 2004, ApJ, 603, 198Gueth, F. & Guilloteau, S. 1999, A&A, 343, 571Guilloteau, S., Bachiller, R., Fuente, A., & Lucas, R. 1992, A&A, 265, L49Gusdorf, A., Pineau Des Forˆets, G., Cabrit, S., & Flower, D. R. 2008a, A&A,490, 695Gusdorf, A., Pineau Des Forˆets, G., Cabrit, S., & Flower, D. R. 2008b, A&A,490, 695Gustafsson, M., Ravkilde, T., Kristensen, L. E., et al. 2010, A&A, 513, A5 + Hartigan, P., Morse, J. A., & Raymond, J. 1994, ApJ, 436, 125Herbig, G. H. 1998, ApJ, 497, 736Hern´andez, J., Calvet, N., Hartmann, L., et al. 2005, AJ, 129, 856Hirano, N., Liu, S.-Y., Shang, H., et al. 2006, ApJ, 636, L141Hollenbach, D. & McKee, C. F. 1989, ApJ, 342, 306Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 18Kaufman, M. J., Wolfire, M. G., Hollenbach, D. J., & Luhman, M. L. 1999, ApJ,527, 795Kristensen, L. E., Ravkilde, T. L., Field, D., Lemaire, J. L., & Pineau Des Forˆets,G. 2007, A&A, 469, 561Lada, C. J., Muench, A. A., Luhman, K. L., et al. 2006, AJ, 131, 1574Lee, C.-F., Ho, P. T. P., Palau, A., et al. 2007, ApJ, 670, 1188McCaughrean, M. J., Rayner, J. T., & Zinnecker, H. 1994, ApJ, 436, L189Nisini, B., Bacciotti, F., Giannini, T., et al. 2005, A&A, 441, 159Nisini, B., Caratti o Garatti, A., Giannini, T., & Lorenzetti, D. 2002a, A&A, 393,1035Nisini, B., Codella, C., Giannini, T., & Richer, J. S. 2002b, A&A, 395, L25O’Connell, B., Smith, M. D., Froebrich, D., Davis, C. J., & Eisl¨o ffff