We report on our observing campaign of the compact binary merger GW190814, detected by the Advanced LIGO and Advanced Virgo detectors on August 14th, 2019. This signal has the best localisation of any observed gravitational wave (GW) source, with a 90% probability area of 18.5 deg$^2$, and an estimated distance of ~ 240 Mpc. We obtained wide-field observations with the Deca-Degree Optical Transient Imager (DDOTI) covering 88% of the probability area down to a limiting magnitude of $w$ = 19.9 AB. Nearby galaxies within the high probability region were targeted with the Lowell Discovery Telescope (LDT), whereas promising candidate counterparts were characterized through multi-colour photometry with the Reionization and Transients InfraRed (RATIR) and spectroscopy with the Gran Telescopio de Canarias (GTC). We use our optical and near-infrared limits in conjunction with the upper limits obtained by the community to constrain the possible electromagnetic counterparts associated with the merger. A gamma-ray burst seen along its jet's axis is disfavoured by the multi-wavelength dataset, whereas the presence of a burst seen at larger viewing angles is not well constrained. Although our observations are not sensitive to a kilonova similar to AT2017gfo, we can rule out high-mass (> 0.1 M$_{\odot}$) fast-moving (mean velocity >= 0.3c) wind ejecta for a possible kilonova associated with this merger.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

The Fermi-LAT collaboration presented the second gamma-ray burst (GRB) catalog covering its first 10 years of operations. A significant fraction of afterglow-phase light curves in this catalog cannot be explained by the closure relations of the standard synchrotron forward-shock model, suggesting that there could be an important contribution from another process. In view of the above, we derive the synchrotron self-Compton (SSC) light curves from the reverse shock in the thick- and thin-shell regime for a uniform-density medium. We show that this emission could explain the GeV flares exhibited in some LAT light curves. Additionally, we demonstrate that the passage of the forward shock synchrotron cooling break through the LAT band from jets expanding in a uniform-density environment may be responsible for the late time ($\approx10^2$ s) steepening of LAT GRB afterglow light curves. As a particular case, we model the LAT light curve of GRB 160509A that exhibited a GeV flare together with a break in the long-lasting emission, and also two very high energy photons with energies of 51.9 and 41.5 GeV observed 76.5 and 242 s after the onset of the burst, respectively. Constraining the microphysical parameters and the circumburst density from the afterglow observations, we show that the GeV flare is consistent with a SSC reverse-shock model, the break in the long-lasting emission with the passage of the synchrotron cooling break through the Fermi-LAT band and the very energetic photons with SSC emission from the forward shock when the outflow carries a significant magnetic field ($R_{\rm B} \simeq 30$) and it decelerates in a uniform-density medium with a very low density ($n=4.554^{+1.128}_{-1.121}\times 10^{-4}\,{\rm cm^{-3}}$).

We present the afterglow light curves produced by the deceleration of the non-relativistic ejecta mass in a stratified circumstellar medium with a density profile $n(r)\propto r^{-k}$ with $k=0$, $1$, $1.5$, $2$ and $2.5$. Once the ejecta mass is launched with equivalent kinetic energy parametrized by $E(>\beta)\propto \beta^{-\alpha}$ (where beta is the ejecta velocity) and propagates into the surrounding circumstellar medium, it first moves with constant velocity (the free-coasting phase), and later it decelerates (the Sedov-Taylor expansion). We present the predicted synchrotron and synchrotron-self Compton light curves during the free-coasting phase, and the subsequent Sedov-Taylor expansion. In particular cases, we show the corresponding light curves generated by the deceleration of several ejecta masses with different velocities launched during the coalescence of binary compact objects and the core-collapse of dying massive stars which will contribute at distinct timescales, frequencies, and intensities. We find that before the radiation escapes from the kilonova (KN) and produces a peak in the light curve, the KN afterglow ejecta is very weak and the corresponding ejecta velocity moves as $\propto t^{-\frac{2}{\alpha+1}}$. Finally, using the multi-wavelength observations and upper limits collected by a large campaign of orbiting satellites and ground telescopes, we constrain the parameter space of both the KN afterglow in GW170817 and the possibly generated KN afterglow in S190814bv. Further observations on timescales of years post-merger are needed to derive tighter constraints.

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.

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.

