Search

Energetic GeV photons expected from the closest and the most energetic Gamma-ray bursts (GRBs) provide an unique opportunity to study the very-high-energy emission as well as the possible correlations with lower energy bands in realistic GRB afterglow models. In the standard GRB afterglow model, the relativistic homogeneous shock is usually considered to be fully adiabatic, however, it could be partially radiative. Based on the external forward-shock scenario in both stellar wind and constant-density medium. We present a radiative-adiabatic analytical model of the synchrotron self-Compton (SSC) and synchrotron processes considering an electron energy distribution with a power-law index of 1 < p < 2 and 2 $\leq$ p. We show that the SSC scenario plays a relevant role in the radiative parameter $\epsilon$, leading to a prolonged evolution during the slow cooling regime. In a particular case, we derive the Fermi/LAT light curves together with the photons with energies $\geq$ 100 MeV in a sample of nine bursts from the second Fermi/LAT GRB catalog that exhibited temporal and spectral indices with $\geq$ 1.5 and $\approx$ 2, respectively. These events can hardly be described with closure relations of the standard synchrotron afterglow model, and also exhibit energetic photons above the synchrotron limit. We have modeled the multi-wavelength observations of our sample to constrain the microphysical parameters, the circumburst density, the bulk Lorentz factor and the mechanism responsible for explaining the energetic GeV photons., Comment: Accepted in MNRAS. 26 pages, 7 figures, 6 tables

Gamma-ray bursts (GRBs) are one of the most exciting sources that offer valuable opportunities for investigating the evolution of energy fraction given to magnetic fields and particles through microphysical parameters during relativistic shocks. The delayed onset of GeV-TeV radiation from bursts detected by the \textit{Fermi} Large Area Telescope (\textit{Fermi}-LAT) and Cherenkov Telescopes provide crucial information in favor of the external-shock model. Derivation of the closure relations (CRs) and the light curves in external shocks requires knowledge of GRB afterglow physics. In this manuscript, we derive the CRs and light curves in a stratified medium with variations of microphysical parameters of the synchrotron and SSC afterglow model radiated by an electron distribution with a hard and soft spectral index. Using Markov Chain Monte Carlo simulations, we apply the current model to investigate the evolution of the spectral and temporal indexes of those GRBs reported in the Second Gamma-ray Burst Catalog (2FLGC), which comprises 29 bursts with photon energies above 10 GeV and of those bursts (GRB 180720B, 190114C, 190829A and 221009A) with energetic photons above 100 GeV, which can hardly be modeled with the CRs of the standard synchrotron scenario. The analysis shows that i) the most likely afterglow model using synchrotron and SSC emission on the 2FLGC corresponds to the constant-density scenario, and ii) variations of spectral (temporal) index keeping the temporal (spectral) index constant could be associated with the evolution of microphysical parameters, as exhibited in GRB 190829A and GRB 221009A., Comment: 34 pages, 11 Figures, 9 Tables, Accepted in MNRAS

The second {\itshape Fermi}/LAT gamma-ray burst (GRB) catalog (2FLGC) spanning the first decade of operations by the LAT collaboration was recently released. The closure relations of the synchrotron forward shock (FS) model are not able to reproduce a sizeable portion of the afterglow-phase light curves in this collection, indicating that there may be a large contribution from some other mechanism. Recently, synchrotron self-Compton (SSC) light curves from the reverse shock (RS) regions were derived in the thick- and thin-shell regime for a constant-density medium, and it was demonstrated that analytical light curves could explain the~GeV flare observed in several bursts from 2FLGC, including GRB 160509A. Here, we generalise the SSC RS scenario from the constant density to a stratified medium, and show that this contribution helps to describe the early light curves exhibited in some {\itshape Fermi}/LAT-detected bursts. As a particular case, we model a sample of eight bursts that exhibited a short-lasting emission with the synchrotron and SSC model from FS and RS regions, evolving in a stellar-wind environment, constraining the microphysical parameters, the circumburst density, the bulk Lorentz factor, and the fraction of shock-accelerated electrons. We demonstrate that the highest-energy photons can only be described by the SSC from the forward-shock region., Comment: 35 pages, 12 figures, 4 tables, accepted in mnras

Gamma-ray bursts (GRBs) are one of the most promising transient events for studying multi-wavelength observations in extreme conditions. Observation of GeV photons from bursts would provide crucial information on GRB physics, including the off-axis emission. The Second Gamma-ray Burst Catalog (2FLGC) was announced by the Fermi Large Area Telescope (Fermi-LAT) Collaboration. This catalog includes 29 bursts with photon energy higher than 10 GeV. While the synchrotron forward-shock model has well explained the afterglow data of GRBs, photon energies greater than 10 GeV are very difficult to interpret within this framework. To study the spectral and temporal indices of those bursts described in 2FLGC, Fraija et al. (2022a) proposed the closure relations (CRs) of the synchrotron self-Compton (SSC) emitted from an on-axis jet which decelerates in stellar-wind and the constant-density medium. In this paper, we extend the CRs of the SSC afterglow from an on-axis scenario to an off-axis, including the synchrotron afterglow radiation that seems off-axis. In order to investigate the spectral and temporal index evolution of those bursts reported in 2FLGC, we consider the hydrodynamical evolution with energy injection in the adiabatic and radiative regime for an electron distribution with a spectral index of $1

The Galactic gamma-ray diffuse emission (GDE) is emitted by cosmic rays (CRs), ultra-relativistic protons and electrons, interacting with gas and electromagnetic radiation fields in the interstellar medium. Here we present the analysis of TeV diffuse emission from a region of the Galactic Plane over the range in longitude of $l\in[43^\circ,73^\circ]$, using data collected with the High Altitude Water Cherenkov (HAWC) detector. Spectral, longitudinal and latitudinal distributions of the TeV diffuse emission are shown. The radiation spectrum is compatible with the spectrum of the emission arising from a CR population with an "index" similar to that of the observed CRs. When comparing with the \texttt{DRAGON} \textit{base model}, the HAWC GDE flux is higher by about a factor of two. Unresolved sources such as pulsar wind nebulae and TeV halos could explain the excess emission. Finally, deviations of the Galactic CR flux from the locally measured CR flux may additionally explain the difference between the predicted and measured diffuse fluxes.

