A&A 466, 565-577 (2007)
DOI: 10.1051/0004-6361:20066734
M. D. Smith1,2 - B. O'Connell1,3 - C. J. Davis4
1 - Armagh Observatory, College Hill, Armagh BT61 9DG, Northern Ireland, UK
2 - Centre for Astrophysics & Planetary Science, The University of Kent, Canterbury CT2 7NH, UK
3 - Physics Department, Trinity College Dublin, College Green, Dublin 2, Ireland
4 - Joint Astronomy Centre, 660 N.A'ohoku Place, University Park, Hilo, Hawaii 96720, USA
Received 10 November 2006 / Accepted 30 January 2007
Abstract
The near-infrared twin jets emanating from the HH 212-mm protostar are remarkable
for their symmetry. By performing integral field spectroscopy with the UIST imager-spectrometer
on UKIRT, we investigate the chains of bright knots and arcs within the jets to
gain insight into the underlying physics and dynamics. We obtain numerous images
associated with line emission from vibrationally-excited molecular hydrogen and
the [Fe II] line at 1.64 m. This allows us to study the spatial variation in
excitation of the inner knots and outer bow-shaped objects.
We find that the excitation properties are consistent with outward-moving bow shocks
close to the plane of the sky. However, there is a gradient in excitation transverse to the jet axis
across the inner knots on the scale of 0.1 arcsec. This C-shaped inner symmetry
suggests a transverse source motion rather than precession, possibly related to the jet
bending and the transverse gradient in radial velocity. Moreover, the bow models predict that the
iron emission should peak further ahead of the molecular emission than
actually observed. This leads us to propose that each inner knot consists of two closely-spaced asymmetric bows, as found for the outer bows which clearly occur in distinct
pairs, well-separated in a lower density environment. The weak inter-knot emission
may then be generated within oblique shock waves resulting from
the deflection of fluid across asymmetric bow flanks.
Key words: ISM: jets and outflows - circumstellar matter - infrared: ISM - Herbig-Haro objects - stars: formation
The near-infrared outflow HH 212 is powered by a protostar, HH 212-mm, concealed within a cool
condensation, IRAS 05413-0104. The K-band emission is dominated by lines of molecular hydrogen with the 1-0 S(1) transition at 2.12 m producing spectacular images (Zinnecker et al. 1998).
The highly collimated outflow contains two prominent inner knots which appeared point-like in the initial low resolution observations, at first suggesting an interpretation as a young binary system
(Reipurth 1989; Zinnecker et al. 1992). However, deeper observations revealed that the two bright points are the first pair in a series of knots and bow shocks which delineate the outflow
(Zinnecker et al. 1998,1996). Here, we present images of the bright features taken in a range of near-infrared lines and consider how these structures are formed. The entire H2 outflow and the features we study are displayed in Fig. 1.
![]() |
Figure 1:
The entire HH 212 outflow at 2.122 ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
IRAS 05413-0104 is located in the L1630 molecular cloud about 90
north-east of the Horsehead nebula at a distance of
460 pc. The protostar has been detected at 25, 60 and 100
m but is too cool and embedded to be detected as yet at 12
m (Beichman et al. 1986). It is a low-mass protostar with a bolometric luminosity estimated to be about 8
and a bolometric temperature of
56 K (Froebrich et al. 2003). The
ratio of submillimetre to bolometric luminosity is
about 0.009, which bestows the Class 0 classification (e.g. Froebrich 2005).
The outflow itself covers an extent of 240
,
or 0.54 pc. The total H2 1-0 S(1) luminosity is 6.9
10
(not corrected for extinction)
(Stanke 2000). A series of inner knots, which are
regularly spaced with inter-knot distances of
or
2000 AU, extend on both sides of the bipolar H2 outflow. Emission is also detected in the inter-knot spaces: the knots appear to be connected by a thin stream of emission
(McCaughrean et al. 2002; Zinnecker et al. 1998). Both series of knots
terminate in small bow shocks at a distance of about 13 000 AU from
the central position (NB1 & SB1). These bow shocks are composed of two separate
arcs nested together. Whereas the brightness of the knots decreases
with distance from the source (McCaughrean et al. 2002; Zinnecker et al. 1998), the bow shocks break the trend and
appear relatively bright again. This could partly be due to a reduction in the extinction as the outflow breaks free of its maternal
envelope. Other significant features are as follows.
The Integral Field Unit (IFU) described below delivers near infrared
images at various wavelengths simultaneously. Consequently,
changes in atmospheric transmission or seeing have no influence on
relative knot morphologies or line ratios, and accurate
subtraction of the continuum is possible. Moreover, the relative positions of
the emission peaks for various transition lines of H2 and
[Fe II] (at 1.644 m) can be determined without the
need for registration using field stars. These observations and the
data reduction techniques are described in Sect. 2 and
the results presented in Sect. 3. In Sect. 4,
the results are analysed and a bow shock model is presented.
The United Kingdom Infra-Red Telescope (UKIRT) 1.5 m Imager
Spectrometer UIST is designed to switch quickly between imaging
and spectroscopy modes. It is equipped with a 1024
1024 InSb array. In spectroscopy mode, the internal optics yield 0.12
per pixel. UIST is capable of spectroscopy over a two dimensional field of view using the Integral Field Unit (IFU). The IFU covers a 3.3
6.0
(rotatable) sky area with
a plate scale of 0.24
0.12
.
An image-slicing mirror divides the incident image into 18 adjacent 0.24
6.0
"slitlets'', of which 14 are
usable. The emission through each slitlet is directed onto a dispersing grism and fed through the UIST optical system as if it were a single long-slit, the 14 spectral images from the useable slices
being projected onto the array simultaneously. The resulting spectral
image thus contains the spatial and spectral information from each
slitlet, and can be reformatted to construct
3.3
6.0
images of the observed object,
one image for each resolution element of the grism. The data reduction
process yields a 3D data cube which consists of a stack of 2D images
spanning the full 1024-pixel wavelength range available from the grism.
The inner knots (NK 1 and SK 1) and bow shocks (NB 1/NB 2 and
SB 1/SB 2) of HH 212 were observed with the IFU on November 19 and 20, 2002 and November 22 and December 6, 2003, respectively. The 2002
observations were acquired during instrument commissioning; the 2003
data were obtained in service mode. On each occasion, the long
(6.0
) axis of the IFU was aligned with the jet axis, using a position angle of 22
E of N. The HK grism was used, which gives a wavelength range of 1.4-2.5
m and wavelength
increment of 1.07
10
m. The exposure time per
frame was 120 s with 10 object frames and 10 sky frames taken
for each knot and bow shock. Therefore, the total on-source exposure
time was 1200 s for each object.
Halogen flat field and Argon arc exposures were used to calibrate the data; the flats were used to measure the pixel-to-pixel response across the array and to map the extent and location of the spectrum from each slice on the array. The arc exposures were used to wavelength calibrate each spectrum. The bright stars HD 38529 (G4V, K = 4.41) and HIP 29487 (A2V, K = 6.60) were then used to correct for telluric absorption and to flux calibrate the 2002 and 2003 data, respectively.
Data reduction and flux calibration were performed using the orac-DR
pipeline software specifically written for the UIST IFU,
supplemented by standard STARLINK routines. The data reduction
process leading from raw frames to a flux calibrated 3D data cube is
described in detail on the UKIRT IFU
webpage. The
seeing throughout the observations was
0.5
.
![]() |
Figure 2:
Spectra of the NK 1 knot ( upper) & the SK 1 knot ( lower panel) between 1.5 ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 3:
Extracted line emission images for NK 1.
The peak flux positions were estimated by Gaussian fitting of the image peaks and are indicated
by the crosses. X-axis offsets increase toward the west; Y-axis offsets increase roughly northward along the jet axis (see Fig. 1). A continuum emission image, from data between 1.7 ![]() ![]() |
Open with DEXTER |
Spectra and images were extracted from the data cubes containing the NK 1, SK 1, NB 1/NB 2 and SB 1/SB 2 information (Figs. 2-7). Note that the corresponding Tables referred to below are presented as online materials only.
In the NK 1 and SK 1 spectra (Fig. 2)
the continuum flux is close to zero across the entire
spectral range and the telluric features have largely
disappeared, indicating that the atmospheric absorption correction and
flux calibration were relatively accurate. Nevertheless,
transmission is poor between 1.8 and
2.0
m and
beyond
2.4
m. In these regions, line flux calibrations
are therefore still likely to be inaccurate. In order to estimate the
error in the flux measurements, the same sky region was used to
measure the aperture integrated flux across a broad range of
wavelengths. The standard deviation of these values was then used as
the average 1
error.
The extracted images of NK 1 and SK 1 are displayed in Figs. 3 and 4. These show a compact, knot-like structure. However, weak emission also extends to the north and south of both knots. This inter-knot emission can be seen in deep VLT H2 1-0 S(1) images (McCaughrean et al. 2002). The origin of this emission is mysterious. Note that no continuum emission is detected in either knot or stream and no continuum subtraction was performed.
![]() |
Figure 4: Extracted line emission images for SK 1. See Fig. 3 for details. |
Open with DEXTER |
The photometric results for the inner bows NB 1/NB 2 and SB 1/SB 2 are presented in Tables 2 and 3. The features NB 1 and NB 2 are distinguishable in our extracted images (Fig. 6) so their flux measurements are listed separately. The features SB 1 and SB 2 are not separated, however, so their flux measurements are presented together. The fluxes of the prominent emission lines in our extracted spectra are considerably weaker than for the inner knots, the emission being distributed over a larger area. Spectra of the bows are presented in Fig. 5.
![]() |
Figure 5:
Spectra of the NB 1 bow ( top panel) and SB 1 and SB 2 ( lower panel) between 1.4 and 2.5 ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
Extracted line emission images of NB 1 and NB 2.
NB 2 is the northern, upper-most feature.
Where possible, the positions of the peaks were determined via Gaussian fitting and are indicated by crosses (see Fig. 3 for further details).
A continuum emission image, from data between 1.7 ![]() ![]() |
Open with DEXTER |
![]() |
Figure 7: Extracted line emission images of SB 1/SB 2 (see Fig. 3 for details). These two faint knots, labelled in Fig. 1, are essentially unresolved here. The large faint cross in each panel indicates the peak of the 1-0 S(1) line. |
Open with DEXTER |
The 1-0 S(1) S/N ratios for NB 1 and SB 1 are 23 and
14, respectively (determined from the peak flux measurements). Because of the weaker fluxes, images from the higher
excitation emission lines were difficult or impossible to extract.
In our H2 images (Figs. 6 and 7) the NB 2 bow structure is unresolved whereas NB 1
displays an eastern wing in the 1-0 lines. The SB 1/SB 2 pair
displays wing emission typical of bow shocks with the leading bow SB2
responsible for the emission south of the peak.
Note that [Fe II] emission is not detected in either bow. This is not strange since the expected [Fe II] flux level from a bow shock is usually much lower than the level of the 1-0 S(1) flux, as also found specifically for HH 212 by Caratti o Garatti et al. (2006). As with the inner knots, no continuum emission is detected.
![]() |
Figure 8: Column density ratio (CDR) diagrams for NK 1 ( left panels) and SK 1 ( right panels). The upper panels display data assuming no extinction whereas in the lower panels the extinction has been adjusted to minimise the difference between the 1-0 S-branch and Q-branch lines originating from the same upper energy level. H2 v = 1-0 transitions are represented by squares, 2-1 transitions by crosses and 3-2 by triangles. The faint squares represent the 1-0 Q-branch measurements. The CDRs predicted by the best fit J-type bow shock model are also shown; the solid line represents the first vibrational level, dotted represents the second, and the dot-dashed line shows the third vibrational level. |
Open with DEXTER |
The molecular excitation of the inner knots are displayed in the Column Density Ratio diagrams of Fig. 8. We first derived the required column of H2 in the upper energy level, kTj, of each observed line. However, to permit the relevant error bars to be seen, we then divide the columns by the columns predicted by a uniform medium of temperature 2000 K, taking local thermodynamic equilibrium with an ortho-para ratio of 3. Hence, a horizontal line on a CDR diagram would indicate a temperature of 2000 K. The effect of added extinction (lower panels) is to provide an improved overall alignment, supporting the interpretation of collisional excitation.
The figures indicate that the measured fluxes for the inner knots and bows of HH 212 are consistent with collisional excitation in shock heated gas. Ortho-para ratios determined from the 1-0 S sequence are close to the thermodynamic equilibrium value of 3 for the inner knots. However, the 1-0 S(0) values for the bows are a factor of 2 below the expected value. This is a consequence of poor sky subtraction near the corresponding wavelength, as evident from the spectra (Fig. 5).
The vibrational excitation is here measured by the integrated flux
ratio Rv = (2-1)/(1-0) S(1) and yields the vibrational
excitation temperature 5594 K
(Eislöffel et al. 2000). The highest excitation is found in the
inner knots, 2000 AU from the source where R = 0.10
0.02 for
NK 1 and 0.16
0.02 for SK 1 which correspond to excitation
temperatures of 2167
170 K and 2586
170 K,
respectively.
Lower excitation conditions are found further downwind in each jet, at
a distance of 15 000 AU from the source where the NB 1/2
and SB 1/2 bow shocks are located. The excitation ratio measured for
NB 1 is 0.05
0.04 (
= 1698
326 K). For SB 1 a ratio of 0.06
0.03 (
= 1807
277 K) is found. These results are not in complete agreement with the excitation temperatures derived from deep K-band spectroscopy at the NTT by
Tedds et al. (2002). They found high excitation temperatures
for both the inner knots as well as the bow shocks (
2300 K). However, when we consider only the peak flux in each case,
i.e. by extracting spectra just from the emission peak in our data
cubes, rather than from the entire knot, equally high excitation
temperatures are found (although the accompanying errors are
relatively large) for the inner knots and bows. This suggests that the
IFU observations are more sensitive to the lower excitation extended
emission which extends away from the bows. Note that the spectra presented by
Tedds et al. (2002) were obtained with a long-slit, which
did not pass through the lateral wings of each knot/bow shock.
For SK 1, the (CDR) is not linearly related to Tj due to the
high 3-2 S(3) column. Therefore, the excited gas is not of constant
temperature. A curved surface of varying temperature is consistent
with this result, i.e. as one would expect from excitation in a bow shock,
rather than a planar shock. SK 1 does indeed have a "triangular'',
bow-shaped morphology. For NK 1, however, a linear fit is
clearly possible in Fig. 8. Deeper observations of emission lines originating
from the higher vibrational levels are required to better constrain the models.
The K-band extinction for each knot and bow is determined using the
ratio RK = (1-0 Q(3))/(1-0 S(1)) as described in
O'Connell et al. (2005). The differential extinction between the
two lines is then given by
.
From this, we
estimate AK values of 1.2
0.2 mag and 1.4
0.2 mag
for NK 1 and SK 1. Much lower extinctions of AK < 0.5 mag are
found for NB 1 and SB 1, supporting the idea that they have broken
out of the dense protostellar envelope and are now propagating through
a medium of reduced density.
CDR diagrams for NB 1, NB 2 and SB 1/SB 2 have also been obtained but are not presented. The lack of detections of high excitation transitions makes it impossible to identify a temperature dependence. The CDRs of all the Q-branch lines relative to the 1-0 S-Branch lines constrain the extinction. An AK extinction close to zero produces the best results, confirming that the foreground column density is low.
The low excitation is consistent with the lack of detectable optical emission from the outer bows. This requires these bows to be weak shocks in which the molecules remain intact even at the apex of the bow, as would appear to be the case from the images (McCaughrean et al. 2002; Zinnecker et al. 1998).
An advantage of the IFU is that multi-wavelength information is gathered simultaneously. Positional variations in the molecular excitation can therefore be accurately determined without the misalignment uncertainty produced by image co-registration. The left panels of Fig. 9 display the NK 1 and SK 1 1-0 S(1) images overlaid with [Fe II] flux contours. Whereas the shift in the peak position is clearly detected between the 1-0 S(1) and [Fe II], these images show that the rest of the emission is also shifted; the bulk of the [Fe II] emission lies ahead of the 1-0 S(1) emission (in the jet flow direction), as would be expected if they were forward moving bow shocks with high-excitation tips.
Also shown are maps of the (2-1)/(1-0) S(1) ratio with overplotted contours of the 1-0 S(1) (middle panels) and [Fe II] (right panels) emission. For NK 1, no variation can be detected in the excitation map but the SK 1 ratio map shows a clear increase in the vibrational excitation towards the south with the [Fe II] shifted in the same direction. Therefore, SK 1 closely corresponds to the general structure expected from a bow shock with excitation increasing towards the bow apex.
The peak flux positions are given as offsets relative to the peak of
the H2 1-0 S(1) line in Table 1. These were
determined by fitting a 2D Gaussian
using a STARLINK object detection routine. This method thus
determines the peak under the condition that the intrinsic spatial profiles are
Gaussian. The detection threshold was set in order to include only the
higher value pixels (typically about 15-25) which are near the peak
position. The pixel scale is 0.24
in the x-direction and 0.12
in the y-direction. However, the errors in the peak locations are considerably less than the pixel size and are related to the seeing. It is straightforward to show that the theoretical error
on the peak position through the Gaussian fitting is approximately
given by
where
is the standard deviation of
the Gaussian and S/N, the signal to noise ratio, is the ratio of
peak flux to the noise level for each line transition. Therefore, the
high S/N ratios for the 1-0 S(1) lines of NK 1 and SK 1 of
49 and
110, respectively, lead to small resultant
errors, as indicated in Figs. 10 and 11.
![]() |
Figure 9: Left: H2 1-0 S(1) image of NK 1 ( upper panels) and SK 1 ( lower panels) overlaid with [Fe II] contours. Middle and right: H2 (2-1)/(1-0) S(1) ratio image of NK 1 ( upper) and SK 1 ( lower) overlaid with H2 1-0 S(1) and [Fe II] contours. |
Open with DEXTER |
![]() |
Figure 10:
The x ( upper panel) and y ( lower panel) peak flux positions relative to the 1-0 S(1) line are plotted against the upper energy level of each transition for
NK 1. Crosses represent the v = 1-0 transitions, squares the v = 2-1, and the triangle represents the [Fe II] 4D7/2 - 4F9/2 transition at 1.644 ![]() |
Open with DEXTER |
![]() |
Figure 11: The x ( upper panel) and y( lower panel) peak flux positions relative to the 1-0 S(1) line are plotted against the upper energy of each transition for SK 1. The symbols are as in Fig. 10 and the diamond represents the 3-2 S(3) peak position. In SK 1, the excitation increases from east to west ( upper panel - like NK 1) but also from north to south ( lower panel) along the jet axis. Peak postions were derived from the convolved data. |
Open with DEXTER |
The absolute positional errors are significantly larger, being
related to the telescope pointing. However, the integral field technique may
yield much smaller relative errors through the above
method. Moreover, a limit to the accuracy can be estimated from the
data by comparing the peak locations from lines originating from the
same upper energy level yet well separated in wavelength. On physical
principles, the line pairs possess identical flux distributions,
differing only in relative flux, possibly varying if there are sharp
spatial gradients in extinction. Most significantly, for SK1, the peak
locations from each of the five pairs agree to within 0.02
except for the y-offsets from the 1-0 S(4) and Q(6) transitions where the agreement is still within 2
(see Fig. 11).
For NK1, there are more significant discrepancies evident in the peak
locations of some of the weaker para 1-0 H2 line pairs
(Fig. 10, the difference reaching 0.1
in the y-offset). Given the accuracy of the SK 1 line peaks, this
suggests that the NK 1 discrepancies root from an intrinsically
non-Gaussian (asymmetric) central region to the line profile. Then, a convolution with the seeing and signal-to-noise would tend to remove
the skewness, shifting the peak position. This is indeed supported by
the line images of NK 1 (Fig. 3) which display a steeper
fall-off in flux towards the north-west.
The peak position offsets are plotted against the upper level
temperature of each transition in Figs. 10 and 11 for NK 1 and SK 1, respectively. Again, note that we assume Gaussian fits. The error is then related to the seeing of 0.5
and the signal-to-noise.
For NK 1 no trend can be seen along the outflow axis (in the y-direction). However, a clear trend is found perpendicular to the flow axis (in the x-direction). This unexpected result for an advancing bow shock indicates that the excitation increases slightly from east to west (left to right in all the figures presented here).
For SK 1 the higher excitation lines are shifted in the y-direction towards the south, away from the central source. Interestingly, in the x-direction the same trend is measured as in NK 1, with increasing excitation from east to west.
The offsets measured for the bows NB and SB show no trends (not displayed here). The detected emission is too weak to locate the higher excitation line emission peaks. More sensitive observations are required to complete such an analysis.
We wish to explain the spatial variation of the gas excitation evident as a) the offset of the [Fe II] emission ahead of the H2 peaks in the inner knots NK1 and SK1, b) a shift in the higher excitation
H2 emission of SK1 towards the south and c) westward shifts in the higher excitation
H2 emission lines across both inner knots. The variation corresponds to a change in gas temperature of order of 1000 K across a distance of 0.1
or
7
1014 cm. We consider three possibilities below.
The cooling length across a radiative shock, in which the
shock front and cooling layer are distinct, is very short, typically
under 1013 cm. Alternatively, with ambipolar diffusion,
the length scale of a C-type shock is
10
15/(ni cm-3) cm where ni is the ion number density (see
Smith & Brand 1990a). Note, however, that the section of the
shock across which the gas remains vibrationally excited is a moderate
fraction of this since the gas is no longer excited once accelerated
(Smith & Brand 1990b). Nevertheless, provided ni is less than
1 cm-3, the shock thickness could contribute to the
excitation gradient.
A second possibility is that the excitation varies smoothly along the surface of a curved bow shock. Since the shocks are clearly curved into bow-like shapes, this almost certainly is important. Modelling which includes the above two effects is carried out below.
Thirdly, two adjacent shocks of different strength may contribute. The shocks could be associated with forward and reverse shocks, the latter possibly taking the form of a Mach disk. Given the smooth variations across the inner knots, the two cooling zones would have to be only partially resolved. Although this cannot be excluded, the presence of emission from the iron ion in the inner knots, which is also displaced towards the location of high molecular excitation, suggests that the temperature variations occur across the surface of a bow shock.
Possible reasons for the displacement of higher excitation emission to the west in both knots are as follows.
![]() |
Figure 12:
The x ( upper panel) and y ( lower panel) offsets relative to
1-0 S(1) position predicted by the model. Note that
6.7 ![]() ![]() |
Open with DEXTER |
We attempt to model the inner knots NK1 and SK1 using various bow shock
configurations. We employ a steady state model which assumes a fixed
geometry. The three dimensional curved bow surface is described by
![]() |
(1) |
Table 4:
Model parameters derived to fit the inner knots NK 1 and SK 1 with a J-type shock.
is taken as the hydrogen nucleon density.
Table 5: Observed and predicted bow shock luminosities and line ratios.
NK 1 and SK 1 are similar in terms of their luminosities and
measured flux ratios. A consistent model is found to be a J-type
bow shock for which sufficient ions are assumed so that the
magnetic field is frozen to the neutral gas. Molecular hydrogen
dissociates when the shock speed exceeds 24 km s-1. The
model parameters are listed in Table 4. The model
1-0 S(1) luminosity and predicted line ratios are given in
Table 5 alongside the measured values for NK 1 and SK 1.
We also present the model images for SK 1 in Fig. 13 since this knot has a slightly higher flux and signal-to-noise ratio (S/N = 110). The derived shape of the model bow, with extended molecular flanks, is necessary to explain the low excitation (which is generally very high for planar J-shocks), as shown by the fits in Fig. 8. On the other hand, the quite low bow speed of 36 km s-1 (relative to the upstream jet flow) maintains the compact shape of the features (whereas fast J-shocks generate long tailed emission line regions).
The alternative C-type model generally consists of an apex region which is J-type
and bow flanks which are curved C-shocks in which the molecules are not dissociated.
We cannot exclude this model with the present data. However, we find that a paraboloidal bow of speed 65 km s-1 emits excessive [Fe II] emission from the bow cap. A more acute bow with shape s = 1.5, as also advanced for S233 N6 (Khanzadyan et al. 2004), is consistent with the HH 212 knot fluxes.
In order to introduce asymmetry into the modelled structure the
magnetic field was orientated at an angle of -54
relative to
the direction of motion of the bow (the geometry is described in
O'Connell et al. (2004)). Note that a high Alfvén speed is
necessary for the magnetic field to provide the asymmetry in a J-type
bow (Table 4).
In order to match the seeing, the model images were convolved with a Gaussian of FWHM = 3.5
1015 cm which is equivalent to the seeing throughout the observations of
0.5
at the adopted distance of 460 pc. The greatest similarity to the
observed bow appearance and line ratios was found for a bow of
velocity 36 km s-1 and pre-shock hydrogen nucleic density of
9.0
104 cm-3. The resulting bow is displayed in
Fig. 13, which shows the predicted appearance in the
1-0 S(1) and 2-1 S(1) lines of H2 as well as the [Fe II] emission which is generated from near the bow apex. Both the
unconvolved and convolved structures are shown. Figure 13
also shows the positions of the predicted emission peaks. The peaks in the model data
were determined using the same method as was used to find
the observed peak positions described above.
In Fig. 12, the predicted x and y peak positions from the SK 1 model bow relative to the 1-0 S(1) peak are plotted against the upper energy level for a number of emission line images (most of which were detected for SK 1). The trend in the y-direction bears a close resemblance to the observed offset positions for SK 1 in Fig. 11. Higher excitation line peaks are located closer to the front of the bow where the strongest excitation conditions are experienced. The introduction of an oblique magnetic field direction induces an asymmetry in the bow appearance. However, the offset trend in the x-direction does not match the observed trend. Firstly, the trend takes place on a much smaller scale than the observed shifts in positions and, secondly, the [Fe II] peak position is offset by a significant amount relative to the H2 lines, unlike the observed offset.
Figure 14 shows an image of the model 2-1 S(1)/1-0 S(1) vibrational excitation ratio. The excitation is seen to steadily increase towards the bow front. For comparison with Fig. 9 the predicted [Fe II]] emission is indicated by the overplotted contours. The similarities between the observed appearance and luminosities of the knots with the predictions of the best-fit model lead us to conclude that the bright structures aligning the inner jet regions of HH 212 are consistent with bow shocks. However, the asymmetry is probably due to some other effect rather than an oblique magnetic field (e.g. rotation or density gradient). Higher resolution observations should reveal the underlying curved structures which trace the outflowing jet material.
The appearance of the double bow shocks NB 1/NB 2 and SB 1/SB 2 is puzzling. Close examination of the jet knots in the high resolution image presented in McCaughrean et al. (2002) reveals that some of the knots can also be resolved into two components.
Variations in jet velocity can steepen and give rise to a series of double shocks which are occasionally visible as double bows in synthetic images of the 1-0 S(1) emission line derived from numerical simulations (Völker et al. 1999). Double shocks may also be due to the presence of a forward ambient shock front and a reverse Mach disk. In both cases, the separation between the two shock components is reduced by the strong cooling of the shocked gas. The clear separation seen in the outer bows would then be due to the lower density medium in these regions.
![]() |
Figure 13:
A J-type bow shock model for SK 1. The x and y axis units
represent 3.9 ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Such a configuration also provides an explanation for the inter-knot separation as illustrated in Fig. 15. Whereas the forward shock results in the entrainment of ambient material which is deflected outward, a curved reverse shock has the effect of entraining and refocusing the jet material into a more collimated flow. At high densities the forward and reverse shocks come into close contact and the reverse shock is forced to adopt a curved shape to match the forward bow. The reverse shock focuses and propels the jet material into a thin stream which approaches the forward shock. The high momentum jet penetrates through the forward shock at or close to the apex forming a continuous stream of inter-knot emission due to shock heating where the shearing of flow velocities occurs. This effect is also observed in numerical simulations, leading to a thin nose of material in advance of the bow shock which then appears as two shoulders (Rosen & Smith 2004).
![]() |
Figure 14: The model (2-1 S(1)/1-0) S(1) ratio image is shown with overplotted contours corresponding to the model [Fe II] emission for comparison with the equivalent SK1 data extracted from Fig. 9. |
Open with DEXTER |
Integral field spectroscopy permits simultaneous imaging over a wide range in wavelength. It can thus be used to reliably trace the excitation conditions across moderately-extended Herbig-Haro objects without concern about the influence of changing atmospheric transmission or seeing between image frames. The main findings for HH 212 are as follows.
![]() |
Figure 15: An explanation for the inter-knot emission is provided by the the combined effect of a forward shock and a curved reverse shock. Whereas the ambient material is entrained and propelled outward by the forward shock, the jet shock entrains the jet material into a more collimated and focused flow which, under certain conditions such as high density, can breach ahead of the forward bow. As shown in the lower sketch, the knots can be asymmetric with oblique shocks exciting gas in the inter-knot region. |
Open with DEXTER |
We have briefly considered mechanisms that could bring about a transverse distribution in excitation. Possibilities include (1) a relative drift of inner cloud material from west to east, (2) jet abrasion at the edges of an asymmetric outflow tunnel, (3) bow shock asymmetries due to the orientation of the magnetic field lines, and (4) a shear in the jet velocity. We find that option (3) is incapable of producing a significant transverse gradient.
A further indication of complex dynamics is the observed variations in the radial velocity
across the knots (Davis et al. 2000). These are generally but not
completely consistent with the observed direction of the radial velocity gradient in the
circumstellar disk. As an alternative to rotation in the jet, the variation
of order of a few km s-1, could be due to a shear in the jet velocity. A jet flow at 5
to the plane of the sky, the radial velocities then imply a shear along the
jet axis in the range 10-30 km s-1. This would be sufficient to significantly displace the point of impact of the bow shocks transverse to the jet axis.
Therefore, a more intricate bow model is indicated. We suggest that, similar to the outer bows, the inner knots may also consist of two adjacent but distinct bow shocks. The bows are not perfectly aligned and the reverse inner bow deflects jet gas preferentially towards the west. The existence of oblique shocks also provides an interpretation for the narrow inter-knot H2 emission (see Fig. 15. Higher spatial and velocity resolution in the 1-0 S(1) line would determine if this is the shock configuration.
Acknowledgements
We wish to thank the Department of Culture, Arts and Leisure, Northern Ireland and INTAS grant 4838 for financial aid. The UK Infrared Telescope is operated by the Joint Astronomy Centre on behalf of the UK Particle Physics and Astronomy Research Council. Some of the data reported here were obtained as part of the UKIRT Service Programme. The remaining data were obtained during UIST commissioning: we acknowledge Suzanne Ramsay Howat, Stephen Todd and Sandy Leggett for their support during these observations. We also thank Ignacio Ugarte Urra for assistance with the images.
Table 1: Photometric results for NK 1 (top) and SK 1 (bottom). The upper energy level of the transition is given in the third column. The x and y peak position offsets are relative to the 1-0 S(1) peak position.
Table 2: Photometric results for NB 1 and NB 2. The upper energy level of the transition is given in the third column. The x and y peak position offsets are relative to the 1-0 S(1) peak position.
Table 3: Photometric results for SB 1/SB 2. The upper energy level of the transition is given in the third column. The x and y peak position offsets are relative to the 1--0 S(1) peak position.