A magnetized protostellar jet launched from the innermost disk at the truncation radius

Abstract

Protostellar jets form when part of the accreting material is energetically ejected from the vicinity of protostars. Understanding where and how they are launched and collimated is crucial to determining their role in the accretion process. The HH 211 jet is a highly collimated, magnetized jet associated with a rotating disk around a young protostar. With Atacama Large Millimeter/submillimeter Array, we have resolved a pristine molecular spine at its base down to the disk. This spine has a high velocity of ~ 107 ± 8 km s⁻¹ but a slow rotation with a specific angular momentum of ~ 4 ± 1 au km s⁻¹, suggesting it to be launched at ~ 0.021 ± 0.005 au in the disk, providing the most stringent constraint yet on current magneto-centrifugal theories of jet production. Quantitative modeling supports the interpretation that the molecular spine represents the dense central component of a magnetized radial wind. This wind, launched by magneto-centrifugal force at the innermost edge of the disk—the truncation radius, removes residual angular momentum from the disk, enabling disk material to accrete onto the protostar. The toroidal field strength to collimate the dense spine in the model also agrees with that previously measured in the jet.

Protostellar jets are believed to play an important role in star formation, carrying away angular momentum from the innermost disk region, allowing disk material to fall onto the central protostars. Their launching mechanism is still undetermined because the launching regions are within ~ 0.1 au of the protostars and thus too small to be resolved with current instruments. Two popular jet-launching models, the X-wind model¹ and the disk-wind model², have been proposed to explain where and how they are launched from the disks by magneto-centrifugal force. In the X-wind model, a wide-angle radial wind is launched from the innermost disk region at the truncation radius called the “X-point” and the jet is only the dense central spine of the wind. In the disk-wind model, an extended wind is launched from a range of disk radii and the jet is the innermost wind launched from a range of the innermost disk radii, with the velocity being mainly poloidal. Thus, these two models can be used to account for the jets with different amounts of specific angular momentum and different degree of radial expansion, especially at the jet base. In particular, the X-wind model has been used to explain the jets with small specific angular momentum and large radial expansion^3,4, while the disk-wind model has been used to explain the jets with small to large specific angular momentum^{5, 6, 7–8}. However, current measurements of specific angular momentum of the jets are quite uncertain due to shock contamination and insufficient spatial and velocity resolutions. Besides, accurate measurements of jet velocity and radial expansion are also needed to determine the jet-launching mechanism.

The HH 211 protostellar system is located in Perseus at 321 pc away. It is one of the youngest protostellar systems⁹, with a gravitational collapse age of only 35,000 yrs¹⁰. The central protostar currently has a mass of only 0.06 ± 0.02 ^10,11 and is expected to grow into a star like our Sun⁹. It is associated with an outflow, consisting of a highly collimated jet propagating supersonically along the outflow axis and outflow cavity walls around the jet^{12, 13–14}. The jet lies close to the plane of the sky¹⁵, optimal for measuring radial expansion and jet rotation. More importantly, the jet has been found to be magnetized¹⁶, supporting the role of magnetic field in the launching and collimation. Previous observations with Atacama Large Millimeter/submillimeter Array (ALMA) at 6 au resolution also spatially resolved an accretion disk around the protostar¹¹, connecting the jet launching to the disk accretion.

Figure 1a presents the recent JWST composite map of the HH 211 outflow in H and CO (taken on 2022 August 28) in near infrared at 30−45 au resolution¹⁷, showing the overall structure of the outflow, with a collimated jet surrounded by outflow cavity walls. The outflow is bipolar with a southeastern (SE) component on the approaching (blueshifted) side and a northwestern (NW) component on the receding (redshifted) side^13,18. As shown in Fig. 1b, although the jet is detected in H (red image) at a distance larger than 1000 au away from the protostar, no jet emission is detected in the inner region due to dust extinction by the surrounding dusty envelope.

Fig. 1 [Images not available. See PDF.]

