1LHAASO J0343+5254u: Galactic PeVatron Candidate

• 1LHAASO J0343+5254u is an ultra–high-energy gamma-ray source in the Galactic plane proposed as a Galactic PeVatron candidate due to its >100 TeV photon emission.
• Multi-wavelength studies using XMM-Newton, LOFAR, and near-infrared data reveal conflicting interpretations of an associated extended X-ray source as either a pulsar wind nebula or a merging galaxy cluster.
• Ongoing observational efforts with hard X-ray, infrared spectroscopy, and TeV–PeV monitoring are crucial to resolve the counterpart ambiguity and identify the true cosmic accelerator.

1LHAASO J0343+5254u is an ultra–high-energy (UHE) $\gamma$-ray source in the Galactic plane detected by the Large High Altitude Air Shower Observatory (LHAASO), and is treated in the recent literature as a Galactic PeVatron candidate because LHAASO reports photons above $100\ \mathrm{TeV}$. Its astrophysical interpretation has been shaped by two linked developments: first, the identification of an extended X-ray source within its localization region and the proposal that this source is a pulsar wind nebula (PWN) candidate; second, a multi-wavelength re-analysis arguing that the same X-ray source is more plausibly a massive, merging galaxy cluster seen through heavy Galactic extinction and therefore unrelated to the UHE emission. As a result, the $\gamma$-ray source itself remains without a secure lower-energy counterpart (DiKerby et al., 21 Feb 2025, Edler et al., 3 Mar 2026).

1. Discovery history and catalog context

The immediate historical precursor of 1LHAASO J0343+5254u was the extended source LHAASO J0341+5258, discovered by LHAASO’s KM2A array as a new unidentified $\gamma$-ray source in the Galactic plane with emission up to $\sim 200\ \mathrm{TeV}$. In that earlier analysis, the best-fit position was ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$ and ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$, with an intrinsic Gaussian extension $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$. The sky map above $25\ \mathrm{TeV}$ yielded a pre-trial significance of $8.2\sigma$, reduced to $100\ \mathrm{TeV}$0 under a conservative trial-factor estimate. A single extended Gaussian provided an adequate morphological description at that stage, and the source spectrum between $100\ \mathrm{TeV}$1 and $100\ \mathrm{TeV}$2 was consistent with a power law of photon index $100\ \mathrm{TeV}$3, with integral flux above $100\ \mathrm{TeV}$4 reported as $100\ \mathrm{TeV}$5, about $100\ \mathrm{TeV}$6 of the Crab Nebula flux above the same threshold (Collaboration, 2021).

Subsequent catalog work re-analyzed the region with more data and improved sensitivity and decomposed the original field into multiple components. Within that later framework, 1LHAASO J0343+5254u corresponds to the main UHE component in the region, while the broader area also contains 1LHAASO J0339+5307. The same literature states that 1LHAASO J0343+5254u “also shows UHE emission,” and gives a catalog value $100\ \mathrm{TeV}$7, which is the specific basis for its status as a Galactic PeVatron candidate. The source lies near $100\ \mathrm{TeV}$8, $100\ \mathrm{TeV}$9, corresponding to Galactic coordinates $\gamma$0, $\gamma$1, in the outer Galaxy (DiKerby et al., 21 Feb 2025).

The 2026 multi-wavelength study did not re-analyse the LHAASO data itself. Instead, it treated 1LHAASO J0343+5254u as the PeV/TeV $\gamma$2-ray context within which a putative X-ray counterpart had been claimed, and then reassessed that counterpart using radio, X-ray, and near-infrared observations (Edler et al., 3 Mar 2026).

2. The proposed X-ray counterpart and the PWN interpretation

