Search

Context. High-energy (HE) and very high-energy (VHE) gamma-ray observations from the Galactic center (GC) detected extended emission correlated with the morphology of the central molecular zone (CMZ). Emission in both bands is expected to be produced by hadronic interaction between cosmic rays (CRs) and ambient gas. Aims. We examine if our three previously proposed scenarios for the CR sources and dynamics, which are consistent with the VHE gamma-ray data (1–100 TeV), also match the HE gamma-ray observations (10–300 GeV). Additionally, we analyze the effect of the isotropic Galactic CR "sea" inside the CMZ. Methods. We generated synthetic gamma-ray maps considering a simplified isotropic diffusion, but more realistic dynamics with two diffusion zones (in and out of the CMZ) and polar advection, for mono-energetic particles of 3 TeV. Additionally, we considered two gas distributions for the CMZ (with and without an inner cavity), and CR populations injected from the clusters of young massive stars (the Arches Cluster, the Quintuplet Cluster, and the nuclear star cluster), plus the supernova Sgr A East. Results. Only the combination of more realistic CR dynamics, the CMZ with an inner cavity, CR injection from all proposed sources, and a CR sea similar to that observed in the Solar System reproduced the current HE and VHE gamma-ray detection from the CMZ and was consistent with the observed gamma-rays from Sagittarius A*. Conclusions. The HE and VHE gamma-rays observations of the GC can be reproduced by a unified model for the CRs. [ABSTRACT FROM AUTHOR]

Context. Estimation of the number of cosmic-ray-induced background events is a challenging task for Imaging Atmospheric Cherenkov Telescopes (IACTs). Most approaches rely on a model of the background signal derived from archival observations, which is then normalised to the region of interest (ROI) and respective observation conditions using emission-free regions in the observation. However, this is disadvantageous for the analysis of large, extended γ-ray structures, where no sufficient source-free region can be found. Aims. We aim to address this issue by estimating the normalisation of a three-dimensional background model template from separate, matched observations of emission-free sky regions. As a result, the need for an emission-free region in the field of view of the observation becomes unnecessary. Methods. To this end, we implemented an algorithm to identify observation pairs with the most closely matching observation conditions. We used the open-source analysis package Gammapy to estimate the background rate, facilitating seamless adaptation of the framework to many γ-ray detection facilities. We employed public data from the High Energy Stereoscopic System (H.E.S.S.) to validate this methodology. Results. The analysis demonstrates that employing a background rate estimated through this run-matching approach yields results consistent with those obtained using the standard application of the background model template. Furthermore, we confirm the compatibility of the source parameters obtained through this approach with previous publications, and present an analysis employing the background model template approach, along with an estimation of the statistical and systematic uncertainties introduced by this method. [ABSTRACT FROM AUTHOR]

Context. There is growing evidence from gamma-ray observations at high and very high energies that particle escape is a key aspect shaping the morphological properties of pulsar wind nebulae (PWNe) at various evolutionary stages. Aims. We aim to provide a simple model for the gamma-ray emission from these objects including the transport of particles across the different components of the system. We applied it to sources HESS J1809−193 and HESS J1825−137. Methods. We developed a multi-zone framework applicable to dynamically young PWNe, taking into account the diffusive escape of relativistic electron-positron pairs out of the nebula into the parent supernova remnant (SNR) and their confinement downstream of the magnetic barrier of the forward shock until an eventual release into the surrounding interstellar medium (ISM). Results. For a wide range of turbulence properties in the nebula, the GeV–TeV inverse-Compton radiation from pairs that escaped into the remnant can be a significant if not dominant contribution to the emission from the system. It may dominate the pion-decay radiation from cosmic rays accelerated at the forward shock and advected downstream of it. In the TeV–PeV range, the contribution from particles escaped into the ISM can exceed by far that of the SNR+PWN components. Applied to HESS J1809−193 and HESS J1825−137, we found that spatially extended GeV–TeV emission components can be accounted for mostly from particles escaped into the ISM, while morphologically more compact components above 50 − 100 TeV are ascribed to the PWNe. In these two cases, the model suggests high turbulence in the nebula and a forward shock accelerating cosmic rays up to ∼100 TeV at most. Conclusions. The model provides the temporal and spectral properties of the flux of particles originally energized by the pulsar wind and ultimately released in the ISM. It can be used to constrain the transport of particles in the vicinity of pulsar-PWN-SNR systems from broadband gamma-ray observations, or in studies of the contribution of pulsar-related systems to the local positron flux. [ABSTRACT FROM AUTHOR]

While the data analysis of γ-ray telescopes has now become more robust, some signals may be misinterpretations of a time-variable foreground emission from the Solar system, induced by low-energy cosmic-ray interactions with asteroids. Our goal is to provide emission templates for this time-variable diffuse γ-ray foreground by considering the populations of Main Belt Asteroids, Jovian and Neptunian Trojans, Kuiper Belt Objects, and the Oort Cloud. We model the spatial distribution of all known asteroids by performing 3D-fits to determine their density profiles and calculate their appearances by line-of-sight integrations. Because Earth and the asteroids are moving with respect to each other, we obtain diffuse emission templates varying on time-scales of days to decades. We find that the temporal variability can lead to flux enhancements that may mimic emission features unless properly taken into account. This variation is further enhanced by the Solar cycle as the cosmic-ray spectrum is attenuated by the Solar modulation potential, leading to a relative flux increase of the outer asteroids. The cumulative effect of the time-dependent emission is illustrated for the case of the '511 keV OSSE fountain', and for emission features near the Galactic Centre, both being possible misinterpretations of the Solar system albedo. We recommend that γ-ray data analyses should always take into account the possibility of a time-variable foreground. Due to the ecliptic overlap with the Galactic plane, the Galactic emission is expected to be weaker by 0.1–20 per cent, depending on time (relative planetary motion), energy, and Solar cycle, which has consequences for the interpretation of dark matter annihilation cross sections, cosmic-ray spectra and amplitudes, as well as nucleosynthesis yields and related parameters. [ABSTRACT FROM AUTHOR]