JWST composite image and ALMA SiO and CO images of the HH 211 outflow. (a) The JWST composite image: blue for F335M (several H lines and continuum nebulosity), green for F460M (4.69 μm H line, CO rotational-vibrational band head emission, and continuum nebulosity), and red for F470N (4.69 μm H line) filters. The cross marks the position of the protostar located at ICRS and ¹¹. (b) Zooms in to the inner region in the JWST image. (c) The ALMA SiO image (green) with the JWST image. (d) The ALMA CO image (green) with the JWST image. The CO map is obtained with to km s⁻¹ for the SE component and 19 to 32 km s⁻¹ for the NW component to avoid any significant contamination of the shells, where is the radial velocity relative to the systemic velocity of 9.2 km s⁻¹. The insert shows the innermost part of the CO jet with the continuum map of the disk (gray image) reported before¹¹. The resolutions are 10 au for SiO and CO in panels c and d, and 7.5 au for CO in the insert.

In the following, we present our ALMA observations in SiO J = 8-7 and CO J = 3-2 at submillimeter wavelengths at a few times higher resolution of 10 au (i.e., ) (taken on 2021 August 29 and September 27, 29) to study this inner region. In addition, a velocity resolution of 0.42 km s⁻¹ is used to study its kinematics. Without suffering from dust extinction, the jet is clearly detected in SiO (green in Fig. 1c) and CO (green in Fig. 1d) at submillimeter wavelengths down to the protostar connecting to that detected in H . Here the CO jet is from the emission at high velocity (as discussed in Methods). The previously identified reflection-symmetric wiggle^19,20 can also be seen down to the protostar in the SiO and CO jets.

The SiO jet consists of a chain of knots with a (transverse) width growing with the distance from the protostar. These knots are roughly equally spaced and associated with a broad range of velocities, and appear as bow-like structures at low velocity (see Methods), supporting that they trace the internal shocks produced by a semi-periodical variation in jet velocity^{21, 22–23}. The CO jet shows a similar knotty structure to the SiO jet. Nonetheless, within 600 au of the protostar, the CO jet is narrower and passes through the center of the SiO jet, and thus can be considered as the central spine of the jet. As shown in the insert in Fig. 1d at 7.5 au resolution, the CO jet can be seen down to the center of the disk (gray image) to within 10 au of the protostar, with the innermost part further in blocked by the disk. More importantly, the CO jet there is roughly continuous and well collimated with a roughly constant radius (i.e., a half of the transverse Gaussian FWHM) estimated to be 13 au, much smaller than those measured further away in atomic lines²⁴.

Figure 2 shows the intensity maps and position-velocity (PV) diagrams of the SiO jet and CO jet within 200 au of the protostar to study the kinematics at the base. The first SiO knots E1 and W1 (as marked by the red line segments) are located at 100 au in the SE jet component and au in the NW jet component, respectively. They are associated with a broad range of velocities of 40 km s⁻¹ (see panels b and d), supporting that they are internal shocks formed in the jet. They appear as wide bow structures at low velocities as shown in Figures S1c and S1d. Moreover, a double-peak line profile characteristic of a bow shock^25,26 is also seen for knot W1 (see Fig. 3b), with one peak at high velocity and the other at low velocity. Only a single-peak line profile is seen for knot E1 at high velocity (see Fig. 3a) because the emission at low velocity is much fainter and does not show a clear peak. Therefore, they are likely wide internal bow shocks. However, since their SiO emission is faint in the bow wings, they appear as knots in the total intensity maps shown in Fig. 2, tracing the bow tips. Near the protostar, the SiO emission traces the V-shaped base of the shells around the jet (see Figs. S1c and S1d for clearer structures). They are associated with a broad range of velocities of 50 km s⁻¹, and thus whatever is driving them must have a significant radial velocity component at the base.

Fig. 2 [Images not available. See PDF.]

Intensity maps and PV diagrams of the SiO and CO jets at the base. The jet is rotated to be aligned with the vertical direction. Red image and red contours are for SiO, while green image and gray image with gray contours are for CO. The first SiO knots are marked with E1 and W1 on the SE and NW side, respectively. Contours start from 3 with a step of 3 for CO and a step of 4 for SiO, where K. Black vertical dashed lines mark the systemic velocity. The velocity resolution in the PV diagrams is 0.42 km s⁻¹.

Fig. 3 [Images not available. See PDF.]

Line profiles of the SiO and CO emission towards knots E1 and W1. The line profiles are obtained by integrating the emission over a circular region with a diameter of 30 au. The profiles have been smoothed with a hanning factor of 5.

