Astrophysical Thermal Jet Classification

• Thermal jet classification is a framework that categorizes astrophysical jets by distinguishing thermal free–free and non-thermal synchrotron emissions and linking them to underlying jet parameters.
• It employs multi-color blackbody modeling combined with non-thermal power-law components to derive key metrics such as magnetization, dimensionless entropy, and launch radius.
• The approach is applied to GRB and protostellar jets, enabling differentiation between pure fireball, hybrid, and mixed emission regimes.

Thermal jet classification encompasses the empirical and physically motivated categorization of astrophysical jets according to the dominant emission processes—thermal versus non-thermal—and the underlying jet structure and central engine parameters. In high-energy transients such as gamma-ray bursts (GRBs) and ionized outflows in star-forming regions (e.g., protostellar jets), rigorous model fitting and inversion techniques have supported the development of structured schemes distinguishing "fireballs", mildly magnetized hybrid jets, and mixed thermal/non-thermal outflows. These frameworks integrate spectral diagnostics, radiative-transfer modeling, and inversion of observables (temperature, flux, spectral index) to constrain fundamental parameters such as magnetization ($\sigma_0$), dimensionless entropy ($\eta$), launch radius ($r_0$), and relativistic electron fraction ($\eta_e^{\rm rel}$).

1. Physical Basis of Thermal and Non-Thermal Emission

Jet classification leverages fundamental differences in emission mechanisms:

• Thermal Free–Free (Bremsstrahlung) Emission: In highly ionized plasma environments, the spontaneous emissivity at frequency $\nu$ is given by $$j_{\rm ff}(\nu) = 6.8 \times 10^{-38}\, n_e^2\, T^{-1/2}\, \exp(-h\nu/kT)\,\bar{g}_{\rm ff}(\nu, T)$$, where $n_e$ is electron density, $T$ temperature, and $\bar{g}_{\rm ff}$ is the Gaunt factor, encapsulating quantum effects. Opacity follows $$\alpha_{\rm ff}(\nu) = 0.08235\, T^{-1.35}\, n_e^2\, \nu^{-2.1}\,L$$.
• Non-Thermal Synchrotron Emission: Relativistic electrons in magnetic fields emit per $$j_{\rm syn}(\nu)\propto N_{\rm rel}\, B^{(p+1)/2}\, \nu^{-(p-1)/2}$$, with $N(E)dE=N_k\,E^{-p}\,dE$ as the electron energy distribution, $B$ the field strength, and resulting spectral index $\alpha_{\rm syn} = -(p-1)/2$. Synchrotron self-absorption is included by $$\alpha_{\rm syn}(\nu)\propto N_{\rm rel}\, B^{(p+2)/2}\, \nu^{-(p+4)/2}$$ (Mohan et al., 2023).

In GRB jets, photospheric emission is empirically described using single or multiple blackbody components (mBB), while non-thermal contributions are fit by power laws (PL) or cut-off power laws (CPL) (Song et al., 2023).

2. Empirical Model Fitting and Inversion Methodologies

The principal schemes for classifying thermal jets employ multi-component spectral fits and "top-down" diagnostic inversions:

• Multi-Color Blackbody (mBB) Modeling: The photon flux is modeled as a superposition of blackbodies over a temperature range $[T_{\min}, T_{\max}]$, weighted by $F(T)\propto T^q$ or $dF/dT$ with index $m$—encoding the steepness/curvature of the distribution:

$$F(T) = F_{\rm max}\, (T/T_{\rm max})^q,\,\, T_{\min}\leq T\leq T_{\max}$$

or

$$\frac{dF}{dT} = \frac{m+1}{(T_{\rm max}/T_{\min})^{m+1}-1}\, F_{\rm mBB}\, T_{\min}\, \left(\frac{T}{T_{\min}}\right)^m$$

• Combined mBB + Non-Thermal Component: In cases where mBB broadening fails, a PL or CPL component $$N_{\rm NT}(E)=K_{\rm NT}\,E^\alpha\,\exp[-E/E_{\rm c}]^\beta$$ is added. The strength and slope of the PL tail (typically $\alpha \lesssim -1.7$ and up to 30% of $F_{\rm mBB}$) and fit improvement ($\Delta\mathrm{BIC}\gtrsim 6$) serve as diagnostics for genuine non-thermal processes (Song et al., 2023).
• Top-Down Magnetization Inversion: The observed $T_{\rm max}$ and $F_{\rm max}$ are inverted via semi-analytic prescriptions to yield engine parameters:
  • Launch radius $r_0$ from acceleration laws
  • Dimensionless entropy $\eta = L_m / (\dot{m}c^2)$
  • Magnetization $\sigma_0$ using $$1+\sigma_0 \propto (kT_{\rm obs}/F_{\rm BB})^{4/3} L_w r_0^{2/3}$$

Solutions to $\{r_{\rm ph}, T_{\rm obs}, F_{\rm BB}\}$ for measured data constrain jet regimes and discriminate between hot fireball ($\sigma_0 \approx 0$) and hybrid ($\sigma_0 \sim 1$) (Song et al., 2023).

In protostellar jets, two-component radio spectrum models fit thermal bremsstrahlung and non-thermal synchrotron to multi-epoch, multi-frequency data, leveraging $\chi^2$ minimization across a coarse parameter grid ($n_e$, $p$, $\eta_e^{\rm rel}$, $B$, $\delta\theta$) (Mohan et al., 2023).

