Search

Radio galaxies have long been considered as potential sources of ultra-high-energy cosmic rays (UHECRs). Recent analyses of the UHECR spectrum, composition, and arrival directions indicate that the nearest radio galaxy, Centaurus A, could be linked to the reported dipole anisotropy, though the mechanism underlying the acceleration remains elusive. In this Letter, we explore UHECR acceleration in the kiloparsec-scale jets of radio galaxies, exemplified by Centaurus A. Using high-resolution relativistic magneto-hydrodynamic and test-particle simulations without sub-grid physics, we investigate the acceleration of the highest-energy particles in the turbulent sheath of a fast-moving jet. Our findings demonstrate that acceleration close to the maximum theoretical expectation is possible. When extrapolated to nearby radio galaxies, our results suggest that the kiloparsec-scale jets of Centaurus A could account for the dipole anisotropy in UHECRs, while more potent Fanaroff-Riley type II radio galaxies may account for the observed UHECR spectrum with a rigidity cutoff at a few Exavolts., Comment: 11 pages, 9 figures, 1 table, 3 appendices, accepted to ApJL

Recent cosmic ray (CR) measurements have revealed unexpected anomalies in secondary CRs, namely deviations from the predictions of the so-called standard Galactic CR paradigm regarding the composition and energy spectra of the products of interactions of primary (accelerated) CRs with interstellar gas: (i) antiparticles (positrons and antiprotons), (ii) light elements of the (Li, Be, B) group, and (iii) diffuse gamma rays. We argue that the new measurements can still be explained within the standard CR paradigm but with an additional assumption that CRs spend a significant part of their lifetime near their formation sites. The latter can be realized if CRs propagate more slowly in these localized regions than in the interstellar medium (ISM). Postulating that CRs accumulate on average energy-independent "grammage" of $0.7 \ \rm g/cm^2$ near the major contributors to galactic CRs, one can explain self-consistently the new measurements of the B/C ratio by DAMPE and the diffuse ultra-high-energy gamma-rays by LHAASO, involving a minimal number of model parameters: the energy-dependent "grammage" in the interstellar medium $\rm \lambda \approx 8 (E/10 \ GeV)^{-0.55}~\rm g/cm^{2}$ and the average CR acceleration (sourcee) spectrum, $\rm Q(E) \propto E^{-2.3}$.

Large-scale magnetic fields are observed off the midplanes of disk galaxies, indicating that they harbour magnetised halos. These halos are crucial to studies of galaxy evolution, galactic-scale outflows, and feedback from star formation activity. Identifying the magnetised halo of the Milky Way is challenging because of the potential contamination from foreground emission arising in local spiral arms. Additionally, it is unclear how our magnetic halo is influenced by recently revealed large-scale structures such as the X-ray emitting eROSITA Bubbles, which, according to previous simulations, might be transient structures powered by the Galactic Center or the Galaxy's star-forming ring. Here we report the identification of several kpc-scale magnetised structures based on their polarized radio emission and their gamma-ray counterparts, which can be interpreted as the radiation of relativistic electrons. These non-thermal structures extend far above and below the Galactic plane and are spatially coincident with the thermal X-ray emission from the eROSITA Bubbles. The morphological consistency of these structures suggests a common origin, which can be sustained by Galactic outflows driven by the active star-forming regions located at 3-5 kpc from the Galactic Centre. These results reveal how X-ray-emitting and magnetised halos of spiral galaxies can be related to intense star formation activities and suggest that the X-shaped coherent magnetic structures observed in their halos can stem from galaxy outflows., Comment: Initially submitted on March 2nd, 2024

Cosmic rays are energetic nuclei that permeate the entire Galactic disk. Their existence requires the presence of powerful particle accelerators. While Galactic supernova explosions may supply the required energy, there is growing evidence that they cannot explain all of the observed properties of cosmic rays, such as their maximum particle energy and isotopic composition. Among Galactic objects, winds from stellar clusters meet the energetic requirement and provide a suitable environment for particle acceleration. The recent detection of some of these objects in gamma rays confirms that they indeed harbor high-energy particles.However, as most supernovae explode inside stellar clusters, it is difficult to distinguish the contribution of winds to particle acceleration. Here we report the detection of young star clusters in the nearby Vela molecular ridge star forming region. The young age of the systems guarantees an unbiased estimate of the stellar CR luminosity free from any supernova or pulsar contamination and allows us to draw conclusions on the acceleration efficiency and the total power supplied by these objects. We demonstrate that much more than 1% of the wind mechanical power is converted into CRs and consequently conclude that a small but non-negligible fraction ~ 1-10% of the CR population is contributed by stellar clusters., Comment: Accepted Manuscript version

