Mixed Jet-Cocoon Region in Astrophysics

• The mixed jet–cocoon region is defined as the finite-width interface where jet material and cocoon (or ambient) material interpenetrate, leading to significant mixing and instability.
• Pressure balance and various instabilities—such as Kelvin–Helmholtz and kink—shape its structure, influencing jet collimation, confinement, and radiation processes.
• Observational and simulation studies in GRBs, AGN, and microquasars highlight the region's role in jet propagation, energy dissipation, and diverse emission signatures.

The mixed jet–cocoon (MJC) region is the finite-width zone in a jet–cocoon system where jet material and cocoon or ambient material are no longer cleanly separable by an idealized contact surface. In magnetically confined gamma-ray-burst (GRB) jets, it is the part of the outflow above the collimation transition, at $z>z^\star$, where a once high-$\sigma$ jet is deflected and confined by the hot cocoon, becomes unstable, and is transformed into hot, low-$\sigma$, mixed jet–cocoon material (Levinson et al., 2012). In structured GRB and merger-jet models, it is the shear layer between an ultra-relativistic core and a broader cocoon, typically occupying $r_0<r<r_2$, where Kelvin–Helmholtz turbulence, magnetic entrainment, or shear particle acceleration operate (Garcia-Garcia et al., 2024, Wang et al., 16 Feb 2025). In AGN, microquasar, and collisionless jet studies, the same term denotes the turbulent or magnetized boundary layer surrounding shocked jet plasma and extending back from the head along the cocoon walls, or the contact layer that replaces a hydrodynamic discontinuity (Walg et al., 2013, Matsumoto et al., 2019, Bosch-Ramon et al., 2011, Dieckmann et al., 2018). Across these usages, the MJC region is the locus where confinement, instability, entrainment, and emission couple jet propagation to cocoon evolution.

1. Definitions and placement within jet–cocoon topology

The literature does not assign a single universal geometry to the MJC region. Instead, the term is attached to related structures that all occupy the dynamical interface between a fast jet and its cocoon or ambient surroundings.

Context	MJC placement	Characteristic description
Magnetically confined GRB jet	$z>z^\star$, below the head	deflected, cocoon-confined, low-$\sigma$ mixed flow
Structured GRB or merger jet	$r_0<r<r_2$	shear layer between core and cocoon
AGN or microquasar cocoon	from head back along cocoon walls, or between recollimation shock and head	turbulent backflow and entrainment zone
Collisionless pair jet	few $c/\omega_p$ thick contact layer	magnetic piston with electron interpenetration

In the analytic model of Levinson & Begelman, the MJC region is the portion of the outflow above the collimation transition and prior to reaching the head. There the jet cross-section is cylindrical, with radius $R_j(z)$ fixed by pressure balance at the jet boundary, and the cocoon consists of an inner light, shocked-jet layer and an outer shocked ambient layer, both at roughly uniform pressure $p_c(z)$ (Levinson et al., 2012).

In the semi-analytical model for a structured jet in a magnetized medium, the MJC region is the shear layer at the interface between the relativistic jet sheath and the surrounding shocked cocoon gas. The outflow is partitioned into an inner spine, a jet sheath, an MJC layer, a hot cocoon, and the unshocked ambient medium. In this definition, jet plasma and shocked ambient ejecta are Kelvin–Helmholtz unstable, and a fraction $\sigma$0 of the tangled external magnetic field is entrained into the cocoon (Garcia-Garcia et al., 2024).

In the AGN simulations of radial jet stratification, the MJC region is the layer of shocked plasma in which jet material and shocked ambient medium interpenetrate and exchange momentum and energy. It is identified with passive tracers and occupies the zone where both jet and ambient tracers coexist. In those simulations, the mixing layer runs from the jet head back along the cocoon walls, extending radially from $\sigma$1 to $\sigma$2, and axially from the Mach disk down to $\sigma$3 (Walg et al., 2013).

In high-mass microquasar studies, the MJC region is the volume downstream of the last strong recollimation shock and upstream of the jet head, where shocked jet plasma backflows, mixes turbulently with shocked ambient material, and inflates a hot overpressured cocoon. In collisionless pair-jet simulations, the corresponding structure is a transition zone across which positive charges remain largely segregated by a magnetic-pressure barrier, while ambient and jet electrons interpenetrate and rapidly mix; this layer explicitly replaces the hydrodynamic contact discontinuity (Bosch-Ramon et al., 2011, Dieckmann et al., 2018).

2. Confinement, pressure balance, and geometric scalings

The geometrical form of the MJC region is set by pressure balance. In the magnetically dominated jet model of Levinson & Begelman, the boundary condition is

$\sigma$4

For an equilibrium, rigidly rotating jet,

