Quasi-Universal Structured Jet Model

• The QUSJ model is a structured jet framework that describes gamma-ray burst jets with a universal core-plus-wing design, where the core produces high luminosity and the wings rapidly decline.
• Analytical models, hydrodynamic simulations, and Bayesian analyses support its use in unifying the observed diversity in long and short GRB afterglow light curves and spectral properties.
• The model attributes variations in observed luminosity to differences in viewing angles rather than intrinsic engine activity, offering a consistent explanation for the GRB luminosity function.

The Quasi-Universal Structured Jet (QUSJ) model is a physically motivated framework describing gamma-ray burst (GRB) jets—incorporating both long (LGRBs) and short (SGRBs) bursts—as a one-parameter family of jets with a universal core-plus-wing angular structure. The diversity of observed isotropic-equivalent luminosities, afterglow light curves, and spectral properties among GRBs is largely attributed to variations in the viewing angle with respect to the jet axis, rather than intrinsic differences in central engine activity or progenitor system characteristics. Key predictions, supporting evidence, and population-level implications for jet structure in both LGRBs and SGRBs have been substantiated by analytical models, hydrodynamic simulations, and hierarchical Bayesian analyses of GRB catalogs and multi-messenger samples including gravitational wave counterparts (Xie et al., 2020, Salafia et al., 2023, Salafia et al., 2019, Salafia et al., 2015, Pracchia et al., 7 Jan 2026).

1. Theoretical Basis and Motivation

The QUSJ model postulates that all GRB jets exhibit an axisymmetric angular structure characterized by a narrowly collimated “core” of quasi-uniform luminosity and Lorentz factor, surrounded by steeply declining “wings” resulting from jet–envelope interaction and jet–cocoon dynamics. For LGRBs, core energy release originates from newly-born millisecond magnetars (surface field $B_p \sim 10^{15}$ G, radius $R \sim 10^6$ cm, initial spin $\Omega_0$), leading to a total rotational energy $E_{\rm rot} \sim 2 \times 10^{52}$ erg. Magnetic dipole spin-down powers an ultra-relativistic jet, while the remainder of spin-down luminosity contributes to a mildly relativistic, long-lived X-ray wind. Jet propagation through stellar or merger ejecta establishes the angular structure: a bright core collimated by confining material, and a cocoon or sheath with a steep angular energy decay (Xie et al., 2020, Salafia et al., 2019, Pracchia et al., 7 Jan 2026).

Observed diversity in GRB luminosities, including the continuum from typical high-luminosity GRBs to low-luminosity events and off-axis counterparts (e.g., GRB 170817A/GW170817), is interpreted as a geometric effect of varying observer viewing angles rather than requiring fundamentally different progenitors or energy budgets (Salafia et al., 2023, Salafia et al., 2015, Salafia et al., 2019).

2. Mathematical Formulation of Jet Structure

The angular structure is parameterized by the isotropic-equivalent luminosity $L(\theta)$ or energy per solid angle $\epsilon(\theta)$. For a power-law profile:

$$L(\theta) = \begin{cases} L_c, & \theta \leq \theta_c \ L_c \left( \frac{\theta}{\theta_c} \right)^{-k}, & \theta > \theta_c \end{cases}$$

where $L_c$ is the core luminosity, $\theta_c$ is the core half-opening angle, and $k$ is the power-law decay index (Xie et al., 2020, Salafia et al., 2023, Salafia et al., 2015). Gaussian profiles and smoothly-broken power-law representations are also employed:

$$\ell(\theta_v) = \left[1 + \left(\frac{\theta_v}{\theta_c}\right)^s \right]^{-\alpha_L/s} \left[1 + \left(\frac{\theta_v}{\theta_w}\right)^s \right]^{-(\beta_L - \alpha_L)/s}$$

with $s$ controlling the smoothness of the break, and $\ell(\theta_v)$ representing the dimensionless angular luminosity profile as a function of viewing angle $\theta_v$ (Salafia et al., 2023, Pracchia et al., 7 Jan 2026).

For SGRBs and LGRBs, typical parameter ranges derived from fits to GRB populations and afterglow observations are:

• Core angle $\theta_c \sim 2^\circ$–$3^\circ$
• Decay index $k$ or $\alpha_L \sim 4$–$5$ (power-law wing)
• On-axis core luminosity $\log_{10} L_c\ (\textrm{erg\,s}^{-1}) \sim 52.5$–$52.7$
• Exponential/Gaussian width $\theta_c \approx 3^\circ$ (Xie et al., 2020, Salafia et al., 2023, Salafia et al., 2015, Salafia et al., 2019, Pracchia et al., 7 Jan 2026)

The distribution of $L_c$ among events follows $p(L_c) \propto L_c^{-A}$ above a minimum scale, typically with $A \sim 2.7$–$3.2$, and a log-normal or power-law cutoff at low luminosity.

3. Physical Origins from Central Engines and Jet–Ejecta Interaction

For magnetar-driven GRBs, strong magnetic dipole spin-down ($L_{\rm sd}(t)$) creates an initial Poynting-flux–dominated outflow, a fraction of which collimates through the fueling and breakout phase to form the “core” of the jet. The remainder contributes to a mildly relativistic cocoon or wind observed as a plateau in X-ray afterglow light curves (Xie et al., 2020). Hydrodynamic modeling and numerical simulations demonstrate that collimation at jet breakout leads to a robust, almost progenitor-independent core+wing structure for both collapsar and neutron-star merger progenitors (Salafia et al., 2019).