On the other hand, the CO jet inner to the first SiO knots traces the pristine jet without any clear sign of shocks (i.e., a localized broad range of velocities) in the PV diagrams. The NW jet component has a mean radial velocity of 22.6 km s⁻¹ (see Methods). It has a velocity range increasing from the first knot W1 to the base, as outlined by the green curves, indicating that the radial component of the jet velocity increases towards the base. This behavior has been seen before in HH 212 and can be attributed to a radial expansion of the jet at the base⁴. Around the mean jet radial velocity, the peak emission also shows a sinusoidal PV structure that could be due to a semi-periodical variation in jet velocity with a wavelength of 50-60 au. On the other hand, the SE jet component has a mean radial velocity of −20 km s⁻¹ (see Methods). Although no clear PV structure can be identified around that velocity, the velocity range of the jet also appears to increase towards the base from the first knot E1. At knots E1 and W1, the CO emission is also associated with a range of velocities around the mean jet radial velocity, but not as broad as that seen in SiO. However, it is difficult to determine if the CO emission there also shows a double-peak line profile (see Fig. 3), because the emission at low velocity is contaminated by the internal shells and resolved out by the interferometry. Therefore, the CO emission there could trace the postshock gas or weakly shocked jet material.

Jet rotation can be measured from the PV diagrams cut across the jet axis. Figure 4 shows such PV diagrams constructed from the CO jet roughly midway between the protostar and the first SiO knots, i.e., 50 au in the SE jet component and −64 au in the NW jet component (with the cuts indicated with the white lines in Fig. 2a,c, respectively), where the jet is not significantly affected by internal shocks and contaminated by shell emission. To increase the sensitivity, we first convolved the map of the jet with a beam of 32 au 6.4 au ( ) elongated along the jet axis with position angle (P.A.) 116 and then produced the PV diagrams. As shown in Fig. 4b, the NW jet component can be seen around the mean jet radial velocity of 22.6 km s⁻¹ and be roughly separated from the shells. A linear velocity gradient with au (km s⁻¹) (as indicated by the cyan line) centered at 23 km s⁻¹ from 20.8 to 25.4 km s⁻¹ with the redshifted velocity in the northeast (NE) and blueshifted velocity in the southwest (SW) is seen across the jet axis. At each velocity channel, the white open box and the error bar (smaller than the box) indicate the peak position and the associated uncertainty obtained from a Gaussian fit to the intensity profile (see Table 1). The velocity gradient was then obtained by performing a linear fit to the peak positions, taking into account their uncertainties, with a reduced of 0.7. As shown in Fig. 5, this velocity gradient, although very small, can still be seen in the maps at the two velocity ends of the velocity gradient along the cut (marked by the black line). The red and blue curves show the peak positions of the redshifted and blueshifted emission, respectively, derived from the Gaussian fitting. The blueshifted emission beyond −90 au is excluded because it is contaminated by shell emission, forming a V-shaped structure that opens away from the central protostar. As can be seen, the position shift between the redshifted and blueshifted emission decreases towards the protostar and becomes unresolved closer in. This behavior has been seen before in the spinning jet HH 212⁴ because the velocity range of the jet increases with the decreasing distance toward the protostar due to a radial expansion, as found earlier. The velocity gradient has the same velocity sense as that seen in the disk and envelope¹⁰, and is thus likely arising from rotation in the jet. It is detected with a velocity range 4.6 km s⁻¹ across the cut position, thus the jet is estimated to have a specific angular momentum au km s⁻¹ about the jet axis. On the blueshifted side, no clear jet emission can be identified around the mean jet radial velocity of −20 km s⁻¹ likely because the jet and shell merge in the velocity (see Fig. 4a).

Fig. 4 [Images not available. See PDF.]

PV diagrams in CO cut across the jet axis to search for jet rotation. (a) The PV diagram for the SE jet component cut at 50 au and (b) The PV diagram for the NW jet component cut at -64 au. In (b), the cyan line marks the linear velocity gradient seen across the jet axis in the NW jet component, with the redshifted velocity in the northeast (NE) and blueshifted velocity in the southwest (SW). The contours start from 25 K with a step of 10 K. At each velocity channel, the white open box and the error bar (smaller than the box) indicate the peak position and the associated uncertainty obtained from a Gaussian fit to the intensity profile (see Table 1). The velocity gradient was then obtained by performing a linear fit to the peak positions, taking into account their uncertainties.

Fig. 5 [Images not available. See PDF.]