We present a search for intraday transient γ-ray signals using 15.4 yr of the Fermi Large Area Telescope data. The search is based on a recently developed variable-size sliding-time-window (VSSTW) analysis and is aimed at studying variable γ-ray emission from gamma-ray bursts (GRBs) and the Sun. We refined the algorithm for searches for transient sources in order to solve the search problem within a reasonable amount of CPU time. These refinements allowed us to increase the number of gamma-ray bursts, solar flares, and quiescent solar events detected with the VSSTW technique by several times compared to the previous VSSTW search. The current search revealed a new γ-ray signal recorded with Fermi-LAT on 2018 January 12. This signal is probably from a GRB and deserves an exploration of the existing archival multiwavelength observations in order to identify it in an unambiguous way. We also report a γ-ray signal from the solar flare on 2023 December 31 that occurred during the 25th solar cycle. [ABSTRACT FROM AUTHOR]

NGC 1068 is a nearby, widely studied Seyfert II galaxy presenting radio, infrared, X-ray, and γ-ray emission, along with strong evidence for high-energy neutrino emission. Recently, the evidence for neutrino emission was explained in a multimessenger model, whereby the neutrinos originate from the corona of the active galactic nucleus. In this environment, γ-rays are strongly absorbed, so that an additional contribution is necessary, for instance, from the circumnuclear starburst ring. In this work, we discuss whether the radio jet can be an alternative source of the γ-rays between about 0.1 and 100 GeV, as observed by Fermi-LAT. In particular, we include both leptonic and hadronic processes, namely, accounting for inverse Compton emission and signatures from pp as well as pγ interactions. In order to constrain our calculations, we used VLBA and ALMA observations of the radio knot structures, which are spatially resolved at different distances from the supermassive black hole. Our results show that the best leptonic scenario for the prediction of the Fermi-LAT data is provided by the radio knot closest to the central engine. For that to be the case, a magnetic field strength of ∼1 mG is needed as well as a strong spectral softening of the relativistic electron distribution at (1 − 10) GeV. However, we show that neither such a weak magnetic field strength, nor such a strong softening is expected for that knot. A possible explanation for the ∼10 GeV γ-rays could potentially be provided by hadronic pion production in case of a gas density ≳104 cm−3. Nonetheless, this process is not found to contribute significantly to the low-energy end of the Fermi-LAT range. We conclude that the emission sites in the jet are not sufficient to explain the γ-rays across the whole Fermi-LAT energy band. [ABSTRACT FROM AUTHOR]

Context. The presence of clouds during observations with Imaging Atmospheric Cherenkov Telescopes can strongly affect the performance of the instrument due to additional absorption of light and scattering of light beyond the field of view of the instrument. If not corrected for, the presence of clouds leads to increased systematic errors in the results. Aims. One approach to correct for the effects of clouds is to include clouds in Monte Carlo simulations to produce models for primary particle classification, and energy and direction estimation. However, this method is challenging due to the dynamic nature of cloudy conditions and requires extensive computational resources. The second approach focuses on correcting the data themselves for cloud effects, which allows the use of standard simulations. However, existing corrections often prioritise the limitation of systematic errors without optimising overall performance. By correcting the data already at the image level, it is possible to improve event reconstruction without the need for specialised simulations. Methods. We introduce a novel analysis method based on a geometrical model that can correct the data already at the image level given a vertical transmission profile of a cloud. Using Monte Carlo simulations of an array of four of the Large-Sized Telescopes of the Cherenkov Telescope Array, we investigated the effect of the correction on the image parameters and the performance of the system. We compared the data correction at the camera level with the use of dedicated simulations for clouds with different transmissions and heights. Results. The proposed method efficiently corrects the extinction of light in clouds, eliminating the need for dedicated simulations. Evaluation using Monte Carlo simulations demonstrates improved gamma-ray event reconstruction and overall system performance. [ABSTRACT FROM AUTHOR]

Aims. We present the first multiwavelength study of Mrk 501 that contains simultaneous very-high-energy (VHE) γ-ray observations and X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). Methods. We used radio-to-VHE data from a multiwavelength campaign carried out between March 1, 2022, and July 19, 2022 (MJD 59639 to MJD 59779). The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and several other instruments that cover the optical and radio bands to complement the IXPE pointings. We characterized the dynamics of the broadband emission around the X-ray polarization measurements through its multiband fractional variability and correlations, and compared changes observed in the polarization degree to changes seen in the broadband emission using a multi-zone leptonic scenario. Results. During the IXPE pointings, the VHE state is close to the average behavior, with a 0.2–1 TeV flux of 20%–50% of the emission of the Crab Nebula. Additionally, it shows low variability and a hint of correlation between VHE γ-rays and X-rays. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings, taken in March 2022 (MJD 59646 to 59648 and MJD 59665 to 59667), with a synchrotron peak frequency > 1 keV. For the third IXPE pointing, in July 2022 (MJD 59769 to 59772), the synchrotron peak shifts toward lower energies and the optical/X-ray polarization degrees drop. All three IXPE epochs show an atypically low Compton dominance in the γ-rays. The X-ray polarization is systematically higher than at lower energies, suggesting an energy stratification of the jet. While during the IXPE epochs the polarization angles in the X-ray, optical, and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. Such results further support the hypothesis of an energy-stratified jet. We modeled broadband spectra taken simultaneous to the IXPE pointings, assuming a compact zone that dominates in the X-rays and the VHE band, and an extended zone stretching farther downstream in the jet that dominates the emission at lower energies. NuSTAR data allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change in the magnetization and/or the emission region size, which directly connects the shift in the synchrotron peak to lower energies with the drop in the polarization degree. [ABSTRACT FROM AUTHOR]

