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.

The detection of the first electromagnetic counterpart to the binary neutron star (BNS) merger remnant GW170817 established the connection between short $\gamma$-ray bursts and BNS mergers. It also confirmed the forging of heavy elements in the ejecta (a so-called kilonova) via the r-process nucleosynthesis. The appearance of non-thermal radio and X-ray emission, as well as the brightening, which lasted more than 100 days, were somewhat unexpected. Current theoretical models attempt to explain this temporal behavior as either originating from a relativistic off-axis jet or a kilonova-like outflow. In either scenario, there is some ambiguity regarding how much energy is transported in the non-thermal electrons versus the magnetic field of the emission region. Combining the VLA (radio) and Chandra (X-ray) measurements with observations in the GeV-TeV domain can help break this ambiguity, almost independently of the assumed origin of the emission. Here we report for the first time on deep H.E.S.S. observations of GW170817 / GRB 170817A between 124 and 272 days after the BNS merger with the full H.E.S.S. array of telescopes, as well as on an updated analysis of the prompt (<5 days) observations with the upgraded H.E.S.S. phase-I telescopes. We discuss implications of the H.E.S.S. measurement for the magnetic field in the context of different source scenarios.

We report on the detection of very-high-energy (VHE; $E > 100$ GeV) $\gamma$-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multi-wavelength observations with Fermi/LAT, XRT and UVOT on board the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum towards the VHE $\gamma$-ray regime, we deduce a 95% confidence level upper limit on the unknown redshift of KUV 00311-1938 of z < 0.98, and of PKS 1440-389 of z < 0.53. When combined with previous spectroscopy results the redshift of KUV 00311-1938 is constrained to $0.51 \leq z < 0.98$ and for PKS 1440-389 to $0.14 \lessapprox z < 0.53$.

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.

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.

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.

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.

List of proceedings from the HAWC Collaboration presented at the 36th International Cosmic Ray Conference, 24 July - 1 August 2019, Madison, Wisconsin, USA.

We report on broadband observations of CTA 102 ($z=1.037$) during the active states in 2016-2017. In the $\gamma$-ray band, Fermi LAT observed several prominent flares which followed a harder-when-brighter behavior: the hardest photon index $\Gamma=1.61\pm 0.10$ being unusual for FSRQs. The peak $\gamma$-ray flux above 100 MeV $(3.55\pm0.55)\times10^{-5}\:{\rm photon\:cm^{-2}\:s^{-1}}$ observed on MJD 57738.47 within 4.31 minutes, corresponds to an isotropic $\gamma$-ray luminosity of $L_{\gamma}=3.25\times10^{50}\:{\rm erg\:s^{-1}}$, comparable with the highest values observed from blazars so far. The analyses of the Swift UVOT/XRT data show an increase in the UV/optical and X-ray bands which is contemporaneous with the bright $\gamma$-ray periods. The X-ray spectrum observed by Swift XRT and NuSTAR during the $\gamma$-ray flaring period is characterized by a hard photon index of $\sim1.30$. The shortest e-folding time was $4.08\pm1.44$ hours, suggesting a very compact emission region $R\leq\delta\times2.16\times10^{14}$ cm. We modeled the spectral energy distribution of CTA 102 in several periods (having different properties in UV/optical, X-ray and $\gamma$-ray bands) assuming a compact blob inside and outside the BLR. We found that the high-energy data are better described when the infrared thermal radiation of the dusty torus is considered. In the flaring periods when the correlation between the $\gamma$-ray and UV/optical/X-ray bands is lacking, the $\gamma$-ray emission can be produced from the interaction of fresh electrons in a different blob, which does not make a dominant contribution at lower energies.

We report on a detailed investigation of the $\gamma$-ray emission from 26 non-blazar AGNs based on the Fermi LAT data accumulated for 7 years. The photon index of non-blazar AGNs changes in the range of 1.84-2.86 and the flux varies from a few times $10^{-9} photon\: cm^{-2} s^{-1}$ to $10^{-7} photon\: cm^{-2}s^{-1}$. Over long time periods, power-law provides an adequate description of the $\gamma$-ray spectra of almost all sources. Significant curvature is observed in the $\gamma$-ray spectra of NGC 1275, NGC 6251, SBS 0846+513 and PMN J0948+0022 and their spectra are better described by log-parabola or power-law with exponential cut-off models. The $\gamma$-ray spectra of PKS 0625-25 and 3C 380 show a possible deviation from a simple power-law shape, indicating a spectral cutoff around the observed photon energy of $E_{cut}=131.2\pm88.04$ GeV and $E_{cut}=55.57\pm50.74$ GeV, respectively. Our analysis confirms the previous finding of an unusual spectral turnover in the $\gamma$-ray spectrum of Cen A: the photon index changes from $2.75\pm0.02$ to $2.31\pm0.1$ at $2.35\pm0.08$ GeV. In the $\Gamma-L_{\gamma}$ plane, the luminosity of non-blazar AGNs is spread in the range of $10^{41}-10^{47}\: erg\: s^{-1}$, where the lowest luminosity have FRI radio galaxies (but typically appear with a harder photon index) and the highest- SSRQs/NLSY1s (with softer photon indexes). We confirm the previously reported short-timescale flux variability of NGC 1275 and 3C 120. The $\gamma$-ray emission from NLSY1s, 1H 0323+342, SBS 0846+513 and PMN J0948+0022, is variable, showing flares in short scales sometimes accompanied by a moderate hardening of their spectra (e.g., on MJD 56146.8 the $\gamma$-ray photon index of SBS 0846+513 was $1.73\pm0.14$). 3C 111, Cen A core, 3C 207, 3C 275.1, 3C 380, 4C+39.23B, PKS 1502+036 and PKS 2004-447 show a long-timescale flux variability in the $\gamma$-ray band.