We present a systematic search for short-duration gamma-ray bursts (GRBs) in the local Universe based on 14 years of observations with the Neil Gehrels Swift Observatory. We cross-correlate the GRB positions with the GLADE catalogue of nearby galaxies, and find no event at a distance $\lesssim$100 Mpc and four plausible candidates in the range 100 Mpc$\lesssim$$D$$\lesssim$200 Mpc. Although affected by low statistics, this number is higher than the one expected for chance alignments to random galaxies, and possibly suggests a physical association between these bursts and nearby galaxies. By assuming a local origin, we use these events to constrain the range of properties for X-ray counterparts of neutron star mergers. Optical upper limits place tight constraints on the onset of a blue kilonova, and imply either low masses ($\lesssim10^{-3}\,M_{\odot}$) of lanthanide-poor ejecta or unfavorable orientations ($\theta_{obs}\gtrsim$30 deg). Finally, we derive that the all-sky rate of detectable short GRBs within 200 Mpc is $1.3^{+1.7}_{-0.8}$ yr$^{-1}$ (68\% confidence interval), and discuss the implications for the GRB outflow structure. If these candidates are instead of cosmological origin, we set a upper limit of $\lesssim$2.0 yr$^{-1}$ (90\% confidence interval) to the rate of nearby events detectable with operating gamma-ray observatories, such as Swift and Fermi.

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.

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

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.

The production of both gravitational waves and short gamma-ray bursts (sGRBs) is widely associated with the merger of compact objects. Several studies have modelled the evolution of the electromagnetic emission using the synchrotron emission produced by the deceleration of both a relativistic top-hat jet seen off-axis, and a wide-angle quasi-spherical outflow (both using numerical studies). In this study we present an analytical model of the synchrotron and synchrotron self-Compton (SSC) emission for an off-axis top-hat jet and a quasi-spherical outflow. We calculate the light curves obtained from an analytic model in which the synchrotron and SSC emission (in the fast- or slow-cooling regime) of an off-axis top-hat jet and a quasi-spherical outflow are decelerated in either a homogeneous or a wind-like circumburst medium. We show that the synchrotron emission of the quasi-spherical outflow is stronger than that of the off-axis jet during the first $\sim$ 10 - 20 days, and weaker during the next $\gtrsim$ 80 days. Moreover, we show that if the off-axis jet is decelerated in a wind-like medium, then the SSC emission is very likely to be detected. Applying a MCMC code to our model (for synchrotron emission only), we find the best-fit values for the radio, optical and X-ray emission of GRB 170817A which are in accordance with observations. For GRB 170817A, we find using our model that the synchrotron emission generated by the quasi-spherical outflow and off-axis top-hat jet increase as $F_\nu\propto t^{\alpha}$ with $\alpha\lesssim 0.8$ and $\alpha>3$, respectively. Finally, we obtain the correspondent SSC light curves which are in accordance with the very-high-energy gamma-ray upper limits derived with the GeV - TeV observatories.

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.

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

GRB 190114C, a long and luminous burst, was detected by several satellites and ground-based telescopes from radio wavelengths to GeV gamma-rays. In the GeV gamma-rays, the Fermi LAT detected 48 photons above 1 GeV during the first hundred seconds after the trigger time, and the MAGIC telescopes observed for more than one thousand seconds very-high-energy (VHE) emission above 300 GeV. Previous analysis of the multi-wavelength observations showed that although these are consistent with the synchrotron forward-shock model that evolves from a stratified stellar-wind to homogeneous ISM-like medium, photons above few GeVs can hardly be interpreted in the synchrotron framework. In the context of the synchrotron forward-shock model, we derive the light curves and spectra of the synchrotron self-Compton (SSC) model in the stratified and homogeneous medium. In particular, we study the evolution of these light curves during the stratified-to-homogeneous afterglow transition. Using the best-fit parameters reported for GRB 190114C we interpret the photons beyond the synchrotron limit in the SSC framework and model its spectral energy distribution. We conclude that low-redshift GRBs described under a favourable set of parameters as found in the early afterglow of GRB 190114C could be detected at hundreds of GeVs, and also afterglow transitions would allow that VHE emission could be observed for longer periods.

