An extreme particle accelerator in the Galactic plane: HESS J1826$-$130

The unidentified very-high-energy (VHE; E $>$ 0.1 TeV) $\gamma$-ray source, HESS J1826$-$130, was discovered with the High Energy Stereoscopic System (HESS) in the Galactic plane. The analysis of 215 h of HESS data has revealed a steady $\gamma$-ray flux from HESS J1826$-$130, which appears extended with a half-width of 0.21$^{\circ}$ $\pm$ 0.02$^{\circ}_{\text{stat}}$ $\pm$ 0.05$^{\circ}_{\text{sys}}$. The source spectrum is best fit with either a power-law function with a spectral index $\Gamma$ = 1.78 $\pm$ 0.10$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$ and an exponential cut-off at 15.2$^{+5.5}_{-3.2}$ TeV, or a broken power-law with $\Gamma_{1}$ = 1.96 $\pm$ 0.06$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$, $\Gamma_{2}$ = 3.59 $\pm$ 0.69$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$ for energies below and above $E_{\rm{br}}$ = 11.2 $\pm$ 2.7 TeV, respectively. The VHE flux from HESS J1826$-$130 is contaminated by the extended emission of the bright, nearby pulsar wind nebula (PWN), HESS J1825$-$137, particularly at the low end of the energy spectrum. Leptonic scenarios for the origin of HESS J1826$-$130 VHE emission related to PSR J1826$-$1256 are confronted by our spectral and morphological analysis. In a hadronic framework, taking into account the properties of dense gas regions surrounding HESS J1826$-$130, the source spectrum would imply an astrophysical object capable of accelerating the parent particle population up to $\gtrsim$200 TeV. Our results are also discussed in a multiwavelength context, accounting for both the presence of nearby supernova remnants (SNRs), molecular clouds, and counterparts detected in radio, X-rays, and TeV energies.


Introduction
The High Energy Stereoscopic System (HESS) Galactic Plane Survey (HGPS) revealed the presence of 16 new γ-ray sources in the very high energy (VHE; E > 0.1 TeV) domain (H.E.S.S. Collaboration et al. 2018). Among them, HESS J1826−130 is located in a region that is exceptionally rich in VHE γray sources, also encompassing the nearby pulsar wind nebula (PWN) HESS J1825−137 and the γ-ray binary system, LS 5039. HESS J1826−130 belongs to a growing class of VHE γ-ray sources affected by source confusion, caused by the influence of luminous neighboring γ-ray emitters, similar to the case of HESS J1641−463 (Abramowski et al. 2014). HESS J1826−130 is located just ∼ 0.7 • to the north of HESS J1825−137 (H.E.S.S. Collaboration et al. 2019a), a bright PWN whose extended emission prevented an earlier discovery of HESS J1826−130 at VHEs.
While HESS J1826−130 was clearly distinguishable from HESS J1825−137 in the HGPS, a deep study accounting for its morphological and spectral properties beyond the results obtained with the standard HGPS analysis tools was still lacking. In this paper, the results obtained from the analysis of 215 hours of VHE γ-ray data recorded with the HESS telescopes is presented, yielding an accurate morphological description and a detailed spectral characterization of the source, including an energydependent emission study of the region. The results reported here on HESS J1826−130 are discussed in a multiwavelength context, discussing in particular its possible association with the spatially coincident supernova remnant, G18.45-0.42, and the Eel Nebula (Roberts et al. 2007), as well as with the VHE source 2HWC J1825−134 (Abeysekara et al. 2017), which displays emission up to energies of around 100 TeV.
Furthermore, HESS J1826−130 is surrounded by dense molecular clouds that provide support for a hadronic emission scenario. A recent investigation of the diffuse VHE emission around the Galactic center (GC) region demonstrated the presence of PeV particles within the central 10 pc of the Galaxy (HESS Collaboration et al. 2016). This is the first robust detection of a VHE cosmic-ray (CR) accelerator which operates as a PeVatron. The MAGIC Collaboration recently published (MAGIC Collaboration et al. 2020) the results of observations of the diffuse γ-ray emission in the vicinity of the GC. However, the GC alone would not be able to sustain the CR population close to PeV energies today unless the GC PeVatron was more powerful in the past. This opens up the possibility of other active PeV accelerators in the Galaxy that contribute to the observed Galactic CR flux around the spectral feature known as the knee, that is, a potentially new CR source population in the Galaxy. HESS J1826−130 may belong to such a new Galactic CR source population.
In Sect. 2, we present a summary of the multiwavelength (MWL) information available in the field of view (FoV) around HESS J1826−130. The HESS observations and data analysis are reported in Sect. 3, including morphological and spectral results, as well as a study of the contamination from the nearby PWN HESS J1825−137. Results are discussed in Sect. 4, with concluding remarks in Sect. 5.