Context. Important information on the evolution of the jet can be obtained by comparing the physical state of the plasma at its propagation through the broad-line region (where the jet is most likely formed) into the intergalactic medium, where it starts to significantly decelerate. Aims. We compare the constraints on the physical parameters in the innermost ($\leq$ pc) and outer ($\geq$ kpc) regions of the 3C 120 jet by means of a detailed multiwavelength analysis and theoretical modeling of their broadband spectra. Methods.The data collected by Fermi LAT, Swift and Chandra are analyzed together and the spectral energy distributions are modeled using a leptonic synchrotron and inverse Compton model, taking into account the seed photons originating inside and outside of the jet. The model parameters are estimated using the MCMC method. Results. The $\gamma$-ray flux from the inner jet of 3C 120 was characterized by rapid variation from MJD 56900 to MJD 57300. Two strong flares were observed on April 24, 2015 when, within 19.0 minutes and 3.15 hours the flux was as high as $(7.46\pm1.56)\times10^{-6}photon\:cm^{-2}\:s^{-1}$ and $(4.71\pm0.92)\times10^{-6}photon\:cm^{-2}\:s^{-1}$ respectively. The broadband emission in the quiet and flaring states can be described as SSC emission while IC scattering of dusty torus photons cannot be excluded for the flaring states. The X-ray emission from the knots can be well reproduced by IC scattering of CMB photons only if the jet is highly relativistic (since even when $\delta=10$ still $U_{\rm e}/U_B\geq80$). These extreme requirements can be somewhat softened assuming the X-rays are from the synchrotron emission of a second population of very-high-energy electrons. Conclusions. We found that the jet power estimated at two scales is consistent, suggesting that the jet does not suffer severe dissipation, it simply becomes radiatively inefficient.

We report on a detailed analysis of the $\gamma$-ray light curve of NGC 1275 using the Fermi large area telescope data accumulated in 2008-2017. Major $\gamma$-ray flares were observed in October 2015 and December 2016/January 2017 when the source reached a daily peak flux of $(2.21\pm0.26)\times10^{-6}\:{\rm photon\:cm^{-2}\:s^{-1}}$, achieving a flux of $(3.48\pm0.87)\times10^{-6}\:{\rm photon\:cm^{-2}\:s^{-1}}$ within 3 hours, which corresponds to an apparent isotropic $\gamma$-ray luminosity of $\simeq3.84\times10^{45}\:{\rm erg\:s^{-1}}$. The most rapid flare had e-folding time as short as $1.21\pm0.22$ hours which had never been previously observed for any radio galaxy in $\gamma$-ray band. Also $\gamma$-ray spectral changes were observed during these flares: in the flux versus photon index plane the spectral evolution follows correspondingly a counter clockwise and a clockwise loop inferred from the light curve generated by an adaptive binning method. On December 30, 2016 and January 01, 2017 the X-ray photon index softened ($\Gamma_{\rm X}\simeq 1.75-1.77$) and the flux increased nearly $\sim3$ times as compared with the quiet state. The observed hour-scale variability suggests a very compact emission region ($R_\gamma\leq5.22\times10^{14}\:(\delta/4)\:{\rm cm}$) implying that the observed emission is most likely produced in the subparsec-scale jet if the entire jet width is responsible for the emission. During the active periods the $\gamma$-ray photon index hardened, shifting the peak of the high energy spectral component to $>{\rm GeV}$, making it difficult to explain the observed X-ray and $\gamma$-ray data in the standard one-zone synchrotron self-Compton model.

We report the analysis of Fermi Large Area Telescope data from five years of observations of the broad line radio galaxy 3C 120. The accumulation of larger data set results in the detection of high-energy $\gamma$-rays up to 10 GeV, with a detection significance of about $8.7\sigma$. A power-law spectrum with a photon index of $2.72\pm0.1$ and integrated flux of $F_{\gamma}=(2.35\pm0.5)\times10^{-8}\:\mathrm{photon\:cm}^{-2}s^{-1}$ above 100 MeV well describe the data averaged over five year observations. The variability analysis of the light curve with 180-, and 365- day bins reveals flux increase (nearly twice from its average level) during the last year of observation. This variability on month timescales indicates the compactness of the emitting region. The $\gamma$-ray spectrum can be described as synchrotron self-Compton (SSC) emission from the electron population producing the radio-to-X-ray emission in the jet. The required electron energy density exceeds the one of magnetic field only by a factor of 2 meaning no significant deviation from equipartition.