With almost 15 yr of data taken by the Large Area Telescope (LAT) onboard the Fermi satellite we discovered an extended source of GeV emission in the region of the very high energy (TeV) source 1LHAASO J1945+2424. This TeV source is more extended than the LAT source. The spectrum of the GeV emission is hard (with a photon spectral index ∼1.5) and connects smoothly with that of the TeV source, indicating a likely common origin. In order to explain the origin of the γ-rays we explore scenarios that are typically used for supernova remnants (SNRs) and pulsar wind nebulae (PWNe). For an SNR with a single particle population, a leptonic particle distribution in the form of a broken power law with a break energy of ∼3.7 TeV explains the spectra well, while in the hadronic scenario a simple power law with a hard spectral index of ∼1.64 is necessary. In the PWN scenario, reasonable parameters are obtained for a source age of 10 kyr and current pulsar spin-down luminosity of ∼1034 erg s−1. [ABSTRACT FROM AUTHOR]

Based on the expected population of core collapse supernova remnants and the huge number of detected pulsars in the Galaxy, still representing only a fraction of the real population, pulsar wind nebulae are likely to constitute one of the largest classes of extended Galactic sources in many energy bands. For simple evolutionary reasons, the majority of the population is made of evolved systems, whose detection and identification are complicated by their reduced luminosity, the possible lack of X-ray emission (that fades progressively away with the age of the pulsar), and by their modified morphology with respect to young systems. Nevertheless they have gained renewed attention in recent years, following the detection of misaligned X-ray tails protruding from an increasing number of nebulae created by fast moving pulsars, and of extended TeV halos surrounding aged systems. Both these features are clear signs of an efficient escape of particles, with energy close to the maximum acceleration limit of the pulsar. Here we discuss the properties of those evolved systems and what we have understood about the process of particle escape, and the formation of observed features. [ABSTRACT FROM AUTHOR]

Context. The number of satellites launched into Earth's orbit has almost tripled in the last three years (to over 4000) due to the increasing commercialisation of space. Multiple satellite constellations, consisting of over 400 000 individual satellites, have either been partially launched or are proposed for launch in the near future. Many of these satellites are highly reflective, resulting in a high optical brightness that affects ground-based astronomical observations. Despite this caveat, the potential effect of these satellites on gamma-ray-observing Imaging Atmospheric Cherenkov Telescopes (IACTs) has largely been assumed to be negligible due to their nanosecond-scale integration times. However, this assumption has not been verified to date. Aims. As IACTs are sensitive to optical wavelength light, we aim to identify satellite trails in data taken by the High Energy Stereoscopic System (H.E.S.S.) IACT array. In particular, this study is aimed at quantifying the potential effects on data quality and extensive air shower event classification and reconstruction. Methods. Using night sky background measurements from H.E.S.S., we determined which observation times and pointing directions are affected most by these satellite trails. We then evaluated their impact on the standard Hillas parameter variables used for event analysis. Results. Due to the brightest trails, false trigger events can occur, however, for most modern analyses, the effect on astronomical results will be minimal. We observe a mild increase in the rate of trail detections over time (approximately doubling in three years), which is partially correlated with the number of satellite launches. Overall, the fraction of H.E.S.S. data affected (~0.2% of dark time observations) is currently minimal. We note that these trails could still have a non-negligible effect on future Cherenkov Telescope Array observations if advanced analysis techniques designed to lower the energy threshold of the instrument are applied. [ABSTRACT FROM AUTHOR]

Aims. Since it launched on 22 March 2021, the 1U-sized CubeSat GRBAlpha operates and collects scientific data on high-energy transients, making it the smallest astrophysical space observatory to date. GRBAlpha is an in-orbit demonstration of a gamma-ray burst (GRB) detector concept suitably small to fit into a standard 1U volume. As was demonstrated in a companion paper, GRBAlpha adds significant value to the scientific community with accurate characterization of bright GRBs, including the recent outstanding event of GRB 221009A. Methods. The GRB detector is a 75 × 75 × 5 mm CsI(Tl) scintillator wrapped in a reflective foil (ESR) read out by an array of SiPM detectors, multi-pixel photon counters by Hamamatsu, driven by two separate redundant units. To further protect the scintillator block from sunlight and protect the SiPM detectors from particle radiation, we applied a multi-layer structure of Tedlar wrapping, anodized aluminium casing, and a lead-alloy shielding on one edge of the assembly. The setup allows observations of gamma radiation within the energy range of 70–890 keV with an energy resolution of ~30%. Results. Here, we summarize the system design of the GRBAlpha mission, including the electronics and software components of the detector, some aspects of the platform, and the current semi-autonomous operations. In addition, details are given about the raw data products and telemetry in order to encourage the community to expand the receiver network for our initiatives with GRBAlpha and related experiments. [ABSTRACT FROM AUTHOR]

Neutrino emission from-ray bursts (GRBs) has been sought for a long time, and stringent limits on the most accredited GRB emission models have been obtained from IceCube. Multiwavelength GRB observations of the last decades have improved our knowledge of the GRB emission parameters, such as the Lorentz factor and the luminosity, which can vary from one GRB to another by several orders of magnitude. Empirical correlations among such parameters have been identified during the prompt phase, with direct implications on GRB models. In this work, we use the PSLab open-access code, developed for IceCube data analyses, to search for individual neutrino emission from the prompt and afterglow phases of selected GRBs, and for stacking emission from the ensemble of such GRBs. For the afterglow phase, we focus on GRBs with X-ray flares and plateaus in particular. While past stacking searches assumed the same GRB fluence at Earth, we present a stacking scheme based on physically motivated GRB weights. Moreover, we conceive a new methodology for the prompt phase that uses the empirical correlations to infer the GRB luminosity and Lorentz factor, when redshift measurements are not available. We do not observe any significant neutrino excess. Hence, we set constraints on the GRB neutrino fluxes and on relevant GRB parameters, including the magnetic field in the jet. Notably, the baryon loading is found to be less than ten for typical GRB jets. [ABSTRACT FROM AUTHOR]