$\sigma$5

with $\sigma$6. If $\sigma$7, then $\sigma$8. In the same framework, the cocoon pressure is written as

$\sigma$9

and, in the strong reverse-shock regime,

$\sigma$0

with $\sigma$1 and $\sigma$2 determined by whether the head is subrelativistic or relativistic (Levinson et al., 2012).

In the magnetized-medium semi-analytic treatment, pressure balance is imposed simultaneously at the jet head, at the jet–cocoon interface, and at the cocoon–medium interface: $\sigma$3

$\sigma$4

$\sigma$5

These relations determine $\sigma$6, $\sigma$7, $\sigma$8, and $\sigma$9, and thereby define the angular width of the mixed layer,

$r_0<r<r_2$0

and its physical thickness,

$r_0<r<r_2$1

For fixed jet luminosity and density profile, $r_0<r<r_2$2 decreases with increasing $r_0<r<r_2$3 and increases with lower $r_0<r<r_2$4, where $r_0<r<r_2$5 is the ambient magnetization (Garcia-Garcia et al., 2024).

Large-scale jet systems admit analogous geometric parameterizations. In microquasar jet–cocoon systems, the fast jet of radius $r_0<r<r_2$6 is surrounded by a cocoon of radius $r_0<r<r_2$7, and the relevant MJC structure is a thin shear layer of thickness $r_0<r<r_2$8, with $r_0<r<r_2$9. In restarted 3C 84, the mini-cocoon inferred from 5 GHz space-VLBI appears as a low-surface-brightness envelope roughly $z>z^\star$0 pc across, with east and west limb-brightened shells of thickness $z>z^\star$1 pc and $z>z^\star$2 pc, respectively; the nearly constant cocoon pressure was proposed as a natural explanation for the almost cylindrical jet profile seen at 22 GHz (Zhang et al., 25 Jun 2025, Savolainen et al., 2021).

3. Instabilities, entrainment, and composition of the mixed layer

The MJC region is produced and maintained by instability. In magnetically confined GRB jets, the dominant destruction mechanism of the ordered toroidal field is a current-driven kink instability, with comoving growth time

$z>z^\star$3

and lab-frame growth length

$z>z^\star$4

Because $z>z^\star$5 whenever the jet is well collimated, the instability has time to grow before the fluid reaches the head. The result is essentially complete destruction of the toroidal field, conversion of Poynting flux into heat, and entrainment of cocoon plasma, so that the flow in the MJC region becomes low-$z>z^\star$6 mixed material (Levinson et al., 2012).

In structured jets propagating through a magnetized medium, Kelvin–Helmholtz shear drives turbulent entrainment of ambient gas and tangled magnetic flux into the cocoon, while magnetic tension and pressure oppose further shear mixing and provide additional lateral support. This competition is explicitly controlled by the parameter $z>z^\star$7, identified with the fraction of ambient magnetic-energy density allowed to penetrate the cocoon (Garcia-Garcia et al., 2024).

In relativistic hydrodynamic jet propagation, the MJC layer is driven by several instabilities at once. Matsumoto & Masada identify oscillation-induced Rayleigh–Taylor instability as the generic outcome when the effective inertia of the jet exceeds that of the cocoon during radial oscillations, with the effective inertia ratio

$z>z^\star$8

They further discuss Richtmyer–Meshkov, Kelvin–Helmholtz, and centrifugal instabilities as concurrent contributors. Their 3D simulations show finger-like protrusions of jet fluid penetrating into the cocoon by $z>z^\star$9, cocoon fluid intruding to similar depth, a mixed-mass fraction of order $\sigma$0 in the annulus around the jet core, and an MJC thickness $\sigma$1 that grows with propagation distance (Matsumoto et al., 2019).

The AGN spine–sheath simulations quantify mixing by the “absolute mixing” and “mass-weighted mixing” scalars

$\sigma$2

with $\sigma$3. In those models, the mixing-layer thickness grows approximately linearly with distance downstream,

$\sigma$4

and the extent and efficiency of mixing depend sensitively on whether the jet is homogeneous, isothermal spine–sheath, or piecewise isochoric spine–sheath (Walg et al., 2013).

Collisionless pair jets replace hydrodynamic entrainment by kinetic separation. The “magnetic piston,” a narrow band of strong $\sigma$5 with thickness $\sigma$6, repels ambient protons but does not impede electrons. Positive charges therefore remain largely segregated, whereas ambient and jet electrons leak across the piston on scales $\sigma$7 and become indistinguishable. In this sense, the MJC is a few electron skin depths thick and is defined by mixed electrons rather than mixed ions (Dieckmann et al., 2018).

4. Regimes set by mixing, magnetization, and external stratification

