---
title: Pulsar Halo Model Overview
url: https://www.emergentmind.com/topics/pulsar-halo-model
type: topic
---

# Pulsar Halo Model Overview

Pulsar halo model denotes a family of astrophysical formulations for describing spatially extended environments associated with pulsars. In the high-energy literature, a pulsar halo is the extended nonthermal region outside the compact pulsar wind nebula (PWN) in which escaped relativistic \(e^\pm\) propagate through the surrounding medium and produce diffuse emission, principally TeV \(\gamma\) rays through inverse Compton scattering; in a separate population-synthesis usage, the term can also refer to canonical radio pulsars that have migrated into the Galactic halo after natal kicks [2207.04011], [1802.04690]. Most current usage concerns TeV halos, because they couple pulsar spin-down, pair injection, near-source transport, multiwavelength radiation, and the local positron problem into a single observable framework.

## 1. Terminology, source class, and physical setting

In the TeV-halo sense, the model distinguishes halos from PWNe both geometrically and dynamically. PWNe are the regions in which pulsar-driven outflows still dominate the local plasma and magnetic environment, whereas pulsar halos are the larger regions in which escaped pairs propagate in the ambient medium and radiate there [2207.04011]. The same review organizes this evolution into a three-stage sequence and places the canonical halo systems in “Stage 3,” where the pulsar is outside, or dynamically decoupled from, the parent supernova remnant and the escaped pairs diffuse in the interstellar medium.

The observational prototypes are Geminga, Monogem / PSR B0656+14, and PSR J0622+3749. Their TeV halos extend to scales of at least \(\sim 20\) pc and are much larger than the X-ray PWNe, which in Geminga are only \(\sim 0.2\)–\(0.3\) pc [2207.04011]. This scale separation is one of the principal empirical arguments that the emission is not simply from the compact nebula. Later work extends the same phenomenology to HESS J1831\(-\)098, which is well described by a pulsar-halo model although a simple 2D Gaussian cannot be rejected, and to PSR B1055\(-\)52, for which H.E.S.S. reports the first detection of extended TeV emission around the pulsar and interprets it as the third detected system of its kind [2510.03183], [2510.02802].

A separate and older population-synthesis use of “pulsar halo model” concerns ordinary, non-recycled radio pulsars in the gaseous circumgalactic medium. In that formulation, the halo is operationally defined as the region outside the Galactic disc and within the virial radius, with YMW16 used to identify the effective edge of the Galaxy along a given line of sight [1802.04690]. The physical motivation is different, but the common theme is migration of pulsars or pulsar-powered particles into environments beyond the compact nebula and thin disc.

## 2. Injection, spin-down evolution, and source prescriptions