The Second Gamma-ray Burst Catalog (2FLGC) was announced by the Fermi Large Area Telescope (Fermi-LAT) Collaboration. It includes 29 bursts with photon energy higher than 10 GeV. Gamma-ray burst (GRB) afterglow observations have been adequately explained by the classic synchrotron forward-shock model, however, photon energies greater than 10 GeV from these transient events are challenging, if not impossible, to characterize using this afterglow model. Recently, the closure relations (CRs) of the synchrotron self-Compton (SSC) forward-shock model evolving in a stellar wind and homogeneous medium was presented to analyze the evolution of the spectral and temporal indexes of those bursts reported in 2FLGC. In this work, we provide the CRs of the same afterglow model, but evolving in an intermediate density profile ($\propto {\rm r^{-k}}$) with ${\rm 0\leq k \leq2.5}$, taking into account the adiabatic/radiative regime and with/without energy injection for any value of the electron spectral index. The results show that the current model accounts for a considerable subset of GRBs that cannot be interpreted in either stellar-wind or homogeneous afterglow SSC model. The analysis indicates that the best-stratified scenario is most consistent with ${\rm k=0.5}$ for no-energy injection and ${\rm k=2.5}$ for energy injection., Comment: 13 pages, 3 figures

Blazars have been pointed out as promising high-energy (HE) neutrinos sources, although the mechanism is still under debate. The blazars with a hard-TeV spectrum, which leptonic models can hardly explain, can be successfully interpreted in the hadronic scenarios. Recently, Aguilar et al. proposed a lepto-hadronic two-zone model to explain the multi-wavelength observations of the six best-known extreme BL Lacs and showed that the hadronic component could mainly interpret very-high-energy (VHE) emission. In this work, we apply this hadronic model to describe the VHE gamma-ray fluxes of 14 extreme BL Lacs and estimate the respective HE neutrino flux from charge-pion decay products. Finally, we compare our result with the diffuse flux observed by the IceCube telescope, showing that the neutrino fluxes from these objects are negligible., Comment: Accepted in Journal of High Energy Astrophysics

While the dominant radiation mechanism gamma-ray bursts (GRBs) remains a question of debate, synchrotron emission is one of the foremost candidates to describe the multi-wavelength afterglow observations. As such, it is expected that GRBs should present some degree of polarization across their evolution - presenting a feasible means of probing these bursts' energetic and angular properties. Although obtaining polarization data is difficult due to the inherent complexities regarding GRB observations, advances are being made, and theoretical modeling of synchrotron polarization is now more relevant than ever. In this manuscript, we present the polarization for a fiduciary model where the synchrotron forward-shock emission evolving in the radiative-adiabatic regime is described by a radially stratified off-axis outflow. This is parameterized with a power-law velocity distribution and decelerated in a constant-density and wind-like external environment. We apply this theoretical polarization model for selected bursts presenting evidence of off-axis afterglow emission, including the nearest orphan GRB candidates observed by the Neil Gehrels Swift Observatory and a few Gravitational Wave (GWs) events that could generate electromagnetic emission. In the case of GRB 170817A, we require the available polarimetric upper limits in radio wavelengths to constrain its magnetic field geometry., Comment: In submission. 18 pages, 7 figures, 3 tables

As we further our studies on Gamma-ray bursts (GRBs), both on theoretical models and observational tools, more and more options begin to open for exploration of its physical properties. As transient events primarily dominated by synchrotron radiation, it is expected that the synchrotron photons emitted by GRBs should present some degree of polarization throughout the evolution of the burst. Whereas observing this polarization can still be challenging due to the constraints on observational tools, especially for short GRBs, it is paramount that the groundwork is laid for the day we have abundant data. In this work, we present a polarization model linked with an off-axis spreading top-hat jet synchrotron scenario in a stratified environment with a density profile $n(r)\propto r^ {-k}$. We present this model's expected temporal polarization evolution for a realistic set of afterglow parameters constrained within the values observed in the GRB literature for four degrees of stratification $k=0,1,1.5 {\rm \, and\,} 2$ and two magnetic field configurations with high extreme anisotropy. We apply this model and predict polarization from a set of GRBs exhibiting off-axis afterglow emission. In particular, for GRB 170817A, we use the available polarimetric upper limits to rule out the possibility of a extremely anisotropic configuration for the magnetic field., Comment: 17 pages, 3 tables, 6 figures. Accepted for publication in the Astrophysical Journal (ApJ)

The GRB 190829A has been widely studied due to its nature and the high energy emission presented. Due to the detection of a very-high-energy component by the High Energy Stereoscopic System and the event's atypically middling luminosity, it has been categorized in a select, limited group of bursts bordering classic GRBs and nearby sub-energetic events. Given the range of models utilized to adequately characterize the afterglow of this burst, it has proven challenging to identify the most probable explanation. Nevertheless, the detection of polarization data provided by the MASTER collaboration has added a new aspect to GRB 190829A that permits us to attempt to explore this degeneracy. In this paper, we present a polarization model coupled with a synchrotron forward-shock model -- a component in all models used to describe GRB 190829A's afterglow -- in order to fit the polarization's temporal evolution with the existing upper limits ($\Pi < 6\%$). We find that the polarization generated from an on-axis emission is favored for strongly anisotropic magnetic field ratios, while an off-axis scenario cannot be fully ruled out when a more isotropic framework is taken into account., Comment: 16 pages, 7 figures, 1 table