The MJC region is not a single dynamical state; it spans several limiting regimes. In long GRBs inside a massive progenitor star, the head motion is initially sub-relativistic and the jet remains strongly confined, with $\sigma$8, all the way to the stellar surface. Because $\sigma$9 in the collimation zone, field destruction can occur deep inside the star. Levinson & Begelman give a quantitative example in which $r_0<r<r_2$0 and $r_0<r<r_2$1 cm yield confinement out to $r_0<r<r_2$2, with $r_0<r<r_2$3 at breakout (Levinson et al., 2012).

A complementary breakout framework parameterizes the degree of stellar–jet mixing with a global parameter $r_0<r<r_2$4, defined as the fraction of shocked stellar mass entrained into the jet cocoon. In that model, no mixing, $r_0<r<r_2$5, gives a pure relativistic fireball and predicts a quasi-thermal pulse at $r_0<r<r_2$6 keV with $r_0<r<r_2$7, $r_0<r<r_2$8 s, and opening angle $r_0<r<r_2$9 rad; full mixing, $c/\omega_p$0, gives Newtonian expansion and optical/UV cooling-envelope emission; partial mixing, $c/\omega_p$1, yields a broad-angle, mildly relativistic outflow with $c/\omega_p$2 and a multi-timescale UV, optical, and X-ray signature (Nakar et al., 2016).

In merger-like or compact-object environments threaded by tangled magnetic fields, the same symbol $c/\omega_p$3 controls magnetic entrainment rather than baryon loading. Low mixing, $c/\omega_p$4, produces a slower-broader jet with a broader and more energetic cocoon, and a thick MJC layer with $c/\omega_p$5. High mixing, $c/\omega_p$6, boosts cocoon pressure, yields a faster-narrower jet with a narrow and less-energetic cocoon, and makes the MJC layer thin, $c/\omega_p$7. In the limit $c/\omega_p$8 G and $c/\omega_p$9, magnetic pressure dominates, the jet may be “laser-colimated,” $R_j(z)$0, and $R_j(z)$1, so the MJC layer collapses (Garcia-Garcia et al., 2024).

On larger scales, recollimation and head propagation regulate analogous mixed regions. In high-mass microquasars, the MJC region forms between the last strong recollimation shock and the reverse shock at the head; at $R_j(z)$2 yr, the first strong recollimation shock occurs at $R_j(z)$3 cm, additional pinching shocks appear at $R_j(z)$4 cm and $R_j(z)$5 cm, and the reverse shock lies at $R_j(z)$6 cm. In AGN cocoon simulations, the mixing layer extends roughly $R_j(z)$7 kpc behind the head, with morphology that depends on whether the jet is homogeneous or radially structured (Bosch-Ramon et al., 2011, Walg et al., 2013).

5. Radiation physics and observational manifestations

Because the MJC region is both dissipative and mixed, it is frequently invoked as an emission site. In structured GRB jets, Wang, Huang, and Liang model the MJC as a shear layer with exponential velocity profile

$R_j(z)$8

or equivalently

$R_j(z)$9

and derive a shear-acceleration timescale

$p_c(z)$0

For $p_c(z)$1, $p_c(z)$2 G, and $p_c(z)$3 cm, they obtain $p_c(z)$4. Fits to GRBs 090926A, 131108A, and 160509A give $p_c(z)$5 G and $p_c(z)$6, with the MJC $p_c(z)$7 component dominating in the keV–MeV band and the combined MJC plus core model reproducing Band or Band-Cut spectra (Wang et al., 2024).

The same series applies the MJC to off-axis and optical-flash phenomenology. For GRB 170817A, adopting $p_c(z)$8, $p_c(z)$9, $\sigma$00 G, $\sigma$01, and viewing angle $\sigma$02 rad yields an MJC synchrotron peak in the infrared and an SSC peak at hundreds of keV; in this off-axis case, the prompt spectrum is dominated by the MJC component because the core is strongly Doppler suppressed (Wang et al., 16 Feb 2025). For GRBs 990123, 080319B, and 130427A, fitted MJC parameters $\sigma$03 and $\sigma$04 G give $\sigma$05, so that MJC synchrotron reproduces bright prompt optical flashes while MJC-SSC produces X-ray excess and MeV–GeV flashes (Wang et al., 9 Aug 2025). A later extension distinguishes weak-scattering MJC emission, in which electrons become quasi-thermal with radial temperature profile

$\sigma$06

from strong-scattering MJC emission, in which shear acceleration creates a broader non-thermal hump. In the weak-scattering fit to GRB 090902B, $\sigma$07 keV and $\sigma$08 reproduce the narrow quasi-thermal peak (Wang et al., 24 May 2026).