Key dynamical dependencies include:

• $E_0$ (core energy) controlled by jet luminosity and breakout duration
• $\theta_c$ set by pressure and collimation balance at breakout
• Wing slope $k$ determined by the duration of the post-breakout engine activity relative to the expansion timescale (Salafia et al., 2019)
• Lorentz factor profile: $$\Gamma(\theta) = \Gamma_0 [1 + (\theta/\theta_c)^2]^{-m}$$ with $m \sim 1/3$–$1/2$ at large $\theta$ (Salafia et al., 2019)

These physical mechanisms provide natural explanations for the uniformity of angular structure across a range of progenitor and engine parameters.

4. Population Synthesis and Luminosity Function Predictions

The QUSJ hypothesis successfully accounts for the observed luminosity function (LF) of both LGRBs and SGRBs, fitting a broken power-law shape spanning $L \sim 10^{46}$–$10^{54}$ erg s$^{-1}$, with the low-luminosity end populated by off-axis (large-$\theta_v$) events. Monte Carlo and Bayesian hierarchical inference frameworks derive best-fit parameter ranges, including viewing angles, jet energies, and core sizes, by comparing simulated and observed distributions (Xie et al., 2020, Salafia et al., 2023, Pracchia et al., 7 Jan 2026).

Joint GW and electromagnetic (EM) analyses utilizing the QUSJ structure predict rates for coincident GW-SGRB detections. For the O4 observing run, the expected joint GW+SGRB detection rates range from 0.2 to 1.3 yr$^{-1}$ depending on sample completeness (Salafia et al., 2023).

Selection effects and flux thresholds are rigorously incorporated in population analyses to avoid biases in jet-structure parameter inference (Pracchia et al., 7 Jan 2026). The mapping from the universal angular structure to the apparent luminosity function is a fundamental component of interpreting the rates and energetics of observed GRBs.

5. Observational Consequences and Correlations

A direct outcome of the QUSJ model is the prediction of a correlation between prompt gamma-ray luminosity and the plateau X-ray wind luminosity (e.g., $L_w \propto L_j^{0.9}$), stemming from their mutual decline with observer viewing angle. The appearance of an intrinsic Yonetoku relation $E_p \propto L^{0.4}$ is a consequence of the co-varying angular profiles of both $L(\theta)$ and $E_p(\theta)$; however, selection effects steepen the observed correlation relative to the intrinsic one (Salafia et al., 2023, Salafia et al., 2015, Pracchia et al., 7 Jan 2026).

In afterglow modeling, the break time and shape of multiwavelength light curves depend critically on viewing angle relative to $\theta_c$, unifying observed diversity within a single structural paradigm (Salafia et al., 2015).

X-ray surveys, such as those anticipated with the Einstein Probe, are predicted to observe a substantial fraction of plateau emissions arising from off-axis (faint) or wind components, further testing the structure model (Xie et al., 2020).

6. Robustness, Limitations, and Future Directions

Assumptions underlying the QUSJ model include minimal intrinsic variability in engine parameters and progenitor profiles, leading to a “quasi-universal” shape. While current models employ narrow log-normal distributions for initial jet luminosities and core angles, broader progenitor diversity is likely to widen the observed distributions but does not erase the core+wing structure (Salafia et al., 2019). Semi-analytic treatments reproduce key hydrodynamic features within $\sim$20% of relativistic simulations, though refinements are needed for highly magnetized outflows and richer ejecta environments (Salafia et al., 2019).

Empirical validation across both LGRB and SGRB populations, the unification of LL-GRBs and classical events, and agreement with GW-inferred merger rates collectively support the QUSJ hypothesis (Xie et al., 2020, Salafia et al., 2023, Pracchia et al., 7 Jan 2026). As the number of multi-messenger (GW+EM) and off-axis detections increases, constraints on wing steepness, cocoon emission, and viewing-angle effects will further sharpen the model's applicability and expose any remaining deviations from universality.

7. Comparative Overview of Jet Structure Models

Model Type	Profile Form	Core Angle ($\theta_c$)	Power-Law Index ($k$/$\alpha_L$)
Uniform (Top-hat)	flat up to $\theta_c$, zero outside	$2^\circ$–$5^\circ$	--
Power-Law Structure	$L(\theta) \propto \theta^{-k}$	$2^\circ$–$3^\circ$	$k \sim 4$–$8$
Gaussian Structure	$L(\theta) \propto \exp[-(\theta/\theta_c)^2]$	$3^\circ$	Apparent $k \sim 3$–$4$
QUSJ (Broken PL)	Core+wing PL, possible wide-angle break	$2^\circ$–$3^\circ$	$\alpha_L \sim 4.7$

Both steep power-law and smooth (e.g., Gaussian) intrinsic profiles yield nearly indistinguishable, broken power-law “apparent” profiles once relativistic beaming and integration over emission patches are included, providing a unified explanation for the observed GRB luminosity functions and afterglow phenomenology (Salafia et al., 2015). The QUSJ model thus subsumes a range of structured jet prescriptions into a comprehensive, observationally anchored framework.