HESS J1857+026: Extended Galactic TeV Source

• HESS J1857+026 is an extended Galactic TeV gamma-ray source, initially unidentified and now linked to an evolved pulsar wind nebula powered by PSR J1856+0245.
• Its broadband spectrum spans GeV to multi-TeV energies, exhibiting hard emission with power-law indices and a spectral turnover near 100 GeV, as confirmed by multiple instruments.
• The source shows energy-dependent, asymmetric morphology and suppressed particle diffusion, indicating complex PWN evolution and environmental interactions in the Galactic plane.

Searching arXiv for papers on HESS J1857+026 to ground the article in the published literature. HESS J1857+026 is an extended Galactic very-high-energy $\gamma$-ray source discovered in the H.E.S.S. Galactic plane survey and long regarded as a prototypical ambiguous TeV emitter: it is bright, morphologically complex, and located near the energetic pulsar PSR J1856+0245, yet for much of its observational history it lacked a confirmed X-ray nebula, supernova-remnant shell, or other secure lower-energy counterpart. Subsequent H.E.S.S., Fermi-LAT, MAGIC, VERITAS, radio, and multiwavelength studies have progressively reframed it from an unidentified TeV source into a system usually interpreted as being dominated by an evolved pulsar wind nebula (PWN), while retaining unresolved questions about source multiplicity, environmental structure, and a possible additional hadronic component at lower $\gamma$-ray energies (Collaboration et al., 2014, Chen, 19 Dec 2025, Eagle et al., 9 Jun 2026).

1. Discovery, localization, and astrophysical setting

HESS J1857+026 was discovered by H.E.S.S. as part of its Galactic plane survey as an extended TeV source without a clear lower-energy counterpart. In the original H.E.S.S. analysis summarized by later work, the source had an intrinsic extension of $$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$, an inclination of $(-3 \pm 49_{\rm stat})^\circ$ counter-clockwise from the RA axis, and already showed a “tail-like” feature extending northward. That northern feature immediately raised the possibility that the source was either strongly asymmetric or already blended with an additional weaker TeV component (Collaboration et al., 2014).

The nearby radio pulsar PSR J1856+0245 rapidly became central to interpretation of the source. Its measured parameters include a spin period of about $81$ ms, a characteristic age of $21$ kyr, and a spin-down luminosity of $4.6\times10^{36}\ {\rm erg\,s^{-1}}$, making it energetically capable of powering a TeV nebula. Early distance estimates from dispersion measure placed it at roughly $9$ kpc, but that estimate was explicitly described as uncertain by factors of $2\!-\!3$; later discussions also cited $\sim 6.3$ kpc from the YMW16 electron-density model, while radio environmental work argued that $\gamma$0 kpc is preferred if the pulsar is physically associated with the surrounding H I cavity-like structure. This distance uncertainty propagates directly into luminosity, size, and efficiency estimates, and remains one of the source’s basic systematic uncertainties (Rousseau et al., 2012, Petriella et al., 2021).

From the outset, HESS J1857+026 was therefore defined by a combination of favorable PWN energetics and poor multiwavelength identifiability. That duality explains why the source has remained a recurrent test case for the distinction between relic PWNe, composite $\gamma$1-ray regions, and larger-scale environmental accelerators.

2. Broadband $\gamma$2-ray spectrum

The source is now measured across the GeV-to-multi-TeV domain, but the spectral interpretation has evolved instrument by instrument. An early Fermi-LAT analysis of 36 months of data found no $\gamma$3-ray pulsations from PSR J1856+0245, but did detect significant GeV emission coincident with HESS J1857+026. That GeV source was fitted by a simple power law with spectral index $\gamma$4, integrated energy flux $\gamma$5, and luminosity $\gamma$6, implying a $\gamma$7-ray efficiency of $\gamma$8 for $\gamma$9 (Rousseau et al., 2012).

MAGIC then extended the VHE spectrum down to $$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$0 GeV and measured one of the source’s canonical spectral results: $$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$1 with

$$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$2

$$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$3

over $$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$4 GeV to $$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$5 TeV, with fit quality $$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$6. That measurement was explicitly described as bridging the gap between the GeV emission measured by Fermi-LAT and the multi-TeV emission measured by H.E.S.S., while providing continuous coverage of the spectral turnover close to $$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$7 GeV (Collaboration et al., 2014).

Representative measurements are summarized below (Rousseau et al., 2012, Collaboration et al., 2014, Chen, 19 Dec 2025, Eagle et al., 9 Jun 2026).