Outside the GRB context, observational signatures emphasize boundaries and cocoons rather than prompt spectra. In AGN spine–sheath jets, the MJC layer is associated with slower mixed plasma, $\sigma$09, enhanced magnetic fields, limb-brightening in radio maps, spectral flattening, and hot-spot morphologies that differ between homogeneous and structured jets (Walg et al., 2013). In 3C 84, RadioAstron 5 GHz imaging revealed low-intensity emission from a cocoon-like structure around the restarted jet, with measured cocoon flux $\sigma$10 Jy at 5 GHz, $\sigma$11, minimum energy $\sigma$12 erg, and cocoon pressure $\sigma$13; roughly half of the jet power is inferred to go into the mini-cocoon enthalpy (Savolainen et al., 2021). In M87, Dainotti et al. report a soft X-ray dip approximately $\sigma$14 in projected size along the north side of the inner jet, with a $\sigma$15 surface-brightness depression and a soft-band significance of $\sigma$16; if interpreted as a cosmic-ray cocoon, the required energy is $\sigma$17 erg, but no firm discovery is claimed (Dainotti et al., 2012).

Microquasar studies extend MJC emission to high energies. In jet-termination models, the cocoon is treated as a homogeneous one-zone emitter with $\sigma$18, $\sigma$19 G in one representative case, radio fluxes of $\sigma$20 mJy at 5 GHz, X-ray fluxes of $\sigma$21 in 1–10 keV, and VHE $\sigma$22-ray fluxes of $\sigma$23 above 100 GeV (Bosch-Ramon et al., 2011). A separate shear-acceleration model for large-scale microquasar MJC layers predicts hard CR spectra up to $\sigma$24, of order a few PeV, and extended $\sigma$25-ray lobes up to $\sigma$26, with a comparable neutrino flux if $\sigma$27 interactions dominate (Zhang et al., 25 Jun 2025).

6. Conceptual scope, controversies, and open problems

A central conceptual issue is that “mixed jet–cocoon region” is not a single invariant object across the literature. In some works it is a post-collimation, low-$\sigma$28 channel produced by current-driven destruction of ordered magnetic field; in others it is a geometrically thin shear layer with explicit angular width $\sigma$29; in others it is a turbulent backflow zone, a mini-cocoon, a cosmic-ray cavity, or a collisionless magnetic piston (Levinson et al., 2012, Garcia-Garcia et al., 2024, Savolainen et al., 2021, Dieckmann et al., 2018). This suggests that the term is best understood as a family resemblance across jet–cocoon interface problems rather than as a single standardized fluid region.

The largest quantitative uncertainty is the degree of mixing. In breakout cocoon models, the absence of a wide-angle thermal $\sigma$30-ray pulse expected for $\sigma$31 is taken to indicate that mixing at some level must occur; partial and full mixing lead to qualitatively different terminal Lorentz factors and transient signatures (Nakar et al., 2016). In magnetized-medium semi-analytic models, $\sigma$32 controls how much of the ambient tangled field is entrained into the cocoon; two-dimensional hydrodynamical simulations were used to constrain this parameter, and the reported global collimation trend depends on magnetic-energy density more than on whether the field is tangled or poloidal (Garcia-Garcia et al., 2024).

Another open issue is observational ambiguity. The M87 soft X-ray dip is morphologically consistent with a dynamic cosmic-ray cocoon, but the highly non-uniform X-ray background prevents a definitive identification, and no clear radio-synchrotron sheath or enhanced boundary polarization coincident with the dip was found (Dainotti et al., 2012). By contrast, the 3C 84 mini-cocoon is supported by space-VLBI imaging, regularized-maximum-likelihood reconstructions, and closure-phase tests, but its particle content, encoded in the parameter $\sigma$33, remains constrained only indirectly (Savolainen et al., 2021).

Finally, dimensionality and microphysics matter. Matsumoto & Masada show that strictly 2D axisymmetric runs do not develop the non-axisymmetric modes responsible for MJC growth in 3D, implying that Rayleigh–Taylor, Richtmyer–Meshkov, centrifugal, and Kelvin–Helmholtz contributions cannot be inferred from axisymmetric calculations alone (Matsumoto et al., 2019). In collisionless pair jets, the MJC cannot be identified with a classical hydrodynamic contact discontinuity because electron mixing and proton segregation are controlled by a magnetic-pressure barrier rather than by collisional fluid separation (Dieckmann et al., 2018). In GRB prompt-emission models, the distinction between weak-scattering thermalized electrons and strong-scattering shear-accelerated electrons changes the MJC spectrum from a narrow quasi-thermal component to a broad non-thermal hump (Wang et al., 24 May 2026). The persistent open problem is therefore not merely where the MJC region is located, but which instability, transport regime, and composition define it in a given source class.