Multiwavelength environment
As mentioned before, HESS J1826−130 is located in a region with several VHE γ-ray sources. At VHEs, the TeV sky around the source is dominated by one of the brightest and largest PWNs detected at VHEs, HESS J1825−137. This latter source is thought to be powered by PSR J1826−1334 (also known as PSR B1823−13), a powerful pulsar with a spin-down luminosity of ∼ 3 × 10 36 erg s −1 , located about 13' from the peak position of HESS J1825−137 in the VHE band. At GeV energies, both the spectral and the morphological properties of HESS J1825−137 resemble and smoothly connect with those observed at VHEs, as reported in Grondin et al. (2011) using the Fermi-LAT.
The first evidence for γ-ray emission originating from the north of PSR J1826−1334 was found by the detection of 3EG J1826−1302 (also called GeV J1825−1310) in 1999 with the EGRET telescope on-board the Compton Gamma Ray Observatory (Hartman et al. 1999). Follow-up X-ray observations of the same region led to the detection of a diffuse X-ray source, AX J1826.1−1300, in the 2−10 keV energy band (Roberts et al. 2001a). This X-ray source is characterized by two emission peaks located on top of a diffuse X-ray nebula. Further observations of AXJ1826.1−1300 in 2007 with the Chandra X-ray Observatory resolved the source in more detail (Roberts et al. (2007)). The southern peak was found to be likely associated with a stellar cluster, whereas the northern peak was resolved and associated with a PWN. Due to a 4 trail of hard X-ray emission connected to this new nebula, the authors called it the Eel Nebula. Roberts et al. (2007) were the first to introduce a new separate VHE γ-ray source north of HESS J1825−137, which they named HESS J1826−131 and associated with the Eel Nebula. Until that time, the remaining VHE emission to the north of HESS J1825−137 was thought to be connected to the bright source itself.
The Fermi/LAT detected the radio-quiet γ-ray pulsar PSR J1826−1256 (Abdo et al. 2009), with a spin down luminos-ityĖ = 3.6 × 10 36 erg s −1 and a characteristic age τ c = 14.4 kyr, whose position is consistent with the Eel Nebula, representing a significant step forward in our understanding of the various sources contributing to the γ-ray emission from this region. In 2017, HESS J1826−130 was reported as a new source within the HGPS. The HGPS used a dedicated source-identification algorithm based on a multi-Gaussian component fit to a given region of interest covered by the survey (H.E.S.S. Collaboration et al. 2018). In particular for HESS J1826−130, the spatially extended VHE γ-ray emission region matched the position of 3EG J1826−1302, also overlapping with the diffuse X-ray nebula observed with ASCA (which displays a radius of about 15 ; see Fig. 11 in (Roberts et al. 2001b)) as well as the Eel Nebula (with an extension of about 4 (Roberts et al. 2007)). The latest results at VHEs of the region include the significant detection of an emission region, 2HWC J1825−134, included in the 2 nd HAWC catalog (Abeysekara et al. 2017). Due to the large spatial extension of the HAWC emission region and the relatively low angular separation between HESS J1826−130 and HESS J1825−137, the two sources were not distinguishable in this catalog. Further studies, however, using a much larger dataset and employing updated analysis methods, provide a statistical improvement at the level of ∼ 5σ when accounting for the presence of the two HESS VHE sources.
Within the surroundings of HESS J1826−130, there are also three supernova remnants (SNRs): G18.1−0.1, G18.6−0.2, and G18.45−0.42, which are detected in the radio band (Brogan et al. 2006;Anderson et al. 2017). Detailed studies of the gas dynamics and densities in the direction of the region north of HESS J1825−137 were carried out by Voisin et al. (2016), following the discovery of a molecular cloud north of PSR J1826−1334 by Lemiere et al. (2005). These studies revealed the presence of dense molecular clouds with particle density values up to n H 2~7 × 10 2 cm −3 . The bulk of molecular gas in the vicinity of HESS J1826−130 is located at V LSR = 45-60 km/s, corresponding to a kinematic distance of 4.0 kpc, and V LSR = 60-80 km/s, corresponding to a kinematic distance of 4.6 kpc. The former gas distance estimate is consistent with the dispersion measurement of PSR J1826−1334. Voisin et al. (2016) concluded that due to the general spatial match between the molecular gas and the TeV emission, HESS J1826−130 might be explained by a hadronic emission scenario, with the progenitor SNR of HESS J1825−137 being a suitable local CR source candidate. However, a leptonic scenario connected to the Eel Nebula could not be ruled out. Furthermore, the two SNRs, G18.1−0.1 and G18.6−0.2, were considered unlikely to be directly related to the TeV excess due to their offset positions and their small angular diameters.

