Gaussian Structured Jet

• Gaussian structured jet is defined by an exponential decline in kinetic energy and Lorentz factor with increasing polar angle, offering a continuous, core-to-wing profile.
• It accounts for a wide range of GRB observations, including prompt emission variations and orphan afterglows, by linking angular structure to beaming effects.
• The model’s analytic simplicity aids population synthesis studies, though simulations suggest a more complex two-component structure may be required at larger angles.

A Gaussian structured jet describes a relativistic outflow in which the kinetic energy and initial Lorentz factor decrease smoothly with polar angle from a central axis, following a Gaussian profile. Originating in gamma-ray burst (GRB) research, this parameterization contrasts with traditional “top-hat” (uniform) or power-law jet models by prescribing a continuous, exponential angular structure. The Gaussian jet profile has played a central role in interpreting diverse multiwavelength light curves from GRBs, especially those seen significantly off-axis or lacking a prompt gamma-ray counterpart.

1. Mathematical Formulation and Physical Parameters

The Gaussian jet structure is defined for polar angle $\theta$ from the jet axis by:

$$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$

$$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$

where $\varepsilon(\theta)$ is the kinetic energy per solid angle, $\Gamma(\theta)$ the bulk Lorentz factor, $\varepsilon_0$ and $\Gamma_0$ the on-axis (“core”) values, $\theta_c$ the Gaussian core width, and $\theta_w$ a truncation (“wing”) angle beyond which the jet is assumed to cut off or sharply decline (Cunningham et al., 2020, Xu et al., 2023, Saleem, 2019).

The total (beaming-corrected) jet energy is

$$E_K = 2\pi \int_0^{\theta_w} \varepsilon(\theta)\sin\theta\,d\theta$$

This parameterization translates into the isotropic-equivalent energy as a function of viewing angle $$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$0, $$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$1, and also governs the observed Lorentz beaming and arrival time structure of the GRB afterglow.

2. Jet Structure, Apparent Isotropic Energy, and Viewing Angle Dependence

While the intrinsic Gaussian profile features a sharp exponential drop-off, the apparent structure as seen by an observer at inclination $$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$2 is regulated by both the local emissivity and Doppler boosting. For $$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$3, the observer measures a nearly flat (“core”) isotropic-equivalent energy. For $$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$4, $$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$5 transitions to a steep power-law decline due to the combined effects of off-axis local energy suppression and reduced beaming (Salafia et al., 2015):

$$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$6

with $$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$7. This mapping is critical for interpreting both observed luminosity functions and the diversity of prompt/afterglow behaviors (Salafia et al., 2015, Guo et al., 2020).

3. Multiwavelength Afterglow and Prompt Emission Implications

The dynamical and radiative consequences of a Gaussian jet profile manifest in several observational domains:

• Prompt Emission: For viewing angles within the core ($$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$8), the event appears as a classical, high-$$\varepsilon(\theta) = \varepsilon_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$9 GRB. At larger $$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$0, the exponential suppression produces low-luminosity ($$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$1 erg/s) or “failed” GRBs (Guo et al., 2020). The peak photon energy shifts into X-ray or even lower energies, linking X-ray flashes to structured jet off-axis viewing.
• Afterglow: Off-axis afterglows exhibit slowly rising light curves—often peaking at late times compared to the population norm. The post-peak decay is generally smoother than the abrupt “jet break” of a top-hat, and can show plateaued or even double-peaked features at intermediate angles, especially if the Gaussian wing is energetic (Lamb et al., 2017, Obayashi et al., 2023, Huang et al., 2019).
• Orphan Afterglows: Gaussian jets naturally predict a population of “orphan" afterglows—optical/radio transients without prompt gamma emission—arising from low-$$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$2 wings or far off-axis lines of sight. These are now actively identified in synoptic surveys (Xu et al., 2023).

4. Population Synthesis, GRB Rates, and Multi-messenger Observability

Models incorporating Gaussian structured jets have been extensively applied to the luminosity distributions and detectability of both short and long GRBs, incorporating cosmic star-formation, compact binary merger rates, and delay-time distributions (Guo et al., 2020, Luo et al., 2022, Saleem, 2019). For GW/GRB joint detections, the structure imposes a strong restriction: prompt gamma-ray detection is suppressed beyond $$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$3, sharply reducing the rate of multi-messenger associations at large binary inclinations (Saleem, 2019). For instance, adopting $$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$4 as measured for GW170817, the joint GW+GRB detection rate drops precipitously for $$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$5.

The Gaussian structured jet model addresses the wide luminosity range of observed short GRBs with a single underlying angular structure, providing a natural explanation for events like GRB 170817A (extremely faint, observed at $$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$6 off-axis) (Guo et al., 2020).

5. Comparative Model Assessment, Hydrodynamics, and Physical Origin