3. Spectral-Index Diagnostics and Classification Criteria

The spectral index $$\alpha(\nu_1,\nu_2) = \log[S(\nu_2)/S(\nu_1)] / \log(\nu_2/\nu_1)$$ distinguishes emission mechanisms:

Jet Type	Typical $\alpha$	Relativistic Fraction $\eta_e^{\rm rel}$	Polarization Signature
Thermal	$\alpha \gtrsim 0$	$\lesssim 10^{-7}$	None
Non-Thermal	$\alpha \lesssim -0.15$	$10^{-7}$–$10^{-4}$	$P > 1\%$
Hybrid/Mixed	$-0.15 < \alpha < +0.11$	$10^{-7}$–$10^{-4}$	Partial ($P < 1\%$) or subtle

A plausible implication is that such combined criteria facilitate robust segregation of outflow types even in ambiguous cases bridging the thermal/non-thermal divide (Mohan et al., 2023).

For GRB jets, further quantitative thresholds are employed:

• If $(1+\sigma_0)$ fits unity (within uncertainties): fireball (thermal-dominated)
• If $(1+\sigma_0)\gtrsim 1.2$ in any time bin: hybrid (mild magnetization)
• Spectral curvature index $m>0.2$ and no PL: thermal; $m<0$ or required PL: hybrid (Song et al., 2023).

4. Application to Exemplar Systems

GRB Case Studies

• GRB 210121A (Pure Fireball): $m\approx 0.5$, $kT_{\rm max}\approx 400$ keV, $F_{\rm max}\approx 1.6\times10^{-5}$ erg/cm$^2$/s, $(1+\sigma_0)=1.0$–$1.2$, no PL required. Classified as pure hot fireball, spectrum consistent with non-dissipative photospheric emission (Song et al., 2023).
• GRB 210610B (Mildly Magnetized Hybrid): $m\approx 0.8$, $kT_{\rm max}\approx 133$ keV, $F_{\rm max}\approx 3.0\times10^{-7}$, PL tail with $\alpha\sim-1.7$, $(1+\sigma_0) = 1.2$–$2.0$, $\eta\sim400$–$800$. PL statistically required ($\Delta$BIC$\,{>}\,6$) (Song et al., 2023).
• GRB 221022B (Thermal to Non-Thermal Transition): Early bins $m\approx 0.6$, $kT_{\rm max}\approx 98$ keV, $F_{\rm max}\approx 6.3\times10^{-7}$, no PL. Later bins: PL emerges, non-thermal component dominates; $(1+\sigma_0)\approx 1.0\pm 0.2$ initially (Song et al., 2023). This suggests transition via internal shock as $\Gamma(t)$ increases.

Protostellar Jet (HH 80–81)

• Electron densities: $200$–$3\times10^5$ cm$^{-3}$
• Relativistic fraction: $10^{-7}$–$10^{-4}$
• Spectral indices: $-0.52 \leq \alpha \leq +0.13$
  • IRAS 18162–2048 (central): thermal ($\alpha \gtrsim +0.1$)
  • HH 80, HH 81: mixed or non-thermal ($\alpha < -0.15$ for efficient particle acceleration) (Mohan et al., 2023).

A plausible implication is that such parameter mapping enables differentiation between central thermal beams and peripheral synchrotron-emitting shocked shells.

5. Generalized Classification Frameworks

Structured schemes classify jets by key observables and inferred parameters:

• Thermal (Fireball, Central Source Dominated): High base density, thick emission region, negligible relativistic electrons, rising or flat radio spectrum, lack of polarization.
• Non-Thermal (Shocked Layer, Hybrid Outflows): Lower density, thin layer, efficient acceleration, pronounced synchrotron contribution, negative spectral index, detected polarization (Mohan et al., 2023).
• Hybrid/Mixed Cases: Intermediate values; both emission mechanisms significant, partial polarization, spectral curvature.

Multi-frequency coverage (0.3–15 GHz), angular resolution ($\lesssim1''$), polarization measurements, and optical line diagnostics are essential for practical classification. Modeling with free–free and synchrotron emission allows inference of $\eta_e^{\rm rel}$, $p$, and $B$ (Mohan et al., 2023).

6. Limitations and Prospective Developments

Current techniques—including multi-color blackbody (mBB), combination models (mBB+NT), and "top-down" inversions—robustly separate extreme cases (pure fireball vs. hybrid jets). However, sensitivity to magnetization is limited when $\sigma_0 \lesssim 2$, as uncertainties in $r_0$, $L_w$, and binning propagate. Distinction between $\sigma_0 = 0.1$ and $\sigma_0 = 1$ remains imprecise (Song et al., 2023).

Additional diagnostics, notably polarization measurements (to probe Poynting dominance), very-high-energy light curves, thermal precursor components, and sophisticated radiative transfer modeling (accounting for subphotospheric dissipation), are recommended to resolve ambiguities and refine classification. Multi-wavelength afterglow studies and numerical simulations of structured jets represent further pathways for breaking degeneracies and solidifying jet categorization.

In sum, thermal jet classification synthesizes empirical spectral analysis and physical inversion for coherent discrimination of jet types. While proficient at identifying clear cases, its capacity for resolving transitional regimes is limited, motivating future research integrating polarimetric, temporal, and radiative-transfer data (Song et al., 2023, Mohan et al., 2023).