HESS J1702-420 is a peculiar TeV complex with a morphology changing from a diffuse (HESS J1702-420B source) at $\lesssim 2$ TeV to point-like (HESS J1702-420A) at $\gtrsim 10$ TeV energies. The morphology and the spectral properties of HESS J1702-420 could be understood in terms of a (diffusive) hadronic or leptonic models in which the observed TeV emission arises correpondingly from proton-proton or IC-radiation of relativistic particles present in the region. In this work we perform searches of the X-ray counterpart of HESS J1702-420B source originated from the synchrotron emission of the primary or secondary relativistic electrons produced within leptonic or hadronic models. Such an emission can be extended and remain beyond the detection capabilities of a narrow-FoV instruments such as XMM-Newton. We utilise the publicly available first 6-months eROSITA dataset (DR1) fully covering selected for the analysis region of $> 5^\circ$-radius around HESS J1702-420. We discuss biases connected to variable plasma temperature/neutral hydrogen column density in the region and present results based on background modelling approach. The performed analysis does not allow us to detect the extended X-ray counterpart of HESS J1702-420 of $0.07^\circ - 3^\circ$-radii sizes. The derived upper limits are significantly higher than the expected hadronic model flux of the X-ray counterpart. For the leptonic model the derived limits indicate the magnetic field in the region $B\lesssim 2\mu$G. We argue, that the further advances in the diffuse X-ray counterpart searches could be achieved either with next generation missions or Msec-long observational campaigns with currently operating instruments., Comment: accepted to A&A

The brightest gamma ray burst (GRB) ever observed, GRB221009A, produced a surprisingly large flux of gamma rays with multi-TeV energies, which are expected to be absorbed in interactions with extragalactic background light (EBL). If the highest energy gamma rays were produced at the source, their spectral shape would have to exhibit a nonphysical spike even for the lowest levels of EBL. We show that, for widely accepted models of EBL, the data can be explained by secondary gamma rays produced in cosmic ray interactions along the line of sight, as long as the extragalactic magnetic fields (EGMFs) are $10^{-16}$G or smaller, assuming 1 Mpc correlation length. Our interpretation supports the widely held expectation that GRB jets can accelerate cosmic rays to energies as high as 10 EeV and above, and it has implications for understanding the magnitudes of EGMFs., Comment: 6 pages, 3 figures

Ultrarelativistic particles called cosmic rays permeate the Milky Way, propagating through the Galactic turbulent magnetic fields. The mechanisms under which these particles increase their energy can be reasonably described by current theories of acceleration and propagation of cosmic rays. There are, however, still many open questions as to how to reach petaelectronvolt (PeV) energies, the maximum energy believed to be attained in our Galaxy, and in which astrophysical sources (dubbed {\it PeVatrons}) this ultra-high energy acceleration happens. In this article, we describe the theoretical conditions for plasma acceleration to these energies, and the Galactic sources in which these conditions are possible. These theoretical predictions are then confronted with the latest experimental results, summarising the state-of-the-art of our current knowledge of PeVatrons. We finally describe the prospects to keep advancing the understanding of these elusive objects, still unidentified more than one hundred years after the discovery of cosmic rays., Comment: Published in Nature Astronomy, Volume 8, pages 425-431 (2024)

HESS J1641-463 is an unidentified gamma-ray source with a hard TeV gamma-ray spectrum, and thus it has been proposed to be a possible candidate for cosmic ray (CR) accelerators up to PeV energies (a PeVatron candidate). The source spatially coincides with the radio supernova remnant (SNR) G338.5+0.1, but has not yet been fully explored in the X-ray band. We analyzed newly taken NuSTAR data, pointing at HESS J1641-463, with 82 ks effective exposure time. There is no apparent X-ray counterpart of HESS J1641-463, while nearby stellar cluster, Mercer 81, and stray-light X-rays are detected. Combined with the archival Chandra data, partially covering the source, we derived an upper limit of $\sim 6\times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in 2-10 keV ($\sim 3\times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in 10-20 keV). If the gamma-ray emission is originated from decay of $\pi^0$ mesons produced in interactions between CR protons and ambient materials, secondary electrons in the proton-proton interactions can potentially emit synchrotron photons in the X-ray band, which can be tested by our X-ray observations. Although the obtained X-ray upper limits cannot place a constraint on the primary proton spectrum, it will be possible with a future hard X-ray mission., Comment: 15 pages, 8 figures, 4 tables, accepted for publication in ApJ