A major step in the source’s interpretation came from deep XMM-Newton observations in February 2024, which reported a previously unknown extended X-ray source inside the LHAASO region: XMMU J034124.2+525720. The X-ray centroid was measured at $\gamma$3, $\gamma$4, and the main extraction region was a circle of radius $\gamma$5. Surface-brightness profiles showed that the emission is intrinsically extended rather than PSF-broadened, reaching $\gamma$6 to the south and dropping more steeply to the north. The source was therefore described as a slightly asymmetric nebula, extended to $\gamma$7, elongated roughly southward (DiKerby et al., 21 Feb 2025).

Its global X-ray spectrum was fitted with an absorbed power law. For the combined exposure of $\gamma$8, the best-fit parameters were $\gamma$9, photon index $\gamma$0, and absorbed flux $\gamma$1 in $\gamma$2–$\gamma$3, with $\gamma$4 for $\gamma$5 degrees of freedom. No variability was found across the individual observations. Spatially resolved spectroscopy further showed radial softening: $\gamma$6 in the core, $\gamma$7 in the inner annulus, and $\gamma$8 in the outer annulus. No X-ray pulsations were detected from the innermost $\gamma$9 region, although the PN mode was not optimized for timing (DiKerby et al., 21 Feb 2025).

On that basis, the source was interpreted as a candidate PWN associated with 1LHAASO J0343+5254u. The argument rested on three features emphasized in the original report: extended non-thermal X-ray emission, a hard featureless spectrum, and significant spectral softening with radius. The source was also presented as morphologically analogous to other X-ray PWNe associated with TeV sources, with the X-ray nebula more compact than the associated TeV emission and offset within the $\sim 200\ \mathrm{TeV}$0-ray contours. Under that interpretation, the absence of a detected pulsar was not taken as decisive, because some TeV PWNe are powered by faint or unidentified pulsars (DiKerby et al., 21 Feb 2025).

3. Radio, X-ray, and NIR re-interpretation as an obscured merging cluster

A later multi-wavelength study revisited XMMU J034124.2+525720 using new LOFAR continuum radio imaging at $\sim 200\ \mathrm{TeV}$1 and $\sim 200\ \mathrm{TeV}$2, alternative X-ray spectral modeling, and archival near-infrared data. That analysis identified several radio structures whose morphologies and spectra were described as suggestive of a radio halo, a radio relic, and tailed radio galaxies, all of which are typically found in merging galaxy clusters rather than PWNe (Edler et al., 3 Mar 2026).

The most prominent radio feature is an arc-like northern source, designated the Arc, lying north of the X-ray emission and aligned with the major axis of the X-ray structure. Its integrated flux densities were measured as $\sim 200\ \mathrm{TeV}$3 and $\sim 200\ \mathrm{TeV}$4, with spectral index $\sim 200\ \mathrm{TeV}$5 for $\sim 200\ \mathrm{TeV}$6. Two tailed sources were also identified: Tail N with $\sim 200\ \mathrm{TeV}$7, $\sim 200\ \mathrm{TeV}$8, and $\sim 200\ \mathrm{TeV}$9; and Tail S with ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$0, ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$1, and ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$2. After subtraction of compact sources and lower-resolution imaging, faint diffuse central emission coincident with the X-ray contours was detected with ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$3, ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$4, and ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$5. The same study also searched for circularly polarized radio emission as a pulsar diagnostic and found no Stokes ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$6 emission above ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$7 (Edler et al., 3 Mar 2026).

The X-ray spectrum was then re-modeled with explicit background treatment. A single absorbed power law again gave an acceptable fit with ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$8, ${\rm R.A.} = 55.34^\circ \pm 0.11^\circ$9, and ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$0. However, a thermal APEC model for the intracluster medium (ICM) yielded a slightly better fit. With metallicity fixed to ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$1, the best-fit parameters were ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$2 and ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$3, with ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$4. Allowing the metallicity to vary produced ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$5, ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$6, ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$7, and ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$8. The reported preference for the thermal interpretation over the power law was ${\rm Dec.} = 52.97^\circ \pm 0.07^\circ$9, corresponding to $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$0 in a likelihood-ratio test (Edler et al., 3 Mar 2026).