Instrument/band	Principal spectral result	Citation set
Fermi-LAT, $$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$8–$$(0.11 \pm 0.08_{\rm stat})^\circ \times (0.08 \pm 0.03_{\rm stat})^\circ$$9 GeV	Power law, $(-3 \pm 49_{\rm stat})^\circ$0; $(-3 \pm 49_{\rm stat})^\circ$1	Fermi-LAT 2012
H.E.S.S., $(-3 \pm 49_{\rm stat})^\circ$2–$(-3 \pm 49_{\rm stat})^\circ$3 TeV	Single power law, photon index $(-3 \pm 49_{\rm stat})^\circ$4	H.E.S.S. result as summarized later
MAGIC, $(-3 \pm 49_{\rm stat})^\circ$5 GeV–$(-3 \pm 49_{\rm stat})^\circ$6 TeV	Power law, $(-3 \pm 49_{\rm stat})^\circ$7, $(-3 \pm 49_{\rm stat})^\circ$8	MAGIC 2014
VERITAS, $(-3 \pm 49_{\rm stat})^\circ$9–$81$0 TeV	Power law, $81$1, $81$2; upper limit beyond $81$3 TeV	VERITAS 2025
HAWC synthesis, $81$4–$81$5 TeV	Power law with exponential cutoff, $81$6 TeV for Gaussian spatial model	Multi-instrument 2026

Later VERITAS work confirmed that the source is well described by a simple power law below $81$7 TeV, with $81$8 and $81$9, and set an upper limit beyond $21$0 TeV. A subsequent multi-instrument synthesis incorporating HAWC found clear spectral curvature at higher energies and favored a power law with exponential cutoff, with $21$1 near $21$2–$21$3 TeV. Taken together, these results imply that the source is hard in the GeV band, remains bright through the TeV regime, and is not adequately characterized by a single instrument-dependent shorthand such as “featureless TeV power law” once the highest-energy data are included.

3. Morphology and its energy dependence

Morphology is the source’s defining complication. MAGIC’s detailed analysis separated the field into two estimated-energy ranges and showed that the sky map is qualitatively energy dependent. In $21$4, HESS J1857+026 appears as a single extended emitter with centroid at

$21$5

$21$6

and intrinsic extension

$21$7

Above $21$8 TeV, however, the emission resolves into two statistically significant spatially distinct components, named MAGIC J1857.2+0263 and MAGIC J1857.6+0297, with significances of $21$9 and $4.6\times10^{36}\ {\rm erg\,s^{-1}}$0, respectively. The southwestern component associated with the original source core is elongated, with intrinsic extensions $4.6\times10^{36}\ {\rm erg\,s^{-1}}$1 and $4.6\times10^{36}\ {\rm erg\,s^{-1}}$2 along the major and minor axes and major-axis inclination $4.6\times10^{36}\ {\rm erg\,s^{-1}}$3 counter-clockwise from the RA axis, while the northern component was fitted by a circular 2D Gaussian and reported as compatible with a point source (Collaboration et al., 2014).

This result was already foreshadowed in an earlier MAGIC conference analysis, which found one broad component in $4.6\times10^{36}\ {\rm erg\,s^{-1}}$4 with intrinsic Gaussian width $4.6\times10^{36}\ {\rm erg\,s^{-1}}$5, but two significant structures above $4.6\times10^{36}\ {\rm erg\,s^{-1}}$6 TeV. That early interpretation emphasized a tension with a simple one-zone PWN picture, because the northern high-energy tail lacked a low-energy counterpart (Klepser et al., 2011).

The GeV morphology was also revised substantially with longer Fermi-LAT exposure. A 2023 analysis using more than 13 years of LAT data reported that the GeV emission in the region is composed of two extended components. The hard component, Src T, has index $4.6\times10^{36}\ {\rm erg\,s^{-1}}$7 in $4.6\times10^{36}\ {\rm erg\,s^{-1}}$8–$4.6\times10^{36}\ {\rm erg\,s^{-1}}$9 GeV, is spatially coincident with HESS J1857+026, and has a 68% containment radius that decreases from $9$0 below $9$1 GeV to $9$2 above $9$3 GeV. The same work identified two softer GeV structures in the broader field, Src A and 4FGL J1857.9+0313e, each extended and spatially associated with molecular clumps rather than with the hard TeV core (Guo et al., 2023).

VERITAS confirmed the extended nature of HESS J1857+026 but adopted a deliberately simpler baseline description. In a symmetric-Gaussian fit, the centroid was

$9$4

with Gaussian extension

$9$5

Its significance maps in $9$6–$9$7 TeV and $9$8–$9$9 TeV again showed that PSR J1856+0245 is offset from the center of the VHE emission and that a northern component becomes visible at higher energies, but the authors explicitly left open whether that northern structure is a distinct source or instead part of the same source with energy-dependent extent (Chen, 19 Dec 2025).

The combined morphological literature therefore supports two statements simultaneously. First, HESS J1857+026 is unquestionably extended. Second, its spatial structure depends on energy and analysis choice. This suggests that any single-Gaussian representation is an effective summary rather than a complete physical description.

