Herbig-Haro Object

A Herbig-Haro object (abbreviated HH object) is a small, transient patch of bright emission nebulosity associated with a newborn Class 0, Class I, or Class II T Tauri star, formed when narrow, highly collimated bipolar jets of partially ionized gas ejected from the star–disk boundary at supersonic velocities (v_jet ≈ 100 to 600 km s⁻¹) collide with ambient molecular cloud material or slower previously ejected wind shells. Discovered independently between 1946 and 1952 by American astronomer George Howard Herbig and Mexican astronomer Guillermo Haro while surveying dark lanes near NGC 1999 in the Orion Nebula complex (identifying the prototype knot systems HH 1, HH 2, and HH 3), Herbig-Haro objects were initially thought to be sites of in situ stellar birth until Richard Schwartz (1975) demonstrated that their optical spectra match radiative shock waves rather than photoionized H II regions. Lasting only 10³ to 10⁵ years and spanning 0.01 to several parsecs in length, Herbig-Haro outflows record the history of episodic mass accretion onto young stars and regulate star formation efficiency in open clusters and Bok globules by injecting supersonic momentum and turbulence into the parent molecular cloud.

Why does an object collapsing under gravity to form a star simultaneously launch powerful supersonic outflows? As gas spirals inward through a protoplanetary disk toward a central protostar (M_⋆ ≈ 0.1–2 M_☉), it must shed nearly all of its specific angular momentum j = (G M_⋆ r)^1/2 to accrete onto the stellar surface. Across all young stellar objects, the mass-outflow rate carried by the bipolar Herbig-Haro jet (Ṁ_jet) is tightly coupled to the disk mass-accretion rate (Ṁ_acc) at an empirical ratio of:

Ṁ_jet / Ṁ_acc ≈ 0.05 to 0.15

Two complementary magnetohydrodynamic (MHD) magnetocentrifugal launching mechanisms convert gravitational accretion energy and rotational angular momentum into collimated poloidal outflow along open magnetic field lines:

1. The Disk-Wind Model (Blandford & Payne, 1982; Pudritz & Norman, 1983): Open poloidal magnetic field lines (B_p) anchored across a range of radii (r₀ ≈ 0.1 to 5 astronomical units) in the ionized surface of a Keplerian disk rotate with the local angular velocity Ω₀ = (G M_⋆ / r₀³)^1/2. When the poloidal field lines are inclined outward by more than 30° from the vertical rotation axis (60° from the disk midplane), centrifugal force along the rigid rotating field line overcomes stellar gravity, flinging ionized gas outward like a bead on a spinning wire. Beyond the Alfvén surface (where outflow speed equals the Alfvén wave speed v_A = B_p / (4πρ)^1/2), differential rotation winds the magnetic field into a tight helical toroidal magnetic hoop (B_φ) whose inward Lorentz pinch force (J × B) collimates the outflow into a narrow cylindrical jet.
2. The X-Wind Model (Shu et al., 1994): A kilogauss stellar magnetosphere truncates the inner protoplanetary disk at the corotation radius R_X ≈ 3–5 R_⋆ (~0.03–0.05 au), where the Keplerian orbital period of the inner disk matches the stellar rotation period. While roughly two-thirds of the disk gas at R_X is funneled along closed stellar field lines onto high-latitude accretion hotspots on the T Tauri star, the remaining one-third is launched centrifugally along open field lines from the narrow X-region, carrying away the excess angular momentum that would otherwise spin the contracting star past break-up velocity.

In both mechanisms, the terminal jet velocity scales directly with the local gravitational escape velocity at the launching footpoint r₀:

v_jet ≈ (2 r_A / r₀)^1/2 v_esc(r₀) ≈ 150 to 500 km s⁻¹

Consequently, Herbig-Haro jets exhibit an "onion-skin" kinematic velocity structure: a fast, hot axial spine (v ≈ 200–400 km s⁻¹) launched from the innermost sub-au region, sheathed by slower, wider-angle molecular disk winds (v ≈ 10–50 km s⁻¹) launched from 1 to 20 au.

Although a Herbig-Haro jet emerges from the star–disk interface at v_jet ≈ 300 km s⁻¹, the internal sound speed of its T ≈ 10⁴ K gas is only c_s ≈ 10 km s⁻¹, yielding an internal Mach number of ℳ = v_jet / c_s ≈ 20 to 40. Where the supersonic jet head impinges on the stationary ambient molecular cloud—or where a newly ejected high-velocity pulse (v₂) overtakes a slower earlier pulse (v₁ < v₂) inside the jet beam—a double-shock working surface forms:

1. The Forward Bow Shock and Reverse Mach Disk

Every Herbig-Haro working surface consists of two opposing shock fronts separated by a dense layer of compressed post-shock gas:

• The Forward Bow Shock: Accelerates and sweeps up stationary ambient interstellar gas (or slower upstream wind material), curving backward into a parabolic or conical bow-shock wings geometry around the jet head.
• The Reverse Shock (Mach Disk): A planar shock front facing upstream into the jet beam that decelerates the incoming supersonic jet flow from v_jet down to the working-surface advance speed v_ws = v_jet / [1 + (ρ_amb / ρ_jet)^1/2].