Hydrodynamic simulations of jet propagation through stellar (collapsar) or merger ejecta find that the post-breakout angular jet structure often differs from an ideal Gaussian. Specifically, simulations reveal a two-component structure: a flat core up to $$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$7, followed by a jet-cocoon interface exhibiting a power-law decline in energy with angle, $$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$8, with $$\Gamma(\theta) = \Gamma_0 \exp\left[-\frac{\theta^2}{2\theta_c^2}\right] \quad (0 \leq \theta \leq \theta_w)$$9 determined by the degree of Rayleigh–Taylor mixing (Gottlieb et al., 2020). The exponential Gaussian form is not favored at large angles; fitting the simulation output with a Gaussian underpredicts the energy in the intermediate and wing regions by factors of 10–100 once $\varepsilon(\theta)$0, and is excluded at $\varepsilon(\theta)$1 by binned $\varepsilon(\theta)$2 and K–S tests (Gottlieb et al., 2020).

Nevertheless, the Gaussian parameterization remains widely used for phenomenological fitting and population studies due to its analytic simplicity and strong predictive power for the core-dominated regime (Cunningham et al., 2020, Guo et al., 2020, Luo et al., 2022). In some cases, complex angular structures arising from precessing narrow top-hat jets are shown to produce effectively Gaussian envelopes due to superposition over many precessional periods (Huang et al., 2019).

6. Applications to Specific Events and Survey Prospects

The Gaussian structured jet model has successfully explained multi-band afterglow light curves of several prominent GRBs:

• GRB 160625B: Detailed MCMC modeling demonstrates a strong ($\varepsilon(\theta)$3) preference for a Gaussian jet over a top-hat, with a core angle $\varepsilon(\theta)$4 and inferred energies up to $\varepsilon(\theta)$5 erg when allowing for a non-unity electron participation fraction $\varepsilon(\theta)$6 (Cunningham et al., 2020).
• AT2021any: An “orphan" optical afterglow candidate, best explained as a failed GRB viewed on-axis from a Gaussian jet with $\varepsilon(\theta)$7, $\varepsilon(\theta)$8, reflecting high baryon loading and suppressed prompt emission (Xu et al., 2023).
• GRB 221009A: The unprecedented TeV afterglow is well modeled by a Gaussian jet ($\varepsilon(\theta)$9, $\Gamma(\theta)$0 erg), reproducing multi-band light curves and naturally explaining the absence of detectable high-energy neutrinos, even for off-axis viewing (Mondal et al., 17 Nov 2025, Mondal et al., 4 Dec 2025).

Broadband fitting, including wind-driven circumburst media, shows that the Gaussian jet profile remains compatible with observed features when accounting for Klein–Nishina effects and external attenuation (Mondal et al., 4 Dec 2025). Survey simulations indicate that only $\Gamma(\theta)$110% of GRB afterglows exceed next-generation Cherenkov Telescope Array (CTA) sensitivity in a wind environment, predominantly for near-core views and high kinetic energies.

7. Statistical Testing, Population Constraints, and Model Limitations

Global statistical comparisons indicate that a universal Gaussian jet profile cannot simultaneously account for all properties of the observed SGRB population—especially when coupled with a simple Gaussian delay-time distribution. Only when paired with a lognormal delay and, in some cases, an additional cocoon component, does the Gaussian profile marginally pass joint KS goodness-of-fit tests (Luo et al., 2022). The two-component (core + cocoon) or core + power-law wing models are generally favored by hydrodynamics and population synthesis (Gottlieb et al., 2020, Luo et al., 2022, Salafia et al., 2015).

The exponential wings of the Gaussian profile are often too faint at large angles compared to simulation results, which show an extended, energetic tail. For many practical observational analyses—particularly for core-aligned or modestly off-axis events—the Gaussian structured jet is a robust framework, but at large angles or for detailed energetic reconstructions, deviations must be considered (Gottlieb et al., 2020, Salafia et al., 2015).

---

In summary, the Gaussian structured jet provides a powerful, mathematically tractable parameterization for relativistic jet outflows in GRB studies, explaining a wide diversity of afterglow behaviors and luminosity distributions—especially at modest viewing angles from the jet axis. However, hydrodynamic modeling and population-level statistics indicate that real jets may require a more complex, two-component angular structure, especially to account for energy distribution at large angles and implications for multi-messenger detection rates. The Gaussian structured jet remains central to ongoing phenomenological analyses, survey forecasting, and the interpretation of high-cadence, multi-messenger observations (Cunningham et al., 2020, Guo et al., 2020, Mondal et al., 17 Nov 2025, Mondal et al., 4 Dec 2025, Luo et al., 2022, Gottlieb et al., 2020, Salafia et al., 2015, Saleem, 2019, Xu et al., 2023, Obayashi et al., 2023).