Near-infrared analysis using the UKIDSS Galactic Plane Survey strengthened the cluster interpretation. In a target region of radius $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$1 around the X-ray peak, the study found a $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$2 overdensity of cataloged galaxies relative to a local background annulus. Restricting to red galaxies with extinction-corrected $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$3 color in $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$4 increased the significance to $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$5: the expected number was $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$6, while $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$7 were observed. Candidate host galaxies for the tailed radio sources were among these red near-infrared sources. On this basis, the study favored interpreting XMMU J034124.2+525720 as a massive, merging galaxy cluster in a highly extinct region of the Galactic plane (Edler et al., 3 Mar 2026).

4. Spatial coincidence and physical disassociation

The counterpart controversy is driven by projected coincidence. The XMM-Newton pointing discussed in the multi-wavelength re-analysis was centered at Galactic coordinates $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$8 and covered both 1LHAASO J0339+5307 and 1LHAASO J0343+5254u. The extended X-ray source XMMU J034124.2+525720 was described there as “positionally coincident with the lower-energy component of 1LHAASO J0343+5254u,” and the LOFAR images showed it lying within the plotted LHAASO $$\sigma_{\rm ext} = 0.29^\circ \pm 0.06^\circ_{\rm stat} \pm 0.02^\circ_{\rm sys}$$9 containment region (Edler et al., 3 Mar 2026).

The later study nonetheless argued that the X-ray/radio/NIR source is probably unrelated to the UHE $25\ \mathrm{TeV}$0-ray emitter. The central reason is that, if the thermal ICM interpretation is correct, the object is extragalactic at $25\ \mathrm{TeV}$1–$25\ \mathrm{TeV}$2. The same work explicitly noted that extragalactic PeV $25\ \mathrm{TeV}$3 rays at such redshift would be strongly attenuated by pair production on background photon fields, and further pointed out that LHAASO has not detected even the nearest rich galaxy clusters. Within that framework, the galaxy cluster is taken to be an intervening or background object that happens to lie in projection within the LHAASO localization region, whereas 1LHAASO J0343+5254u remains likely associated with a yet-to-be-discovered Galactic source (Edler et al., 3 Mar 2026).

This has become the main objective correction to the earlier PWN-based interpretation. The X-ray source is not established as the counterpart of 1LHAASO J0343+5254u; rather, the current literature presents two incompatible identifications for the same extended X-ray object, with the more recent study favoring the cluster explanation and explicitly concluding that the counterpart of the $25\ \mathrm{TeV}$4-ray source remains unknown (DiKerby et al., 21 Feb 2025, Edler et al., 3 Mar 2026).

5. Emission scenarios advanced for the region

Before the cluster reinterpretation, the PWN hypothesis was developed quantitatively as a leptonic model for the X-ray and $25\ \mathrm{TeV}$5-ray emission. Using XMM-Newton data, LHAASO $25\ \mathrm{TeV}$6-ray measurements, Fermi-LAT upper limits at the X-ray source position, radio non-detections, and Swift-BAT hard-X-ray limits, the 2025 study modeled a relativistic electron population with

$25\ \mathrm{TeV}$7

taking $25\ \mathrm{TeV}$8, $25\ \mathrm{TeV}$9, and $8.2\sigma$0. The best-fit model gave $8.2\sigma$1, $8.2\sigma$2, magnetic field $8.2\sigma$3, IR scaling factor $8.2\sigma$4, normalization $8.2\sigma$5, and total electron energy $8.2\sigma$6. A crucial feature of that model was that solar-neighborhood radiation fields were insufficient: with $8.2\sigma$7, inverse-Compton emission could account for only $8.2\sigma$8 of the observed LHAASO flux, whereas elevated IR photon fields made a fully leptonic explanation viable (DiKerby et al., 21 Feb 2025).

The same paper also stressed environmental ambiguity. Fermi-LAT data showed a strong compact GeV source, 4FGL J0340.4+5302, offset from the X-ray nebula and consistent with a point source, while no significant GeV emission was detected at the X-ray PWN position. Under an assumed $8.2\sigma$9 power law at that position, the reported upper limit was