CO Channel maps at the two velocity ends of the linear velocity gradient across the jet axis. The black line indicates the cut position. Blue contours for the velocity km s⁻¹ while red contours for the velocity km s⁻¹. Red contours start from 20 K with a step of 10 K. Blue contours start from 22 K with a step of 11 K. The thick contours indicate the highest intensity levels and highlight the velocity gradient. The red and blue curves show the peak positions of the redshifted and blueshifted emission, respectively, derived from the Gaussian fitting. The ellipse shows the beam size.

Table 1

Best-fit parameters for the centroids in the PV diagram using a Gaussian Profile with .

Velocity	Amplitude (A)	Peak Position ( )	FWHM
(km/s)	(K)	(au)	(au)
21.01	66.93 ± 4.04	− 1.07 ± 0.84	28.35 ± 1.98
21.43	70.86 ± 4.22	− 2.07 ± 0.80	27.25 ± 1.87
21.86	75.32 ± 3.25	− 1.42 ± 0.59	27.68 ± 1.38
22.28	77.74 ± 2.75	− 0.48 ± 0.47	26.95 ± 1.10
22.70	75.48 ± 2.85	0.63 ± 0.49	26.41 ± 1.15
23.13	73.69 ± 2.76	0.59 ± 0.46	25.09 ± 1.08
23.55	72.20 ± 1.98	0.62 ± 0.33	24.50 ± 0.77
23.97	71.13 ± 1.89	1.01 ± 0.33	25.58 ± 0.78
24.40	69.96 ± 2.17	1.02 ± 0.37	24.07 ± 0.86
24.82	63.24 ± 2.75	1.61 ± 0.43	20.37 ± 1.02
25.24	51.09 ± 3.71	1.82 ± 0.70	19.71 ± 1.65

Independent of the uncertain launching conditions, for a jet launched by magneto-centrifugal force from a rotating disk, the launching radius at the foot point of the jet in the disk can be derived assuming a conservation of energy and angular momentum along the field line⁶. In particular, it can be derived from the observed velocity ( ) and specific angular momentum ( ) of the jet at a large distance from the launching point and the mass of the protostar ( ). With ^10,11, km s⁻¹ (see Methods), and au km s⁻¹, the launching radius is estimated to be au, where the specific angular momentum au km s⁻¹. Here, Equation 2 in Lee et al.²⁷, which was derived for a small value as estimated here, was used to estimate the value. Since the protostar has a theoretical radius of 0.01 au²⁸, which is half of the launching radius, the jet is likely launched at the innermost edge of the disk truncated by the stellar field¹. The magnetic lever arm parameter, which determines the efficiency of acceleration and transport of angular momentum, is . Therefore, the Alfvén radius within which most of the acceleration is achieved is au.

We can further constrain the jet-launching models by comparing our results to the models more quantitatively. Since the CO jet in HH 211 is launched from the innermost edge of the disk, we adopt the asymptotic X-wind model²⁹ for which analytic solutions exist and 3.7. We defer comparing our results with the disk-wind model, because it is not clear if a large widening factor of 620 to produce the observed jet radius is feasible for a small ⁸.

Figure 6 shows the density, velocity, and streamlines of the wind in the X-wind model. To mimic the observations, we also include a dusty disk to block the innermost jet and impose an outer boundary of the wind defined by the observed CO shell (see Methods) as outlined by the red curve. As can be seen, the model launches a radial wind with a dense spine along the outflow axis. The tenuous wide-angle radial wind component is expected to produce the observed broad range of velocities for the V-shaped SiO shells at the base and the wide bow structures³⁰ for the first knots W1 and E1, although simulations are needed to confirm these. This model is self-similar and scaled with the mass of the protostar , the radius of the launching point , and the mass-loss rate of the wind . To achieve the observed jet velocity, it requires 0.017 au, similar to the estimated value. Assuming that the accretion rate in the disk is smaller than the infall rate in the envelope, which is ¹⁰ , we have . Here we set .

Fig. 6 [Images not available. See PDF.]

The X-wind model for the jet. Left panel shows the density distribution (color image), velocity field (vectors), and streamlines (black curves) of the wind in our model. Thick red curves mark the CO shell for the outer boundary of the wind. Black hashed region shows the dusty disk. Right panel shows the model in 3D created with VisIt. It also shows the magnetic field at the base.