HESS observations and data analysis results
The High Energy Stereoscopic System is an array of five imaging atmospheric Cherenkov telescopes located in the Khomas Highland of Namibia, 1800 m above sea level. HESS in phase I comprised four 12 m diameter telescopes that have been fully operational since 2004. A fifth telescope (CT5), with a larger mirror diameter of 28 m and newly designed camera (Bolmont et al. 2014), was added to the center of the array and has been operational since September 2012. The HESS phase I array configuration is sensitive to γ-ray energies between 100 GeV and several tens of TeV, while the energy threshold of the array was lowered to 10's of GeV with the addition of CT5. The VHE HESS data presented in this paper were recorded with the HESS phase I array configuration, which can measure extensive air showers with an angular resolution better than 0.1 • and an average energy resolution of 15% for an energy of 1 TeV (Aharonian et al. 2006c).
The VHE γ-ray observations of the FoV around HESS J1826−130 were carried out between 2004 and 2015.
Article number, page 3 of 9 A&A proofs: manuscript no. HESSJ1826-130_An_extreme_particle_accelerator_in_the_Galactic_plane_arxiv About 140 h of observations were recorded with the HESS phase I array configuration (between 2004 and 2012), while ∼95 h were obtained with the HESS phase II array configuration in a split observation mode without participation of the CT5 telescope. These data sets provide a total acceptance-corrected live-time of 215 h of HESS data after the application of quality selection criteria (Aharonian et al. 2006c).
The data have been analyzed with the HESS Analysis Package (HAP) for shower reconstruction, and a multivariate analysis technique (Ohm et al. 2009) was applied to provide an improved distinction between hadron and γ-ray events. A cross-check analysis was performed using an independent calibration and analysis method (de Naurois & Rolland 2009), which gives compatible results with the main analysis.

Detection and morphological analysis
In order to achieve an improved angular resolution and reduce strong contamination coming from the nearby bright source, HESS J1825−137, the source position and morphology of HESS J1826−130 were obtained with a cut configuration optimized for high energies that requires a minimum of 160 photoelectrons per image. The cosmic-ray background level was estimated using the ring background model (Berge et al. 2007). Using the full 215 h dataset, HESS J1826−130 was detected with a statistical significance of 22.6σ 1 , determined following Equation (17) in Li & Ma (1983). Figure Fig. 2, with reduced contamination from HESS J1825−137, potentially due to its energy-dependent morphology and hinting at a harder VHE spectrum.
The best-fit position and extension of the source have been determined by convolving 2D Gaussian models with the HESS PSF and fitting this to the excess event distribution (for E > 0.4 TeV). The HESS PSF was derived from simulations and modeled as a weighted sum of three 2D Gaussian functions for the morphology analysis. A log-likelihood ratio test (LLRT) was used to select the best morphology model that represents the HESS data. The extension (radius) of the source is estimated to be 0.21 • ± 0.02 • stat ± 0.05 • sys , which is significantly larger than the HESS PSF.