TeV-halo models generally treat the pulsar as a continuous point-like source of relativistic \(e^\pm\). For Geminga, one representative source term is written as
\[
Q(E,\mathbi{r},t)=\left\{ \begin{aligned}
& Q_0(E/E_{100})^{-p}\,\delta(\mathbi{r}-\mathbi{r}_s)\,[(t_s+t_{\rm sd})/(t+t_{\rm sd})]^2\,, & t>0 \\
& 0\,, & t<0
\end{aligned} \right.
\]
with \(t_s=342\ {\rm kyr}\), \(t_{\rm sd}=10\ {\rm kyr}\), and \(p=2.24\) in the Geminga analysis of superdiffusion [2101.01438]. The same work emphasizes that \(\sim 100\) TeV electrons cool in about \(10\) kyr, much shorter than the pulsar age, so the current TeV halo is dominated by recently injected particles.

Other formulations retain continuous spin-down injection but replace the phenomenological power law by a spectrum derived from acceleration at the pulsar-wind termination shock. In the DSA-based Geminga model, the halo source term is
\[
Q(E,\mathbf{r},t) = 
\begin{cases}
q(E)\delta(\mathbf{r}-\mathbf{r_{s}})\dfrac{(1 + t/t_{sd})^{-2}}{(1 + t_s/t_{sd})^{-2}}, & t \ge 0 \\
0, & t < 0
\end{cases}
\]
with \(q(E)=p^2f_0(p)\) supplied by the shock solution, \(t_{sd}=10\ \mathrm{kyr}\), \(t_s=342\ \mathrm{kyr}\), and present spin-down power \(L = 3.26 \times 10^{34}\ \mathrm{erg\,s^{-1}}\) [2507.10276]. In this class of models, the injection cutoff is not merely imposed phenomenologically; it is controlled by the finite-upstream DSA solution and by the diffusion coefficient in the acceleration zone.

Many multiwavelength halo calculations adopt a power law with exponential cutoff,
\[
Q(E)=Q_0 E^{-\gamma} \exp(-E/E_c),
\]
normalized by an efficiency \(\eta\) that converts a fraction of the pulsar rotational energy or spin-down power into escaping pairs. Representative values in the literature surveyed here include \(\gamma=1.85\), \(E_c=200\) TeV, and \(\eta=0.12\) for a benchmark Geminga X-ray/TeV model [2512.13846], and \(\gamma=1.4\), \(E_c=200\) TeV, and \(\eta=0.2\) for the LHAASO J0621+3755 / PSR J0622+3749 halo model [2504.02185]. Taken together, these works suggest that continuous spin-down-modulated injection is standard, while the detailed spectral form remains a major source of model dependence.

## 3. Transport prescriptions

The common transport backbone is the diffusion-loss equation. In isotropic form it is written as
\[
\frac{\partial N(E_e, \mathbi{r}, t)}{\partial t}
=
\nabla \cdot[D(E_e)\nabla N(E_e, \mathbi{r}, t)]
+
\frac{\partial[b(E_e)N(E_e, \mathbi{r}, t)]}{\partial E_e}
+
Q(E_e, \mathbi{r}, t),
\]
or equivalently, in the Geminga superdiffusion study,
\[
\frac{\partial N(E, \mathbi{r}, t)}{\partial t}
=
-D(E,\alpha)(-\Delta)^{\alpha/2}N(E,\mathbi{r},t)
+
\frac{\partial[b(E)N(E,\mathbi{r},t)]}{\partial E}
+
Q(E,\mathbi{r},t),
\]
where \(\alpha=2\) recovers normal diffusion and \(\alpha<2\) gives a Lévy-flight-like propagator with heavy tails [2101.01438]. In that application, the HAWC surface-brightness profile of Geminga yields a monotonic increase of \(\chi^2\) as \(\alpha\) decreases from \(2\) to \(1\), and models with \(\alpha<1.32\) or \(<1.40\), depending on the angular range used, are disfavored at \(95\%\) confidence level.

The standard isotropic slow-diffusion picture is often embedded in a two-zone construction. A generic form is
\[
D(E,r)=
\left\{
\begin{array}{ll}
D_0(E/100\,{\rm TeV})^{\delta_0}, & r<r_0 \\
D_{\rm ISM}(E/100\,{\rm TeV})^{\delta_{\rm ISM}}, & r\ge r_0 ,
\end{array}
\right.
\]
or, in source-specific studies, a piecewise constant inner suppressed region plus normal Galactic diffusion outside it [2207.04011]. This formulation is used to reconcile compact TeV morphologies with eventual positron escape to Earth, and it reappears in Monogem modeling, in candidate-halo analyses, and in CTA detectability forecasts [2406.13426], [2310.01016].

Two alternatives are treated explicitly in the literature. The first is isotropic unsuppressed diffusion with a ballistic-to-diffusive transition. In the review formalism this is described with a generalized Jüttner propagator, finite mean free path \(\lambda(E)=3D(E)/c\), and asymptotic profiles \(n(r)=Q/4\pi r^2c\) in the ballistic regime and \(n(r)=Q/4\pi Dr\) in the diffusive regime [2207.04011]. The second is anisotropic diffusion,
\[
\frac{\partial N_e}{\partial t}
=
\frac{1}{r}\frac{\partial}{\partial r}\left(rD_{\perp}\frac{\partial N_e}{\partial r} \right)
+
D_{\parallel}\frac{\partial^2N_e}{\partial z^2}
-
\frac{\partial}{\partial E_e}\left(\dot{E}_eN_e\right)
+
Q(E_e,t)\delta(r)\delta(z),
\]
with \(D_\perp \sim M_A^4 D_\parallel\), so compact projected halos can arise from slow perpendicular transport and favorable viewing geometry rather than from a genuinely small scalar diffusion coefficient [2207.04011].

The resulting controversy is not merely formal. The three frameworks—suppressed isotropic diffusion, ballistic-corrected standard diffusion, and anisotropic diffusion—can all be tuned to reproduce current TeV surface-brightness profiles, but they imply very different local transport coefficients, magnetic geometries, and positron escape properties [2401.17982].

## 4. Radiation, morphology, and representative systems

The emitted TeV \(\gamma\) rays are modeled as inverse-Compton emission from the transported \(e^\pm\), with synchrotron emission from the same particles providing the X-ray counterpart. A useful pair of scaling relations quoted in the review literature is
\[
E_\gamma\approx(E_e/m_ec^2)^2\epsilon=25(E_e/100{\rm TeV})^2\ {\rm TeV},
\]
for IC on the CMB, and
\[
E_{\rm syn}\simeq 2(E_e/100{\rm TeV})^2(B/3\mu {\rm G})\ {\rm keV},
\]
for synchrotron emission. The characteristic halo size is correspondingly estimated as
\[
r_{\rm diff}=2\sqrt{D(E_e)\min[t_{\rm c},t_{\rm age}]},
\]
which makes energy-dependent morphology a central observable [2207.04011].

Current modeling is strongly driven by multiwavelength tension. In Geminga, homogeneous isotropic slow diffusion with \(B=3\,\mu\)G tends to overpredict the X-ray flux by about an order of magnitude, which has motivated either \(B\lesssim 1\,\mu\)G, radial inhomogeneity in the diffusion coefficient, or anisotropic transport [2207.04011]. The dedicated X-ray diagnostic study argues that the three transport models predict distinct X-ray surface-brightness profiles even when their TeV profiles are similar, and identifies Monogem and PSR J0622+3749 as especially promising systems for large-field-of-view X-ray tests [2401.17982].

Several source-specific analyses now provide concrete parameter benchmarks:

| System | Key model result | Source |
|---|---|---|
| Geminga | \(\alpha<1.32\) or \(<1.40\) disfavored at \(95\%\) CL in superdiffusion fits; normal diffusion gives the best fit | [2101.01438] |
| Geminga | \(D_{100}=4.4^{+0.66}_{-0.57}\times10^{27}\ \mathrm{cm^2\,s^{-1}}\); fit requires a rapid increase of \(D(E)\) above \(\sim 100\) TeV | [2507.10276] |
| Monogem / PSR B0656+14 | favored two-zone GALPROP models use \(D_{\rm SDZ}(4\,{\rm GV})=1.3\times10^{26}\ \mathrm{cm^2\,s^{-1}}\), present-day SDZ sizes \((30,50)\) or \((50,70)\) pc, and predict a positron contribution smaller than \(10\%\) | [2406.13426] |
| LHAASO J0621+3755 / PSR J0622+3749 | benchmark parameters \(\gamma=1.4\), \(E_c=200\) TeV, \(D_0=2\times10^{25}\ \mathrm{cm^2\,s^{-1}}\), \(r_b\sim30\) pc, \(B\lesssim1\,\mu\mathrm G\) | [2504.02185] |
| HESS J1831\(-\)098 | halo model is plausible, with \(\Delta AIC \sim 2.1\) relative to a Gaussian\(\times\)power-law model | [2510.03183] |
| PSR B1055\(-\)52 | H.E.S.S. reports \((2.05 \pm 0.32_{\rm stat} \pm 0.72_{\rm sys})^\circ\) extension and \(D_0(100\,{\rm TeV})\) in the range \(0.32\)–\(1.28\times10^{27}\ \mathrm{cm^2\,s^{-1}}\), far below the canonical ISM value | [2510.02802] |

The broad implication is that suppressed diffusion around pulsars is repeatedly inferred, but the spectral shape and the size of the suppression region are not universal. LHAASO J0621+3755 is explicitly argued to require a harder injected spectrum and a smaller suppression zone than Geminga [2504.02185].

## 5. Population synthesis, diffuse emission, and the Galactic-halo usage

Population-level pulsar-halo models extrapolate the source-class picture from nearby prototypes to the Milky Way. One influential application argues that unresolved TeV halos explain the Milagro TeV excess: with a lepton injection spectrum \(\propto E^{-1.7}e^{-E/100\,{\rm TeV}}\) and \(10\%\) conversion of spin-down power into \(e^\pm\) above 1 GeV, the cumulative halo component overtakes hadronic diffuse emission above \(\sim 500\) GeV and exceeds the standard background by about \(\sim 3\) at \(3.5\) TeV and \(\sim 8\) at \(15\) TeV [1707.01905].

A later joint population synthesis of SNRs, PWNe, and halos embeds halos into an explicit PWN-to-halo evolutionary sequence in which halo injection begins when the pulsar exits the nebula. Under the assumption that all middle-aged pulsars in the bow-shock phase develop halos with Geminga-like or Monogem-like transport, the predicted number of detectable halos ranges from \(30\%\) to \(80\%\) of the number of detectable PWNe in H.E.S.S. and HAWC-like surveys, while CTA Galactic Plane Survey forecasts reach about \(170\) PWNe and up to \(100\) halos [2207.11178]. A complementary CTA study, using full spatial-spectral likelihoods, finds that under optimistic Geminga-like assumptions the Galactic Plane Survey could detect about \(300\) halos, identify more than \(70\) through significant energy-dependent morphology, and characterize up to \(30\) in enough detail to study diffusion conditions [2310.01016].

This extrapolation is contested. A phenomenological two-zone study of Geminga and B0656+14 argues that strong suppression by \(2\)–\(3\) orders of magnitude at \(\sim 100\) TeV is required by HAWC and Fermi-LAT data, but that the suppression extent can be as small as \(30\) pc for both objects [2206.11803]. The same paper concludes that if all nearby middle-aged pulsars develop such halos, their combined positron flux, including Geminga, would saturate the \(>100\) GeV AMS-02 measurement for injection efficiencies much smaller than those inferred for the canonical halos. It therefore regards it as simpler to assume that most middle-aged pulsars do not develop halos, with a possible occurrence rate of only \(\sim 5\)–\(10\%\) [2206.11803].

A distinct, non-TeV usage of “pulsar halo model” appears in radio pulsar population synthesis. In that framework, about \(33\%\) of all pulsars beaming toward Earth are in the Galactic halo if magnetic-field and inclination-angle evolution are ignored, but the fraction drops to \(\sim 1.5\%\) if both evolve as a falling exponential [1802.04690]. The same study finds that current surveys are sensitive to only \(\sim 0.5\%\) of halo pulsars beaming toward us in the favorable non-evolving scenario, and that time-domain periodicity searches are more effective than single-pulse or Fourier-domain searches.

## 6. Open problems, observational tests, and computational frameworks

The leading unresolved problem is the origin of the transport environment itself. The review literature discusses self-generated turbulence, Bell instability, SNR-induced turbulence, ballistic-to-diffusive transition physics, and anisotropic diffusion in an ordered field, but does not identify a single accepted mechanism [2207.04011]. This uncertainty propagates directly into the positron problem: in one-zone slow diffusion Geminga contributes very little at Earth, whereas mildly superdiffusive transport or finite-size two-zone models can produce a much larger local positron flux [2101.01438], [2207.04011].

Current observational leverage is limited by morphology and background systematics. In the DSA-based Geminga test, only \(D_{100}\) is used as a hard fit constraint, the HAWC morphology bins are broad, and the strongest future discriminator is a rapidly increasing \(D(E)\) above \(100\) TeV [2507.10276]. In the X-ray band, wide-field imaging has become critical because narrow-field telescopes undersample degree-scale halos. The SRG/ART-XC Geminga study derives a \(68\%\)-confidence upper limit \(B\le 6\,\mu\)G within a benchmark suppressed-diffusion model and finds, through simulation, that a 20-day exposure would yield a \(68\%\) probability of detecting the halo for a \(3\,\mu\)G field [2512.13846].

Large-field-of-view X-ray observations are also proposed as the clearest route to transport-model discrimination. The dedicated diagnostic study concludes that eROSITA and Einstein Probe, with sufficiently long exposures, can distinguish the isotropic suppressed-diffusion, ballistic-corrected, and anisotropic diffusion models through their X-ray surface-brightness profiles, especially for Monogem and PSR J0622+3749 [2401.17982]. This suggests that the next decisive advance may come not from another TeV radial profile alone, but from combined TeV-plus-X-ray morphology.

The computational side of the subject has correspondingly broadened. PHECT, introduced as a lightweight computation tool for pulsar-halo emission, implements `normal`, `normEx`, `normJ`, `superdiff`, `2zoneNumr`, `2zSNR`, and `aniso` transport models, with finite-volume discretizations designed to remain stable on non-uniform grids and in the presence of discontinuous coefficients [2508.13667]. Its existence reflects a general trend: pulsar halo modeling has moved from single-source phenomenology toward a modular numerical discipline in which alternative transport hypotheses can be confronted with increasingly precise \(\gamma\)-ray and X-ray data.

Taken together, these developments define the contemporary pulsar halo model as a multi-regime framework rather than a single equation set. Its stable elements are continuous pulsar-powered \(e^\pm\) injection, radiative cooling by IC and synchrotron, and degree-scale emission outside the compact PWN. Its unsettled elements are the transport operator, the size and physical origin of the slow-diffusion region, the role of magnetic anisotropy, and the population-level commonness of the phenomenon.

Source: https://www.emergentmind.com/topics/pulsar-halo-model