Context. Young massive star clusters (YMCs) have increasingly become the focus of discussions on the origin of galactic cosmic rays (CRs). The proposition that CRs are accelerated inside superbubbles (SBs) blown by the strong winds of these clusters avoids issues faced by the standard paradigm of acceleration at supernova remnant shocks. Aims. We provide an interpretation of the latest TeV γ-ray observations of the region around the YMC Westerlund 1 taken with the High Energy Stereoscopic System (H.E.S.S.) in terms of diffusive shock acceleration at the cluster wind termination shock, taking the spectrum and morphology of the emission into account. As Westerlund 1 is a prototypical example of a YMC, this study is relevant to the general question about the role of YMCs for the Galactic CR population. Methods. We generated model γ-ray spectra, characterised particle propagation inside the SB based on the advection, diffusion, and cooling timescales, and constrained key parameters of the system. We considered hadronic emission from proton-proton interaction and subsequent pion decay and leptonic emission from inverse Compton scattering on all relevant photon fields, including the cosmic microwave background, diffuse and dust-scattered starlight, and the photon field of Westerlund 1 itself. The effect of the magnetic field on cooling and propagation is discussed. Klein-Nishina effects are found to be important in determining the spectral evolution of the electron population. Results. A leptonic origin of the bulk of the observed γ-rays is preferable. The model is energetically plausible, consistent with the presence of a strong shock, and allows for the observed energy-independent morphology. The hadronic model faces two main issues: confinement of particles to the emission region, and an unrealistic energy requirement. [ABSTRACT FROM AUTHOR]

Context. Particle-accelerating colliding-wind binaries (PACWBs) are systems that are formed by two massive and hot stars and produce nonthermal radiation. The key elements of these systems are fast winds and the shocks that they create when they collide. Binaries with nonaccreting young pulsars have also been detected as nonthermal emitters, again as a consequence of the wind–wind interaction. Black holes might produce nonthermal radiation by this mechanism if they accrete at super-Eddington rates. In such cases, the disk is expected to launch a radiation-driven wind, and if this wind has an equatorial component, it can collide with the companion star yielding a PACWB. These systems are supercritical colliding wind binaries. Aims. We aim to characterize the particle acceleration and nonthermal radiation produced by the collision of winds in binary systems composed of a superaccreting black hole and an early-type star. Methods. We estimated the terminal velocity of the disk-driven wind by calculating the spatial distribution of the radiation fields and their effect on disk particles. We then found the location of the wind collision region and calculated the timescales of energy gain and losses of relativistic particles undergoing diffusive particle acceleration. With this information, we were able to compute the associated spectral energy distribution of the radiation. We calculated a number of specific models with different parameters to explore this scenario. Results. We find that the interaction of winds can produce nonthermal emission from radio up to tens of GeV, with luminosities in the range of ∼1033–1035 erg s−1, which for the most part are contributed by electron synchrotron and inverse Compton radiation. Conclusions. We conclude that supercritical colliding wind binaries, such as some ultraluminous X-ray sources and some Galactic X-ray binaries, are capable of accelerating cosmic rays and producing nonthermal electromagnetic emission from radio to γ-rays, in addition to the thermal components. [ABSTRACT FROM AUTHOR]

Context. Certain types of supernova remnants (SNRs) in our Galaxy are assumed to be PeVatrons, capable of accelerating cosmic rays (CRs) to ~ PeV energies. However, conclusive observational evidence for this has not yet been found. The SNR G106.3+2.7, detected at 1–100 TeV energies by different γ-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape, which can be divided into a head and a tail region with different physical conditions. However, in which region the 100 TeV emission is produced has not yet been identified because of the limited position accuracy and/or angular resolution of existing observational data. Additionally, it remains unclear as to whether the origin of the γ-ray emission is leptonic or hadronic. Aims. With the better angular resolution provided by new MAGIC data compared to earlier γ-ray datasets, we aim to reveal the acceleration site of PeV particles and the emission mechanism by resolving the SNR G106.3+2.7 with 0.1° resolution at TeV energies. Methods. We observed the SNR G106.3+2.7 using the MAGIC telescopes for 121.7 h in total – after quality cuts – between May 2017 and August 2019. The analysis energy threshold is ~0.2 TeV, and the angular resolution is 0.07−0.1°. We examined the γ-ray spectra of different parts of the emission, whilst benefitting from the unprecedented statistics and angular resolution at these energies provided by our new data. We also used measurements at other wavelengths such as radio, X-rays, GeV γ-rays, and 10 TeV γ-rays to model the emission mechanism precisely. Results. We detect extended γ-ray emission spatially coincident with the radio continuum emission at the head and tail of SNR G106.3+2.7. The fact that we detect a significant γ-ray emission with energies above 6.0 TeV from only the tail region suggests that the emissions above 10 TeV detected with air shower experiments (Milagro, HAWC, Tibet ASγ and LHAASO) are emitted only from the SNR tail. Under this assumption, the multi-wavelength spectrum of the head region can be explained with either hadronic or leptonic models, while the leptonic model for the tail region is in contradiction with the emission above 10 TeV and X-rays. In contrast, the hadronic model could reproduce the observed spectrum at the tail by assuming a proton spectrum with a cutoff energy of ~1 PeV for that region. Such high-energy emission in this middle-aged SNR (4−10 kyr) can be explained by considering a scenario where protons escaping from the SNR in the past interact with surrounding dense gases at present. Conclusions. The γ-ray emission region detected with the MAGIC telescopes in the SNR G106.3+2.7 is extended and spatially coincident with the radio continuum morphology. The multi-wavelength spectrum of the emission from the tail region suggests proton acceleration up to ~PeV, while the emission mechanism of the head region could either be hadronic or leptonic. [ABSTRACT FROM AUTHOR]