Many gamma-ray bursts (GRBs) have been observed from radio wavelengths, and a few at very-high energies (VHEs, > 100GeV). The HAWC gamma-ray observatory is well suited to study transient phenomena at VHEs due to its large field of view and duty cycle. These features allow for searches of VHE emission and can probe different model assumptions of duration and spectra. In this paper, we use data collected by HAWC between December 2014 and May 2020 to search for emission in the energy range from 80 to 800 GeV coming from a sample 47 short GRBs that triggered the Fermi, Swift and Konus satellites during this period. This analysis is optimized to search for delayed and extended VHE emission within the first 20 s of each burst. We find no evidence of VHE emission, either simultaneous or delayed, with respect to the prompt emission. Upper limits (90% confidence level) derived on the GRB fluence are used to constrain the synchrotron self-Compton forward-shock model. Constraints for the interstellar density as low as $10^{-2}$ cm$^{-3}$ are obtained when assuming z=0.3 for bursts with the highest keV-fluences such as GRB 170206A and GRB 181222841. Such a low density makes observing VHE emission mainly from the fast cooling regime challenging., Comment: 20 pages, 8 figures

Sub-relativistic materials launched during the merger of binary compact objects and the core-collapse of massive stars acquire velocity structures when expanding in a stratified environment. The remnant (either a spinning magnetized neutron star (NS) or a central black hole) from the compact-object or core-collapse could additionally inject energy into the afterglow via spin-down luminosity or/and by accreting fall-back material, producing a refreshed shock, modifying the dynamics, and leading to rich radiation signatures at distinct timescales and energy bands with contrasting intensities. We derive the synchrotron light curves evolving in a stratified environment when a power-law velocity distribution parametrizes the energy of the shock, and the remnant continuously injects energy into the blastwave. As the most relevant case, we describe the latest multi-wavelength afterglow observations ($\gtrsim 900$ days) of the GW170817/GRB 170817A event via a synchrotron afterglow model with energy injection of a sub-relativistic material. The features of the remnant and the synchrotron emission of the sub-relativistic material are consistent with a spinning magnetized NS and the faster "blue" kilonova afterglow, respectively. Using the multi-band observations of some short-bursts with evidence of kilonova, we provide constraints on the expected afterglow emission., Comment: 36 pages, 14 figures, accepted for publication in ApJ

Gamma-ray bursts (GRBs) are fascinating extragalactic objects. They represent a fantastic opportunity to investigate unique properties not exhibited in other sources. Multi-wavelength afterglow observations from some short- and long-duration GRBs reveal an atypical long-lasting emission that evolves differently from the canonical afterglow light curves favoring the off-axis emission. We present an analytical synchrotron afterglow scenario, and the hydrodynamical evolution of an off-axis top-hat jet decelerated in a stratified surrounding environment. The analytical synchrotron afterglow model is shown during the coasting, deceleration (off- and on-axis emission), and the post-jet-break decay phases, and the hydrodynamical evolution is computed by numerical simulations showing the time evolution of the Doppler factor, the half-opening angle, the bulk Lorentz factor, and the deceleration radius. We show that numerical simulations are in good agreement with those derived with our analytical approach. We apply the current synchrotron model and describe successfully the delayed non-thermal emission observed in a sample of long and short GRBs with evidence of off-axis emission. Furthermore, we provide constraints on the possible afterglow emission by requiring the multi-wavelength upper limits derived for the closest Swift-detected GRBs and promising gravitational-wave events., Comment: 36 pages, 16 figures, accepted for publication in ApJ

The emission of the so-called extreme BL Lacs poses challenges to the particle acceleration models. The hardness of their spectrum, $\lesssim 2$, in the high-energy band demands unusual parameters using the standard one-zone synchrotron self-Compton (SSC) model with a deficient magnetized plasma. Some authors use either two-zone or hadronic/lepto-hadronic models to relax these atypical values. In this work, we present a lepto-hadronic two-zone model to explain the multi-wavelength observations of the six best-known \textit{extreme} BL Lacs. The very-high-energy gamma-ray observations are described by the photo-hadronic processes in a blob close to the AGN core and by SSC and external inverse Compton-processes in an outer blob. The photo-hadronic interactions occur when accelerated protons in the inner blob interact with annihilation line photons from a sub-relativistic pair plasma. The X-ray observations are described by synchrotron radiation from the outer blob. The parameter values found from the description of the spectral energy distribution for each object with our phenomenological model are similar to each other, and lie in the typical range reported in BL Lacs., Comment: Accepted in Monthly Notices of the Royal Astronomical Society

We investigate the evolution of the afterglow produced by the deceleration of the non-relativistic material due to its surroundings. The ejecta mass is launched into the circumstellar medium with equivalent kinetic energy expressed as a power-law velocity distribution $E\propto (\Gamma\beta)^{-\alpha}$. The density profile of this medium follows a power law $n(r)\propto r^{-k}$ with $k$ the stratification parameter, which accounts for the usual cases of a constant medium ($k=0$) and a wind-like medium ($k=2$). A long-lasting central engine, which injects energy into the ejected material as ($E\propto t^{1-q}$) was also assumed. With our model, we show the predicted light curves associated with this emission for different sets of initial conditions and notice the effect of the variation of these parameters on the frequencies, timescales and intensities. The results are discussed in the Kilonova scenario., Comment: Presented at the 37th International Cosmic Ray Conference (ICRC2021), Berlin, Germany

In light of the most recent observations of late afterglows produced by the merger of compact objects or by the core-collapse of massive dying stars, we research the evolution of the afterglow produced by an off-axis top-hat jet and its interaction with a surrounding medium. The medium is parametrized by a power law distribution of the form $n(r)\propto r^{-k}$ is the stratification parameter and contains the development when the surrounding density is constant ($k=0$) or wind-like ($k=2$). We develop an analytical synchrotron forward-shock model when the outflow is viewed off-axis, and it is decelerated by a stratified medium. Using the X-ray data points collected by a large campaign of orbiting satellites and ground telescopes, we have managed to apply our model and fit the X-ray spectrum of the GRB afterglow associated to SN 2020bvc with conventional parameters. Our model predicts that its circumburst medium is parametrized by a power law with stratification parameter $k=1.5$., Comment: Presented at the 37th International Cosmic Ray Conference (ICRC2021), Berlin, Germany