De-excitation gamma-ray lines, produced by nuclei colliding with protons, provide information about astrophysical environments where particles have kinetic energies of $10-100$ MeV per nucleon. In general, such environments can be categorized into two types: the interaction between non-thermal MeV cosmic rays and ambient gas, and the other is thermal plasma with a temperature above a few MeV. In this paper, we focus on the latter type and investigate the production of de-excitation gamma-ray lines in very hot thermal plasma, especially the dependence of the line profile on the plasma temperature. We have calculated the line profile of prompt gamma rays from $^{12}$C and $^{16}$O and found that when nuclei have a higher temperature than protons, gamma-ray line profiles can have a complex shape unique to each nucleus species. This is caused by anisotropic gamma-ray emission in the nucleus rest frame. We propose that the spectroscopy of nuclear de-excitation gamma-ray lines may enable to probe energy distribution in very hot astrophysical plasmas. This diagnostics can be a new and powerful technique to investigate the physical state of a two-temperature accretion flows onto a black hole, especially the energy distributions of the protons and nuclei, which are difficult to access for any other diagnostics., Comment: 12 pages, Accepted for publication in MNRAS

We explore the implications of the light curve of the early TeV gamma-ray afterglow of GRB221009A reported by the LHAASO collaboration. We show that the reported offset of the reference time, $T_*$, allows the determination of the relativistic jet activation time, which occurs approximately $200\,\mathrm{s}$ after the GBM trigger time and closely precedes the moment at which GBM was saturated. We find that while the LHAASO data do not exclude the homogeneous circumburst medium scenario, the progenitor wind scenario looks preferable, finding excellent agreement with the expected size of the stellar bubble. We conclude that the initial growth of the light curve is dominated by processes internal to the jet or by gamma-gamma attenuation on the photons emitted during the prompt phase. Namely, either the activation of the acceleration process or the decrease of internal gamma-gamma absorption can naturally explain the initial rapid flux increase. The subsequent slow flux growth phase observed up to $T_*+18\,\mathrm{s}$ is explained by the build-up of the synchrotron radiation -- the target for inverse Compton scattering, which is also supported by a softer TeV spectrum measured during this period. The duration of this phase allows an almost parameter-independent determination of the jet's initial Lorentz factor, $\Gamma_0\approx600$, and magnetic field strength, $B'\sim0.3\,\mathrm{G}$. These values appear to match well those previously revealed through spectral modeling of the GRB emission., Comment: 20 pages, 6 figures, submitted

HESS J1702-420 is a multi-TeV gamma-ray source with an unusual energy-dependent morphology. The recent H.E.S.S. observations suggest that the emission is well described by a combination of point-like HESS J1702-420A (dominating at highest energies, $\gtrsim$ 30 TeV ) and diffuse ($\sim$ 0.3$^\circ$) HESS J1702-420B (dominating below $\lesssim$ 5TeV) sources with very hard (${\Gamma} \sim 1.5$) and soft (${\Gamma}$ ~2.6) power-law spectra, respectively. Here we propose a model which postulates that the proton accelerator is located at the position of HESS J1702-420A and is embedded into a dense molecular cloud that coincides with HESS J1702-420B. In the proposed model, the VHE radiation of HESS J1702-420 is explained by the pion-decay emission from the continuously injected relativistic protons propagating through a dense cloud. The energy-dependent morphology is defined by the diffusive nature of the low-energy protons propagation, transiting sharply to (quasi) ballistic propagation at higher energies. Adopting strong energy dependence of the diffusion coefficient, $D \propto E^\beta$ with $\beta \geq 1$, we argue that HESS J1702-420 as the system of two gamma-ray sources is the result of the propagation effect. Protons injected by a single accelerator at the rate $Q_0 \simeq 10^{38} \, (n_0/100 \, \rm cm^{-3})^{-1}\, (d/ \, 0.25\,kpc)^{-1} \rm erg/s$ can reasonably reproduce the morphology and fluxes of two gamma-ray components., Comment: accepted to ApJ