Spectral analysis
The green circular region with a radius of 0.22 • centered at the best-fit position of HESS J1826−130 (see Fig. 1) was used as the integration region for measuring the differential VHE γ-ray spectrum of the source. The reflected background model (Berge  Table 1. Analysis statistics of the HESS J1826−130 region using different energy thresholds (E th ) corresponding to the VHE γ-ray excess profiles shown in Fig. 2. The events below the energy threshold shown in the E th column were not taken into account in the analysis. The number of events within the source region (N On ), within the ring region used for estimating the background level (N Off ), the alpha factor (α), the corresponding number of excess events, and the Li&Ma statistical significance (Li & Ma 1983) (Piron et al. 2001). Several spectral models were used to fit the data, including a simple power-law (PL) model, a broken power-law (BPL) model, and a power-law model with an exponential cut-off (ECPL; see Table 2). The spectral fit is performed between 0.42 TeV and 56.2 TeV. The best-fit model is found to be the ECPL, with a p-value = 0.41 and a likelihood ratio test (LLRT) against the PL model providing a statistical improvement at the ∼ 3.2σ level. The differential VHE emission is thus best represented by dN/    0.69 stat ± 2.02 sys ) × 10 −13 cm −2 s −1 TeV −1 , Γ = 1.78 ± 0.10 stat ± 0.20 sys and a cut-off energy of E c = 15.2 +5.5 −3.2 TeV as shown in Fig. 3. The integral flux level above 1 TeV is Φ(> 1 TeV) = (10.2 ± 0.5 stat ± 2.0 sys ) × 10 −13 cm −2 s −1 and corresponds to 4.9% of the Crab Nebula flux at the same energies 2 . A 2σ lower limit on the cut-off energy of the source is derived at ∼ 10 TeV. A fit to a BPL model gives a break energy at 11.2 ± 2.7 TeV, which is compatible with the cut-off energy obtained from the ECPL model. The comparison between the ECPL and the BPL model does not show any preference from obtained p-values. Both models clearly indicate a spectral steepening of the spectrum above 10 TeV. The integral flux is found to be constant within the HESS dataset at different timescales, from run-wise (∼30 min) to year-wise light curves.
The integration region used to measure the VHE spectrum of HESS J1826−130 is strongly contaminated, especially at lower energies, by the relatively softer spectrum of HESS J1825−137. Therefore, the obtained spectrum is affected by this contamination. Taking into account the relatively softer contribution from the contaminating source, the intrinsic spectrum of HESS J1826−130 is expected to be even harder with respect to the spectral results discussed above.

Contamination study
The estimation of intrinsic spectra for sources strongly contaminated by neighboring bright and extended sources presents a nontrivial challenge. Although morphology studies in different energy bands can provide a powerful tool for new discoveries, such a method does not help in providing a quantitative characterization of the contamination, as systematic uncertainties in the flux estimation are expected to dominate. In general, contamination effects can be neglected if they are at the 10% level or lower, since this is, for example, the typical HESS systematic error on the source flux. For contamination exceeding 10%, in particular for spectral studies, a dedicated analysis is needed.
Performing contamination studies for HESS J1826−130 requires knowledge about the energy-dependent morphology of HESS J1825−137, which is highly asymmetrical. The contamination was estimated in this case using the morphology model that represents the HESS data best fit. The model includes a 2D Gaussian function accounting for HESS J1826−130 along with other 2D Gaussian functions representing the other sources, HESS J1825−137 and the gamma-ray binary system LS 5039, present in the FoV. Taking the ratio of the total number of excess counts within the integration region (Fig. 1, green circle) used for measuring HESS J1826−130's spectrum between the model components can give a quantitative estimation of the contamination. By using such an approach, the contamination from HESS J1825−137 is estimated to be ∼40% below 1.5 TeV and ∼20% above 1.5 TeV, while it is strongly reduced (<10%) at energies above 2.0 TeV. Given the relatively large angular distance of LS 5039 to the integration region (around 2.0 • ) encompassing HESS J1826−130, we find that the contamination from the binary system is negligible.
The spectral contamination from HESS J1825−137 has been estimated based on the fact that as the separation between contaminated and contaminating source increases, one expects the effect of contamination to decrease. To test the biasing impact of the contamination, the HESS J1826−130 integration region was split into a near and far half-circle, referred to as the south and north semi-circles shown in Fig. 4. Spectra from two control regions at both sides of the source position were also extracted in order to determine whether or not the contribution from HESS J1825−137 Article number, page 5 of 9 A&A proofs: manuscript no. HESSJ1826-130_An_extreme_particle_accelerator_in_the_Galactic_plane_arxiv Table 3. Comparison of ECPL spectral models for the HESS J1826−130 integration region (Total), plus the north and south semi-circle regions. The integral flux (in % Crab Nebula) is given for the same energy range (above 1 TeV).
One expects that the contamination in the north semi-circle region should be less with respect to the south semi-circle region, and the spectrum should reflect the intrinsic spectral properties of HESS J1826−130. The spectral fit results, assuming an ECPL model and obtained from the north and south semi-circle regions, are given in Table 3 along with the total results. As expected, the north semi-circle spectrum is slightly harder than those of the total and south semi-circle regions. Although the contamination from HESS J1825−137 could potentially affect the spectrum of HESS J1826−130, our studies demonstrate that the effects on the spectral parameters are at the level of our analysis systematic uncertainties. The intrinsic spectrum of HESS J1826−130, corrected for this contamination, is found to be well described by an ECPL model where Φ 0 = 9.2 × 10 −13 cm −2 s −1 TeV −1 , Γ = 1.7 and a cut-off energy around E c = 16 TeV.
A further investigation of the intrinsic spectrum of HESS J1826−130 was performed using a 3D cube-analysis approach 3 . With such a template-based analysis method, uncontaminated spectra can be derived, since all sources located in the FoV are simultaneously taken into account. The 3D analysis results published in Ziegler (2018) are found to be compatible with the results presented here.