Figure 7 shows the model results. We only model the NW jet component because the kinematics of the SE jet component is contaminated by the shells. For comparison, we also show the observed intensity of the jet in Fig. 7a. As can be seen, this model can readily produce the key features in the observations, e.g., a highly collimated jet structure with roughly similar intensity (see Fig. 7b), a velocity range increasing toward the base (see Fig. 7c), and a linear velocity gradient across the jet axis (see Fig. 7d). Nonetheless, comparing to the observed jet, the jet in the model is narrower and has a broader velocity range away from the base. This implies that the streamlines of the wind in the model are required to be more collimated away from the base with less expansion velocity perpendicular to the jet axis. Further work is needed to determine if this can be achieved by fine-tuning the values and the functional forms of the specific angular momentum and ratio of magnetic field to mass-flux carried on a streamline in the X-wind model²⁹.

Fig. 7 [Images not available. See PDF.]

Model results in comparison with the observations in CO for the NW jet component at the base. (a,b) The observed and model intensity maps of the jet, respectively. The first SiO knot is marked with W1. (c) The PV diagram along the jet axis, as shown in Fig. 2d. (d) The PV diagram along the jet axis, as shown in Fig. 4b. Red contours are from the model.

Previously, SiO line polarization was detected in this jet at a distance of 300−500 au and the field strength was estimated to be 15 mG assuming a density of ¹⁶ . In this model, the toroidal field is dominant over the poloidal field at large distance and has a field strength G at those distances. The cylindrical radius of the wind at the assumed is 10 au (see Fig. 6) and thus the field strength is 10 mG, also roughly in agreement with the observed value.

In summary, the X-wind model can provide a reasonable match to the observations of the CO jet at the base, supporting that the jet is only the dense central spine of the wide-angle radial wind. The wind removes residual angular momentum from the disk, enabling disk material to accrete onto the protostar. This model has also been used to reproduce the radial expansion and rotation of the knots detected in SiO at the base of the HH 212 jet^4,27. However, the SiO knots can be affected by shocks. In addition, SiO emission traces only the densest spine of the wind due to its high critical density. On the other hand, the CO jet base here in HH 211 appears not to be affected by shocks. It also appears wider due to the lower critical density of the CO emission, and thus covering more streamlines of the wind, providing stronger constraint on the model, especially on the collimation. On the other hand, since we defer comparing our results with the disk-wind model due to a lack of appropriate analytical formulation available, further work is needed to check whether the disk-wind model could also provide a possible interpretation.

Acknowledgements

We thank the referees for their constructive comments that improve the paper. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2019.1.00570.S and 2015.1.00024.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSTC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. C.-F.L. acknowledges grants from the National Science and Technology Council of Taiwan (110-2112-M-001-021-MY3 and 112-2112-M-001-039-MY3) and the Academia Sinica (Investigator Award AS-IA-108-M01).

Author contributions

C.-F. Lee and K.S Jhan led the project. C.-F. Lee led the analysis, and discussion, and drafted the manuscript. K.S. Jhan and A. Moraghan commented on the manuscript and participated in the discussion.

Data availability

This paper makes use of the following ALMA data: ADS/JAO.ALMA#2019.1.00570.S and 2015.1.00024.S. The data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request.

Declarations

Competing interests

The authors declare no competing interests.

Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1038/s41598-025-11602-w.

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Shu, F. H., Najita, J. R., Shang, H., & Li, Z.-Y. X-Winds Theory and Observations. Protostars and Planets IV, 789 (2000).

2. Konigl, A. & Pudritz, R. E. Disk winds and the accretion-outflow connection. In Protostars and Plane IV (eds. Mannings, V., Boss, A. P. & S) 759 (2000). https://doi.org/10.48550/arXiv.astro-ph/9903168