Context. The Crab Nebula emits exceptionally bright non-thermal radiation across the entire wavelength range from the radio to the most energetic photons. So far, the underlying physical model of a relativistic wind from the pulsar terminating in a hydrodynamic standing shock has remained fairly unchanged since the early 1970s when it was first introduced. One of the predictions of this model is an increase in the toroidal magnetic field downstream from the shock where the flow velocity drops quickly with increasing distance until it reaches its asymptotic value, matching the expansion velocity of the nebula. Aims. The magnetic field strength in the nebula is poorly known. Using the recent measurements of the spatial extension and improved spectroscopy of the gamma-ray nebula, it has become –for the first time – feasible to determine in a robust way both the strength as well as the radial dependence of the magnetic field in the downstream flow. Methods. In this work, we introduce a detailed radiative model which was used to calculate the emission from non-thermal electrons (synchrotron and inverse Compton) as well as from thermal dust present in the Crab Nebula in a self-consistent way to compare it quantitatively with observational data. Special care was given to the radial dependence of the electron and seed field density. Results. The radiative model was used to estimate the parameters related to the electron populations responsible for radio and optical/X-ray synchrotron emission. In this context, the mass of cold and warm dust was determined. A combined fit based upon a χ2 minimisation successfully reproduced the complete data set used. For the best-fitting model, the energy density of the magnetic field dominates over the particle energy density up to a distance of ≈1.3 rs (rs: distance of the termination shock from the pulsar). The very high energy (VHE: E > 100 GeV) and ultra-high energy (UHE: E > 100 TeV) gamma-ray spectra set the strongest constraints on the radial dependence of the magnetic field, favouring a model where B(r) = (264 ± 9) μG(r/rs)−0.51 ± 0.03. For a collection of VHE measurements during epochs of higher hard X-ray emission, a significantly different solution B(r) = (167 ± 5) μG(r/rs)−0.29(+0.03, −0.06) is found. Conclusions. The high energy (HE: E > 100 MeV) and VHE gamma-ray observations of the Crab Nebula lift the degeneracy of the synchrotron emission between particle and magnetic field energy density. The reconstructed magnetic field and its radial dependence indicates a ratio of Poynting to kinetic energy flux σ ≈ 0.1 at the termination shock, which is ≈30 times larger than estimated up to now. Consequently, the confinement of the nebula would require additional mechanisms to slow the flow down through, for example, excitation of small-scale turbulence with possible dissipation of the magnetic field. [ABSTRACT FROM AUTHOR]

Pulsar wind nebula (PWN) HESS J1825-137 is one of the largest and brightest gamma-ray sources in Milky Way. It was found that the morphology of the nebula strongly depends on energy and the photon spectrum varies obviously with distance from the pulsar in the >0.1 TeV energy band. In the frame of a spatially-dependent particle transport model, we investigate the TeV gamma-ray observed features of HESS J1825-137 by using the Markov chain Monte Carlo method. Our results are as follows: (i) the observed features can be reasonably explained by the competition between the transport (advection and diffusion) and energy cooling losses (adiabatic, synchrotron, and inverse Compton) of the high energy electrons, and that a parameter ξ ≡ τp/τC is introduced to describe this competition, which is the function of electron energy and position within the nebula, where τ p and τ C are total transport and energy cooling timescales, respectively; (ii) for the high energy electrons (E e > 0.1 TeV) within the nebula, the particle diffusion dominates over its advection, and the radiative cooling mostly dominates over adiabatic cooling; and (iii) the change of the radial extent with photon energy E γ can be expressed as |$R\propto E^{-\beta }_\gamma$|⁠ , where the index β is a function of E γ and the Klein−Nishina effect has a very important role for the variation of β with E γ. From the modelling results, the initial rotation period, braking index, and age of the pulsar are predicted to be 12.16 ms, 1.79, and 44 kyr, respectively. [ABSTRACT FROM AUTHOR]

Cosmological models and their corresponding parameters are widely debated because of the current discrepancy between the results of the Hubble constant, H 0, obtained by SNe Ia, and the Planck data from the cosmic microwave background radiation. Thus, considering high redshift probes like gamma-ray bursts (GRBs) is a necessary step. However, using GRB correlations between their physical features to infer cosmological parameters is difficult because GRB luminosities span several orders of magnitude. In our work, we use a three-dimensional relation between the peak prompt luminosity, the rest-frame time at the end of the X-ray plateau, and its corresponding luminosity in X-rays: the so-called 3D Dainotti fundamental plane relation. We correct this relation by considering the selection and evolutionary effects with a reliable statistical method, obtaining a lower central value for the intrinsic scatter, σ int = 0.18 ± 0.07 (47.1 per cent) compared to previous results, when we adopt a particular set of GRBs with well-defined morphological features, called the platinum sample. We have used the GRB fundamental plane relation alone with both Gaussian and uniform priors on cosmological parameters and in combination with SNe Ia and BAO measurements to infer cosmological parameters like H 0, the matter density in the universe (ΩM), and the dark energy parameter w for a w CDM model. Our results are consistent with the parameters given by the Lambda cold dark matter model but with the advantage of using cosmological probes detected up to z  = 5, much larger than the one observed for the furthest SNe Ia. [ABSTRACT FROM AUTHOR]

MAXI J1820+070 is a low-mass X-ray binary with a black hole (BH) as a compact object. This binary underwent an exceptionally bright X-ray outburst from 2018 March to October, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 h of observations of the MAXI J1820+070 outburst with the H.E.S.S. MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi -LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy (HE) gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi -LAT upper limits, and closer to them for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ∼500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential HE and very-HE gamma-ray emitting region should be located at a distance from the BH ranging between 1011 and 1013 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA. [ABSTRACT FROM AUTHOR]