Crucially, for internal working surfaces generated by episodic accretion velocity variability (Δv = v₂ − v₁ ≈ 30 to 80 km s⁻¹), although the knot moves across the sky with a high bulk proper-motion velocity of 200 to 400 km s⁻¹, the effective shock jump velocity entering the Mach disk and bow shock is only Δv ≈ 30 to 80 km s⁻¹. According to the Rankine–Hugoniot jump conditions for a strong shock (γ = 5/3), the immediate post-shock temperature is:

T_shock = (3 μ m_H / 16 k_B) Δv² ≈ 1.4 × 10⁵ K × (Δv / 100 km s⁻¹)²

Rather than heating the gas to X-ray temperatures of 10⁷ K (as would occur in a 400 km s⁻¹ stationary wall collision) and dissociating all molecules, internal velocity pulses produce moderate post-shock temperatures of T_shock ≈ 2 × 10⁴ to 10⁵ K that cool rapidly via optical forbidden lines and near-infrared H₂ emission (radiative shocks).

2. Spectroscopic Signature of Herbig-Haro Shock Cooling

Unlike photoionized Strömgren spheres or planetary nebulae (where ultraviolet stellar photons maintain a high ionization state dominated by [O III]), the post-shock cooling zone behind a Herbig-Haro shock recombines rapidly while remaining warm (T_e ≈ 5,000–10,000 K), producing an unmistakable low-excitation forbidden-line spectrum:

• Optical Low-Ionization Forbidden Lines: Intense emission from neutral and singly ionized species with low ionization potentials, specifically [S II] λλ671.6, 673.1 nm, [O I] λλ630.0, 636.3 nm, [N II] λλ654.8, 658.3 nm, [Fe II], and Balmer Hα 656.3 nm. On standard Baldwin–Phillips–Terlevich (BPT) diagnostic diagrams, HH objects occupy the shock-excited regime with high [S II]/Hα ≳ 0.5 and [O I]/Hα ≳ 0.2.
• Bow-Shock Apex vs. Wing Excitation Stratification: At the stagnation apex of a bow shock (where the normal shock velocity v_⊥ = v_ws is highest, ~80–150 km s⁻¹), higher-excitation [O III] λ500.7 nm emission appears, whereas along the oblique bow-shock wings (where only the perpendicular velocity component v_⊥ = v_ws sin θ is thermalized), lower-velocity J-type and C-type (continuous magnetic) shocks excite ro-vibrational H₂ 2.12 μm (v = 1–0 S(1)) and far-infrared CO and H₂O cooling lines.

High-angular-resolution observations across optical (Hubble Space Telescope), near- and mid-infrared (James Webb Space Telescope), and millimeter/submillimeter (ALMA) wavebands have established a set of benchmark Herbig-Haro systems that trace jet launching and shock physics across the full evolutionary sequence from deeply embedded Class 0 protostars to optically visible Classical T Tauri stars:

When time-series imaging with the Hubble Space Telescope, the James Webb Space Telescope (JWST), and wide-field optical CCD cameras traced Herbig-Haro knot chains across entire molecular clouds (such as HH 34 and HH 111), astronomers discovered that individual sub-parsec HH jets are merely the innermost, brightest segments of parsec-scale giant Herbig-Haro flows extending 1 to 10 parsecs from their driving protostars. Because the spacing between successive bow-shock working surfaces along a giant HH flow increases systematically from Δx ≈ 100–300 au near the source (corresponding to dynamical intervals Δt = Δx / v_jet ≈ 5–20 yr) to Δx ≈ 0.01–0.05 pc farther out (Δt ≈ 300–1,500 yr) and ~0.5–1 pc between giant terminal bow shocks (Δt ≈ 5,000–10,000 yr), Herbig-Haro flows act as fossil tape recordings of episodic disk accretion instabilities—linking short-term magnetospheric flares to decadal EX Lupi bursts and millennial FU Orionis eruptions.

As these supersonic jets propagate through their parent molecular cloud, prompt bow-shock entrainment and turbulent mixing along the jet shear layer sweep up 1 to 20 M_☉ of cold ambient CO and H₂ gas into wide bipolar molecular outflows moving at 10 to 30 km s⁻¹. By carving bipolar cavities through the infalling envelope and replenishing supersonic magnetohydrodynamic turbulence across parsec scales, Herbig-Haro jets limit the final stellar mass of their driving protostar and prevent runaway global free-fall collapse of the surrounding molecular cloud.

1. Herbig, G. H. (1951). "The spectra of two nebulous objects near NGC 1999." The Astrophysical Journal, 113, 697–699. https://doi.org/10.1086/145440
2. Haro, G. (1952). "Herbig's nebulous objects near NGC 1999." The Astrophysical Journal, 115, 572–573. https://doi.org/10.1086/145576
3. Blandford, R. D., & Payne, D. G. (1982). "Hydromagnetic flows from accretion discs and the production of radio jets." Monthly Notices of the Royal Astronomical Society, 199(4), 883–903. https://doi.org/10.1093/mnras/199.4.883
4. Reipurth, B., & Bally, J. (2001). "Herbig-Haro flows: Probes of early stellar evolution." Annual Review of Astronomy and Astrophysics, 39, 403–455. https://doi.org/10.1146/annurev.astro.39.1.403
5. Frank, A., Ray, T. P., Cabrit, S., et al. (2014). "Jets and outflows from star to cloud: Observations confront theory." In Protostars and Planets VI (pp. 451–474). University of Arizona Press. https://doi.org/10.2458/azu_uapress_9780816531240-ch020