Relativistic electrons are an essential component in many astrophysical sources, and their radiation may dominate the high-energy bands. Inverse Compton (IC) emission is the radiation mechanism that plays the most important role in these bands. The basic properties of IC, such as the total and differential cross sections, have long been studied; the properties of the IC emission depend strongly not only on the emitting electron distribution but also on the properties of the target photons. This complicates the phenomenological studies of sources, where target photons are supplied from a broad radiation component. We study the spectral properties of IC emission generated by a power-law distribution of electrons on a power-law distribution of target photons. We approximate the resulting spectrum by a broken-power-law distribution and show that there can be up to three physically motivated spectral breaks. If the target photon spectrum extends to sufficiently low energies, $\varepsilon_{\mathrm{min}}< m_e^2c^4/E_{\mathrm{max}}$ ($m_e$ and $c$ are electron mass and speed of light, respectively; $\varepsilon_{\mathrm{min}}$ and $E_{\mathrm{max}}$ are the minimum/maximum energies of target photons and electrons, respectively), then the high energy part of the IC component has a spectral slope typical for the Thomson regime with an abrupt cutoff close to $E_{\mathrm{max}}$. The spectra typical for the Klein-Nishina regime are formed above $m_e^2c^4/\varepsilon_{\mathrm{min}}$. If the spectrum of target photons features a cooling break, i.e., a change of the photon index by $0.5$ at $\varepsilon_{\mathrm{br}}$, then the transition to the Klein-Nishina regime proceeds through an intermediate change of the photon index by $0.5$ at $m_e^2c^4/\varepsilon_{\mathrm{br}}$., Comment: 12 pages, 7 figures

MeV nuclear de-excitation lines serve as a unique tool to study low-energy cosmic rays (CRs), containing both spectral and elemental information of the interacting material. In this paper, we estimated the possible nuclear de-excitation lines from the young supernova remnant Cassiopeia A. Given different CR spectral shapes and interacting materials, we found the predicted fluxes of strong narrow line emissions from the remnant are highly model-dependent, ranging from about $1\times10^{-10}\,{\rm \,cm^{-2}\,s^{-1}}$ to $1\times10^{-6}\, {\rm \,cm^{-2}\,s^{-1}}$ for the 4.44 MeV narrow line and from about $4\times10^{-11}\,{\rm \,cm^{-2}\,s^{-1}}$ to $2\times10^{-7}{\rm \,cm^{-2}\,s^{-1}}$ for the 6.13 MeV narrow line, respectively. Based on the new estimation, we also discussed the detection probability of these line emissions against the MeV diffuse Galactic background under different assumptions of instrument response functions., Comment: 6 pages, 6 figures, accepted for publication in MNRAS

Electron Compton scattering of target photons into the gamma-ray energy band (inverse Compton scattering --IC--) is commonly expected to dominate the very high energy spectra in gamma-ray bursts especially during the afterglow phase. For sufficiently large center-of-mass energies in these collisions, the effect of the electron recoil starts reducing the scattering cross section (the Klein-Nishina regime). The IC spectra generated in the Klein-Nishina regime is softer and has a smaller flux level compared to the synchrotron spectra produced by the same electrons. The detection of afterglow emission from nearby GRB 190819A in the very high energy (VHE) domain with H.E.S.S. has revealed an unexpected feature: the slope of the VHE spectrum matches well the slope of the X-ray spectra, despite expectations that for the IC production process, the impact of the Klein-Nishina effect should be strong. The multi-wavelength spectral energy distribution appears to be inconsistent with predictions of one-zone synchrotron-self-Compton models. We study the possible impact of two-zone configuration on the properties of IC emission when the magnetic field strength differs considerably between the two zones. Synchrotron photons from the strong magnetic field zone provide the dominant target for cooling of the electrons in the weak magnetic field zone, which results in a formation of hard electron distribution and consequently of a hard IC emission. We show that the two-zone model can provide a good description of the X-ray XRT and VHE H.E.S.S. data., Comment: 13 pages, 7 figures, ApJ submitted

