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 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 - 1015.

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.

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.

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.

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.

Simulations of dark matter show a discrepancy between the expected number of Galactic dark matter sub-halos and how many have been optically observed. Some of these unseen satellites may exist as dark dwarf galaxies: sub-halos like dwarf galaxies with no luminous counterpart. Assuming WIMP dark matter, it may be possible to detect these unseen sub-halos from gamma-ray signals originating from dark matter annihilation. The High Altitude Water Cherenkov Observatory (HAWC) is a very high energy (500 GeV to 100 TeV) gamma ray detector with a wide field-of-view and near continuous duty cycle, making HAWC ideal for unbiased sky surveys. We perform such a search for gamma ray signals from dark dwarfs in the Milky Way halo. We perform a targeted search of HAWC gamma-ray sources which have no known association with lower-energy counterparts, based on an unbiased search of the entire sky. With no sources found to strongly prefer dark matter models, we calculate the ability of HAWC to observe dark dwarfs. We also compute the HAWC sensitivity to potential future detections for a given model of dark matter substructure.

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.

The Andromeda Galaxy (M31) is a nearby ($\sim$780 kpc) galaxy similar to our own Milky Way. Observational evidence suggests that it resides in a large halo of dark matter (DM), making it a good target for DM searches. We present a search for gamma rays from M31 using 1017 days of data from the High Altitude Water Cherenkov (HAWC) Observatory. With its wide field of view and constant monitoring, HAWC is well-suited to search for DM in extended targets like M31. No DM annihilation or decay signal was detected for DM masses from 1 to 100 TeV in the $b\bar{b}$, $t\bar{t}$, $\tau^{+}\tau^{-}$, $\mu^{+}\mu^{-}$, and $W^{+}W^{-}$ channels. Therefore we present limits on those processes. Our limits nicely complement the existing body of DM limits from other targets and instruments. Specifically the DM decay limits from our benchmark model are the most constraining for DM masses from 25 TeV to 100 TeV in the $b\bar{b}, t\bar{t}$ and $\mu^{+}\mu{-}$ channels. In addition to DM-specific limits, we also calculate general gamma-ray flux limits for M31 in 5 energy bins from 1 TeV to 100 TeV.

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}$.

SS 433 is a binary system containing a supergiant star that is overflowing its Roche lobe with matter accreting onto a compact object (either a black hole or neutron star). Two jets of ionized matter with a bulk velocity of $\sim0.26c$ extend from the binary, perpendicular to the line of sight, and terminate inside W50, a supernova remnant that is being distorted by the jets. SS 433 differs from other microquasars in that the accretion is believed to be super-Eddington, and the luminosity of the system is $\sim10^{40}$ erg s$^{-1}$. The lobes of W50 in which the jets terminate, about 40 pc from the central source, are expected to accelerate charged particles, and indeed radio and X-ray emission consistent with electron synchrotron emission in a magnetic field have been observed. At higher energies (>100 GeV), the particle fluxes of $\gamma$ rays from X-ray hotspots around SS 433 have been reported as flux upper limits. In this energy regime, it has been unclear whether the emission is dominated by electrons that are interacting with photons from the cosmic microwave background through inverse-Compton scattering or by protons interacting with the ambient gas. Here we report TeV $\gamma$-ray observations of the SS 433/W50 system where the lobes are spatially resolved. The TeV emission is localized to structures in the lobes, far from the center of the system where the jets are formed. We have measured photon energies of at least 25 TeV, and these are certainly not Doppler boosted, because of the viewing geometry. We conclude that the emission from radio to TeV energies is consistent with a single population of electrons with energies extending to at least hundreds of TeV in a magnetic field of $\sim16$~micro-Gauss.

Pierre Auger observatory reported the distribution of arrival directions of the highest energy cosmic rays. These events were collected in 10 years of operations with declinations between -90$^\circ$ and +45$^\circ$. The IceCube neutrino telescope reported the detection of 82 extraterrestrial neutrinos in the High-Energy Starting Events catalog. The highest-energy neutrino event (IC35) reported in this catalog had an energy of 2004$^{+236}_{-262}$ TeV and its reconstruction was centered at RA=$208.4^\circ$ and DEC=$-55.8^\circ$ (J2000). Being Centaurus A the nearest radio-loud active Galactic nucleus and one of the potential candidates for accelerating cosmic rays up to $\sim\, 10^{20}$ eV, we show that the ultra-high-energy cosmic rays (UHECRs) around the direction of Centaurus A (15$^\circ$ radius) could be accelerated inside the giant lobes. Studying the composition of UHECRs through the photo-disintegration processes, the maximum energy that these can reach in the giant lobes and the average deflecting angles that nuclei undergo due to Galactic and extragalactic magnetic fields, it is shown that the most promise candidates of UHECR composition are light nuclei such as carbon/nitrogen nuclei. Considering the interactions of the relativistic carbon nuclei with photon fields (extragalactic background light and synchrotron) and materials (inside and outside the lobes), we do not find enough evidence to associate the IC35 event with the UHECRs reported by Pierre Auger observatory.

An indirect measurement of the antiproton flux in cosmic rays is possible as the particles undergo deflection by the geomagnetic field. This effect can be measured by studying the deficit in the flux, or shadow, created by the Moon as it absorbs cosmic rays that are headed towards the Earth. The shadow is displaced from the actual position of the Moon due to geomagnetic deflection, which is a function of the energy and charge of the cosmic rays. The displacement provides a natural tool for momentum/charge discrimination that can be used to study the composition of cosmic rays. Using 33 months of data comprising more than 80 billion cosmic rays measured by the High Altitude Water Cherenkov (HAWC) observatory, we have analyzed the Moon shadow to search for TeV antiprotons in cosmic rays. We present our first upper limits on the $\bar{p}/p$ fraction, which in the absence of any direct measurements, provide the tightest available constraints of $\sim1\%$ on the antiproton fraction for energies between 1 and 10 TeV.

The results of 8-year R-band photopolarimetric data of blazar Mrk 421 collected from February 2008 to May 2016 are presented, along with extensive multiwavelength observations covering from radio to TeV $\gamma$-rays around the flares observed in May 2008, March 2010, and April 2013. The most important results are found in 2013 when the source displayed in the R-band a very high brightness state of $11.29\pm 0.03$ mag ($93.60\pm 1.53$ mJy) on April 10th and a polarization degree of ($11.00\pm0.44$) on May 13th. The analysis of the optical data shows that the polarization variability is due to the superposition of two polarized components that might be produced in two distinct emitting regions. An intranight photopolarimetric variability study carried out over 7 nights after the 2013 April maximum found flux and polarization variations on the nights of April 14, 15, 16 and 19. In addition, the flux shows a minimum variability timescale of $\Delta\,t=$2.34$\pm$0.12 hours, and that the polarization degree presented variations of $\sim$ (1 - 2) in a timescale of $\Delta\,t\sim$ minutes. Also, a detailed analysis of the intranight data shows a coherence length of the large-scale magnetic field of $l_B\simeq$ 0.3 pc which is the same order of magnitude to the distance traveled by the relativistic shocks. This result suggests that there is a connection between the intranight polarimetric variations and spatial changes of the magnetic field. Analysis of the complete R-band data along with the historical optical light curve found for this object shows that Mrk 421 varies with a period of 16.26 $\pm$ 1.78 years.