Discussion
The study of VHE γ-ray sources displaying spectra extending above several tens of TeV is required to understand the origin of the highest energy CRs close to the knee in the CR spectrum. These studies are particularly promising in cases in which dense gas regions are coincident with the source of interest. In the case of HESS J1826−130, the spatial coincidence of a dense molecular hydrogen region found along the line of sight (see Voisin et al. (2016)), suggests that the radiation could actually be produced by protons with energies of several hundred TeV colliding with this gas. Figure 5 displays the HESS excess map obtained at energies above 20 TeV together with column density profiles derived from Nobeyama 12 CO(1-0) (Umemoto et al. 2017) and SGPS HI data (McClure-Griffiths et al. 2005), for two different gas velocity ranges, 45-60 km/s (corresponding to a kinematic distance of 4.0 kpc) and 60-80 km/s (corresponding to a kinematic distance of 4.6 kpc). The scaling X-factors, X HI = 1.8 × 10 18 cm −2 (K km s −1 ) −1 and X CO = 2.0 × 10 20 cm −2 (K km s −1 ) −1 , were taken from Dickey & Lockman (1990) A. W. et al. (2004), respectively. The Nobeyama 12 CO(1-0) data was smoothed to match the SGPS beam size of 2'. The CO and HI distribution display two local maxima that do not coincide with the E > 20 TeV peak, although they still provide target densities at the level of ∼ 5 × 10 2 cm −3 at the position of HESS J1826−130.
The best-fit parent proton population in this scenario would have a spectral index of ∼1.7 and a cut-off energy around 200 TeV 4 . The relative hardness of the photon spectrum could be due to the highest energy protons diffusing into dense clouds, while lower energy protons might still be confined within the accelerating source, or efficiently excluded from the clouds if the diffusion coefficient inside the clouds is suppressed (Gabici et al. 2009). A hadronic scenario in which runaway protons accelerated by one of the coincident or nearby SNRs emit TeV photons when interacting with the dense ambient gas found along the line of sight could explain the hard spectrum of HESS J1826−130. We note, however, that the filled centers of SNR G18.5-0.04 and SNR G18.1-00.1 are possibly too old to be able to accelerate protons up to hundreds of TeV. The γ-ray luminosity, L γ , of the source is 8 × 10 33 erg s −1 for a distance of 4 kpc, corresponding to the estimated distances of the SNRs G18.45−0.42 (Karpova et al. 2019), G018.1−00.1, or G018.6−00.2 5 . This would translate into an energy output in accelerated protons W pp = L γ × t pp , of 6 × 10 49 (n/1 cm −3 ) −1 erg, where the gas density value is n = 5 × 10 2 cm −3 and t pp is the typical timescale for gamma-rays to be produced through neutral pion decay following the collision of relativistic protons.
The relatively hard spectrum found at VHE gamma rays could also be produced in a leptonic scenario by an electron population with a spectral index close to 2.0 and a cut-off at around 70 TeV. These electrons could be accelerated by the pulsar PSR J1826-1256 and up-scatter cosmic microwave background (CMB) and IR photons to produce the observed TeV emission. PSR J1826−1256 is powering the Eel Nebula (PWN G018. 5-0.4, Roberts et al. 2007), which displays an elongated X-ray trail of shocked material associated with its fast proper motion through the surrounding medium. An estimate for the magnetic field inside the Eel Nebula of ∼20-30 µG is provided in Keogh (2010). Recently, Duvidovich et al. (2019) made use of XMM-Newton observations of the region around PSR J1826-1256 and discussed such a possibility. Based on the spectral softening found along the nebula, they considered X-rays produced by electrons with energies E e 150 TeV radiating synchrotron emission in a magnetic field B > 2µG. In a parallel study, Karpova et al. (2019) also reported on the X-ray analysis of this very same XMM-Newton dataset, complemented with additional ∼ 90 ks of Chandra observations of the region. These data support the association of HESS J1826−130 with the Eel Nebula, although somewhat harder spectral indices close to the PWN are derived (Γ X ∼ 1.2), followed by softer emission, with Γ X ∼ 2.5 at larger distances along the PWN tail.
The excess map reported in Fig. 5 at E > 20 TeV traces emission from the vicinity of PSR J1826−1256. TeV photons could be produced by these synchrotron X-ray-emitting electrons through inverse Compton scattering off the CMB and IR photon fields. For the latter, an energy density at a distance of 4 kpc ∼0.23 eV cm −3 was assumed (see e.g., Vernetto & Lipari 2016; we note, however, that the distance to PSR J1826−1256 is not very well constrained; e.g., Wang 2011 favored a distance of about 1.2 kpc, whereas Karpova et al. 2019 found a distance of ∼ 3.5 kpc based on an empirical interstellar absorption-distance correlation). The CMB photon field provides u CMB ∼0.26 eV cm −3 . Electron Lorentz factors of γ e ∼ 1.3 × 10 8 and γ e ∼ 3.6 × 10 7 would be required to reach the VHE γ-ray fluxes reported here for CMB and IR seed photons, respectively, in agreement with the estimates in Duvidovich et al. (2019).
The complex region around HESS J1826−130 has also been observed at TeV energies with the HAWC telescopes. A spatially coincident source, 2HWC J1825−13, is included in the 2 nd HAWC catalog (Abeysekara et al. 2017). Furthermore, the HAWC collaboration recently reported the detection of E > 100 TeV emission from eHWC J1825−134 (HAWC Collaboration et al. 2019). Its position and extension are marked with a white dotted circle in Fig. 5. With a half-width of ∼ 0.36 • , the source encompasses a local peak of the 12 CO distribution. It is also compatible with the location of the E > 20 TeV γ-rays from HESS J1825−137, and adjacent to HESS J1826−130 emission E > 20 TeV. On spectral grounds, both HESS J1825−137 and HESS J1826−130 display statistically preferred cut-offs and/or curvature in their spectral analysis (see Sect. 3.2 and H.E.S.S. Col- 5 We note that Voisin et al. (2016) discuss the possibility of a larger distance to G018.6−00.2 of 11.4 kpc, disfavoring an association of HESS J1826−130 with G018.6−00.2.

Conclusions
Finally, HESS J1826−130 is a particle accelerator in the Galactic plane capable of delivering TeV emission above several tens of TeV. The deep observation campaign on the source reported in this paper provides several constraints to its spectral and morphological properties, although its origin is still uncertain. The presence of nearby SNRs found in the region and/or another yet unresolved accelerator, together with the gas density properties of the surroundings, could make hadronic interactions responsible for the observed emission. If protons are delivered with energies up to hundreds of TeV, HESS J1826−130 could point to the existence of a distinct population of CR sources that, along with the GC, significantly contribute to the Galactic CR flux around the knee feature. In a leptonic scenario, HESS J1826−130 has been proposed as the counterpart of the Eel Nebula powered by the pulsar PSR J1826−1256. Our observations reveal VHE emission above several tens of TeV that comes from the vicinity of the pulsar. Regarding the association of HESS J1826−130 with the recently discovered 100 TeV emitter, 2HWC J1825−13, it is worth noting that the spectrum reported in Fig. 3 extends at least up to ∼ 43 TeV. Therefore, it could be possible that a fraction of the emission reported from eHWC J1825−134 could be partially contributed by HESS J1826−130 (see Salesa Greus & Casanova 2019). On the other hand, HESS J1825−137 also displays VHE emission above tens of TeV, albeit displaying a