3. Lee, C-F; Hirano, N; Palau, A; Ho, PTP; Bourke, TL; Zhang, Q; Shang, H. Rotation and outflow motions in the very low-mass class 0 protostellar system HH 211 at subarcsecond resolution. ApJ; 2009; 699, 1584.2009ApJ..699.1584L1:CAS:528:DC%2BD1MXpslOnsro%3D [DOI: https://dx.doi.org/10.1088/0004-637X/699/2/1584]

4. Lee, C-F; Li, Z-Y; Shang, H; Hirano, N. Magnetocentrifugal origin for protostellar jets validated through detection of radial flow at the jet base. ApJL; 2022; 927, L27.2022ApJ..927L.27L [DOI: https://dx.doi.org/10.3847/2041-8213/ac59c0]

5. Bacciotti, F; Ray, TP; Mundt, R; Eislöffel, J; Solf, J. Hubble space telescope/STIS spectroscopy of the optical outflow from DG Tauri: Indications for rotation in the initial jet channel. ApJ; 2002; 576, 222.2002ApJ..576.222B [DOI: https://dx.doi.org/10.1086/341725]

6. Anderson, JM; Li, Z-Y; Krasnopolsky, R; Blandford, RD. Locating the launching region of T Tauri winds: The case of DG Tauri. ApJL; 2003; 590, L107.2003ApJ..590L.107A [DOI: https://dx.doi.org/10.1086/376824]

7. Coffey, D; Bacciotti, F; Ray, TP; Eislöffel, J; Woitas, J. Further indications of jet rotation in new ultraviolet and optical hubble space telescope STIS spectra. ApJ; 2007; 663, 350.2007ApJ..663.350C1:CAS:528:DC%2BD2sXosVanuro%3D [DOI: https://dx.doi.org/10.1086/518100]

8. Tabone, B; Cabrit, S; Pineau des Forêts, G; Ferreira, J; Gusdorf, A; Podio, L; Bianchi, E; Chapillon, E; Codella, C; Gueth, F. Constraining MHD disk winds with ALMA Apparent rotation signatures and application to HH212. A&A; 2020; 640, A82.2020A&A..640A.82T1:CAS:528:DC%2BB3cXitVKlsLfO [DOI: https://dx.doi.org/10.1051/0004-6361/201834377]

9. Froebrich, D. Which are the youngest protostars? Determining properties of confirmed and candidate class 0 sources by broadband photometry. ApJS; 2005; 156, 169.2005ApJS.156.169F [DOI: https://dx.doi.org/10.1086/426441]

10. Lee, C-F; Kwon, W; Jhan, K-S; Hirano, N; Hwang, H-C; Lai, S-P; Ching, T-C; Rao, R; Ho, PTP. A pseudodisk threaded with a toroidal and pinched poloidal magnetic field morphology in the HH 211 protostellar system. ApJ; 2019; 879, 101.2019ApJ..879.101L1:CAS:528:DC%2BC1MXit1KhsrjK [DOI: https://dx.doi.org/10.3847/1538-4357/ab2458]

11. Lee, C-F; Jhan, K-S; Moraghan, A. First detection of a linear structure in the midplane of the young HH 211 protostellar disk: A spiral arm?. ApJL; 2023; 951, L2.2023ApJ..951L..2L [DOI: https://dx.doi.org/10.3847/2041-8213/acdbca]

12. McCaughrean, MJ; Rayner, JT; Zinnecker, H. Discovery of a molecular hydrogen jet near IC 348. ApJL; 1994; 436, L189.1994ApJ..436L.189M1:CAS:528:DyaK2MXjsVCqsbw%3D [DOI: https://dx.doi.org/10.1086/187664]

13. Gueth, F; Guilloteau, S. The jet-driven molecular outflow of HH 211. A&A; 1999; 343, 571.1999A&A..343.571G1:CAS:528:DyaK1MXit1Slsbc%3D

14. Lee, C-F; Ho, PTP; Palau, A; Hirano, N; Bourke, TL; Shang, H; Zhang, Q. Submillimeter arcsecond-resolution mapping of the highly collimated protostellar Jet HH 211. ApJ; 2007; 670, 1188.2007ApJ..670.1188L1:CAS:528:DC%2BD1cXivFemsg%3D%3D [DOI: https://dx.doi.org/10.1086/522333]

15. Jhan, K-S; Lee, C-F. 25 au angular resolution observations of HH 211 with ALMA: jet properties and shock structures in SiO, CO, and SO. ApJ; 2021; 909, 11.2021ApJ..909..11J1:CAS:528:DC%2BB3MXhtlOisrvO [DOI: https://dx.doi.org/10.3847/1538-4357/abd6c5]

16. Lee, C-F; Li, Z-Y; Hirano, N; Shang, H; Ho, PTP; Zhang, Q. ALMA observations of the very young class 0 protostellar system HH211-mms: A 30 au dusty disk with a disk wind traced by SO?. ApJ; 2018; 863, 94.2018ApJ..863..94L1:CAS:528:DC%2BC1cXisVamsrfP [DOI: https://dx.doi.org/10.3847/1538-4357/aad2da]

17. Ray, TP; McCaughrean, MJ; Caratti o Garatti, A et al. Outflows from the youngest stars are mostly molecular. Nature; 2023; 622, 48.2023Natur.622..48R1:CAS:528:DC%2BB3sXitVWitbvN [DOI: https://dx.doi.org/10.1038/s41586-023-06551-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37619607][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10550818]

18. Hirano, N; Liu, S-Y; Shang, H; Ho, PTP; Huang, H-C; Kuan, Y-J; McCaughrean, MJ; Zhang, Q. SiO J = 5–4 in the HH 211 protostellar jet imaged with the submillimeter array. ApJL; 2006; 636, L141.2006ApJ..636L.141H1:CAS:528:DC%2BD28Xhs1Cqtrg%3D [DOI: https://dx.doi.org/10.1086/500201]

19. Lee, C-F; Hasegawa, TI; Hirano, N; Palau, A; Shang, H; Ho, PTP; Zhang, Q. The reflection-symmetric wiggle of the Young protostellar Jet HH 211. ApJ; 2010; 713, 731.2010ApJ..713.731L1:CAS:528:DC%2BC3cXmtlSmu7w%3D [DOI: https://dx.doi.org/10.1088/0004-637X/713/2/731]

20. Moraghan, A; Lee, C-F; Huang, P-S; Vaidya, B. A study of the wiggle morphology of HH 211 through numerical simulations. MNRAS; 2016; 460, 1829.2016MNRAS.460.1829M1:CAS:528:DC%2BC1cXlvFWmurw%3D [DOI: https://dx.doi.org/10.1093/mnras/stw1089]

21. Raga, AC; Canto, J; Binette, L; Calvet, N. Stellar jets with intrinsically variable sources. ApJ; 1990; 364, 601.1990ApJ..364.601R [DOI: https://dx.doi.org/10.1086/169443]

22. Stone, JM; Norman, ML. Numerical simulations of protostellar jets with nonequilibrium cooling II. Models of pulsed jets. ApJ; 1993; 413, 210.1993ApJ..413.210S [DOI: https://dx.doi.org/10.1086/172989]

23. Lee, C-F; Sahai, R. Magnetohydrodynamic models of the bipolar Knotty jet in Henize 2–90. ApJ; 2004; 606, 483.2004ApJ..606.483L [DOI: https://dx.doi.org/10.1086/381677]

24. Caratti o Garatti, A; Ray, TP; Kavanagh, PJ et al. JWST observations of Young protoStars (JOYS): HH211: Textbook case of a protostellar jet and outflow. A&A; 2024; 691, A134.1:CAS:528:DC%2BB2MXhtFWgurk%3D [DOI: https://dx.doi.org/10.1051/0004-6361/202451350]

25. Raga, A. C. & Bohm, K. H. Predicted long-slit, high-resolution emission-line profiles from interstellar bow shocks. II. Arbitrary Bow Shock Orientation. APJ308, 829. https://doi.org/10.1086/164554 (1986).

26. Hartigan, P; Raymond, J; Hartmann, L. Radiative bow shock models of Herbig-Haro objects. ApJ; 1987; 316, 323.1987ApJ..316.323H1:CAS:528:DyaL2sXhvVGlu78%3D [DOI: https://dx.doi.org/10.1086/165204]

27. Lee, C-F; Ho, PTP; Li, Z-Y; Hirano, N; Zhang, Q; Shang, H. A rotating protostellar jet launched from the innermost disk of HH 212. Nat. Astron.; 2017; 1, 0152.2017NatAs..1.152L [DOI: https://dx.doi.org/10.1038/s41550-017-0152]

28. Stahler, SW. Deuterium and the Stellar birthline. ApJ; 1988; 332, 804.1988ApJ..332.804S1:CAS:528:DyaL1cXlslyhtbc%3D [DOI: https://dx.doi.org/10.1086/166694]

29. Shu, FH; Najita, J; Ostriker, EC; Shang, H. Magnetocentrifugally driven flows from young stars and disks. V. Asymptotic collimation into jets. ApJL; 1995; 455, L155.1995ApJ..455L.155S [DOI: https://dx.doi.org/10.1086/309838]

30. Lee, C-F; Stone, JM; Ostriker, EC; Mundy, LG. Hydrodynamic simulations of jetand wind-driven protostellar outflows. ApJ; 2001; 557, 429.2001ApJ..557.429L1:CAS:528:DC%2BD3MXmslKisL4%3D [DOI: https://dx.doi.org/10.1086/321648]