The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy gamma-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their observation strategy, the size of their field of view and their angular resolution, leading to different analysis approaches. Until now, it has been unclear if the results of observations by both types of instruments are consistent: several of the recently discovered HAWC sources have been followed up by IACTs, resulting in a confirmed detection only in a minority of cases. With this paper, we go further and try to resolve the tensions between previous results by performing a new analysis of the H.E.S.S. Galactic plane survey data, applying an analysis technique comparable between H.E.S.S. and HAWC. Events above 1 TeV are selected for both datasets, the point spread function of H.E.S.S. is broadened to approach that of HAWC, and a similar background estimation method is used. This is the first detailed comparison of the Galactic plane observed by both instruments. H.E.S.S. can confirm the gamma-ray emission of four HAWC sources among seven previously undetected by IACTs, while the three others have measured fluxes below the sensitivity of the H.E.S.S. dataset. Remaining differences in the overall gamma-ray flux can be explained by the systematic uncertainties. Therefore, we confirm a consistent view of the gamma-ray sky between WCD and IACT techniques.

Cosmic rays with energies up to a few PeV are known to be accelerated within the Milky Way. Traditionally, it has been presumed that supernova remnants were the main source of very-high-energy cosmic rays but theoretically it is difficult to get protons to PeV energies and observationally there simply is no evidence to support the remnants as sources of hadrons with energies above a few tens of TeV. One possible source of protons with those energies is the Galactic Center region. Here we report observations of 1-100 TeV gamma rays coming from the 'Cygnus Cocoon', which is a superbubble surrounding a region of OB2 massive star formation. These gamma rays are likely produced by 10-1000 TeV freshly accelerated CRs originating from the enclosed star forming region Cygnus OB2. Hitherto it was not known that such regions could accelerate particles to these energies. The measured flux is likely originated by hadronic interactions. The spectral shape and the emission profile of the Cocoon changes from GeV to TeV energies, which reveals the transport of cosmic particles and historical activity in the superbubble.

Microquasars with high-mass companion stars are promising very-high-energy (VHE; 0.1-100 TeV) gamma-ray emitters, but their behaviors above 10 TeV are poorly known. Using the High Altitude Water Cherenkov (HAWC) observatory, we search for excess gamma-ray emission coincident with the positions of known high-mass microquasars (HMMQs). No significant emission is observed for LS 5039, Cygnus X-1, Cygnus X-3, and SS 433 with 1,523 days of HAWC data. We set the most stringent limit above 10 TeV obtained to date on each individual source. Under the assumption that HMMQs produce gamma rays via a common mechanism, we have performed source-stacking searches, considering two different scenarios: I) gamma-ray luminosity is a fraction $\epsilon_\gamma$ of the microquasar jet luminosity, and II) very-high-energy gamma rays are produced by relativistic electrons up-scattering the radiation field of the companion star in a magnetic field $B$. We obtain $\epsilon_\gamma < 5.4\times 10^{-6}$ for scenario I, which tightly constrains models that suggest observable high-energy neutrino emission by HMMQs. In the case of scenario II, the non-detection of VHE gamma rays yields a strong magnetic field, which challenges synchrotron radiation as the dominant mechanism of the microquasar emission between 10 keV and 10 MeV., Comment: Accepted for publication in ApJL

The study of high-energy gamma rays from passive Giant Molecular Clouds (GMCs) in our Galaxy is an indirect way to characterize and probe the paradigm of the "sea" of cosmic rays in distant parts of the Galaxy. By using data from the High Altitude Water Cherenkov (HAWC) observatory, we measure the gamma-ray flux above 1 TeV of a set of these clouds to test the paradigm. We selected high-galactic latitude clouds that are in HAWC's field-of-view and which are within 1~kpc distance from the Sun. We find no significant excess emission in the cloud regions, nor when we perform a stacked log-likelihood analysis of GMCs. Using a Bayesian approach, we calculate 95\% credible intervals upper limits of the gamma-ray flux and estimate limits on the cosmic-ray energy density of these regions. These are the first limits to constrain gamma-ray emission in the multi-TeV energy range ($>$1 TeV) using passive high-galactic latitude GMCs. Assuming that the main gamma-ray production mechanism is due to proton-proton interaction, the upper limits are consistent with a cosmic-ray flux and energy density similar to that measured at Earth., Comment: 6 figures, 6 tables

The highest-energy known gamma-ray sources are all located within 0.5 degrees of extremely powerful pulsars. This raises the question of whether ultra-high-energy (UHE; $>$ 56 TeV) gamma-ray emission is a universal feature expected near pulsars with a high spin-down power. Using four years of data from the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory, we present a joint-likelihood analysis of ten extremely powerful pulsars to search for UHE gamma-ray emission correlated with these locations. We report a significant detection ($>$ 3$\sigma$), indicating that UHE gamma-ray emission is a generic feature of powerful pulsars. We discuss the emission mechanisms of the gamma rays and the implications of this result. The individual environment that each pulsar is found in appears to play a role in the amount of emission., Comment: accepted by ApJL

We report the ground-level detection of a Galactic Cosmic-Ray (GCR) flux enhancement lasting $\sim$ 17 hr and associated with the passage of a magnetic flux rope (MFR) over the Earth. The MFR was associated with a slow Coronal Mass Ejection (CME) caused by the eruption of a filament on 2016 October 9. Due to the quiet conditions during the eruption and the lack of interactions during the interplanetary CME transport to the Earth, the associated MFR preserved its configuration and reached the Earth with a strong magnetic field, low density, and a very low turbulence level compared to the local background, thus generating the ideal conditions to redirect and guide GCRs (in the $\sim$ 8 to 60 GV rigidity range) along the magnetic field of the MFR. An important negative $B_Z$ component inside the MFR caused large disturbances in the geomagnetic field and a relatively strong geomagnetic storm. However, these disturbances are not the main factors behind the GCR enhancement. Instead, we found that the major factor was the alignment between the MFR axis and the asymptotic direction of the observer.