We present observations of the possible short GRB 180418A in $\gamma$-rays, X-rays, and in the optical. Early optical photometry with the TAROT and RATIR instruments show a bright peak ($\approx$ 14.2 AB mag) between $T+28$ and $T+90$ seconds that we interpret as the signature of a reversal shock. Later observations can be modeled by a standard forward shock model and show no evidence of jet break, allowing us to constrain the jet collimation to $\theta_j> 7^\circ$. Using deep late-time optical observations we place an upper limit of $r>24$ AB mag on any underlying host galaxy. The detection of the afterglow in the \textit{Swift} UV filters constrains the GRB redshift to $z<1.3$ and places an upper bound on the $\gamma$-ray isotropic equivalent energy $E_{\rm{\gamma,iso}} < 3 \times 10^{51}$ erg. The properties of this GRB (e.g. duration, hardness ratio, energetic, and environment) lie at the intersection between short and long bursts, and we can not conclusively identify its type. We estimate that the probability that it is drawn from the population of short GRBs is 10\%-30\%.

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.

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

The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail.

We present Karl G. Jansky Very Large Array radio observations of the long gamma-ray burst GRB 161219B ($z=0.147$) spanning $1-37$ GHz. The data exhibit unusual behavior, including sharp spectral peaks and minutes-timescale large-amplitude variability centered at $20$ GHz and spanning the full frequency range. We attribute this behavior to scattering of the radio emission by the turbulent ionized Galactic interstellar medium (ISM), including both diffractive and refractive scintillation. However, the scintillation is much stronger than predicted by a model of the Galactic electron density distribution (NE2001); from the measured variability timescale and decorrelation bandwidth we infer a scattering measure of $SM\approx {(8-70)\times 10^{-4}}$ kpc m$^{-20/3}$ (up to ${25}$ times larger than predicted in NE2001) and a scattering screen distance of $d_{\rm scr}\approx {0.2-3}$ kpc. We infer an emission region size of $\theta_s \approx {0.9-4}$ $\mu$as ($\approx {(1-4)}\times 10^{16}$ cm) at $\approx4$ days, and find that prior to 8 days the source size is an order of magnitude smaller than model predictions for a uniformly illuminated disk or limb-brightened ring, indicating a slightly off-axis viewing angle or significant substructure in the emission region. Simultaneous multi-hour broadband radio observations of future GRB afterglows will allow us to characterize the scintillation more completely, and hence to probe the observer viewing angle, the evolution of the jet Lorentz factor, the structure of the afterglow emission regions, and ISM turbulence at high Galactic latitudes.

The short GRB 170817A, detected by the Fermi Gamma-ray Burst Monitor, orbiting satellites and ground-based telescopes, was the electromagnetic counterpart of a gravitational-wave transient (GW170817) from a binary neutron star merger. After this merger the $\gamma$-ray light curve exhibited a faint peak at $\sim$ 1.7s and the X-ray, optical and radio light curves displayed an extended emission which increased in brightness up to $\sim$ 160 days. In this paper, we show that the X-ray, optical and radio fluxes are consistent with the synchrotron forward-shock model viewed off-axis when the matter in the outflow is parametrized through a power law velocity distribution. We discuss the origin of the $\gamma$-ray peak in terms of internal and external shocks. We show that the $\gamma$-ray flux might be consistent with a synchrotron self-Compton reverse-shock model observed at high latitudes. Comparing the best-fit values obtained after describing the $\gamma$-ray, X-ray, optical and radio fluxes with our model, we find that the afterglow and $\gamma$-ray emission occurred in different regions and also evidence to propose that the progenitor environment was entrained with magnetic fields and therefore, we argue for the presence of the magnetic field amplification in the binary neutron star merger.