The density of cosmic rays inside molecular clouds determines the ionization rate in the dense cores where stars form. It is also one of the drivers of astrochemistry leading to the creation of complex molecules. Through Fermi Large Area Telescope observations of nearby giant molecular clouds, we observed deficits (holes) in the gamma-ray residual map when modelling with the expected gamma-ray diffuse emission from uniform cosmic rays interacting with the molecular content. We propose that the deficit is due to the lack of penetration of the low-energy (sub-GeV to GeV) cosmic rays into denser regions or clumps. This differs from the prevailing view of fast cosmic ray transport in giant molecular clouds where the magnetic turbulence is suppressed by neutral-ion damping, as our results require a slow diffusion inside dense molecular clumps. Through modelling we find that while the shielding is negligible on the cloud scale, it becomes important in the denser, parsec-sized regions where the gravitational collapse is already at play, changing the initial condition of star formation and astrochemistry., Comment: Published in Nature Astronomy

The current generation of very-high-energy $gamma-$ray (VHE; E above 30 GeV) detectors (MAGIC and H.E.S.S.) have recently demonstrated the ability to detect the afterglow emission of GRBs. However, the GRB prompt emission, typically observed in the 10 keV-10 MeV band, has so far remained undetected at higher energies. Here, we investigate the perspectives of multi-messenger observations to detect the prompt emission of short GRBs in VHE. Considering binary neutron star mergers as progenitors of short GRBs, we evaluate the joint detection efficiency of the Cherenkov Telescope Array (CTA) observing in synergy with the third generation of gravitational wave detectors, such as the Einstein Telescope (ET) and Cosmic Explorer (CE). In particular, we evaluate the expected capabilities to detect and localize gravitational wave events in the inspiral phase and to provide an early warning alert able to drive the VHE search. We compute the amount of possible joint detections by considering several observational strategies, and demonstrate that the sensitivities of CTA make the detection of the VHE emission possible even if it is several orders fainter than the one observed at 10 keV-10 MeV. We discuss the results in terms of possible scenarios of production of VHE photons from binary neutron star mergers by considering GRB prompt and afterglow emissions., Comment: Accepted for publication in Astronomy and Astrophysics

X-ray observations of several kiloparsec-scale extragalactic jets favour a synchrotron origin. The short cooling times of the emitting electrons requires distributed acceleration of electrons up to sub-PeV energies. In a previous paper, we found that this can be self-consistently explained by a shear acceleration model, where particles are accelerated to produce power-law spectra with a spectral index being determined mainly by the velocity profile and turbulence spectrum. In this paper, we perform 3D relativistic magneto-hydrodynamic simulations to investigate the formation of a spine-sheath structure and the development of turbulence for a relativistic jet propagating into a static cocoon. We explore different spine velocities and magnetic field profiles with values being chosen to match typical Fanaroff-Riley type I/II jets. We find that in all cases a sheath is generated on the interface of the spine and the cocoon mainly due to the Kelvin-Helmholtz instability. The large scale velocity profile in the sheath is close to linear. Turbulence develops in both the spine and the sheath, with a turbulent velocity spectrum consistent with Kolmogorov-scaling. The implications for shear particle acceleration are explored, with a focus on the particle spectral index., Comment: 9 pages, 11 figures, 1 table, accepted for publication in MNRAS

We critically assess limits on the maximum energy of protons accelerated within superbubbles around massive stellar clusters, considering a number of different scenarios. In particular, we derive under which circumstances acceleration of protons above peta-electronvolt (PeV) energies can be expected. While the external forward shock of the superbubble may account for acceleration of particles up to 100 TeV, internal primary shocks such as supernova remnants expanding in the low density medium or the collective wind termination shock which forms around a young compact cluster provide more favourable channels to accelerate protons up to 1 PeV, and possibly beyond. Under reasonable conditions, clustered supernovae launching powerful shocks into the magnetised wind of a young and compact massive star cluster are found to be the most promising systems to accelerate protons above 10 PeV. On the other hand, stochastic re-acceleration in the strongly turbulent plasma is found to be much less effective than claimed in previous works, with a maximum proton energy of at most a few hundred TeV., Comment: 10 pages, 5 figures

