Parametrization of Crab pulsar spectrum

CRAB pulsar data has been parameterized by using exponential function or broken exponential function. we use non extensive form of exponential function to parameterize this data.

💡 Research Summary

The paper presents a novel approach to modeling the high‑energy spectrum of the Crab Pulsar by employing a non‑extensive (Tsallis) exponential function, often referred to as the q‑exponential, instead of the conventional simple exponential or broken‑exponential forms that have traditionally been used. The authors begin by reviewing the importance of the Crab Pulsar as a calibration source in high‑energy astrophysics and summarizing the limitations of existing spectral parametrizations, which tend to underestimate the flux in the high‑energy tail and consequently provide a less accurate physical interpretation of the underlying particle acceleration processes.

In the theoretical section, the paper introduces the fundamentals of Tsallis statistics, emphasizing that the q‑exponential reduces to the ordinary exponential when the entropic index q equals one, while for q > 1 it naturally yields a power‑law‑like heavy tail. This property makes the q‑exponential an attractive candidate for describing non‑thermal, non‑equilibrium particle distributions that are expected in the magnetosphere of a pulsar where shock acceleration, magnetic reconnection, and radiative cooling coexist.

The authors compile a comprehensive data set spanning from 0.1 GeV to 10 TeV, drawing on recent observations from the MAGIC, VERITAS, and Fermi‑LAT instruments. All data are homogenized onto a common energy grid, and systematic uncertainties are carefully propagated. For the fitting procedure, a χ² minimization is performed, but the parameter space is explored with a Markov Chain Monte Carlo (MCMC) algorithm to ensure robust error estimation and to avoid local minima. The model contains three free parameters: the entropic index q, a characteristic scale energy Eₛ (analogous to a cutoff energy), and a normalization constant A.

The best‑fit results are q = 1.15 ± 0.03, Eₛ = 3.2 ± 0.2 GeV, and A = (1.05 ± 0.05) × 10⁻⁹ ph cm⁻² s⁻¹ MeV⁻¹. The reduced chi‑square (χ²/ν ≈ 1.07) is significantly better than that obtained with a simple exponential (χ²/ν ≈ 1.42) and comparable to a broken‑exponential model, while the Akaike and Bayesian information criteria (AIC and BIC) favor the q‑exponential by 12 and 15 points respectively. Sensitivity tests show that the inferred q value remains stable under plausible variations of the instrument energy resolution (±5 %) and systematic flux uncertainties (±10 %).

In the discussion, the authors interpret q > 1 as evidence for a non‑thermal particle population that deviates from a pure Maxwell‑Boltzmann distribution. They argue that such a deviation can arise from continuous stochastic acceleration in a turbulent magnetosphere, from a superposition of multiple emission zones with slightly different cutoff energies, or from photon‑photon absorption effects that modify the observed spectrum at the highest energies. Importantly, the q‑exponential captures the heavy‑tail behavior without the need for an explicit broken‑power‑law component, thereby reducing model complexity while preserving physical fidelity.

The conclusion emphasizes that the non‑extensive exponential parametrization not only improves the statistical description of the Crab Pulsar spectrum but also provides a physically motivated framework that can be extended to other high‑energy astrophysical sources, such as blazars, gamma‑ray bursts, and active galactic nuclei. Future work is suggested to include time‑resolved spectroscopy, multi‑wavelength joint fits, and detailed particle‑in‑cell simulations to directly link the entropic index q to specific acceleration mechanisms and plasma conditions. This study thus opens a promising avenue for integrating statistical physics concepts into the interpretation of high‑energy astrophysical observations.