The Earth is bombarded by ultra-relativistic particles, known as cosmic rays (CRs). CRs with energies up to a few PeV (=10$^{15}$ eV), the knee in the particle spectrum, are believed to have a Galactic origin. One or more factories of PeV CRs, or PeVatrons, must thus be active within our Galaxy. The direct detection of PeV protons from their sources is not possible since they are deflected in the Galactic magnetic fields. Hundred TeV $\gamma$-rays from decaying $\pi^0$, produced when PeV CRs collide with the ambient gas, can provide the decisive evidence of proton acceleration up to the knee. Here we report the discovery by the High Altitude Water Cherenkov (HAWC) observatory of the $\gamma$-ray source, HAWC~J1825-134, whose energy spectrum extends well beyond 200 TeV without a break or cutoff. The source is found to be coincident with a giant molecular cloud. The ambient gas density is as high as 700 protons/cm$^3$. While the nature of this extreme accelerator remains unclear, CRs accelerated to energies of several PeV colliding with the ambient gas likely produce the observed radiation., Comment: Accepted for publication in ApJL

The detection of high-energy astrophysical neutrinos and ultra-high-energy cosmic rays (UHECRs) provides a new way to explore sources of cosmic rays. One of the highest energy neutrino events detected by IceCube, tagged as IC35, is close to the UHECR anisotropy region detected by Pierre Auger Observatory. The nearby starburst galaxy (SBG), NGC 4945, is close to this anisotropic region and inside the mean angular error of the IC35 event. Considering the hypernovae contribution located in the SB region of NGC 4945, which can accelerate protons up to $\sim 10^{17} \, {\rm eV}$ and inject them into the interstellar medium, we investigate the origin of this event around this starburst galaxy. We show that the interaction of these protons with the SB region's gas density could explain Fermi-LAT gamma-ray and radio observations if the magnetic field's strength in the SB region is the order of $\sim \rm mG$. Our estimated PeV neutrino events, in ten years, for this source is approximately 0.01 ($4\times10^{-4}$) if a proton spectral index of 2.4 (2.7) is considered, which would demonstrate that IC35 is not produced in the central region of this SBG. Additionally, we consider the superwind region of NGC 4945 and show that protons can hardly be accelerated in it up to UHEs., Comment: Accepted in Physical Review D

The discovery of the TeV point source 2HWC J2006+341 was reported in the second HAWC gamma-ray catalog. We present a follow-up study of this source here. The TeV emission is best described by an extended source with a soft spectrum. At GeV energies, an extended source is significantly detected in Fermi-LAT data. The matching locations, sizes and spectra suggest that both gamma-ray detections correspond to the same source. Different scenarios for the origin of the emission are considered and we rule out an association to the pulsar PSR J2004+3429 due to extreme energetics required, if located at a distance of 10.8 kpc., Comment: 12 pages, 2 figures. To appear in ApJL

The High Altitude Water Cherenkov Gamma-Ray Observatory (HAWC) continuously detects TeV photons and particles within its large field-of-view, accumulating every day a deeper exposure of two thirds of the sky. We analyzed 1523~days of HAWC live data acquired over four and a half years, in a follow-up analysis of {138} nearby ($z<0.3$) active galactic nuclei from the {\em Fermi} 3FHL catalog culminating within $40^\circ$ of the zenith at Sierra Negra, the HAWC site. This search for persistent TeV emission used a maximum-likelihood analysis assuming intrinsic power-law spectra attenuated by pair production of gamma-ray photons with the extragalactic background light. HAWC clearly detects persistent emission from Mkn~421 and Mkn~501, the two brightest blazars in the TeV sky, at 65$\sigma$ and 17$\sigma$ level, respectively. {Weaker evidence for long-term emission is found for three other known very-high energy emitters:} the radiogalaxy M87 and the BL Lac objects VER~J0521+211 and 1ES~1215+303, the later two at $z\sim 0.1$. We find evidence for collective emission from the set of 30 previously reported very high-energy sources that excludes Mkn~421 and Mkn~501 with a random probability $\sim 10^{-5}$. Upper limits are presented for the sample under the power-law assumption and in the predefined (0.5-2.0), (2.0-8.0) and (8.0-32.0) TeV energy intervals., Comment: The HAWC survey of Fermi 3FHL AGNs: 30 pages, 13 figures. Submitted to the Astrophysical Journal

The High Altitude Water Cherenkov (HAWC) and IceCube observatories, through the Astrophysical Multimessenger Observatory Network (AMON) framework, have developed a multimessenger joint search for extragalactic astrophysical sources. This analysis looks for sources that emit both cosmic neutrinos and gamma rays that are produced in photo-hadronic or hadronic interactions. The AMON system is running continuously, receiving sub-threshold data (i.e. data that is not suited on its own to do astrophysical searches) from HAWC and IceCube, and combining them in real-time. We present here the analysis algorithm, as well as results from archival data collected between June 2015 and August 2018, with a total live-time of 3.0 years. During this period we found two coincident events that have a false alarm rate (FAR) of $<1$ coincidence per year, consistent with the background expectations. The real-time implementation of the analysis in the AMON system began on November 20th, 2019, and issues alerts to the community through the Gamma-ray Coordinates Network with a FAR threshold of $<4$ coincidences per year., Comment: 14 pages, 5 figures, 3 tables