We present a broadband X-ray study of W50 (`the Manatee nebula'), the complex region powered by the microquasar SS 433, that provides a test-bed for several important astrophysical processes. The W50 nebula, a Galactic PeVatron candidate, is classified as a supernova remnant but has an unusual double-lobed morphology likely associated with the jets from SS 433. Using NuSTAR, XMM-Newton, and Chandra observations of the inner eastern lobe of W50, we have detected hard non-thermal X-ray emission up to $\sim$30 keV, originating from a few-arcminute size knotty region (`Head') located $\lesssim$ 18$^{\prime}$ (29 pc for a distance of 5.5 kpc) east of SS 433, and constrain its photon index to 1.58$\pm$0.05 (0.5-30 keV band). The index gradually steepens eastward out to the radio `ear' where thermal soft X-ray emission with a temperature $kT$$\sim$0.2 keV dominates. The hard X-ray knots mark the location of acceleration sites within the jet and require an equipartition magnetic field of the order of $\gtrsim$12$\mu$G. The unusually hard spectral index from the `Head' region challenges classical particle acceleration processes and points to particle injection and re-acceleration in the sub-relativistic SS 433 jet, as seen in blazars and pulsar wind nebulae., Comment: 19 pages, 9 figures, 2 tables, accepted for publication in the Astrophysical Journal

The Fermi Bubbles are giant, gamma-ray emitting lobes emanating from the nucleus of the Milky Way discovered in ~1-100 GeV data collected by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope. Previous work has revealed substructure within the Fermi Bubbles that has been interpreted as a signature of collimated outflows from the Galaxy's super-massive black hole. Here we show via a spatial template analysis that much of the gamma-ray emission associated to the brightest region of substructure -- the so-called cocoon -- is likely due to the Sagittarius dwarf spheroidal (Sgr dSph) galaxy. This large Milky Way satellite is viewed through the Fermi Bubbles from the position of the Solar System. As a tidally and ram-pressure stripped remnant, the Sgr dSph has no on-going star formation, but we nevertheless demonstrate that the dwarf's millisecond pulsar (MSP) population can plausibly supply the gamma-ray signal that our analysis associates to its stellar template. The measured spectrum is naturally explained by inverse Compton scattering of cosmic microwave background photons by high-energy electron-positron pairs injected by MSPs belonging to the Sgr dSph, combined with these objects' magnetospheric emission. This finding plausibly suggests that MSPs produce significant gamma-ray emission amongst old stellar populations, potentially confounding indirect dark matter searches in regions such as the Galactic Centre, the Andromeda galaxy, and other massive Milky Way dwarf spheroidals., Comment: Updated to match version accepted for publication in Nature Astronomy (2022). 14 pages main text, 3 main figures, 7 extended data figures. For published version of the paper, see https://rdcu.be/cUZ6X

The recent discovery of a new population of ultra-high-energy gamma-ray sources with spectra extending beyond 100 TeV revealed the presence of Galactic PeVatrons - cosmic-ray factories accelerating particles to PeV energies. These sources, except for the one associated with the Crab Nebula, are not yet identified. With an extension of 1 degree or more, most of them contain several potential counterparts, including Supernova Remnants, young stellar clusters and Pulsar Wind Nebulae (PWNe), which can perform as PeVatrons and thus power the surrounding diffuse ultra-high energy gamma-ray structures. In the case of PWNe, gamma rays are produced by electrons, accelerated at the pulsar wind termination shock, through the inverse Compton scattering of 2.7 K CMB radiation. The high conversion efficiency of pulsar rotational power to relativistic electrons, combined with the short cooling timescales, allow gamma-ray luminosities up to the level of $L_\gamma \sim 0.1 \dot{E}$. The pulsar spin-down luminosity, $\dot E$, also determines the absolute maximum energy of individual photons: $E_{\rm \gamma~\rm max}\approx 0.9 \dot E_{36}^{0.65}~~\rm{PeV}$. This fundamental constraint dominates over the condition set by synchrotron energy losses of electrons for young PWNe with typical magnetic field of $\approx$100$\mu$G with $\dot{E} \lesssim 10^{37}\ \rm erg/s$. We discuss the implications of $E_{\rm \gamma~\rm max}$ by comparing it with the highest energy photons reported by LHAASO from a dozen of ultra-high-energy sources. Whenever a PWN origin of the emission is possible, we use the LHAASO measurements to set upper limits on the nebular magnetic field., Comment: accepted for publication in ApJ Letters