We analyze the Sun as a source for the indirect detection of dark matter through a search for gamma rays from the solar disk. Capture of dark matter by elastic interactions with the solar nuclei followed by annihilation to long-lived mediators can produce a detectable gamma-ray flux. We search three years of data from the High Altitude Water Cherenkov Observatory and find no statistically significant detection of TeV gamma-ray emission from the Sun. Using this, we constrain the spin-dependent elastic scattering cross section of dark matter with protons for dark matter masses above 1 TeV, assuming an unstable mediator with a favorable lifetime. The results complement constraints obtained from Fermi-LAT observations of the Sun and together cover WIMP masses between 4 GeV and $10^6$ GeV. The cross section constraints for mediator decays to gamma rays can be as strong as $\sim10^{-45}$ cm$^{-2}$, which is more than four orders of magnitude stronger than current direct-detection experiments for 1 TeV dark matter mass. The cross-section constraints at higher masses are even better, nearly 7 orders of magnitude better than the current direct-detection constraints for 100 TeV dark matter mass. This demonstration of sensitivity encourages detailed development of theoretical models in light of these powerful new constraints.

Steady gamma-ray emission up to at least 200 GeV has been detected from the solar disk in the Fermi-LAT data, with the brightest, hardest emission occurring during solar minimum. The likely cause is hadronic cosmic rays undergoing collisions in the Sun's atmosphere after being redirected from ingoing to outgoing in magnetic fields, though the exact mechanism is not understood. An important new test of the gamma-ray production mechanism will follow from observations at higher energies. Only the High Altitude Water Cherenkov (HAWC) Observatory has the required sensitivity to effectively probe the Sun in the TeV range. Using three years of HAWC data from November 2014 to December 2017, just prior to the solar minimum, we search for 1--100 TeV gamma rays from the solar disk. No evidence of a signal is observed, and we set strong upper limits on the flux at a few $10^{-12}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV. Our limit, which is the most constraining result on TeV gamma rays from the Sun, is $\sim10\%$ of the theoretical maximum flux (based on a model where all incoming cosmic rays produce outgoing photons), which in turn is comparable to the Fermi-LAT data near 100 GeV. The prospects for a first TeV detection of the Sun by HAWC are especially high during solar minimum, which began in early 2018.

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.

The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.

After two years of operation, the High-Altitude Water Cherenkov (HAWC) Observatory has analyzed the TeV cosmic-ray sky over an energy range between $2.0$ and $72.8$ TeV. The HAWC detector is a ground-based air-shower array located at high altitude in the state of Puebla, Mexico. Using 300 light-tight water tanks, it collects the Cherenkov light from the particles of extensive air showers from primary gamma rays and cosmic rays. This detection method allows for uninterrupted observation of the entire overhead sky (2~sr instantaneous, 8.5~sr integrated) in the energy range from a few TeV to hundreds of TeV. Like other detectors in the northern and southern hemisphere, HAWC observes an energy-dependent anisotropy in the arrival direction distribution of cosmic rays. The observed cosmic-ray anisotropy is dominated by a dipole moment with phase $\alpha\approx40^{\circ}$ and amplitude that slowly rises in relative intensity from $8\times10^{-4}$ at 2 TeV to $14\times10^{-4}$ around 30.3 TeV, above which the dipole decreases in strength. A significant large-scale ($>60^{\circ}$ in angular extent) signal is also observed in the quadrupole and octupole moments, and significant small-scale features are also present, with locations and shapes consistent with previous observations. Compared to previous measurements in this energy range, the HAWC cosmic-ray sky maps improve on the energy resolution and fit precision of the anisotropy. These data can be used in an effort to better constrain local cosmic-ray accelerators and the intervening magnetic fields.

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.

We report on the measurement of the all-particle cosmic ray energy spectrum with the High Altitude Water Cherenkov (HAWC) Observatory in the energy range 10 to 500 TeV. HAWC is a ground based air-shower array deployed on the slopes of Volcan Sierra Negra in the state of Puebla, Mexico, and is sensitive to gamma rays and cosmic rays at TeV energies. The data used in this work were taken from 234 days between June 2016 to February 2017. The primary cosmic-ray energy is determined with a maximum likelihood approach using the particle density as a function of distance to the shower core. Introducing quality cuts to isolate events with shower cores landing on the array, the reconstructed energy distribution is unfolded iteratively. The measured all-particle spectrum is consistent with a broken power law with an index of $-2.49\pm0.01$ prior to a break at $(45.7\pm0.1$) TeV, followed by an index of $-2.71\pm0.01$. The spectrum also respresents a single measurement that spans the energy range between direct detection and ground based experiments. As a verification of the detector response, the energy scale and angular resolution are validated by observation of the cosmic ray Moon shadow's dependence on energy.