Based on the expected population of core collapse supernova remnants and the huge number of detected pulsars in the Galaxy, still representing only a fraction of the real population, pulsar wind nebulae are likely to constitute one of the largest classes of extended Galactic sources in many energy bands. For simple evolutionary reasons, the majority of the population is made of evolved systems, whose detection and identification are complicated by their reduced luminosity, the possible lack of X-ray emission (that fades progressively away with the age of the pulsar), and by their modified morphology with respect to young systems. Nevertheless they have gained renewed attention in recent years, following the detection of misaligned X-ray tails protruding from an increasing number of nebulae created by fast moving pulsars, and of extended TeV halos surrounding aged systems. Both these features are clear signs of an efficient escape of particles, with energy close to the maximum acceleration limit of the pulsar. Here we discuss the properties of those evolved systems and what we have understood about the process of particle escape, and the formation of observed features.

Aims. We aim to search for a hidden leptonic accelerator, such as a high-Ė pulsar, associated with the unidentified TeV object HESS J1702−420A. Methods. We carried out a 72 ks X-ray observation with the XMM-Newton satellite and analyzed the resulting data jointly with the publicly available HESS spectral energy distribution (SED) results to derive constraints on the leptonic contribution to the TeV emission of HESS J1702−420A. A set of scripts dedicated to the multi-wavelength modeling of X-ray and γ-ray data, based on Gammapy, Naima, and Xspec, has been developed in the context of this work and made publicly available along with this paper. Results. No object clearly associated with HESS J1702−420A was found in the XMM-Newton data. After excluding the unidentified object Suzaku src B as a possible X-ray counterpart and classifying it as a new cataclysmic variable source candidate, we derived strict upper limits on the level of diffuse X-ray emission in the HESS J1702−420A region: F(2 − 10 keV)≲5.4  × 10−5 keV  cm−2 s−1 at 2σ (≈95.5%) confidence level. A tight constraint on the magnetic field was derived, under a one-zone leptonic scenario, by jointly fitting the XMM-Newton spectra and the HESS SED: B ≲ 1.45 μG at 2σ level. We additionally report the serendipitous discovery of a new extended X-ray source with a hard spectral index of 1.99 ± 0.45, named XMMU J170147.3−421407 which is likely Galactic. Its classification as a high-speed runaway pulsar wind nebula (PWN), possibly associated with HESS J1702−420A, is not obvious but cannot be ruled out either. Conclusions. The hard γ-ray object HESS J1702−420A remains unidentified, but the absence of a clear X-ray counterpart strongly challenges simple leptonic scenarios. The only remaining possible leptonic counterpart for HESS J1702−420A appears to be a newly discovered X-ray source with extended morphology and hard spectral index, which may be a PWN powered by a high-speed runaway pulsar. [ABSTRACT FROM AUTHOR]

The Crab Nebula is the brightest TeV gamma-ray source in the sky and has been used for the past 25 years as a reference source in TeV astronomy, for calibration and verification of new TeV instruments. The High Altitude Water Cherenkov Observatory (HAWC), completed in early 2015, has been used to observe the Crab Nebula at high significance across nearly the full spectrum of energies to which HAWC is sensitive. HAWC is unique for its wide field of view, nearly 2 sr at any instant, and its high-energy reach, up to 100 TeV. HAWC's sensitivity improves with the gamma-ray energy. Above ∼1 TeV the sensitivity is driven by the best background rejection and angular resolution ever achieved for a wide-field ground array. We present a time-integrated analysis of the Crab using 507 live days of HAWC data from 2014 November to 2016 June. The spectrum of the Crab is fit to a function of the form φ(E)= φ0(E/E0)-α-β In(E/E0). The data is well fitted with values of α = 2.63 ±0.03, β = 0.15 ±0.03, and log10(φ0cm2s TeV)=-12.60±0.02 when E 0 is fixed at 7 TeV and the fit applies between 1 and 37 TeV. Study of the systematic errors in this HAWC measurement is discussed and estimated to be ±50% in the photon flux between 1 and 37 TeV. Confirmation of the Crab flux serves to establish the HAWC instrument's sensitivity for surveys of the sky. The HAWC all-sky survey will be the deepest survey of the northern sky ever conducted in the multi-TeV band.

The last measurement of the diffuse emission spectrum of the Milky Way in the megaelectronvolt (MeV) photon energy range was performed by CGRO/COMPTEL more than 20 yr ago. We report a new analysis with the spectrometer SPI aboard INTEGRAL in the band 0.5–8.0 MeV, finally superseding the signal-to-noise ratio (S/N) of the historic observations. This is possible thanks to an elaborate instrumental background model and careful considerations of the selected data, which are strongly affected by solar activity. We base our analysis on energy-dependent spatial template fitting in a region of Δl ×Δb = 95° ×95° around the Galactic centre. Our flux estimates are consistent with COMPTEL measurements and show no 'MeV bump' The spectrum follows a power-law shape with index −1.39 ± 0.09stat ± 0.10syst and an integrated flux of (5.7 ± 0.8stat ± 1.7syst) x 10−8 erg cm−2 s−1 between 0.5 and 8.0 MeV. We find that cosmic-ray electrons and propagation models consistent with the latest Fermi/LAT, Voyager 1, and AMS-02 data are broadly in agreement with the inferred inverse Compton spectral shape. However, a mismatch of a factor of 2–3 in normalisation with respect to baseline expectations may point to enhanced target photon densities and/or electron source spectra in the inner Galaxy, slightly modified diffusion properties, or the presence of an unresolved population of MeV γ-ray sources. [ABSTRACT FROM AUTHOR]

A survey of the inner Galaxy region of Galactic longitude l ∈[+15,+50] and latitude b ∈ [-4, +4] is performed using one-third of the High Altitude Water Cherenkov Observatory, operated during ts construction phase. To address the ambiguities arising from unresolved sources in the data, we use a maximum likelihood technique to identify point source candidates. Ten sources and candidate sources are identified in this analysis. Eight of these are associated with known TeV sources but not all have ifferential fluxes that are compatible with previous measurements. Three sources are detected with significances >5 sigma; after accounting for statistical trials, and are associated with known TeV sources.