4. PSR J1856+0245 and the relic-PWN interpretation

The association with PSR J1856+0245 remains the strongest engine-based interpretation for at least part of the emission. The pulsar’s spin-down luminosity of $2\!-\!3$0 is sufficient to power the observed $2\!-\!3$1-ray output, and the absence of detected $2\!-\!3$2-ray pulsations in Fermi-LAT data made a dominantly nebular interpretation more plausible than a pulsar-dominated magnetospheric one. Broadband modeling of the GeV–TeV spectral energy distribution in the 2012 Fermi-LAT study found that a low-magnetic-field leptonic nebula was favored, with a relativistic Maxwellian plus power-law tail electron distribution fitting the data significantly better than a simple power-law injection. In that modeling, the best-fit final magnetic field was $2\!-\!3$3, the age was $2\!-\!3$4 kyr, and the hard LAT spectrum plus TeV data made it difficult to match the SED with a simple single power-law electron population (Rousseau et al., 2012).

X-ray observations have constrained that PWN interpretation without decisively confirming it. XMM-Newton detected the pulsar as a hard point source with unabsorbed $2\!-\!3$5–$2\!-\!3$6 keV flux $2\!-\!3$7, corresponding to $2\!-\!3$8 at $2\!-\!3$9 kpc, but Chandra found no evidence for a surrounding compact nebula. For an annulus from $\sim 6.3$0 to $\sim 6.3$1 around the pulsar, the compact-nebula upper limit corresponds to $\sim 6.3$2 at $\sim 6.3$3 kpc. The same study concluded that PSR J1856+0245 remains the most viable association for HESS J1857+026 despite the lack of X-ray morphological confirmation (Nice et al., 2013).

MAGIC sharpened the relic-PWN picture for the southwestern TeV component. MAGIC J1857.2+0263 lies close to PSR J1856+0245, is elongated, and when isolated has a size above $\sim 6.3$4 TeV that agrees well with the original H.E.S.S. extension. The pulsar lies near one end of the major axis rather than at the TeV centroid. Assuming $\sim 6.3$5 kpc, the observed elongation corresponds to a physical scale of about $\sim 6.3$6, whereas proper motion over $\sim 6.3$7 kyr with a typical transverse speed of $\sim 6.3$8–$\sim 6.3$9 could at most account for an offset of $\gamma$00. On that basis, the authors suggested that evolutionary effects such as interaction of the PWN with the SNR reverse shock are likely important in shaping the elongated and offset morphology, and explicitly interpreted MAGIC J1857.2+0263 as the relic PWN of PSR J1856+0245 (Collaboration et al., 2014).

Later VERITAS work reinforced the same qualitative picture. The pulsar is clearly offset from the center of the VHE emission, there is still no confirmed X-ray or shell-type counterpart, and the morphology can be interpreted in terms of relic or asymmetric PWN evolution rather than a compact nebula centered on the pulsar. At the same time, that work emphasized that the source remains “mysterious” precisely because the PWN interpretation is strong but not yet uniquely confirmed by lower-energy imaging (Chen, 19 Dec 2025).

5. Environment, alternative counterparts, and competing large-scale interpretations

The principal challenge to a single-source PWN interpretation comes from the surrounding interstellar medium. MAGIC argued that the northern TeV component, MAGIC J1857.6+0297, is likely unrelated to the pulsar/PWN system and may instead be associated with a foreground molecular-cloud/H II-region complex. In archival VLA Galactic Plane Survey 21 cm continuum and H I data, Galactic Ring Survey $\gamma$01 data, and Spitzer GLIMPSE $\gamma$02 images, the relevant gas structures were found mainly at local-standard-of-rest velocities near $\gamma$03, corresponding to kinematic distances of roughly $\gamma$04–$\gamma$05 kpc. The ultra-compact H II region U36.40+0.02 has $\gamma$06 and near kinematic distance $\gamma$07 kpc; a nearby cloud complex includes G036.59-00.06 and G036.74-00.16. In $\gamma$08 channel maps, an incomplete shell-like feature in the velocity range $\gamma$09–$\gamma$10 coincides with MAGIC J1857.6+0297, and in a single channel at $\gamma$11 the feature resembles a possible cavity or wind-blown bubble with inferred radius $\gamma$12. Because the northern TeV peak does not coincide with the densest $\gamma$13 peaks, MAGIC regarded a straightforward hadronic target-cloud interpretation as disfavored and slightly favored a leptonic scenario in which electrons upscatter infrared-to-UV photons from U36.40+0.02 or related stellar sources, while explicitly not excluding a hadronic origin (Collaboration et al., 2014).