$100\ \mathrm{TeV}$00

Nobeyama CO observations further revealed five molecular clouds within the LHAASO region: four at $100\ \mathrm{TeV}$01–$100\ \mathrm{TeV}$02 with $100\ \mathrm{TeV}$03–$100\ \mathrm{TeV}$04, and one associated with an AGB star at $100\ \mathrm{TeV}$05 with $100\ \mathrm{TeV}$06. These were identified as potential hadronic targets if relativistic protons are present (DiKerby et al., 21 Feb 2025).

For the precursor source LHAASO J0341+5258, the discovery paper had already shown that both leptonic and hadronic phenomenologies were plausible. In that earlier modeling, a leptonic interpretation required electrons with cutoff energy $100\ \mathrm{TeV}$07 and magnetic field $100\ \mathrm{TeV}$08 to satisfy the diffuse X-ray upper limit, with total electron energy above $100\ \mathrm{TeV}$09 of order $100\ \mathrm{TeV}$10. A hadronic interpretation used a proton spectrum with index $100\ \mathrm{TeV}$11 and cutoff $100\ \mathrm{TeV}$12, giving

$100\ \mathrm{TeV}$13

Because that analysis was performed before the later decomposition of the field, these parameterizations apply to the broader unresolved region rather than uniquely to 1LHAASO J0343+5254u, but they remain relevant as a template for the kinds of leptonic–hadronic degeneracy present in this sky area (Collaboration, 2021).

A plausible implication is that the source environment permits multiple Galactic interpretations even after the X-ray cluster reinterpretation: a hidden pulsar or TeV halo, a supernova-remnant or superbubble system, or cosmic-ray interactions with nearby molecular material remain compatible with the source’s unresolved counterpart status.

6. Current status and observational priorities

The current state of the literature is therefore bifurcated but not symmetric. The earlier XMM-Newton work presented XMMU J034124.2+525720 as a candidate PWN associated with the Galactic PeVatron 1LHAASO J0343+5254u and showed that a fully leptonic model with elevated IR photon fields can reproduce the observed keV–PeV spectral energy distribution (DiKerby et al., 21 Feb 2025). The later multi-wavelength study argued instead that the same X-ray source is a massive, merging galaxy cluster at $100\ \mathrm{TeV}$14, located behind strong extinction and unrelated to the LHAASO emission, and therefore that the counterpart of 1LHAASO J0343+5254u remains unidentified (Edler et al., 3 Mar 2026).

Several observational tests follow directly from that disagreement. The cluster study stated that future observations in the hard X-ray regime will allow certainty about the nature of XMMU J034124.2+525720 because the thermal and power-law models diverge strongly at higher energies; it also proposed near-infrared spectroscopy of the red galaxies in the vicinity and deep high-resolution X-ray spectroscopy to detect ICM emission lines and refine redshift and metallicity (Edler et al., 3 Mar 2026). The PWN study, by contrast, emphasized NuSTAR for locating the synchrotron cutoff, Chandra for resolving a compact pulsar from the nebula, deep radio pulsar searches, improved CO and HI mapping, and continued TeV–PeV spectroscopy to constrain the cutoff and spatial substructure (DiKerby et al., 21 Feb 2025).

Accordingly, 1LHAASO J0343+5254u is best described at present as a UHE Galactic-plane $100\ \mathrm{TeV}$15-ray source whose $100\ \mathrm{TeV}$16-ray existence is secure, whose projected X-ray field contains a well-studied extended source, but whose true counterpart is unsettled. The strongest recent result is not a positive identification but a negative one: the extended X-ray source once proposed as a PWN is now argued to be an obscured merging galaxy cluster, which, if confirmed, removes it from consideration as the source of the PeV emission and re-opens the search for the Galactic accelerator responsible for 1LHAASO J0343+5254u (Edler et al., 3 Mar 2026).