BL Lacs are one subclass of blazars with highly energetic γ -ray emission, which is strongly boosted by a relativistic beaming effect. The latest catalogue of the 10 years of Fermi/LAT data (Abdollahi et al. 2020) and the γ -ray Doppler factors in Pei et al. (2020) provide us with a large BL Lac sample to study their jet emission morphologies and intrinsic properties. In this paper, we collected a sample of 294 Fermi BL Lacs and probed the correlations between the γ -ray emissions and luminosity distances. Our analyses give following conclusions: (1) the observed γ -ray emissions are really boosted by the γ -ray Doppler factor, and the intrinsic γ -ray emissions are closely correlated with luminosity distances. (2) the morphology of jet emissions for HBLs may be continuous, while that for IBLs may be the case of a moving sphere in the γ -ray bands. [ABSTRACT FROM AUTHOR]

Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV γ-rays and TeV neutrinos on a timescale of several months. We perform the first systematic search for γ-ray emission in Fermi Large Area Telescope data in the energy range from to from the ensemble of 147 SNe Type IIn exploding in a dense CSM. We search for a γ-ray excess at each SNe location in a one-year time window. In order to enhance a possible weak signal, we simultaneously study the closest and optically brightest sources of our sample in a joint-likelihood analysis in three different time windows (1 year, 6 months, and 3 months). For the most promising source of the sample, SN 2010jl (PTF 10aaxf), we repeat the analysis with an extended time window lasting 4.5 years. We do not find a significant excess in γ-rays for any individual source nor for the combined sources and provide model-independent flux upper limits for both cases. In addition, we derive limits on the γ-ray luminosity and the ratio of γ-ray-to-optical luminosity ratio as a function of the index of the proton injection spectrum assuming a generic γ-ray production model. Furthermore, we present detailed flux predictions based on multi-wavelength observations and the corresponding flux upper limit at a 95% confidence level (CL) for the source SN 2010jl (PTF 10aaxf).

The first limits on the prompt emission from the long gamma-ray burst (GRB) 130427A in the >100 GeV energy band are reported. GRB 130427A was the most powerful burst ever detected with a redshift z ≲ 0.5 and featured the longest lasting emission above 100 MeV. The energy spectrum extends at least up to 95 GeV, clearly in the range observable by the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory, a new extensive air shower detector currently under construction in central Mexico. The burst occurred under unfavorable observation conditions, low in the sky and when HAWC was running 10% of the final detector. Based on the observed light curve at MeV-GeV energies, eight different time periods have been searched for prompt and delayed emission from this GRB. In all cases, no statistically significant excess of counts has been found and upper limits have been placed. It is shown that a similar GRB close to zenith would be easily detected by the full HAWC detector, which will be completed soon. The detection rate of the full HAWC detector may be as high as one to two GRBs per year. A detection could provide important information regarding the high energy processes at work and the observation of a possible cut-off beyond the Fermi Large Area Telescope energy range could be the signature of gamma-ray absorption, either in the GRB or along the line of sight due to the extragalactic background light.

In this work, we re-analyze the Fermi-LAT observation of NGC 1275 to search for axion-like particle (ALP) effects and constrain ALP parameters. Instead of fitting the observed spectrum with ALP models, we adopt an alternative method for the analysis of this source which calculates the irregularity of the spectrum. With the newly used method, we find no spectral oscillation for the NGC 1275 and rule out couplings gaγ>3×10−12GeV−1 around ALP mass of ma∼1 neV at 95% confidence level, which is more stringent than the previous results. We also show that the constraints can be further improved by combining the observation of PKS 2155-304. We suggest that with more sources taken into account, we could obtain a much wider exclusion region.

γ -ray is a unique probe for extreme events in the universe. Detecting the γ -ray provides an important opportunity to understand the composition of universe, the evolution of stars, the origin of cosmic rays, etc. γ -ray astronomy involves in various frontier scientific issues, and the observed energy spectrum spans over a wide range from a few hundreds of keV to a few hundreds of TeV. Different γ -ray telescopes are in need for different scientific goals and spectral bands. In this work, 5 kinds of space- and ground-based γ -ray observing techniques were summarized, including the Coded-aperture telescopes, Compton telescopes, Pair-production telescopes, Imaging Atmospheric Cherenkov Telescopes, and Extensive Air Shower Arrays. The progress in γ -ray astronomy in the past 70 years, motivated by the observation capability, was reviewed. Great achievements have been made in the high-energy domain and very-high-energy domain, while because of the limited missions conducted, as well as a lower sensitivity comparing with other domains, discoveries in low- and medium-energy are few, and due to the high observation difficulty, as well as the late start, relevant scientific yields in ultra- and extremely-high energy are limited. Moreover, the future planned missions and capabilities of the γ -ray telescopes and their possible scientific outputs were discussed. Among these missions, low- and medium-energy space telescopes e-ASTROGAM (enhanced-ASTROGAM), AMEGO (All-sky Medium Energy Gamma-ray Observatory), and very-high-energy ground-based arrays LHAASO (Large High Altitude Air Shower Observatory) and CTA (Cherenkov Telescope Array) greatly improve sensitivity than their corresponding last generation, thus expect very likely to further expand our knowledge on the γ -universe. [ABSTRACT FROM AUTHOR]