We present a new catalog of TeV gamma-ray sources using 1523 days of data from the High Altitude Water Cherenkov (HAWC) observatory. The catalog represents the most sensitive survey of the Northern gamma-ray sky at energies above several TeV, with three times the exposure compared to the previous HAWC catalog, 2HWC. We report 65 sources detected at $\geq$ 5 sigma significance, along with the positions and spectral fits for each source. The catalog contains eight sources that have no counterpart in the 2HWC catalog, but are within $1^\circ$ of previously detected TeV emitters, and twenty sources that are more than $1^\circ$ away from any previously detected TeV source. Of these twenty new sources, fourteen have a potential counterpart in the fourth \textit{Fermi} Large Area Telescope catalog of gamma-ray sources. We also explore potential associations of 3HWC sources with pulsars in the ATNF pulsar catalog and supernova remnants in the Galactic supernova remnant catalog., Comment: 20 pages, 11 figures. Matches the version published in The Astrophysical Journal. See accompanying data release at https://data.hawc-observatory.org/datasets/3hwc-survey/index.php

We present the detection of VHE gamma-ray emission above 100 TeV from HAWC J2227+610 with the HAWC observatory. Combining our observations with previously published results by VERITAS, we interpret the gamma-ray emission from HAWC J2227+610 as emission from protons with a lower limit in their cutoff energy of 800 TeV. The most likely source of the protons is the associated supernova remnant G106.3+2.7, making it a good candidate for a Galactic PeVatron. However, a purely leptonic origin of the observed emission cannot be excluded at this time., Comment: Accepted for publication by ApJL

The detection of a long flaring activity from blazar TXS 0506+056 in temporal and spatial coincidence with the energetic neutrino IceCube-170922A provided evidence about the photo-hadronic interactions in this source. However, analysis of the archival neutrino and multi-wavelength data from the direction of this blazar between September 2014 and March 2015 revealed a "neutrino flare" without observing quasi-simultaneous activity in the gamma-ray bands, posing challenges to established models. Electron-positron ($e^\pm$) pairs generated from the accretion disks have been amply proposed as a mechanism of bulk acceleration of sub-relativistic and relativistic jets. These pairs annihilate inside the source producing a line around the electron mass which is expected to be blueshifted in the observed frame (on Earth) and redshifted in the frame of the dissipation region of the jet. The redshifted photons in the dissipation region interact with accelerated protons, producing high-energy neutrinos that contribute significantly to the diffuse neutrino flux in the $\sim$ 10 - 20 TeV energy range in connection with gamma-rays from photo-pion process which can be detected by future MeV orbiting satellites. Based on this phenomenological model we can explain the "neutrino flare" reported in 2014-2015., Comment: 10 pages, 4 figures. Accepted for publication in MNRAS

Very-high-energy (VHE) emission is usually interpreted in the synchrotron-self Compton (SSC) scenario, and expected from the low-redshift and high-luminosity gamma-ray bursts (GRBs), as GRB 180720B and GRB 190114C. Recently, VHE emission was detected by the H.E.S.S. telescopes from one of the closest burst GRB 190829A which was associated with the supernova (SN) 2019oyw. In this paper, we present a temporal and spectral analysis from optical bands to Fermi-LAT energy range over multiple observational periods beginning just after the BAT trigger time and extending for almost three months. We show that the X-ray and optical observations are consistent with synchrotron forward-shock emission evolving between the characteristic and cooling spectral breaks during the early and late afterglow in a uniform-density medium. Modeling the light curves together with its spectral energy distribution, it is shown that the outflow expands with an initial bulk Lorentz factor of $\Gamma\sim 30$, which is high for a low-luminosity GRBs and low for a high-luminosity GRBs. The values of the initial bulk Lorentz factor and the isotropic equivalent energy suggest that GRB 190829A is classified as an intermediate-luminosity burst and consequently, it becomes the first burst of this class in being detected in the VHE gamma-ray band by an imaging atmospheric Cherenkov telescope, and, in turn, the first event without being simultaneously observed by the Fermi-LAT instrument. Analyzing the intermediate-luminosity bursts with $z\lesssim 0.2$ such as GRB 130702A, we show that bursts with intermediate luminosity are potential candidates to be detected in very-high energies., Comment: 26 pages, 7 figures, Submitted to ApJ

Cosmic rays, along with stellar radiation and magnetic fields, are known to make up a significant fraction of the energy density of galaxies such as the Milky Way. When cosmic rays interact in the interstellar medium, they produce gamma-ray emission which provides an important indication of how the cosmic rays propagate. Gamma rays from the Andromeda Galaxy (M31), located 785 kpc away, provide a unique opportunity to study cosmic-ray acceleration and diffusion in a galaxy with a structure and evolution very similar to the Milky Way. Using 33 months of data from the High Altitude Water Cherenkov Observatory, we search for TeV gamma rays from the galactic plane of M31. We also investigate past and present evidence of galactic activity in M31 by searching for Fermi Bubble-like structures above and below the galactic nucleus. No significant gamma-ray emission is observed, so we use the null result to compute upper limits on the energy density of cosmic rays $>10$ TeV in M31. The computed upper limits are approximately ten times higher than expected from the extrapolation of the Fermi LAT results., Comment: Accepted for publication in ApJ

Due to the high energies and long distances to the sources, astrophysical observations provide a unique opportunity to test possible signatures of Lorentz invariance violation (LIV). Superluminal LIV enables the decay of photons at high energy. The High Altitude Water Cherenkov (HAWC) Observatory is among the most sensitive gamma-ray instruments currently operating above 10 TeV. HAWC finds evidence of 100 TeV photon emission from at least four astrophysical sources. These observations exclude, for the strongest of the limits set, the LIV energy scale to $2.2\times10^{31}$ eV, over 1800 times the Planck energy and an improvement of 1 to 2 orders of magnitude over previous limits., Comment: Accepted for publication in Physical Review Letters