Gamma-ray bursts (GRBs) are among the most luminous sources in the universe and the nature of their emission up to very high energy is one of the most important open issue connected with the study of these peculiar events. The High Altitude Water Cherenkov (HAWC) gamma-ray observatory, installed at an altitude of 4100 m a. s. l. in the state of Puebla (Mexico), has completed its second year of full operations. Thanks to its instantaneous field of view of ~2 sr and its high duty cycle ($\ge$ 95%), HAWC is an ideal instrument for the study of transient phenomena such as GRBs. We performed a search for TeV emission delayed with respect to, and of longer duration than the prompt emission observed by satellites. We present here the results obtained by observing at the position of a sample of GRBs detected by the Fermi and Swift satellites from December 2014 to February 2017. The upper limits resulting from this analysis are presented and theoretical implications are discussed.

We present multi-wavelength observations and modeling of the exceptionally bright long $\gamma$-ray burst GRB 160625B. The optical and X-ray data are well-fit by synchrotron emission from a collimated blastwave with an opening angle of $\theta_j\approx 3.6^\circ$ and kinetic energy of $E_K\approx 2\times10^{51}$ erg, propagating into a low density ($n\approx 5\times10^{-5}$ cm$^{-3}$) medium with a uniform profile. The forward shock is sub-dominant in the radio band; instead, the radio emission is dominated by two additional components. The first component is consistent with emission from a reverse shock, indicating an initial Lorentz factor of $\Gamma_0\gtrsim 100$ and an ejecta magnetization of $R_B\approx 1-100$. The second component exhibits peculiar spectral and temporal evolution and is most likely the result of scattering of the radio emission by the turbulent Milky Way interstellar medium (ISM). Such scattering is expected in any sufficiently compact extragalactic source and has been seen in GRBs before, but the large amplitude and long duration of the variability seen here are qualitatively more similar to extreme scattering events previously observed in quasars, rather than normal interstellar scintillation effects. High-cadence, broadband radio observations of future GRBs are needed to fully characterize such effects, which can sensitively probe the properties of the ISM and must be taken into account before variability intrinsic to the GRB can be interpreted correctly.

List of proceedings from the HAWC Collaboration presented at the 35th International Cosmic Ray Conference, 12 July - 20 July 2017, Bexco, Busan, Korea.

The origin of the gamma-ray burst (GRB) prompt emission still defies explanation, in spite of recent progress made, for example, on the occasional presence of a thermal component in the spectrum along with the ubiquitous non-thermal component that is modelled with a Band function. The combination of finite duration and aperiodic modulations make GRBs hard to characterise temporally. Although correlations between GRB luminosity and spectral hardness on one side and time variability on the other side have long been known, the loose and often arbitrary definition of the latter makes the interpretation uncertain. We characterise the temporal variability in an objective way and search for a connection with rest-frame spectral properties for a number of well-observed GRBs. We studied the individual power density spectra (PDS) of 123 long gamma-ray bursts with measured redshift, rest-frame peak energy Ep,i of the time-averaged nuFnu spectrum, and well-constrained PDS slope alpha detected with Swift, Fermi and past spacecraft. The PDS were modelled with a power law either with or without a break adopting a Bayesian Markov chain Monte Carlo technique. We find a highly significant Ep,i-alpha anti-correlation. The null hypothesis probability is ~ 10^-9. In the framework of the internal shock synchrotron model, the Ep,i-alpha anti-correlation can hardly be reconciled with the predicted Ep,i propto Gamma^-2, unless either variable microphysical parameters of the shocks or continual electron acceleration are assumed. Alternatively, in the context of models based on magnetic reconnection, the PDS slope and Ep,i are linked to the ejecta magnetisation at the dissipation site, so that more magnetised outflows would produce more variable GRB light curves at short timescales (<= 1 s), shallower PDS, and higher values of Ep,i.