We present a machine learning (ML) based method for automated detection of Gamma-Ray Burst (GRB) candidate events in the range 60–250 keV from the AstroSat Cadmium Zinc Telluride Imager data. We use density-based spatial clustering to detect excess power and carry out an unsupervised hierarchical clustering across all such events to identify the different light curves present in the data. This representation helps us to understand the instrument's sensitivity to the various GRB populations and identify the major non-astrophysical noise artefacts present in the data. We use Dynamic Time Warping (DTW) to carry out template matching, which ensures the morphological similarity of the detected events with known typical GRB light curves. DTW alleviates the need for a dense template repository often required in matched filtering like searches. The use of a similarity metric facilitates outlier detection suitable for capturing previously unmodelled events. We briefly discuss the characteristics of 35 long GRB candidates detected using the pipeline and show that with minor modifications such as adaptive binning, the method is also sensitive to short GRB events. Augmenting the existing data analysis pipeline with such ML capabilities alleviates the need for extensive manual inspection, enabling quicker response to alerts received from other observatories such as the gravitational-wave detectors. [ABSTRACT FROM AUTHOR]

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) γ -ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE;E >100 MeV) γ -rays, poses intriguing questions on the location of the γ -ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a γ -ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two γ -ray peaks, while no clear connection was observed between the X-ray and γ -ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing γ -ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of

We extend our previous model-independent determination of the intergalactic background light, based purely on galaxy survey data, out to a wavelength of 5 μm. Our approach enables us to constrain the range of photon densities, based on the uncertainties from observationally determined luminosity densities and colors. We further determine a 68% confidence upper and lower limit on the opacity of the universe to γ-rays up to energies of 1.6/(1 + z) TeV. A comparison of our lower limit redshift-dependent opacity curves to the opacity limits derived from the results of both ground-based air Cerenkov telescope and Fermi-LAT observations of PKS 1424+240 allows us to place a new upper limit on the redshift of this source, independent of IBL modeling. © 2014. The American Astronomical Society. All rights reserved.

Observations of occultations of bright γ-ray sources by the Sun may reveal predicted pair halos around blazars and/or new physics, such as, e.g., hypothetical light dark matter particles-axions. We use Fermi Gamma-Ray Space Telescope (Fermi) data to analyze four occultations of blazar 3C 279 by the Sun on October 8 each year from 2008 to 2011. A combined analysis of the observations of these occultations allows a point-like source at the position of 3C 279 to be detected with significance of ≈3σ, but does not reveal any significant excess over the flux expected from the quiescent Sun. The likelihood ratio test rules out complete transparency of the Sun to the blazar γ-ray emission at a 3σ confidence level.

TeV-flaring activity with timescales as short as tens of minutes and an orphan TeV flare have been observed from the blazar Markarian 421 (Mrk 421). The TeV emission from Mrk 421 is believed to be produced by leptonic synchrotron self-Compton (SSC) emission. In this scenario, correlations between the X-ray and the TeV fluxes are expected, TeV orphan flares are hardly explained, and the activity (measured as duty cycle) of the source at TeV energies is expected to be equal to or less than that observed in X-rays if only SSC is considered. To estimate the TeV duty cycle of Mrk 421 and to establish limits on its variability at different timescales, we continuously observed Mrk 421 with the Milagro observatory. Mrk 421 was detected by Milagro with a statistical significance of 7.1 standard deviations between 2005 September 21 and 2008 March 15. The observed spectrum is consistent with previous observations by VERITAS. We estimate the duty cycle of Mrk 421 for energies above 1 TeV for different hypotheses of the baseline flux and for different flare selections and we compared our results with the X-ray duty cycle estimated by Resconi et al. The robustness of the results is discussed. © 2014. The American Astronomical Society. All rights reserved..

We present the nucleosynthesis of magneto-rotational supernovae (MR-SNe) including neutrino-driven and magneto-rotational-driven ejecta based, for the first time, on 2D simulations with accurate neutrino transport. The models analysed here have different rotation and magnetic fields, allowing us to explore the impact of these two key ingredients. The accurate neutrino transport of the simulations is critical to analyse the slightly neutron-rich and proton-rich ejecta that are similar to the, also neutrino-driven, ejecta in standard supernovae. In the model with strong magnetic field, the r-process produces heavy elements up to the third r-process peak (A ∼ 195), in agreement with previous works. This model presents a jet-like explosion with proton-rich jets surrounded by neutron-rich material where the r-process occurs. We have estimated a lower limit for 56Ni of |$2.5\times 10^{-2} \, \mathrm{M}_\odot$|⁠ , which is still well below the expected hypernova value. Longer simulations including the accretion disc evolution are required to get a final prediction. In addition, we have found that the late evolution is critical in a model with weak magnetic field in which late-ejected neutron-rich matter produces elements up to the second r-process peak. Even if we cannot yet provide conclusions for hypernova nucleosynthesis, our results agree with observations of old stars and radioactive isotopes in supernova remnants. This makes MR-SNe a good additional scenario to neutron star mergers for the synthesis of heavy elements and brings us closer to understand their origin and the role of MR-SNe in the early Galaxy nucleosynthesis. [ABSTRACT FROM AUTHOR]

Pulsar wind nebulae (PWNe) have been established as the most populous class of TeV γ-ray emitters. Since launch, the Fermi Large Area Telescope (LAT) has identified five high-energy (100 MeV < E < 100 GeV) γ-ray sources as PWNe and detected a large number of PWN candidates, all powered by young and energetic pulsars. The wealth of multi-wavelength data available and the new results provided by Fermi-LAT give us an opportunity to find new PWNe and to explore the radiative processes taking place in known ones. The TeV γ-ray unidentified (UNID) sources are the best candidates for finding new PWNe. Using 45 months of Fermi-LAT data for energies above 10 GeV, an analysis was performed near the position of 58 TeV PWNe and UNIDs within 5° of the Galactic plane to establish new constraints on PWN properties and find new clues on the nature of UNIDs. Of the 58 sources, 30 were detected, and this work provides their γ-ray fluxes for energies above 10 GeV. The spectral energy distributions and upper limits, in the multi-wavelength context, also provide new information on the source nature and can help distinguish between emission scenarios, i.e., between classification as a pulsar candidate or as a PWN candidate. Six new GeV PWN candidates are described in detail and compared with existing models. A population study of GeV PWN candidates as a function of the pulsar/PWN system characteristics is presented. © 2013. The American Astronomical Society. All rights reserved.