Primordial Black Holes (PBHs) may have been created by density fluctuations in the early Universe and could be as massive as $> 10^9$ solar masses or as small as the Planck mass. It has been postulated that a black hole has a temperature inversely-proportional to its mass and will thermally emit all species of fundamental particles via Hawking Radiation. PBHs with initial masses of $\sim 5 \times 10^{14}$ g (approximately one gigaton) should be expiring today with bursts of high-energy gamma radiation in the GeV--TeV energy range. The High Altitude Water Cherenkov (HAWC) Observatory is sensitive to gamma rays with energies of $\sim$300 GeV to past 100 TeV, which corresponds to the high end of the PBH burst spectrum. With its large instantaneous field-of-view of $\sim 2$ sr and a duty cycle over 95%, the HAWC Observatory is well suited to perform an all-sky search for PBH bursts. We conducted a search using 959 days of HAWC data and exclude the local PBH burst rate density above $3400~\mathrm{pc^{-3}~yr^{-1}}$ at 99% confidence, the strongest limit on the local PBH burst rate density from any existing electromagnetic measurement., Comment: Corresponding authors: K.L. Engel & A. Peisker. 13 pages, 5 figures

We present the first catalog of gamma-ray sources emitting above 56 and 100 TeV with data from the High Altitude Water Cherenkov (HAWC) Observatory, a wide field-of-view observatory capable of detecting gamma rays up to a few hundred TeV. Nine sources are observed above 56 TeV, all of which are likely Galactic in origin. Three sources continue emitting past 100 TeV, making this the highest-energy gamma-ray source catalog to date. We report the integral flux of each of these objects. We also report spectra for three highest-energy sources and discuss the possibility that they are PeVatrons., Comment: Accepted by Physical Review Letters

List of proceedings from the HAWC Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA., Comment: List of proceedings from the HAWC Collaboration presented at ICRC2019. Corrected typos in the index of the previous version. Follow the "HTML" link to access the list

The detection of the gravitational wave GW170817 defined a breakthrough in multi-messenger astronomy. For the first time, a gravitational wave transient detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo interferometer was associated with a faint electromagnetic gamma-ray counterpart reported by the Gamma-ray Burst Monitor (GBM) aboard on the Fermi satellite. GRB 170817A was followed up by an enormous observational campaign covering a large fraction of the electromagnetic spectrum. In this work, we use the data from High Altitude Water Cherenkov (HAWC) gamma-ray observatory to search for very-high-energy (VHE) TeV photons in coincidence with the X-ray emission from GRB 170817A. Since no counts were observed up to $\sim$120 days after the trigger time, we derive and report the corresponding upper limits in the energy range from 1 to 100 TeV. In addition, we extend the analysis to GRBs with similar features proposed by A. von Kienlin., Comment: Presented at the 36th International Cosmic Ray Conference (ICRC 2019)

An exhaustive analysis of 9-year optical R-band photopolarimetric data of the flat-spectrum radio quasar 3C279 from 2008 February 27 to 2017 May 25 is presented, alongside with multiwavelength observing campaigns performed during the flaring activity exhibited in 2009 February/March, 2011 June, 2014 March/April, 2015 June and 2017 February. In the R-band, this source showed the maximum brightness state of $13.68\pm 0.11$ mag ($1.36\pm0.20$ mJy) on 2017 March 02, and the lowest brightness state ever recorded of $18.20\pm 0.87$ mag ($0.16\pm0.03$ mJy) on 2010 June 17. During the entire period of observations, the polarization degree varied between $0.48\pm0.17$% and $31.65\pm0.77$% and the electric vector position angle exhibited large rotations between $82.98^\circ \pm0.92$ and $446.32^\circ \pm1.95$. Optical polarization data show that this source has a stable polarized component that varied from $\sim$6% (before the 2009 flare) to $\sim$13% after the flare. The overall behavior of our polarized variability data supports the scenario of jet precessions as responsible of the observed large rotations of the electric vector position angle. Discrete correlation function analysis show that the lags between gamma-rays and X-rays compared to the optical R-band fluxes are $\Delta t \sim$ 31 d and $1$ d in 2009. Lags were also found among gamma-rays compared with X-rays and radio of $\Delta t \sim$ 30 d and $43$ d in 2011, and among radio and optical-R band of $\Delta t \sim$ 10 d in 2014. A very intense flare in 2017 was observed in optical bands with a dramatic variation in the polarization degree (from $\sim$ 6% to 20%) in 90 days without exhibiting flaring activity in other wavelengths., Comment: 21 pages, 10 figures. Accepted for publication in ApJS

The origin of recently discovered PeV neutrinos is an unsolved problem. In this work we consider a hadronic scenario to produce PeV neutrinos from a region around giant lobes of Centaurus A. Although ultrahigh-energy cosmic rays (UHECRs) are accelerated and confined by giant lobes, they can escape to be later injected in the inter-group medium where galaxies near the giant lobes provides the condition to confine them. UHECRs interact with low-energy photons and protons producing high-energy photons and neutrinos. We found that the IC35 event cannot be generated neither inside the giant lobes nor galaxies close to the lobes of Centaurus A., Comment: To appear in the Proceedings of the XV International Conference on Topics in Astroparticle and Underground Physics (TAUP2017) at Journal of Physics: Conference Series (JPCS)

The detection of gravitational waves together with their electromagnetic counterpart, in the gamma-ray burst GRB 170817A, marked a new era of multi-messenger astronomy. Several theoretical models have been proposed to explain the atypical behavior of this event. Recently, it was shown that the multi-wavelength afterglow of GRB 170817A was consistent with a synchrotron forward-shock model when the outflow was viewed off-axis, decelerated in a uniform medium and parametrized through a power-law velocity distribution. Motivated by the upper limits on the very-high-energy emission, and the stratified medium in the close vicinity of a binary neutron star merger proposed to explain the gamma-ray flux in the short GRB 150101B, we extend the mechanism proposed to explain GRB 170817A to a more general scenario deriving the synchrotron self-Compton (SSC) and synchrotron forward-shock model when the off-axis outflow is decelerated in a uniform and stratified circumburst density. As particular cases, we show that the delayed and long-lasting afterglow emission observed in GRB 080503, GRB140903A, GRB 150101B, and GRB 160821B could be interpreted by a similar scenario to the one used to describe GRB 170817A. In addition, we show that the proposed scenario agrees with the MAGIC, Fermi-LAT and H.E.S.S upper limits on gamma-ray emission from GRB 160821B and GRB 170817A., Comment: 24 pages, 15 figures. Accepted for publication in ApJ