Timing analysis is a powerful tool with which to shed light on the still obscure emission physics and geometry of the prompt emission of GRBs. Fourier power density spectra (PDS) characterise time series as stochastic processes and can be used to search for coherent pulsations and to investigate the dominant variability timescales. Because of the limited duration and of the statistical properties, modelling the PDS of individual GRBs is challenging, and only average PDS of large samples have been discussed in the literature. We characterise the individual PDS of GRBs in terms of a stochastic process, and carry out for the first time a systematic search for periodic signals and for a link between the PDS and other observables. We present a Bayesian procedure that uses a Markov chain Monte Carlo technique and apply it to study 215 bright long GRBs detected with the Swift Burst Alert Telescope from January 2005 to May 2015. The PDS are modelled with a power-law either with or without a break. Two classes of GRBs emerge: with or without a unique dominant timescale. A comparison with active galactic nuclei (AGNs) reveals similar distributions of PDS slopes. Unexpectedly, GRBs with subsecond-dominant timescales and duration longer than a few tens of seconds in the source frame appear to be either very rare or altogether absent. Three GRBs are found with possible evidence for a periodic signal at ~3 sigma (Gaussian) significance, corresponding to a multitrial chance probability of ~1%. Thus, we found no compelling evidence for periodic signals. The analogy between the PDS of GRBs and of AGNs could tentatively indicate similar stochastic processes that rule BH accretion across different BH mass scales and objects. In addition, we find evidence that short dominant timescales and duration are not completely independent of each other, in contrast with commonly accepted paradigms (abridged).

Gamma-ray burst GRB 140430A was detected by the Swift satellite and observed promptly with the imaging polarimeter RINGO3 mounted on the Liverpool Telescope, with observations beginning while the prompt $\gamma$-ray emission was still ongoing. In this paper, we present densely sampled (10-second temporal resolution) early optical light curves in 3 optical bands and limits to the degree of optical polarization. We compare optical, X-ray and gamma-ray properties and present an analysis of the optical emission during a period of high-energy flaring. The complex optical light curve cannot be explained merely with a combination of forward and reverse shock emission from a standard external shock, implying additional contribution of emission from internal shock dissipation. We estimate an upper limit for time averaged optical polarization during the prompt phase to be as low as P < 12% (1$\sigma$). This suggests that the optical flares and early afterglow emission in this GRB are not highly polarized. Alternatively, time averaging could mask the presence of otherwise polarized components of distinct origin at different polarization position angles.

Aims: We characterize a sample of Gamma-Ray Bursts with low luminosity X-ray afterglows (LLA GRBs), and study their properties. Method: We select a sample consisting of the 12\% faintest X-ray afterglows from the total population of long GRBs (lGRBs) with known redshift. We study their intrinsic properties (spectral index, decay index, distance, luminosity, isotropic radiated energy and peak energy) to assess whether they belong to the same population than the brighter afterglow events. Results: We present strong evidences that these events belong to a population of nearby events, different from that of the general population of lGRBs. These events are faint during their prompt phase, and include the few possible outliers of the Amati relation. Out of 14 GRB-SN associations, 9 are in LLA GRB sample, prompting for caution when using SN templates in observational and theoretical models for the general lGRBs population.

Prompt gamma-ray and early X-ray afterglow emission in gamma-ray bursts (GRBs) are characterized by a bursty behavior and are often interspersed with long quiescent times. There is compelling evidence that X-ray flares are linked to prompt gamma-rays. However, the physical mechanism that leads to the complex temporal distribution of gamma-ray pulses and X-ray flares is not understood. Here we show that the waiting time distribution (WTD) of pulses and flares exhibits a power-law tail extending over 4 decades with index ~2 and can be the manifestation of a common time-dependent Poisson process. This result is robust and is obtained on different catalogs. Surprisingly, GRBs with many (>=8) gamma-ray pulses are very unlikely to be accompanied by X-ray flares after the end of the prompt emission (3.1 sigma Gaussian confidence). These results are consistent with a simple interpretation: an hyperaccreting disk breaks up into one or a few groups of fragments, each of which is independently accreted with the same probability per unit time. Prompt gamma-rays and late X-ray flares are nothing but different fragments being accreted at the beginning and at the end, respectively, following the very same stochastic process and likely the same mechanism.

