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.

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.

The High Altitude Water Cherenkov (HAWC) array is sensitive to high energy Cosmic Rays (CR) in the $\sim 10$ to $\sim 200$ TeV energy range, making it possible to construct maps of the so called "Sun Shadow" ($SS$), i. e. of the deficit of CR coming from the direction of the Sun. In this work, we present the variation of the Relative Intensity of the deficit ($SS_{RI}$) for three years of HAWC observations form 2016 to 2018 in which we found a clear decreasing trend of the ($SS_{RI}$) over the studied period, corresponding to the declining phase of the solar cycle 24. By comparing the $SS_{RI}$ with the photospheric magnetic field evolution, we show that there is a linear relationship between the $SS_{RI}$ and the median photospheric magnetic field of the Active Region belt (-40$^\circ \le$ lat $\le$ 40$^\circ$) and a inverse linear relationship with the polar photospheric magnetic field (lat $\ge \pm$ 60$^\circ$). The former relationship is due to the magnetic field causing a deviation of the CR, whereas the latter reflects the change of the heliospheric field topology from multipolar to dipolar configurations. These relationships are valid only when the median magnetic field is lower than 8 G, during the declining and minimum phases of the solar cycle 24. Finally, we show that relativistic charged particles, in the 10 to 200 TeV energy range, are deflected a few degrees.

Strong electric fields in thunderclouds have long been known to accelerate secondary charged particles. We investigate this effect using three years (2015-2017) of data from the scalar system of the High Altitude Water Cherenkov (HAWC) observatory, which is an air shower array deployed 4100 m a.s.l. in central Mexico. The experimental site is frequently affected by strong thunderstorms, and the detector's high altitude, large area, and high sensitivity to cosmic-ray air showers make it ideal for investigating particle acceleration due to the electric fields present inside the thunder storm clouds. In particular, the scaler system of HAWC records the output of each one of the 1200 PMTs as well as the 2, 3, and 4-fold multiplicities (logic AND in a time window of 30 ns) of each water Cherenkov detectors (WCD) with a sampling rate of 40 Hz. Using data from this scaler system, we identify approximately 100 increases in the scaler rate which is in time coincidence with thunderstorms. These events show high cross correlation between the scaler rate and the electric field, hence can be produced by the acceleration of secondary particle by the thunderstorm electric fields. In this work we present the method of identification of these events and their general characteristics.

The variation in atmospheric pressure is due to changes in mass of the air column above, which in turn resembles the density variation of atmosphere and will affect the decay of secondary particles of cosmic rays. The ground based cosmic ray detectors observe pressure dependent variation in their flux. The High Altitude Water Cherenkov (HAWC) gamma ray observatory is a great detector of secondary particles because of its high altitude, high uptime, and large area (including total photo-cathode area), which makes the HAWC scaler system an ideal instrument for solar modulation studies. Although, in order to perform these studies it is necessary to isolate and remove the atmospheric modulations. The observed rate in each PMT has signatures of both the solar and atmospheric modulations, which makes it difficult to measure the pressure coefficient ($\beta_P$ ). The pressure at the HAWC site shows a periodic behavior ($\sim$ 12 hours), which also reflects in the scalar rates. This periodic property was used to isolate the pressure modulation and $\beta_P$ were estimated with accuracy. Since the pressure dependence is a physical phenomenon, the estimated coefficients for PMTs should be identical, any deviation from this can be due to malfunction of the PMT. This make this method a useful tool to identify the malfunctioning PMTs and help us to isolate them from the analysis. In this analysis we are presenting the method of estimation of the pressure coefficients and its usage to correct the HAWC scalar data to make it suitable for the solar modulations studies.

The GRAPES-3 muon telescope in Ooty, India had claimed detection of a 2 hour (h) high-energy ($\sim$20 GeV) burst of galactic cosmic-rays (GCRs) through a $>$50$\sigma$ surge in GeV muons, was caused by reconnection of the interplanetary magnetic field (IMF) in the magnetosphere that led to transient weakening of Earth's magnetic shield. This burst had occurred during a G4-class geomagnetic storm (storm) with a delay of $\frac{1}{2}$h relative to the coronal mass ejection (CME) of 22 June 2015 (Mohanty et al., 2016). However, recently a group interpreted the occurrence of the same burst in a subset of 31 neutron monitors (NMs) to have been the result of an anisotropy in interplanetary space (Evenson et al., 2017) in contrast to the claim in (Mohanty et al., 2016). A new analysis of the GRAPES-3 data with a fine 10.6$^{\circ}$ angular segmentation shows the speculation of interplanetary anisotropy to be incorrect, and offers a possible explanation of the NM observations. The observed 28 minutes (min) delay of the burst relative to the CME can be explained by the movement of the reconnection front from the bow shock to the surface of Earth at an average speed of 35 km/s, much lower than the CME speed of 700 km/s. This measurement may provide a more accurate estimate of the start of the storm.

Aims. Forbush decrease (FD) is a transient decrease followed by a gradual recovery in the observed galactic cosmic ray intensity. We seek to understand the relationship between the FDs and near-Earth interplanetary magnetic field (IMF) enhancements associated with solar coronal mass ejections (CMEs). Methods. We use muon data at cutoff rigidities ranging from 14 to 24 GV from the GRAPES-3 tracking muon telescope to identify FD events. We select those FD events that have a reasonably clean profile, and magnitude > 0.25%. We use IMF data from ACE/WIND spacecrafts. We look for correlations between the FD profile and that of the one hour averaged IMF. We ask if the diffusion of high energy protons into the large scale magnetic field is the cause of the lag observed between the FD and the IMF. Results. The enhancement of the IMF associated with FDs occurs mainly in the shock-sheath region, and the turbulence level in the magnetic field is also enhanced in this region. The observed FD profiles look remarkably similar to the IMF enhancement profiles. The FDs typically lag the IMF enhancement by a few hours. The lag corresponds to the time taken by high energy protons to diffuse into the magnetic field enhancement via cross-field diffusion. Conclusions. Our findings show that high rigidity FDs associated with CMEs are caused primarily by the cumulative diffusion of protons across the magnetic field enhancement in the turbulent sheath region between the shock and the CME.

We examine the propagation of several CMEs with well-observed flux rope signatures in the field of view of the SECCHI coronagraphs aboard the STEREO satellites using the GCS fitting method of Thernisien, Vourlidas \& Howard (2009). We find that the manner in which they propagate is approximately self-similar; i.e., the ratio ($\kappa$) of the flux rope minor radius to its major radius remains approximately constant with time. We use this observation of self-similarity to draw conclusions regarding the local pitch angle ($\gamma$) of the flux rope magnetic field and the misalignment angle ($\chi$) between the current density ${\mathbf J}$ and the magnetic field ${\mathbf B}$. Our results suggest that the magnetic field and current configurations inside flux ropes deviate substantially from a force-free state in typical coronagraph fields of view, validating the idea of CMEs being driven by Lorentz self-forces.