Early and late multiwavelength observations play an important role in determining the nature of the progenitor, circumburst medium, physical processes and emitting regions associated to the spectral and temporal features of bursts. GRB 180720B is a long and powerful burst detected by a large number of observatories in multiwavelenths that range from radio bands to sub-TeV gamma-rays. The simultaneous multiwavelength observations were presented over multiple periods of time beginning just after the trigger time and extending for more than 30 days. The temporal and spectral analysis of Fermi LAT observations suggests that it presents similar characteristics to other bursts detected by this instrument. Coupled with X-ray and optical observations, the standard external-shock model in a homogeneous medium is favored by this analysis. The X-ray flare is consistent with the synchrotron self-Compton (SSC) model from the reverse-shock region evolving in a thin shell and long-lived LAT, X-ray and optical data with the standard synchrotron forward-shock model. The best-fit parameters derived with the Markov chain Monte Carlo simulations indicate that the outflow is endowed with magnetic fields and that the radio observations are in the self-absorption regime. The SSC forward-shock model with our parameters can explain the LAT photons beyond the synchrotron limit as well as the emission recently reported by the HESS Collaboration., Comment: 18 pages, 9 figures. Minor changes. Accepted for publication in ApJ

We present TeV gamma-ray observations of the Crab Nebula, the standard reference source in ground-based gamma-ray astronomy, using data from the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory. In this analysis we use two independent energy-estimation methods that utilize extensive air shower variables such as the core position, shower angle, and shower lateral energy distribution. In contrast, the previously published HAWC energy spectrum roughly estimated the shower energy with only the number of photomultipliers triggered. This new methodology yields a much improved energy resolution over the previous analysis and extends HAWC's ability to accurately measure gamma-ray energies well beyond 100 TeV. The energy spectrum of the Crab Nebula is well fit to a log parabola shape $\left(\frac{dN}{dE} = \phi_0 \left(E/\textrm{7 TeV}\right)^{-\alpha-\beta\ln\left(E/\textrm{7 TeV}\right)}\right)$ with emission up to at least 100 TeV. For the first estimator, a ground parameter that utilizes fits to the lateral distribution function to measure the charge density 40 meters from the shower axis, the best-fit values are $\phi_o$=(2.35$\pm$0.04$^{+0.20}_{-0.21}$)$\times$10$^{-13}$ (TeV cm$^2$ s)$^{-1}$, $\alpha$=2.79$\pm$0.02$^{+0.01}_{-0.03}$, and $\beta$=0.10$\pm$0.01$^{+0.01}_{-0.03}$. For the second estimator, a neural network which uses the charge distribution in annuli around the core and other variables, these values are $\phi_o$=(2.31$\pm$0.02$^{+0.32}_{-0.17}$)$\times$10$^{-13}$ (TeV cm$^2$ s)$^{-1}$, $\alpha$=2.73$\pm$0.02$^{+0.03}_{-0.02}$, and $\beta$=0.06$\pm$0.01$\pm$0.02. The first set of uncertainties are statistical; the second set are systematic. Both methods yield compatible results. These measurements are the highest-energy observation of a gamma-ray source to date., Comment: published in ApJ

Very-high-energy (VHE; $\geq 10$ GeV) photons are expected from the nearest and brightest Gamma-ray bursts (GRBs). VHE photons, at energies higher than 300 GeV, were recently reported by the MAGIC collaboration for this burst. Immediately, GRB 190114C was followed up by a massive observational campaign covering a large fraction of the electromagnetic spectrum. In this paper, we obtain the LAT light curve of GRB 190114C and show that it exhibits similar features to other bright LAT-detected bursts; the first high-energy photon ($\geq$ 100 MeV) is delayed with the onset of the prompt phase and the flux light curve exhibits a long-lived emission (lasting much longer than the prompt phase) and a short-lasting bright peak (located at the beginning of long-lived emission). Analyzing the multi-wavelength observations, we show that the short-lasting LAT and GBM bright peaks are consistent with the synchrotron self-Compton reverse-shock model and the long-lived observations with the standard synchrotron forward-shock model that evolves from a stratified stellar-wind like medium to a uniform ISM-like medium. Given the best-fit values, a bright optical flash produced by synchrotron reverse-shock emission is expected. From our analysis we infer that the high-energy photons are produced in the deceleration phase of the outflow and some additional processes to synchrotron in the forward shocks should be considered to properly describe the LAT photons with energies beyond the synchrotron limit. Moreover, we claim that an outflow endowed with magnetic fields could describe the polarization and properties exhibited in the light curve of GRB 190114C., Comment: 18 Pages, 5 figures. Accepted for publication in ApJ Letter

We present the first full-sky analysis of the cosmic ray arrival direction distribution with data collected by the HAWC and IceCube observatories in the Northern and Southern hemispheres at the same median primary particle energy of 10 TeV. The combined sky map and angular power spectrum largely eliminate biases that result from partial sky coverage and holds a key to probe into the propagation properties of TeV cosmic rays through our local interstellar medium and the interaction between the interstellar and heliospheric magnetic fields. From the map we determine the horizontal dipole components of the anisotropy $\delta_{0h} = 9.16 \times 10^{-4}$ and $\delta_{6h} = 7.25 \times 10^{-4}~(\pm0.04 \times 10^{-4})$. In addition, we infer the direction ($229.2\pm 3.5^\circ$ RA , $11.4\pm 3.0^\circ$ Dec.) of the interstellar magnetic field from the boundary between large scale excess and deficit regions from which we estimate the missing corresponding vertical dipole component of the large scale anisotropy to be $\delta_N \sim -3.97 ^{+1.0}_{-2.0} \times 10^{-4}$., Comment: Accepted by The Astrophysical Journal