A contrasting environmental interpretation emerged from deep radio work. New VLA observations at $\gamma$14 GHz and $\gamma$15 GHz found no diffuse radio counterpart to HESS J1857+026, no radio PWN around PSR J1856+0245, and no shell-type SNR. The $\gamma$16 GHz image showed no evidence of associated emission above the noise level of $\gamma$17, and the $\gamma$18 GHz image likewise showed no extended radio emission from a PWN powered by the pulsar. Using VGPS H I and CO survey data, that study instead identified a cavity-like H I structure in the velocity range $\gamma$19–$\gamma$20, centered at $\gamma$21, $\gamma$22, with semi-axes $\gamma$23, effective radius $\gamma$24 pc at an adopted distance of $\gamma$25 kpc, expansion velocity $\gamma$26, swept-up H I mass $\gamma$27, kinetic energy $\gamma$28, and dynamical age $\gamma$29 yr. On that basis, the authors argued for a superbubble interpretation and concluded that TeV emission from HESS J1857+026 originates in a superbubble, favoring a single $\gamma$30-ray source rather than a superposition of two distinct sources, while allowing PSR J1856+0245 to contribute as a source of $\gamma$31-rays within the bubble (Petriella et al., 2021).

More recent Fermi-LAT work complicates the picture again rather than resolving it into a single alternative. The 2023 LAT analysis separated a hard component spatially coincident with HESS J1857+026 from softer extended GeV components, Src A and 4FGL J1857.9+0313e, both spatially coincident with molecular clumps in the southwest and northeast. Those soft sources have indices $\gamma$32 and $\gamma$33, respectively, and were interpreted as favoring a hadronic origin, with protons possibly accelerated by a hypothetical SNR associated with PSR J1856+0245. In that framework, the hard GeV/TeV core behaves like a PWN, while at least part of the surrounding GeV emission traces cosmic-ray interactions with ambient gas (Guo et al., 2023).

A 2026 multi-instrument interpretation adopted a compatible but more explicitly hierarchical view. It argued that the dominant MeV–TeV emission mechanism is leptonic inverse-Compton emission from a PWN powered by PSR J1856+0245, but stated that the low-energy component below $\gamma$34 GeV remains ambiguous and could be dominated by hadronic emission originating from a supernova remnant. It also emphasized that the GeV/TeV emission lies in a region of low molecular gas density and lacks a strong spatial match between TeV peaks and CO peaks, weakening a hadronic interpretation for the TeV emission itself (Eagle et al., 9 Jun 2026).

6. Modern synthesis, transport physics, and unresolved issues

The most recent literature converges on an evolved-PWN interpretation above $\gamma$35 GeV, but does so through transport arguments rather than through direct X-ray identification. VERITAS constructed a radial profile centered on the best-fit Gaussian centroid and compared it to the analytic halo-like diffusion form

$\gamma$36

obtaining

$\gamma$37

Assuming distance $\gamma$38 kpc, this corresponds to a physical diffusion length

$\gamma$39

and, using $\gamma$40 with $\gamma$41 kyr, a diffusion coefficient

$\gamma$42

That value was described as about an order of magnitude lower than the Galactic average, though higher than the strongly suppressed values discussed for Geminga-like halos, and therefore indicative of a moderately inhibited-diffusion environment around the pulsar (Chen, 19 Dec 2025).

The 2026 multi-instrument study extended that logic with Fermi-LAT, VERITAS, and HAWC radial surface-brightness profiles. Its time-dependent PWN evolutionary fits yielded system ages between $\gamma$43 kyr and magnetic field strengths between $\gamma$44. In the same framework, the diffusion coefficient at $\gamma$45 TeV was inferred to be $\gamma$46 for $\gamma$47 and $\gamma$48 for $\gamma$49, compared with a Galactic ISM value $\gamma$50. The same study therefore concluded that diffusion is suppressed by roughly $\gamma$51–$\gamma$52 orders of magnitude relative to the average ISM, consistent with other evolved TeV PWNe and TeV halos (Eagle et al., 9 Jun 2026).

What remains unresolved is not the existence of a PWN contribution but the completeness of that interpretation. No confirmed X-ray nebula or shell-type SNR has yet been established; the northern TeV component can still be read either as a second source or as part of a single source with energy-dependent structure; and the low-energy GeV field evidently contains softer extended components that are not simply reducible to the hard PWN-like core. A plausible synthesis is therefore a stratified one: HESS J1857+026 is best understood as a hard, extended, evolved PWN powered by PSR J1856+0245 above $\gamma$53 GeV, embedded in a more complex Galactic-plane environment that may contribute additional hadronic emission and may include large-scale neutral-gas structure such as a superbubble. Higher-resolution radio and especially X-ray observations remain critical for discriminating between these possibilities and for determining whether the northern TeV enhancement is physically distinct or part of the same transport-dominated nebular system (Collaboration et al., 2014, Chen, 19 Dec 2025).