Most Gamma-Ray Bursts (GRBs) detected by the Fermi Gamma-ray Space Telescope exhibit a delay of up to about 10 seconds between the trigger time of the hard X-ray signal as measured by the Fermi GBM and the onset of the MeV-GeV counterpart detected by the LAT. This delay may hint at important physics, whether it is due to the intrinsic variability of the inner engine or it is related to quantum dispersion effects in the velocity of light propagation from the sources to the observer. It is critical to have a proper assessment of how these time delays affect the overall properties of the light curves. We cross-correlated the 5 brightest GRBs of the 1st Fermi LAT Catalog by means of the continuous correlation function (CCF) and of the Discrete Correlation Function (DCF). A maximum in the DCF suggests the presence of a time lag between the curves, whose value and uncertainty are estimated through a Gaussian fitting of the DCF profile and light curve simulation via a Monte Carlo approach. The cross-correlation of the observed LAT and GBM light curves yields time lags that are mostly similar to those reported in the literature, but they are formally consistent with zero. The cross-correlation of the simulated light curves yields smaller errors on the time lags and more than one time lag for GRBs 090902B and 090926A; for all 5 GRBs, the time lags are significantly different from zero and consistent with those reported in the literature, when only the secondary maxima are considered for those two GRBs. The DCF method evidences the presence of time lags between the LAT and GBM light curves and underlines their complexity. While this suggests that the delays should be ascribed to intrinsic physical mechanisms, more sensitivity and larger statistics are needed to assess whether time lags are universally present in the early GRB emission and which dynamical time scales they trace.

We report results of a systematic study of the broad--band (2--2000 keV) time--resolved prompt emission spectra of a sample of Gamma-Ray Bursts (GRBs) detected with both the Wide Field Cameras (WFCs) aboard the BeppoSAX satellite and the BATSE experiment aboard CGRO. The main goal of the paper is to test spectral models of the GRB prompt emission that have recently been proposed. In particular, we test the photospheric model proposed by Ryde and Pe'er (2009), i.e., blackbody plus power--law, the addition of a blackbody emission to the Band function in the cases in which this function does not fit the data, and the Comptonization model developed by Titarchuk et al. (2012). By considering the few spectra for which the simple Band function does not provide a fully acceptable fit to the data (Frontera et al. 2012), only in one case we find a statistically significant better fit by adding a blackbody to this function. We confirm the results found fitting the BATSE spectra alone with a blackbody plus a power law. Instead when the BATSE GRB spectra are joined to those obtained with WFCs (2--28 keV), this model becomes unacceptable in most of time intervals in which we subdivide the GRB light curves. We find instead that the Comptonization model is always acceptable, even in the few cases in which the Band function is inconsistent with the data. We discuss the implications of these results.

We searched for periodic and quasiperiodic signal in the prompt emission of a sample of 44 bright short gamma-ray bursts detected with Fermi/GBM, Swift/BAT, and CGRO/BATSE. The aim was to look for the observational signature of quasiperiodic jet precession which is expected from black hole-neutron star mergers, but not from double neutron star systems. Thus, this kind of search holds the key to identify the progenitor systems of short GRBs and, in the wait for gravitational wave detection, represents the only direct way to constrain the progenitors. We tailored our search to the nature of the expected signal by properly stretching the observed light curves by an increasing factor with time, after calibrating the technique on synthetic curves. In none of the GRBs of our sample we found evidence for periodic or quasiperiodic signals. In particular, for the 7 unambiguously short GRBs with best S/N we obtained significant upper limits to the amplitude of the possible oscillations. This result suggests that BH-NS systems do not dominate the population of short GRB progenitors as described by the kinematic model of Stone, Loeb, & Berger (2013).

From past experiments the average power density spectrum (PDS) of GRBs with unknown redshift was found to be modelled from 0.01 to 1 Hz with a power-law, f^(-alpha), with alpha broadly consistent with 5/3. Recent analyses of the Swift/BAT catalogue showed analogous results in the 15-150 keV band. We carried out the same analysis on the bright GRBs detected by BeppoSAX/GRBM and Fermi/GBM. The BeppoSAX/GRBM data, in the energy range 40-700 keV and with 7.8 and 0.5-ms time resolutions, allowed us to explore for the first time the average PDS at very high frequencies (up to 1 kHz) and reveal a break around 1-2 Hz, previously found in CGRO/BATSE data. The Fermi/GBM data, in the energy band 8-1000 keV, allowed us to explore for the first time the average PDS within a broad energy range. Our results confirm and extend the energy dependence of the PDS slope, according to which harder photons have shallower PDS.