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.

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.

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.

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

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.

The early optical emission of gamma-ray bursts gives an opportunity to understand the central engine and first stages of these events. About 30\% of GRBs present flares whose origin is still a subject of discussion. We present optical photometry of GRB 180620A with the COATLI telescope and RATIR instrument. COATLI started to observe from the end of prompt emission at $T+39.3$~s and RATIR from $T+121.4$~s. We supplement the optical data with the X-ray light curve from \emph{Swift}/XRT. %The optical and X-ray light curves show very unusual behavior with features clearly beyond the standard fireball model. We observe an optical flare from $T+110$ to $T+550$~s, with a temporal index decay $\alpha_\mathrm{O,decay}=1.32\pm 0.01$, and a $\Delta t/t=1.63$, which we interpret as the signature of a reverse shock component. After the initial normal decay the light curves show a long plateau from $T+500$ to $T+7800$~s both in X-rays and the optical before decaying again after an achromatic jet break at $T+7800$~s. Fluctuations are seen during the plateau phase in the optical. Adding to the complexity of GRB afterglows, the plateau phase (typically associated with the coasting phase of the jet) is seen in this object after the ``normal'' decay phase (emitted during the deceleration phase of the jet) and the jet break phase occurs directly after the plateau. We suggest that this sequence of events can be explained by a rapid deceleration of the jet with $t_d\lesssim 40$ s due to the high density of the environment ($\approx 100$ cm$^{-3}$) followed by reactivation of the central engine which causes the flare and powers the plateau phase.

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.

The spectral energy distribution (SED) of blazars is generally understood through the standard one-zone synchrotron self-Compton (SSC) model, where a strong correlation between X-ray and TeV gamma-ray fluxes is expected. Recently, for Mrk 421 the correlation was confirmed and interpreted within the SSC model arguing that a description of the correlation up to large flux values considerably restricts the value of the magnetic field under the assumption of an electron spectral index independent of the flux state. In this work, we extend the analysis done on Mrk 421 for other two high-energy peaked BL Lacs: Mrk 501 and 1ES 1959+650.

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.

With an instantaneous field of view of 2 sr and a duty cycle > 95%, the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory is a perfect instrument for monitoring variable TeV sources. Because radio galaxies are a type of Active Galactic Nuclei (AGN) with their jets misaligned with respect to our line of sight, they may help us to probe the physics of very-high-energy (VHE) emission processes in these objects. Three out of four radio galaxies that have been detected at TeV energies by other facilities are located within the field of view of the HAWC Observatory: M87, NGC 1275, and 3C 264. A search for TeV gamma rays at their locations yields no statistically significant excess of counts. We present corresponding upper limits for each radio galaxy and light curves covering 3 years of data taken with HAWC.

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

The HAWC gamma ray observatory is located at the Sierra Negra Volcano in Puebla, Mexico, at an altitude of 4,100 meters. HAWC is a wide field of view array of 300 water Cherenkov detectors that are continuously surveying 2sr of the sky since March 2015. The large collected data sample consisting in more than 3 years makes HAWC an ideal instrument to perform an unbiased monitoring of blazars in the very-high-energy (VHE) emission. This is particular relevant for Mrk 421, one of the closest and brightest blazars in the gamma-ray/X-ray classified as highsynchrotron-peaked BL Lac class. In this work we present light curves for Mrk 421 and the Crab nebula obtained with the first 17 months of data and a new analysis framework. We compare the results with the light curves reported in [1]. The main advantage of the new framework is the capability to derive fluxes for arbitrary timescales. We show that both, previous and present, methods give consistent results

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.

Detections of very-high-energy (VHE; > 100 GeV) photons from Gamma-ray bursts (GRBs) can provide fundamental information on the involved radiative processes, physical composition of the ejecta and acceleration processes. The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is the best gamma-ray instrument to study transient phenomena over a long period of time from 100 GeV to 100 TeV. Its large field of view and duty cycle (> 95%) allow it to search blindly for sources of GRBs and to follow up on external alerts from satellite instruments, such as Fermi and Swift, searching for their VHE counterpart. We present results from the on-going GRB monitoring program and VHE upper limits of the latest interesting low-redshift and/or powerful bursts in HAWC's field of view: GRB 170206A, GRB 170817A, GRB 180720B and GRB 190114C.

We study the spatial correlation of astrophysical neutrinos detected by IceCube with the geometry of the two large globular structures located in the center of our Galaxy, known as Fermi Bubbles (FB). Using the Fermi-LAT data collected during 8 years and the upper limits derived by the High Altitude Water Cherenkov (HAWC) gamma-ray observatory, we use a hybrid (lepto-hadronic) model to investigate a possible correlation with the high-energy neutrinos in the direction of the Fermi Bubbles. We find that these events are possibly not associated with the Northern Bubble but do not dismiss a possible correlation with the Southern globular structure. We expect in the coming years to improve the gamma-ray observations through the Cerenkov Telescope Array (CTA) and the Southern Gamma-Ray Survey Observatory (SGSO) observatories to test a possible hadronic emission with the Southern Bubble.

With the most recent multi-messenger detection, a new branch in modern astronomy has arisen. The GW170817 event together with the short gamma-ray burst GRB 170817A was the first-ever detection of the gravitational waves and an electromagnetic counterpart. These detections encourage us to think that in the following years we will detect a single event through three different channels: including the mentioned above plus neutrinos from multiple astrophysical sources, like those detected from SN1987A. It is believed that short GRBs are originated in the merger of a black-hole (BH) with a neutron star (NS) or NS-NS scenario. Particularly only in the latter case, several simulations suggest that the magnetic field can be amplified up to $10^{16}$ G. Considering this effect over created thermal neutrinos during the initial stage, we could differentiate short GRB progenitors through the neutrino expected flavor ratio and the opacity created by the baryon-loaded winds ejected in each scenario. Moreover, We find that it is more feasible to detect neutrinos from BH-NS than NS-NS systems. Finally, we also estimate the number of neutrino events expected on ground-based detectors, finding that it is possible to detect neutrinos from an energetic enough source located within a nearby vicinity with Hyper-Kamiokande detector.

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\%.

We describe plans for the development of the Southern Wide-field Gamma-ray Observatory (SWGO), a next-generation instrument with sensitivity to the very-high-energy (VHE) band to be constructed in the Southern Hemisphere. SWGO will provide wide-field coverage of a large portion of the southern sky, effectively complementing current and future instruments in the global multi-messenger effort to understand extreme astrophysical phenomena throughout the universe. A detailed description of science topics addressed by SWGO is available in the science case white paper [1]. The development of SWGO will draw on extensive experience within the community in designing, constructing, and successfully operating wide-field instruments using observations of extensive air showers. The detector will consist of a compact inner array of particle detection units surrounded by a sparser outer array. A key advantage of the design of SWGO is that it can be constructed using current, already proven technology. We estimate a construction cost of 54M USD and a cost of 7.5M USD for 5 years of operation, with an anticipated US contribution of 20M USD ensuring that the US will be a driving force for the SWGO effort. The recently formed SWGO collaboration will conduct site selection and detector optimization studies prior to construction, with full operations foreseen to begin in 2026. Throughout this document, references to science white papers submitted to the Astro2020 Decadal Survey with particular relevance to the key science goals of SWGO, which include unveiling Galactic particle accelerators [2-10], exploring the dynamic universe [11-21], and probing physics beyond the Standard Model [22-25], are highlighted in red boldface.

Centaurus B (Cen B) is one of the closest and brightest radio-loud galaxy in the southern sky. This radio galaxy, proposed as a plausible candidate for accelerating ultra-high-energy cosmic rays (UHECRs), is near the highest-energy neutrino event reported (IC35) in the High-Energy Starting Events catalog. Pierre Auger observatory reported the highest energy comic rays during 10 years of collecting data with some of them around this source. In this paper, the analysis of the gamma-ray spectrum and the light curve above 200 MeV is presented with nine years of cumulative Fermi-LAT data around Cen B. Taking into consideration the multi-wavelength observations carried out about this radio galaxy, leptonic and hadronic scenarios are introduced in order to fit the spectral energy distribution, assuming that the gamma-ray flux is produced in a region close to the core and in the extended lobes. Using the best-fit values found, several physics properties of this radio galaxy are derived. Furthermore, a statistical analysis of the cosmic ray distribution around Cen B is performed, finding that this distribution is not different from the background at a level of significance of 5%. Considering the UHECR event associated to this source by Moskalenko et al. and extrapolating its luminosity to low energies, we do not find enough evidence to associate the highest-energy neutrino event (IC35) with this radio galaxy.

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

We study neutrino production, propagation, and oscillations within an extremely magnetized background of finite-temperature nuclear matter. We focus on three particularly interesting cases and identify the astrophysical scenarios where such a signal may be found. The first case involves nuclear matter with electrons, and it is found during the central-engine stage of, both, short and long gamma-ray bursts (GRBs). Thus, for the short GRB case it will also be associated to gravitational-wave events where there is an electromagnetic counterpart (e.g., GW170817). The second and third scenarios involve the presence of strange-quark matter (SQM). The second scenario occurs if SQM can become negatively charged (SQM$^-$; which may only occur at high pressure) and, thus, it is embedded in a positron plasma. The third case may be found at the interphase where SQM transitions from positive (SQM$^+$) to negative; here, positrons and electrons may constantly annihilate and give a distinctive neutrino signature. Therefore, this may also be a signature of the existence of strange stars. Given the wide range of magnetic fields we find in the literature, we also briefly discuss the maximum limit that a stellar mass compact object may posses.

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 outline the science motivation for SGSO, the Southern Gamma-Ray Survey Observatory. SGSO will be a next-generation wide field-of-view gamma-ray survey instrument, sensitive to gamma-rays in the energy range from 100 GeV to hundreds of TeV. Its science topics include unveiling galactic and extragalactic particle accelerators, monitoring the transient sky at very high energies, probing particle physics beyond the Standard Model, and the characterization of the cosmic ray flux. SGSO will consist of an air shower detector array, located in South America. Due to its location and large field of view, SGSO will be complementary to other current and planned gamma-ray observatories such as HAWC, LHAASO, and CTA.

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

We present optical photometry of the afterglow of the long GRB 180205A with the COATLI telescope from 217 seconds to about 5 days after the {\itshape Swift}/BAT trigger. We analyse this photometry in the conjunction with the X-ray light curve from {\itshape Swift}/XRT. The late-time light curves and spectra are consistent with the standard forward-shock scenario. However, the early-time optical and X-ray light curves show non-typical behavior; the optical light curve exhibits a flat plateau while the X-ray light curve shows a flare. We explore several scenarios and conclude that the most likely explanation for the early behavior is late activity of the central engine.

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.

Multi-wavelength campaigns have been carried out to study the correlation between the very high energy (VHE) $\gamma$-ray and the X-ray emissions in blazars but, no conclusive results have been achieved yet. In this paper, we add Milagro data to the existing VHE $\gamma$-ray data from HEGRA-CT1 and Whipple and test the consistency and robustness of the reported correlation between VHE $\gamma$-ray and X-ray fluxes in Mrk 421. We found that at monthly time scale the correlation is robust, consistent between instruments and described as a linear function. Furthermore, most of the fluxes on shorter time scales are consistent with the correlation within 3 $\sigma_A$ even, where $\sigma_A$ is an estimated intrinsic scatter. However, a break-down of the correlation becomes clearly evident at high states of activity with fluxes $\rm \gtrsim 2.5\times 10^{-10}\, cm^{-2}s^{-1}$ at energies above 400 GeV independently of the time scale, observational period or instrument, even for single flares, the X-ray and VHE $\gamma$-ray emissions lie on the correlation until the VHE $\gamma$-ray flux reaches values higher than the one mentioned above. We have interpreted our results within the one-zone synchrotron self-Compton model. We found that describing a single and unique $\gamma$-ray/X-ray correlation strongly narrows the range of possible values of the magnetic field $B$ when a constant value of the spectral index along the correlation is assumed.

Binary neutron star mergers are believed to eject significant masses with a diverse range of velocities. Once these ejected materials begin to be decelerated by a homogeneous medium, relativistic electrons are mainly cooled down by synchrotron radiation, generating a multiwavelength long-lived afterglow. Analytic and numerical methods illustrate that the outermost matter, the merger shock-breakout material, can be parametrized by power-law velocity distributions $\propto \left(\beta_{\rm c}\Gamma \right)^{-\alpha_s}$. Considering that the shock-breakout material is moving on-axis towards the observer and the relativistic jet off-axis, we compute the light curves during the relativistic and the lateral expansion phase. As a particular case, we successfully describe the X-ray, optical and radio light curves alongside the spectral energy distribution from the recently discovered gravitational-wave transient GW170817, when the merger shock-breakout material moves with mildly relativistic velocities near-Newtonian phase and the jet with relativistic velocities. Future electromagnetic counterpart observations of this binary system could be able to evaluate different properties of these light curves.

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.

Recent X-ray observations in the central compact remnant of Puppis A and other young pulsars have provided convincing evidence about the anti-magnetar nature of the so-called central compact objects (CCOs). The measured period derivative, corrected by the proper movement, provides exceptionally low superficial magnetic fields for such sources. Using the dipole radiation canonical model, the superficial magnetic field of the CCO (RX J0822-4300) in Puppis A was estimated to be $B\simeq 2.9 \times 10^{10}$ G. In this paper, we present a numerical model to interpret the unusual magnetic-field behavior in RX J0822-4300. Using the magnetohydrodynamic simulations during the hypercritical accretion phase, we propose that a variety of MeV neutrinos is created given evidence about the submergence of the magnetic field in the pulsar. We investigate the production, propagation, oscillations, and detection of MeV neutrinos from this source. The detection of neutrino events with their flavor ratios would be a compelling signature of the decreasing evolution of the magnetic field not only in RX J0822-4300 but also in CCO candidates.

The measurement of the period and period derivative, and the canonical model of dipole radiation have provided a method to estimate the low superficial magnetic fields in the so-called Central Compact Objects (CCOs). In the present work, a scenario is introduced in order to explain the magnetic behavior of such CCOs. Based on magnetohydrodynamic simulations of the post core-collapse supernova phase during the hypercritical accretion episode, we argue that the magnetic field of a newborn neutron star could have been early buried. During this phase, thermal neutrinos are created mainly by the pair annihilation, plasmon decay, photo-neutrino emission and other processes. We study the dynamics of these neutrinos in this environment and also estimate the number expected of the neutrino events with their flavor ratios on Earth. The neutrino burst is the only viable observable that could provide compelling evidence of the hypercritical phase and therefore, the hidden magnetic field mechanism as the most favorable scenario to explain the anomalous low magnetic fields estimated for CCOs.

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.

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.

Correlations between optical flashes and gamma-ray emissions in gamma-ray bursts have been searched in order to clarify the question whether these emissions occur at internal and/or external shocks. Among the most powerful gamma-ray bursts ever recorded are GRB 080319B and GRB 130427A which at early phase presented bright optical flashes possible correlated with $\gamma$-ray components. Additionally, both bursts were fortuitously located within the field of view of the TeV $\gamma$-ray Milagro and HAWC observatories, and although no statistically significant excess of counts were collected, upper limits were placed on the GeV - TeV emission. Considering the synchrotron self-Compton emission from internal shocks and requiring the GeV-TeV upper limits we found that the optical flashes and the $\gamma$-ray components are produced by different electron populations. Analyzing the optical flashes together the multiwavelength afterglow observation, we found that these flashes can be interpreted in the framework of the synchrotron reverse-shock model when outflows have arbitrary magnetizations.

For the first time, a short gamma-ray burst (GRB) was unambiguously associated with a gravitational wave (GW) observation from a binary neutron star (NS) merger. This allows us to link the details of the central engine properties to GRB emission models. We find that photospheric models (both dissipative and non-dissipative variants) have difficulties accounting for the observations. Internal shocks give the most natural account of the observed peak energy, viewing angle and total energy. We also show that a simple external shock model can reproduce the observed GRB pulse with parameters consistent with those derived from the afterglow modeling. We find a simple cocoon shock breakout model is in mild tension with the observed spectral evolution, however it cannot be excluded based on gamma-ray data alone. Future joint observations of brighter GRBs will pose even tighter constraints on prompt emission models.

The High Altitude Water Cherenkov observatory (HAWC) is an air shower array devised for TeV gamma-ray astronomy. HAWC is located at an altitude of 4100 m a.s.l. in Sierra Negra, Mexico. HAWC consists of 300 Water Cherenkov Detectors, each instrumented with 4 photomultiplier tubes (PMTs). HAWC re-uses the Front-End Boards from the Milagro experiment to receive the PMT signals. These boards are used in combination with Time to Digital Converters (TDCs) to record the time and the amount of light in each PMT hit (light flash). A set of VME TDC modules (128 channels each) is operated in a continuous (dead time free) mode. The TDCs are read out via the VME bus by Single-Board Computers (SBCs), which in turn are connected to a gigabit Ethernet network. The complete system produces ~ 500 MB/s of raw data. A high-throughput data processing system has been designed and built to enable real-time data analysis. The system relies on off-the-shelf hardware components, an open-source software technology for data transfers (ZeroMQ) and a custom software framework for data analysis (AERIE). Multiple trigger and reconstruction algorithms can be combined and run on blocks of data in a parallel fashion, producing a set of output data streams which can be analyzed in real time with minimal latency (< 5 s). This paper provides an overview of the hardware set-up and an in-depth description of the software design, covering both the TDC data acquisition system and the real-time data processing system. The performance of these systems is also discussed.

On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.

The unexpectedly high flux of cosmic ray positrons detected at Earth may originate from nearby astrophysical sources, dark matter, or unknown processes of cosmic-ray secondary production. We report the detection, using the HighAltitude Water Cherenkov Observatory (HAWC), of extended tera-electron volt gamma-ray emission coincident with the locations of two nearby middle-aged pulsars (Geminga and PSR B0656+14). The HAWC observations demonstrate that these pulsars are indeed local sources of accelerated leptons, but the measured tera-electron volt emission profile constrains the diffusion of particles away from these sources to be much slower than previously assumed. We demonstrate that the leptons emitted by these objects are therefore unlikely to be the origin of the excess positrons, which may have a more exotic origin.

The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is a wide field-of-view observatory sensitive to 500 GeV - 100 TeV gamma rays and cosmic rays. With its observations over 2/3 of the sky every day, the HAWC observatory is sensitive to a wide variety of astrophysical sources, including possible gamma rays from dark matter. Dark matter annihilation and decay in the Milky Way Galaxy should produce gamma-ray signals across many degrees on the sky. The HAWC instantaneous field-of-view of 2 sr enables observations of extended regions on the sky, such as those from dark matter in the Galactic halo. Here we show limits on the dark matter annihilation cross-section and decay lifetime from HAWC observations of the Galactic halo with 15 months of data. These are some of the most robust limits on TeV and PeV dark matter, largely insensitive to the dark matter morphology. These limits begin to constrain models in which PeV IceCube neutrinos are explained by dark matter which primarily decays into hadrons.

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.

GRB 110721A was detected by the Gamma-ray Burst Monitor and the Large Area Telescope (LAT) onboard the Fermi satellite and the Gamma-ray Burst Polarimeter onboard the IKAROS solar mission. Previous analysis done of this burst showed: i) a linear polarization signal with position angle stable ($\phi_p= 160^\circ\pm11$) and high degree of $\Pi=84^{+16}_{-28}$, ii) an extreme peak energy of a record-breaking at 15$\pm$2 MeV, and iii) a subdominant prompt thermal component observed right after the onset of this burst. In this paper, the LAT data around the reported position of GRB 110721A are analysed with the most recent software and then, the LAT light curve above 100 MeV was obtained. The LAT light curve is modelled in terms of adiabatic early-afterglow external shocks when the outflow propagates into a stellar wind. Additionally, we discuss the possible origins and also study the implications of the early-afterglow external shocks on the extreme peak energy observed at 15$\pm$2 MeV, the polarization observations and the subdominant prompt thermal component.

The High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory is an extensive air shower array located in Puebla, Mexico. The closest radio galaxy within the HAWC field of view, M87, has been detected in very high energies. In this work we report upper limits on the TeV {\gamma}-ray flux of the radio galaxy M87. At a distance of 16 Mpc, M87 is a supergiant elliptical galaxy located in the Virgo Cluster that has been observed from radio wavelengths to TeV {\gamma}-rays. Although a single-zone synchrotron self-Compton model has been successfully used to explain the spectral energy distribution of this source up to a few GeV, the {\gamma}-ray spectrum at TeV has been interpreted within different theoretical models. We discuss the implications of these upper limits on the photo-hadronic interactions, as well as the number of neutrino events expected in the IceCube neutrino telescope.

GRB 160625B, one of the brightest bursts in recent years, was simultaneously observed by Fermi and Swift satellites, and ground-based optical telescopes in three different events separated by long periods of time. In this paper the non-thermal multiwavelength observations of GRB 160625B are described and a transition phase from wind-type-like medium to interstellar medium between the early (event II) and the late (event III) afterglow is found. The multiwavelength observations of the early afterglow are consistent with the afterglow evolution starting at $\sim$ 150 s in a stellar wind medium whereas the observations of the late afterglow are consistent with the afterglow evolution in interstellar medium (ISM). The wind-to-ISM transition is calculated to be at $\sim 8\times 10^3$ s when the jet has decelerated, at a distance of $\sim$ 1 pc from the progenitor. Using the standard external shock model, the synchrotron and synchrotron self-Compton emission from reverse shock is required to model the GeV $\gamma$-ray and optical observations in the early afterglow, and synchrotron radiation from the adiabatic forward shock to describe the X-ray and optical observations in the late afterglow. The derived values of the magnetization parameter, the slope of the fast decay of the optical flash and the inferred magnetic fields suggest that Poynting flux-dominated jet models with arbitrary magnetization could account for the spectral properties exhibited by GRB 160625B.

The HAWC gamma ray observatory is located at the Sierra Negra Volcano in Puebla, Mexico, at an altitude of 4,100 meters. HAWC is a wide field of view array of 300 water Cherenkov detectors that are continuously surveying ~ 2sr of the sky, operating since March 2015. The large collected data sample allows HAWC to perform an unbiased monitoring of the BL Lac Mrk 421. This is the closest and brightest known extragalactic high-synchrotron-peaked BL Lac in the gamma-ray/X- ray bands and is extensively monitored by the Large Area Telescope (LAT) on-board the Fermi satellite, and the BAT and XRT instruments of the Swift satellite. In this work, we use 25 months of HAWC data together with Swift-XRT data to characterize potential correlations between both wavelengths. This analysis shows that HAWC and Swift-XRT data are correlated even stronger than expected for quasi-simultaneous observations.

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

The High Altitude Water Cherenkov (HAWC) Gamma-ray Observatory is an extensive air shower detector operating in central Mexico, which has recently completed its first two years of full operations. If for a burst like GRB 130427A at a redshift of 0.34 and a high-energy component following a power law with index -1.66, the high-energy component is extended to higher energies with no cut-off other than from extragalactic background light attenuation, HAWC would observe gamma rays with a peak energy of $\sim$300 GeV. This paper reports the results of HAWC observations of 64 gamma-ray bursts (GRBs) detected by $\mathit{Swift}$ and $\mathit{Fermi}$, including three GRBs that were also detected by the Large Area Telescope ($\mathit{Fermi}$-LAT). An ON/OFF analysis method is employed, searching on the time scale given by the observed light curve at keV-MeV energies and also on extended time scales. For all GRBs and time scales, no statistically significant excess of counts is found and upper limits on the number of gamma rays and the gamma-ray flux are calculated. GRB 170206A, the third brightest short GRB detected by the Gamma-ray Burst Monitor on board the $\mathit{Fermi}$ satellite ($\mathit{Fermi}$-GBM) and also detected by the LAT, occurred very close to zenith. The LAT measurements can neither exclude the presence of a synchrotron self-Compton (SSC) component nor constrain its spectrum. Instead, the HAWC upper limits constrain the expected cut-off in an additional high-energy component to be less than $100~\rm{GeV}$ for reasonable assumptions about the energetics and redshift of the burst.

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.

The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is a wide field of view observatory sensitive to 500 GeV - 100 TeV gamma rays and cosmic rays. It can also perform diverse indirect searches for dark matter (DM) annihilation and decay. Among the most promising targets for the indirect detection of dark matter are dwarf spheroidal galaxies. These objects are expected to have few astrophysical sources of gamma rays but high dark matter content, making them ideal candidates for an indirect dark matter detection with gamma rays. Here we present individual limits on the annihilation cross section and decay lifetime for 15 dwarf spheroidal galaxies within the HAWC field-of-view, as well as their combined limit. These are the first limits on the annihilation cross section and decay lifetime using data collected with HAWC.

We present the development of a real-time flare monitor for the High Altitude Water Cherenkov (HAWC) observatory. The flare monitor has been fully operational since 2017 January and is designed to detect very high energy (VHE; $E\gtrsim100$ GeV) transient events from blazars on time scales lasting from 2 minutes to 10 hours in order to facilitate multiwavelength and multimessenger studies. These flares provide information for investigations into the mechanisms that power the blazars' relativistic jets and accelerate particles within them, and they may also serve as probes of the populations of particles and fields in intergalactic space. To date, the detection of blazar flares in the VHE range has relied primarily on pointed observations by imaging atmospheric Cherenkov telescopes. The recently completed HAWC observatory offers the opportunity to study VHE flares in survey mode, scanning 2/3 of the entire sky every day with a field of view of $\sim$1.8 steradians. In this work, we report on the sensitivity of the HAWC real-time flare monitor and demonstrate its capabilities via the detection of three high-confidence VHE events in the blazars Markarian 421 and Markarian 501.

We present a search of very high energy gamma-ray emission from the Northern $\textit{Fermi}$ Bubble region using data collected with the High Altitude Water Cherenkov (HAWC) gamma-ray observatory. The size of the data set is 290 days. No significant excess is observed in the Northern $\textit{Fermi}$ Bubble region, hence upper limits above $1\,\text{TeV}$ are calculated. The upper limits are between $3\times 10^{-7}\,\text{GeV}\, \text{cm}^{-2}\, \text{s}^{-1}\,\text{sr}^{-1}$ and $4\times 10^{-8}\,\text{GeV}\,\text{cm}^{-2}\,\text{s}^{-1}\,\text{sr}^{-1}$. The upper limits disfavor a proton injection spectrum that extends beyond $100\,\text{TeV}$ without being suppressed. They also disfavor a hadronic injection spectrum derived from neutrino measurements.

We present results from daily monitoring of gamma rays in the energy range $\sim0.5$ to $\sim100$ TeV with the first 17 months of data from the High Altitude Water Cherenkov (HAWC) Observatory. Its wide field of view of 2 steradians and duty cycle of $>95$% are unique features compared to other TeV observatories that allow us to observe every source that transits over HAWC for up to $\sim6$ hours each sidereal day. This regular sampling yields unprecedented light curves from unbiased measurements that are independent of seasons or weather conditions. For the Crab Nebula as a reference source we find no variability in the TeV band. Our main focus is the study of the TeV blazars Markarian (Mrk) 421 and Mrk 501. A spectral fit for Mrk 421 yields a power law index $\Gamma=2.21 \pm0.14_{\mathrm{stat}}\pm0.20_{\mathrm{sys}}$ and an exponential cut-off $E_0=5.4 \pm 1.1_{\mathrm{stat}}\pm 1.0_{\mathrm{sys}}$ TeV. For Mrk 501, we find an index $\Gamma=1.60\pm 0.30_{\mathrm{stat}} \pm 0.20_{\mathrm{sys}}$ and exponential cut-off $E_0=5.7\pm 1.6_{\mathrm{stat}} \pm 1.0_{\mathrm{sys}}$ TeV. The light curves for both sources show clear variability and a Bayesian analysis is applied to identify changes between flux states. The highest per-transit fluxes observed from Mrk 421 exceed the Crab Nebula flux by a factor of approximately five. For Mrk 501, several transits show fluxes in excess of three times the Crab Nebula flux. In a comparison to lower energy gamma-ray and X-ray monitoring data with comparable sampling we cannot identify clear counterparts for the most significant flaring features observed by HAWC.

The first limits on the prompt emission from the long gamma-ray burst (GRB) 130427A in the $>100\nobreakspace\rm{GeV}$ energy band are reported. GRB 130427A was the most powerful burst ever detected with a redshift $z\lesssim0.5$ and featured the longest lasting emission above $100\nobreakspace\rm{MeV}$. The energy spectrum extends at least up to $95\nobreakspace\rm{GeV}$, clearly in the range observable by the High Altitude Water Cherenkov (HAWC) Gamma-ray Observatory, a new extensive air shower detector currently under construction in central Mexico. The burst occurred under unfavourable observation conditions, low in the sky and when HAWC was running 10% of the final detector. Based on the observed light curve at MeV-GeV energies, eight different time periods have been searched for prompt and delayed emission from this GRB. In all cases, no statistically significant excess of counts has been found and upper limits have been placed. It is shown that a similar GRB close to zenith would be easily detected by the full HAWC detector, which will be completed soon. The detection rate of the full HAWC detector may be as high as one to two GRBs per year. A detection could provide important information regarding the high energy processes at work and the observation of a possible cut-off beyond the $\mathit{Fermi}$-LAT energy range could be the signature of gamma-ray absorption, either in the GRB or along the line of sight due to the extragalactic background light.

We observed the afterglow of GRB 130427A with the RATIR instrument on the 1.5-m Harold L. Johnson telescope of the Observatorio Astron\'omico Nacional on Sierra San Pedro M\'artir. Our homogenous $griZYJH$ photometry extends from the night of burst to three years later. We fit a model for the afterglow. There is a significant positive residual which matches the behavior of SN 1998bw in the $griZ$ filters; we suggest that this is a photometric signature of the supernova SN 2013cq associated with the GRB. The peak absolute magnitude of the supernova is $M_r=-18.43\pm0.11$.

We present the first catalog of TeV gamma-ray sources realized with the recently completed High Altitude Water Cherenkov Observatory (HAWC). It is the most sensitive wide field-of-view TeV telescope currently in operation, with a 1-year survey sensitivity of ~5-10% of the flux of the Crab Nebula. With an instantaneous field of view >1.5 sr and >90% duty cycle, it continuously surveys and monitors the sky for gamma ray energies between hundreds GeV and tens of TeV. HAWC is located in Mexico at a latitude of 19 degree North and was completed in March 2015. Here, we present the 2HWC catalog, which is the result of the first source search realized with the complete HAWC detector. Realized with 507 days of data and represents the most sensitive TeV survey to date for such a large fraction of the sky. A total of 39 sources were detected, with an expected contamination of 0.5 due to background fluctuation. Out of these sources, 16 are more than one degree away from any previously reported TeV source. The source list, including the position measurement, spectrum measurement, and uncertainties, is reported. Seven of the detected sources may be associated with pulsar wind nebulae, two with supernova remnants, two with blazars, and the remaining 23 have no firm identification yet.

The radio galaxy IC310 located in the Perseus Cluster is one of the brightest objects in the radio and X-ray bands, and one of the closest active galactic nuclei observed in very-high energies. In GeV - TeV $\gamma$-rays, IC310 was detected in low and high flux states by the MAGIC telescopes from October 2009 to February 2010. Taking into account that the spectral energy distribution (SED) up to a few GeV seems to exhibit a double-peak feature and that a single-zone synchrotron self-Compton (SSC) model can explain all of the multiwavelength emission except for the non-simultaneous MAGIC emission, we interpret, in this work, the multifrequency data set of the radio galaxy IC310 in the context of homogeneous hadronic and leptonic models. In the leptonic framework, we present a multi-zone SSC model with two electron populations to explain the whole SED whereas for the hadronic model, we propose that a single-zone SSC model describes the SED up to a few GeVs and neutral pion decay products resulting from p$\gamma$ interactions could describe the TeV - GeV $\gamma$-ray spectra. These interactions occur when Fermi-accelerated protons interact with the seed photons around the SSC peaks. We show that, in the leptonic model the minimum Lorentz factor of second electron population is exceedingly high $\gamma_e\sim10^5$ disfavoring this model, and in the hadronic model the required proton luminosity is not extremely high $\sim 10^{44}$ erg/s, provided that charge neutrality between the number of electrons and protons is given. Correlating the TeV $\gamma$-ray and neutrino spectra through photo-hadronic interactions, we find that the contribution of the emitting region of IC310 to the observed neutrino and ultra-high-energy cosmic ray fluxes are negligible.

Analysis of data from the Fermi Large Area Telescope (LAT) in the region of the supernova remnant G78.2+2.1 reveals an excess at the position of the TeV source VER J2019+407. The GeV source is extended with a hard spectrum and it is likely the counterpart of the TeV source. The spectrum of the LAT emission connects smoothly with that of VER J2019+407, and the overall GeV--TeV spectrum is best described by a broken power-law with indices 1.8 and 2.5 below and above a break energy of 71 GeV. Several scenarios are considered to explain the nature of this unidentified gamma ray source.

The Crab Nebula is the brightest TeV gamma-ray source in the sky and has been used for the past 25 years as a reference source in TeV astronomy, for calibration and verification of new TeV instruments. The High Altitude Water Cherenkov Observatory (HAWC), completed in early 2015, has been used to observe the Crab Nebula at high significance across nearly the full spectrum of energies to which HAWC is sensitive. HAWC is unique for its wide field-of-view, nearly 2 sr at any instant, and its high-energy reach, up to 100 TeV. HAWC's sensitivity improves with the gamma-ray energy. Above $\sim$1 TeV the sensitivity is driven by the best background rejection and angular resolution ever achieved for a wide-field ground array. We present a time-integrated analysis of the Crab using 507 live days of HAWC data from 2014 November to 2016 June. The spectrum of the Crab is fit to a function of the form $\phi(E) = \phi_0 (E/E_{0})^{-\alpha -\beta\cdot{\rm{ln}}(E/E_{0})}$. The data is well-fit with values of $\alpha=2.63\pm0.03$, $\beta=0.15\pm0.03$, and log$_{10}(\phi_0~{\rm{cm}^2}~{\rm{s}}~{\rm{TeV}})=-12.60\pm0.02$ when $E_{0}$ is fixed at 7 TeV and the fit applies between 1 and 37 TeV. Study of the systematic errors in this HAWC measurement is discussed and estimated to be $\pm$50\% in the photon flux between 1 and 37 TeV. Confirmation of the Crab flux serves to establish the HAWC instrument's sensitivity for surveys of the sky. The HAWC survey will exceed sensitivity of current-generation observatories and open a new view of 2/3 of the sky above 10 TeV.

The merger of a binary neutron star system is suggested as the central engine of short gamma-ray bursts (sGRBs). For the merger process, simulations predict that magnetic field is amplified beyond magnetar field strength by Kelvin-Helmholtz instabilities. With the Large Area Telescope (LAT), bursts have been detected that show a temporally extended component in coincidence with a short-lasting peak at the end of the prompt phase. We show that the presence of these LAT components in a sGRB could provide evidence of magnetic field amplification in the neutron star merger.

The bright, short and hard GRB 090510 was detected by all instruments aboard Fermi and Swift satellites. The multiwavelength observations of this burst presented similar features with the Fermi-LAT-detected gamma-ray bursts. In the framework of the external shock model of early afterglow, a leptonic scenario that evolves in a homogeneous medium is proposed to revisit GRB 090510 and explain the multiwavelength light curve observations presented in this burst. These observations are consistent with the evolution of a jet before and after the jet break. The long-lasting LAT, X-ray and optical fluxes are explained in the synchrotron emission from the adiabatic forward shock. Synchrotron self-Compton emission from the reverse shock is consistent with the bright LAT peak provided that progenitor environment is entrained with strong magnetic fields. It could provide compelling evidence of magnetic field amplification in the neutron star merger.

We present magnetohydrodynamical simulations of a strong accretion onto magnetized proto-neutron stars for the Kesteven 79 (Kes 79) scenario. The supernova remnant Kes 79, observed with the Chandra ACIS-I instrument during approximately 8.3 h, is located in the constellation Aquila at a distance of 7.1 kpc in the galactic plane. It is a galactic and a very young object with an estimate age of 6 kyr. The Chandra image has revealed, for the first time, a point-like source at the center of the remnant. The Kes 79 compact remnant belongs to a special class of objects, the so-called Central Compact Objects, which exhibits no evidence for a surrounding pulsar wind nebula. In this work we show that the submergence of the magnetic field during the hypercritical phase can explain such behavior for Kes 79 and others CCOs. The simulations of such regime were carried out with the adaptive-mesh-refinement code FLASH in two spatial dimensions, including radiative loss by neutrinos and an adequate equation of state for such regime. From the simulations, we estimate that the number of thermal neutrinos expected on the Hyper-Kamiokande Experiment is 733$\pm$364. In addition, we compute the flavor ratio on Earth for a progenitor model.

The closest radio galaxies; Centaurus A, M87 and NGC 1275, have been detected from radio wavelengths to TeV $\gamma$-rays, and also studied as high-energy neutrino and ultra-high-energy cosmic ray potential emitters. Their spectral energy distributions show a double-peak feature, which is explained by synchrotron self-Compton model. However, TeV $\gamma$-ray measured spectra could suggest that very-high-energy $\gamma$-rays might have a hadronic origin. We introduce a lepto-hadronic model to describe the broadband spectral energy distribution; from radio to sub GeV photons as synchrotron self-Compton emission and TeV $\gamma$-ray photons as neutral pion decay resulting from p$\gamma$ interactions occurring close to the core. These photo-hadronic interactions take place when Fermi-accelerated protons interact with the seed photons around synchrotron self-Compton peaks. Obtaining a good description of the TeV $\gamma$-ray fluxes, firstly, we compute neutrino fluxes and events expected in IceCube detector and secondly, we estimate ultra-high-energy cosmic ray fluxes and event rate expected in Telescope Array, Pierre Auger and HiRes observatories. Within this scenario we show that the expected high-energy neutrinos cannot explain the astrophysical flux observed by IceCube, and the connection with ultra-high-energy cosmic rays observed by Auger experiment around Centaurus A, might be possible only considering a heavy nuclei composition in the observed events.

Analysis of gamma-ray emission from the supernova remnant G78.2+2.1 ($\gamma$ Cygni) with 7.2 years of cumulative data from the Fermi-LAT telescope shows a distinct hard, bright and extended component to the north of the shell coincident with the known TeV source VER J2019+407. In the GeV-TeV energy range its spectrum is best described by a broken power-law with indices 1.8 below a break energy of 71 GeV and 2.5 above the break. A broadband spectral energy distribution is assembled and different scenarios for the origin of the gamma-rays are explored. Both hadronic and leptonic mechanisms are able to account for the GeV-TeV observations. In the leptonic framework, a superposition of inverse Compton and nonthermal bremsstrahlung emissions is needed whereas the hadronic scenario requires a cosmic ray population described by a broken power-law distribution with a relatively hard spectral index of $\sim1.8$ below a break particle energy of 0.45 TeV. In addition, the neutrino flux expected from cosmic ray interactions is calculated.

Removing outbursts from multiwavelength light curves of the blazar Mrk~421, we construct outburstless time series for this system. A model-independent power spectrum light curve analysis in the optical, hard X-ray and gamma-rays for this outburstless state and also the full light-curves, show clear evidence for a periodicity of ~ 310 days across all wavelengths studied. A subsequent full maximum likelihood analysis fitting an eclipse model confirms a periodicity of 310 \pm 1 days. The power spectrum of the signal in the outburstless state of the source does not follow a flicker noise behaviour and so, the system producing it is not self-organised. This and the fact that the periodicity is better defined in the outburstless state, strongly suggests that it is not produced by any internal physical processes associated to the central engine. The simplest physical mechanism to which this periodicity could be ascribed is a dynamical effect produced by an orbiting supermassive black hole companion eclipsing the central engine. Interestingly, the optimal eclipse model infers a brightness enhancement of (136.4 \pm 20 )%, suggesting an eclipse resulting in a gravitational lensing brightening. Consisting with this interpretation, the eclipse occurs for only ( 9.7 \pm 0.2)% of the orbital period.

One of the most powerful gamma-ray bursts, GRB 130427A was swiftly detected from GeV $\gamma$-rays to optical wavelengths. In the GeV band, the Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope observed the highest-energy photon ever recorded of 95 GeV, and a bright peak in the early phase followed by emission temporally extended for more than 20 hours. In the optical band, a bright flash with a magnitude of $7.03\pm 0.03$ in the time interval from 9.31 s to 19.31 s after the trigger was reported by RAPTOR in r-band. We study the origin of the GeV $\gamma$-ray emission, using the multiwavelength observation detected in X-ray and optical bands. The origin of the temporally extended LAT, X-ray and optical flux is naturally interpreted as synchrotron radiation and the 95-GeV photon and the integral flux upper limits placed by the HAWC observatory are consistent with synchrotron self-Compton from an adiabatic forward shock propagating into the stellar wind of its progenitor. The extreme LAT peak and the bright optical flash are explained through synchrotron self-Compton and synchrotron emission from the reverse shock, respectively, when the ejecta evolves in thick-shell regime and carries a significant magnetic field.

The IceCube Collaboration initially reported the detection of 37 extraterrestrial neutrinos in the TeV - PeV energy range. The reconstructed neutrino events were obtained during three consecutive years of data taking, from 2010 to 2013. Although these events have been discussed to have an extragalactic origin, they have not been correlated to any known source. Recently, the IceCube Collaboration reported a neutrino-induced muon event with energy of $2.6\pm0.3$ PeV which corresponds to the highest event ever detected. Neither the reconstructed direction of this event (J2000.0), detected on June 11 2014 at R.A.=110$^\circ$.34, Dec.=11$^\circ$.48 matches with any familiar source. Long gamma-ray bursts (lGRBs) are usually associated with the core collapse of massive stars leading relativistic-collimated jets inside stars with high-energy neutrino production. These neutrinos have been linked to the 37 events previously detected by IceCube experiment. In this work, we explore the conditions and values of parameters so that the highest neutrino recently detected could be generated by proton-photon and proton-hadron interactions at internal shocks inside lGRB progenitor star and then detected in IceCube experiment. Considering that internal shocks take place in a relativistic collimated jet, whose (half) opening angle is $\theta_0\sim$ 0.1, we found that lGRBs with total luminosity $L\lesssim 10^{48}$ erg/s and internal shocks on the surface of progenitors such as Wolf-Rayet (WR) and blue super giant (BSG) stars favor this multi-PeV neutrino production, although this neutrino could be associated to $L\sim 10^{50.5}$ ($\sim 10^{50}$) erg/s provided that the internal shocks occur at $\sim 10^9$ ($\sim 10^{10.2}$) cm for a WR (BSG).

List of proceedings from the HAWC Collaboration presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

A survey of the inner Galaxy region of Galactic longitude l in [+15, +50] degree and latitude b in [-4,+4] degree is performed using one-third of the High Altitude Water Cherenkov (HAWC) Observatory operated during its construction phase. To address the ambiguities arising from unresolved sources in the data, we use a maximum likelihood technique to identify point source candidates. Ten sources and candidate sources are identified in this analysis. Eight of these are associated with known TeV sources but not all have differential fluxes compatible with previous measurements. Three sources are detected with significances $>5\,\sigma$ after accounting for statistical trials, and are associated with known TeV sources.

The blazar Mrk 421 is one of the closest, brightest and fastest varying source in the extragalactic X-ray/TeV sky. In the last years, many multi-wavelength campaigns have been carried out to study the correlation between the very high energy (VHE) $\gamma$-ray and X-ray fluxes of this source and, although the activity in these two energy ranges seems to be correlated in many observations, no conclusive results have been achieved yet. In this work we present a robust study of the VHE $\gamma$-ray/X-ray correlation of Mrk 421 with data taken with different VHE experiments on different time scales and for different levels of activity of the source, with special focus on the low activity states. In particular, we discuss the robustness of the correlation at the lowest fluxes corresponding to the quiescent state of Mrk 421.

One of the most energetic gamma-ray bursts GRB 110731A, was observed from optical to GeV energy range by Fermi and Swift Observatories, and by the MOA and GROND optical telescopes. The multiwavelength observations over different epochs (from trigger time to more than 800 s) showed that the spectral energy distribution was better fitted by a wind afterglow model. We present a leptonic model based on an early afterglow that evolves in a stellar wind to describe the multiwavelength light curves observations. In particular, the origin of the LAT emission is explained through the superposition of synchrotron radiation from the forward shock and synchrotron self-Compton emission from the reverse shock. The bulk Lorentz factor required in this model is $\Gamma\simeq520$ and the result suggests that the ejecta must be magnetized.

The first one started in 2012 July 16 (MJD 56124) and the second one in 2013 April 9 (MJD 56391). The multiwavelength data analysis shows that the $\gamma$-ray flare observed in 2012 was not detected in the hard-X ray bands. This result is usually interpreted as an "orphan" flare. In 2013, the analysis of the multiwavelength light curves shows that there are two very bright states detected in the optical R-band. The first one in 2013 April 9 (R =11.74 $\pm$ 0.04) and the second one in May 12 (R =11.62 $\pm$ 0.04). Also, high activity states were detected in the soft and hard X-rays. A discrete correlation function analysis of this last flare shows a strong correlation between the GeV $\gamma$-rays and the optical/hard-X ray emission. These results are discussed in terms of the more adequate standard scenarios that could explain the multiwavelength variations displayed by this blazar.

The unusual GRB 110709B triggered Swift/BAT twice, with a time difference of $\sim 11$ minutes. Its light curve presented three noticeable peaks but only two were originally identified. In this work, we describe each peak as due to a different central-engine phase: the first one is the millisecond-protomagnetar stage, the second one is the BH-formation collapse phase and the last one is the Collapsar scenario with a Blandford-Znajek engine. Additionally, we analyze and explain the afterglow phase evoking the standard fireball model. Our model can successfully describe the timescales, fluxes and spectral indices observed for GRB 110709B.

GRB 110709B presented a peculiar three-peak lightcurve; this burst twice triggered the BAT detector onboard Swift. The two triggers were separated by $\sim 10$ minutes. In order to explain such an event, we unify into a single description the millisecond (ms) protomagnetar and the collapsar central-engine models. We find that such a scenario could produce GRBs with three peaks. One for the ms-protomagnetar stage, a second one for the BH-formation event and a third one for the collapsar phase. We show that the three peaks for GRB 110709B originate from different phases of the same collapsing object. We estimate the energies and timescales of the different episodes of this burst using our model and compare with previous results as well as with a reanalysis we perform on the data. We show that not only the light curve, but also the photon index evolution and the delay between the prompt emission and the afterglow of the second central-engine activity phase point towards a model like the one proposed here. We find that, with reasonable assumptions, our model correctly describes the activity in GRB 110709B. We further suggest careful study of future GRBs lightcurves which may help show the validity of our model. If our model is correct, this would be the first time that the formation of a BH from a core-collapse event is observed unimpededly.

Markarian 421 (Mrk 421) is one of the brightest, fastest and closest BL Lac object known. Its very high energy (VHE) spectrum has been successfully modeled with both leptonic and hadronic models and not conclusive results have been achieved yet about the origin of its VHE emission. Here we investigate the possibility that a fraction of the VHE flares of Mrk 421 are due to hadronic processes and calculate the expected neutrino flux associated. We introduce the obtained neutrino flux in a Monte Carlo simulation to see the expectation for a Km$^{3}$ Cherenkov neutrino telescope.

The nearby active galaxy NGC 1275, has widely been detected from radio to gamma rays. Its spectral energy distribution (SED) shows a double-peak feature, which is well explained by synchrotron self-Compton (SSC) model. However, recent TeV detections might suggest that very-high-energy $\gamma$-rays (E$\geq$100 GeV) may not have a leptonic origin. We test a lepto-hadronic model to describe the whole SED through SSC emission and neutral pion decay resulting from p$\gamma$ interactions. Also, we estimate the neutrino events expected in a Km$^3$ Cherenkov telescope.

Cassiopeia A the youngest supernova remnant known in the Milky Way is one of the brightest radio sources in the sky and a unique laboratory for supernova physics. Although its compact remnant was discovered in 1999 by the Chandra X-Ray Observatory, nowadays it is widely accepted that a neutron star lies in the center of this supernova remnant. In addition, new observations suggest that such neutron star with a low magnetic field and evidence of a carbon atmosphere could have suffered a hypercritical accretion phase seconds after the explosion. Considering this hypercritical accretion episode, we compute the neutrino cooling effect, the number of events and neutrino flavor ratios expected on Hyper-Kamiokande Experiment. The neutrino cooling effect (the emissivity and luminosity of neutrinos) is obtained through numerical simulations performed in a customized version of the FLASH code. Based on these simulations, we forecast that the number of events expected on the Hyper-Kamiokande Experiment is around 3195. Similarly, we estimate the neutrino flavor ratios to be detected considering the neutrino effective potential due to the thermal and magnetized plasma and thanks to the envelope of the star. It is worth noting that our estimates correspond to the only trustworthy method for verifying the hypercritical phase and although this episode took place 330 years ago, at present supernova remnants with these similarities might occur thus confirming our predictions for this phase.

TeV $\gamma$-ray detections in flaring states without activity in X-rays from blazars have attracted much attention due to the irregularity of these "orphan" flares. Although the synchrotron self-Compton model has been very successful in explaining the spectral energy distribution and spectral variability of these sources, it has not been able to describe these atypical flaring events. On the other hand, an electron-positron pair plasma at the base of the AGN jet was proposed as the mechanism of bulk acceleration of relativistic outflows. This plasma in quasi-themal equilibrium called Wein fireball emits radiation at MeV-peak energies serving as target of accelerated protons. In this work we describe the "orphan" TeV flares presented in blazars 1ES 1959+650 and Mrk421 assuming geometrical considerations in the jet and evoking the interactions of Fermi-accelerated protons and MeV-peak target photons coming from the Wein fireball. After describing successfully these "orphan" TeV flares, we correlate the TeV $\gamma$-ray, neutrino and UHECR fluxes through p$\gamma$ interactions and calculate the number of high-energy neutrinos and UHECRs expected in IceCube/AMANDA and TA experiment, respectively. In addition, thermal MeV neutrinos produced mainly through electron-positron annihilation at the Wein fireball will be able to propagate through it. By considering two- (solar, atmospheric and accelerator parameters) and three-neutrino mixing, we study the resonant oscillations and estimate the neutrino flavor ratios as well as the number of thermal neutrinos expected on Earth.

Long TeV $\gamma$-ray campaigns have been carried out to study the spectrum, variability and duty cycle of the BL Lac object Markarian 421. These campaigns have given some evidence of the presence of protons in the jet: i) Its spectral energy distribution which shows two main peaks; one at low energies ($\sim$ 1 keV) and the other at high energies (hundreds of GeV), has been described by using synchrotron proton blazar model. ii) The study of the variability at GeV $\gamma$-rays and X-rays has indicated no significant correlation. iii) TeV $\gamma$-ray detections without activity in X-rays, called "orphan flares" have been observed in this object.\\ Recently, The Telescope Array Collaboration reported the arrival of 72 ultra-high-energy cosmic rays with some of them possibly related to the direction of Markarian 421. The IceCube Collaboration reported the detection of 37 extraterrestrial neutrinos in the TeV - PeV energy range collected during three consecutive years. In particular, no neutrino track events were associated with this source. In this paper, we consider the proton-photon interactions to correlate the TeV $\gamma$-ray fluxes reported by long campaigns with the neutrino and ultra-high-energy cosmic ray observations around this blazar. Considering the results reported by The IceCube and Telescope Array Collaborations, we found that only from $\sim$ 25\% to 70\% of TeV fluxes described with a power law function with exponential cutoff can come from the proton-photon interactions.

High-energy neutrinos generated in collimated jets inside the progenitors of gamma-ray bursts (GRBs) have been related with the events detected by IceCube. These neutrinos, produced by hadronic interactions of Fermi-accelerated protons with thermal photons and hadrons in internal shocks, are the only signature when jet has not broken out or failed. Taking into account that the photon field is thermalized at keV energies and the standard assumption that the magnetic field maintains a steady value throughout the shock region (with a width of $10^{10} - 10^{11}$ cm in the observed frame), we study the effect of thermal and magnetized plasma generated in internal shocks on the neutrino oscillations. We calculate the neutrino effective potential generated by this plasma, the effects of the envelope of the star, and the vacuum on the path to Earth. By considering these three effects, the two (solar, atmospheric and accelerator parameters) and three neutrino mixing, we show that although GeV - TeV neutrinos can oscillate resonantly from one flavor to another, a nonsignificant deviation of the standard flavor ratio (1:1:1) could be expected on Earth.

One of the most energetic gamma-ray burst GRB 110731A was observed from optical to GeV energy range. Previous analysis on the prompt phase revealed similarities with the Large Area Telescope (LAT) bursts observed by Fermi: i) a delayed onset of the high-energy emission ($> 100$ MeV), ii) a short-lasting bright peak at later times and iii) a temporally extended component from this phase and lasting hundreds of seconds. Additionally to the prompt phase, multiwavelength observations over different epochs showed that the spectral energy distribution was better fitted by a wind afterglow model. We present a leptonic model based on an early afterglow that evolves in a stellar wind of its progenitor. We apply this model to interpret the temporally extended LAT emission and the brightest LAT peak exhibited by the prompt phase of GRB 110731A. Additionally, using the same set of parameters, we describe the multiwavelength afterglow observations. The origin of the temporally extended LAT, X-ray and optical flux is explained through synchrotron radiation from the forward shock and the brightest LAT peak is described evoking the synchrotron self-Compton emission from the reverse shock. The bulk Lorentz factor required in this model ($\Gamma\simeq520$) lies in the range of values demanded for most LAT-detected gamma-ray bursts. We show that the strength of the magnetic field in the reverse-shock region is $\sim$ 50 times stronger than in the forward-shock region. This result suggests that for GRB 110731A, the central engine is likely entrained with strong magnetic fields.

The High Altitude Water Cherenkov (HAWC) observatory is a wide field-of-view detector sensitive to gamma rays of 100 GeV to a few hundred TeV. Located in central Mexico at 19 degrees North latitude and 4100 m above sea level, HAWC will observe gamma rays and cosmic rays with an array of water Cherenkov detectors. The full HAWC array is scheduled to be operational in Spring 2015. In this paper, we study the HAWC sensitivity to the gamma-ray signatures of high-mass (multi- TeV) dark matter annihilation. The HAWC observatory will be sensitive to diverse searches for dark matter annihilation, including annihilation from extended dark matter sources, the diffuse gamma-ray emission from dark matter annihilation, and gamma-ray emission from non-luminous dark matter subhalos. Here we consider the HAWC sensitivity to a subset of these sources, including dwarf galaxies, the M31 galaxy, the Virgo cluster, and the Galactic center. We simulate the HAWC response to gamma rays from these sources in several well-motivated dark matter annihilation channels. If no gamma-ray excess is observed, we show the limits HAWC can place on the dark matter cross-section from these sources. In particular, in the case of dark matter annihilation into gauge bosons, HAWC will be able to detect a narrow range of dark matter masses to cross-sections below thermal. HAWC should also be sensitive to non-thermal cross-sections for masses up to nearly 1000 TeV. The constraints placed by HAWC on the dark matter cross-section from known sources should be competitive with current limits in the mass range where HAWC has similar sensitivity. HAWC can additionally explore higher dark matter masses than are currently constrained.

Fanaroff-Riley I radiogalaxies have been observed in TeV gamma-rays during the last decades. The origin of the emission processes related with this energy band is still under debate. Here we consider the case of the two closest Fanaroff-Riley I objects: Centaurus A and M87. Their entire broadband spectral energy distributions and variability fluxes show evidences that leptonic models are not sufficient to explain their fluxes above 100 GeV. Indeed, both objects have been imaged by LAT instrument aboard of Fermi telescope with measured spectra well connected with one-zone leptonic models. However, to explain the TeV spectra obtained with campaigns by H.E.S.S., for Centaurus A, and by VERITAS, MAGIC and H.E.S.S. for M87, different emission processes must be introduced. In this work we evoke hadronic scenarios to describe the TeV gamma-ray fluxes observed and to obtain the expected neutrino counterparts for each considered TeV campaign. With the obtained neutrino spectra we calculate, through Monte Carlo simulations, the expected neutrino event rate in a hypothetical Km$^{3}$ neutrino telescope and we compare the results with what has been observed by IceCube experiment up to now.

Very high energy gamma-ray emission of Fanaroff-Riley I objects is not univocally explained by a single emission model. Leptonic models with one and multi-zone emission regions, occurring in the jet of these objects, are usually used to describe the broadband spectral energy distribution. A correlation between the X-ray and TeV emission is naturally expected within leptonic models whereas a lack of correlation between these two observables represents a challenge and favors the hadronic scenarios. This is the case of M87 as we show here by analyzing its TeV and X-ray emission recorded in the last decade. Furthermore, we point out that the spectra obtained by MAGIC, H.E.S.S. and VERITAS telescopes cannot be described with the same leptonic model introduced by the Fermi-LAT collaboration. We introduce hadronic scenarios to explain the TeV gamma-ray fluxes of this radiogalaxy as products of Fermi-accelarated protons interacting with seed photons in the jet or thermal particles in the giant lobes. By fitting this part of spectral energy distribution as pion decay products, we obtain the expected neutrino counterpart and the luminosity of accelerating protons in the jet and/or lobes. With the expected neutrino fluxes we investigate, through Monte Carlo simulations, the possibility to see the signal from M87 with a Km$^{3}$ neutrino telescope, and compare the results with what has been seen by IceCube experiment up to now. Finally we constrain the features of giant lobes through the observations performed at ultra high energies by TA experiment.

The High-Altitude Water Cherenkov (HAWC) Observatory is sensitive to gamma rays and charged cosmic rays at TeV energies. The detector is still under construction, but data acquisition with the partially deployed detector started in 2013. An analysis of the cosmic-ray arrival direction distribution based on $4.9\times 10^{10}$ events recorded between June 2013 and February 2014 shows anisotropy at the $10^{-4}$ level on angular scales of about $10^\circ$. The HAWC cosmic-ray sky map exhibits three regions of significantly enhanced cosmic-ray flux; two of these regions were first reported by the Milagro experiment. A third region coincides with an excess recently reported by the ARGO-YBJ experiment. An angular power spectrum analysis of the sky shows that all terms up to $\ell=15$ contribute significantly to the excesses.

Primordial Black Holes (PBHs) are gravitationally collapsed objects that may have been created by density fluctuations in the early universe and could have arbitrarily small masses down to the Planck scale. Hawking showed that due to quantum effects, a black hole has a temperature inversely proportional to its mass and will emit all species of fundamental particles thermally. PBHs with initial masses of ~5.0 x 10^14 g should be expiring in the present epoch with bursts of high-energy particles, including gamma radiation in the GeV - TeV energy range. The Milagro high energy observatory, which operated from 2000 to 2008, is sensitive to the high end of the PBH evaporation gamma-ray spectrum. Due to its large field-of-view, more than 90% duty cycle and sensitivity up to 100 TeV gamma rays, the Milagro observatory is well suited to perform a search for PBH bursts. Based on a search on the Milagro data, we report new PBH burst rate density upper limits over a range of PBH observation times. In addition, we report the sensitivity of the Milagro successor, the High Altitude Water Cherenkov (HAWC) observatory, to PBH evaporation events.

VAMOS was a prototype detector built in 2011 at an altitude of 4100m a.s.l. in the state of Puebla, Mexico. The aim of VAMOS was to finalize the design, construction techniques and data acquisition system of the HAWC observatory. HAWC is an air-shower array currently under construction at the same site of VAMOS with the purpose to study the TeV sky. The VAMOS setup included six water Cherenkov detectors and two different data acquisition systems. It was in operation between October 2011 and May 2012 with an average live time of 30%. Besides the scientific verification purposes, the eight months of data were used to obtain the results presented in this paper: the detector response to the Forbush decrease of March 2012, and the analysis of possible emission, at energies above 30 GeV, for long gamma-ray bursts GRB111016B and GRB120328B.

The High Altitude Water Cherenkov (HAWC) observatory is a TeV gamma-ray and cosmic-ray detector currently under construction at an altitude of 4100 m close to volcano Sierra Negra in the state of Puebla, Mexico. The HAWC observatory is an extensive air-shower array comprised of 300 optically-isolated water Cherenkov detectors (WCDs). Each WCD contains $\sim$200,000 liters of filtered water and four upward-facing photomultiplier tubes. In Fall 2014, when the HAWC observatory will reach an area of 22,000 m$^{2}$, the sensitivity will be 15 times higher than its predecessor Milagro. Since September 2012, more than 30 WCDs have been instrumented and taking data. This first commissioning phase has been crucial for the verification of the data acquisition and event reconstruction algorithms. Moreover, with the increasing number of instrumented WCDs, it is important to verify the data taken with different configuration geometries. In this work we present a comparison between Monte Carlo simulation and data recorded by the experiment during 24 hours of live time between 14 and 15 April of 2013 when 29 WCDs were active.

The neutrino signal from SN1987A confirmed the core-collapse scenario and the possible formation of a neutron star. Although this compact object has eluded all observations, theoretical and numerical developments have allowed a glimpse of the fate of it. In particular, a hypercritical accretion model has been proposed to forecast the accretion of $\sim 0.15 \:\mathrm{M_{\odot}}$ in two hours and the subsequent submergence of the magnetic field in the newborn neutron star. In this paper, we revisit the Chevalier's model in a numerical framework, focusing on the neutrino cooling effect on the supernova fall-back dynamics. For that, using a customized version of the FLASH code, we carry out numerical simulations of the accretion of matter onto the newborn neutron star in order to estimate the size of the neutrino-sphere, the emissivity and luminosity of neutrinos. As a signature of this phase, we estimate the neutrinos expected on SK neutrino experiment and their flavor ratios. This is academically important because, although currently it was very difficult to detect 1.46 thermal neutrinos and their oscillations, these fingerprints are the only viable and reliable way to confirm the hypercritical phase. Perhaps new techniques for detecting neutrino oscillations arise in the near future allowing us to confirm our estimates.

Neutrons play an important role in the dynamics of gamma-ray bursts. The presence of neutrons in the baryon-loaded fireball is expected. If the neutrons abundance is comparable to that of protons, important features may be observed such as quasi-thermal multi-GeV neutrinos in coincidence with a subphotospheric $\gamma$-ray emission, nucleosynthesis at later times and rebrightening of the afterglow emission. Additionally, thermal MeV neutrinos are created by electron-positron annihilation, electron (positron) capture on protons (neutrons) and nucleonic bremsstrahlung. Although MeV neutrinos are difficult to detect, quasi-thermal GeV neutrinos are expected in km$^3$ detectors and/or DeepCore+IceCube. In this paper, we show that neutrino oscillations have outstanding implications for the dynamics of the fireball evolution and also that they can be detected through their flavor ratio on Earth. For that we derive the resonance and charged-neutrality conditions as well as the neutrino self-energy and effective potential up to order $m_W^{-4}$ at strong, moderate and weak magnetic field approximation to constrain the dynamics of the fireball. We found important implications: i) resonant oscillations are suppressed for high baryon densities as well as neutrons abundance larger than that of protons and ii) the effect of magnetic field is to decrease the proton-to-neutron ratio aside from the number of multi-GeV neutrinos expected in DeepCore detector. Also we estimate the GeV neutrino flavor ratios along the jet and on Earth.

A high-energy component in the radio galaxy Centaurus A was reported after analyzing four years of Fermi data. The spectrum of this component is described by means of a broken power law with a break energy of 4 GeV and, below and above spectral indices of $\alpha_1$=2.74$\pm$0.03 and $\alpha_2$=2.09$\pm$0.20, respectively. Also a faint $\gamma$-ray flux at TeV energies was detected by H.E.S.S.. In this paper we show that the spectrum at GeV-TeV energies is described through synchrotron self-Compton emission up to a few GeV ($\sim$ 4 GeV) and $\pi^0$ decay products up to TeV energies, although the emission of synchrotron radiation by muons could contribute to the spectrum at GeV energies, if they are rapidly accelerated. Muons and $\pi^0$s are generated in the interactions of accelerated protons with two populations of seed photons which were reported by Compton Gamma-Ray Observatory: one population at intermediate state emission with energy peak of 0.15 MeV and another at low state emission with energy peak of 0.59 MeV. In addition, we show that the reported observations of ultra-high-energy cosmic rays and non high-energy neutrino detection around Centaurus A can be explained through these interactions, assuming that proton spectrum is extended up to ultra-high-energies.

Collapsars are the likely progenitors of Long Gamma-Ray Burst (lGRBs). lGRBs have been observed to last for thousands to tens of thousands of seconds, thus making unlikely the neutrino-driven engine as the main mechanism for driving the jets. In this context, the Blandford-Znajek mechanism seems likely to explain the production of rotational-axis directed jets without the need for large accretion rates. These engines, require magnetic fields between $10^{12}$ G $< B < 10^{15}$ G threading the innerdisk, Kerr-BH region to exist. We derive the neutrino self-energy and the effective potential up to O(1/$M_W^4$) in a weakly and highly magnetized GRB fireball flow which is made up of electrons, protons, neutrons and their anti-particles. We consider neutrino energies of 1-100 MeV which are produced during stellar collapse, merger events or in the fireball itself by electron-positron annihilation, inverse beta decay and nucleonic bremsstrahlung processes. Many of these neutrinos propagate through the fireball and may oscillate resonantly. Using two-neutrino mixing we study the possibility of these oscillations.

Strange-quark matter (SQM) is a likely candidate of the ground state of nuclear matter. Along with many other equations of state (EoSs), SQM seemed to be severely constrained by the recent discoveries of the 1.97 $M_\odot$ PSR J1614-2230 and the 2.01 $M_\odot$ PSR J0348+0432. However with new, $O(\alpha_c^2)$, perturbative calculations, SQM seems to be able to accommodate masses as large as $\sim 2.75 M_\odot$. The literature of SQM stars or strange stars includes estimates of internal magnetic fields as large as $10^{20}$ G, which are unlikely to be formed as they would require $\sim 10^{57}$ erg to be produced. Nonetheless, if strange stars may hold magnetar-strength fields ($10^{15}$ G), their internal fields are likely to reach magnetic fields as large as $10^{17}$ G. We consider neutrinos with energies of some MeV and oscillation parameters from solar, atmospheric and accelerator experiments. We study the possibility of resonant oscillation of neutrinos in strange stars.

The Milagro experiment was a TeV gamma-ray observatory designed to continuously monitor the overhead sky in the 0.1-100 TeV energy range. It operated from 2000 and 2008 and was characterized by a large field of view ($\sim$ 2 sr) and a high duty cycle ($\geq$ 90$\%$). Here we report on the long-term monitoring of the blazar Mrk 421 with Milagro over the period from September 21, 2005 to March 15, 2008. We present a study of the TeV variability of the source and provide upper limits for the measured flux for different time scales, ranging from one week up to one year.

The blazar Markarian 421 (Mrk 421) is one of the brightest sources in the extragalactic X-ray/TeV sky. It is also one of the fastest varying TeV $\gamma$-ray sources, showing flaring activity on time scales as short as tens of minutes. To know the level of activity of this source, Tluczykont et al. (2007) calculated the fraction of time spent by Mrk 421 in flaring states with fluxes above 1 Crab at TeV energies (i.e., TeV - duty cycle). Here we present an alternative approach to calculate the TeV duty cycle of Mrk 421 taking advantage of the continuous monitoring of the source by the Milagro observatory. Milagro was a water Cherenkov detector sensitive at energies between 100 GeV and 100 TeV. We present our estimation of the TeV - duty cycle and study its robustness.

Markarian 421 (Mrk 421) is one the brightest and closest (z=0.031) blazars known (de Vaucouleurs et al 1991). It is also one of the fastest varying TeV $\gamma$-ray sources, with a flaring activity on time scales as short as tens of minutes. The activity of Mrk 421 at different frequencies may reflect the radiation mechanisms involved. Tluczykont et al. (2007) estimated the TeV activity of Mrk 421 through calculating the fraction of time spent in flaring states at TeV energies (TeV duty cycle) by using data from several imaging atmospheric Cherenkov telescopes (IACTs). Since IACT observations are biased towards high flux states they overestimated the TeV duty cycle of Mrk 421. Here we propose an alternative approach to calculate the TeV duty cycle of Mrk 421 that takes advantage of the continuous monitoring of the source by the Milagro experiment, a water Cherenkov detector sensitive to primary $\gamma$-rays between 100 GeV and 100 TeV. We present our estimation of the TeV - duty cycle and study its robustness.

TeV flaring activity with time scales as short as tens of minutes and an orphan TeV flare have been observed from the blazar Markarian 421 (Mrk 421). The TeV emission from Mrk 421 is believed to be produced by leptonic synchrotron self-Compton (SSC) emission. In this scenario, correlations between the X-ray and the TeV fluxes are expected, TeV orphan flares are hardly explained and the activity (measured as duty cycle) of the source at TeV energies is expected to be equal or less than that observed in X-rays if only SSC is considered. To estimate the TeV duty cycle of Mrk 421 and to establish limits on its variability at different time scales, we continuously observed Mrk 421 with the Milagro observatory. Mrk 421 was detected by Milagro with a statistical significance of 7.1 standard deviations between 2005 September 21 and 2008 March 15. The observed spectrum is consistent with previous observations by VERITAS. We estimate the duty cycle of Mrk 421 for energies above 1 TeV for different hypothesis of the baseline flux and for different flare selections and we compare our results with the X-ray duty cycle estimated by Resconi et al. 2009. The robustness of the results is discussed.

Population III stars are believed to be rapidly-rotating sources with a mass range of hundreds to thousand of solar masses. Masses larger than 260 M$_\odot$ are expected to collapse resulting in central rotating Kerr black holes with large rotation parameter a. Due to particle-antiparticle asymmetry is small, we use the neutrino effective potential up to order $1/M_W^4$ in a magnetized plasma at the base of the ejecta to constrain the rotation parameter by means of neutrino oscillations. Additionally, we investigate the implications in the magnetic field, temperature and electron asymmetry.

The spectral energy distribution of giant lobes shows one main peak detected by the Wilkinson Microwave Anisotropy Probe at low energy of $10^{-5}$ eV and a faint $\gamma$-ray flux imaged by Fermi Large Area Telescope at energy $\geq$ 100 MeV. On the other hand, Pierre Auger Observatory associated some ultra-high-energy cosmic rays with the direction of Centaurus A and IceCube reported 28 neutrino-induced events in a TeV - PeV energy range, although none of them related with this direction. In this work we describe the spectra for each of the lobes, the main peak with synchrotron radiation, and the high-energy emission with $pp$ interactions. Obtaining a good description of the main peak, we deduce the magnetic fields, electron densities and the age of the lobes. Describing successfully the $\gamma$-ray emission by pp interactions and considering as targets those thermal particles in the lobes with density in the range $10^{-10}$ to $10^{-4}\, {\rm cm}^{-3}$, we calculate the number of ultra-high-energy cosmic rays. Although $\gamma$-spectrum is well described with any density in the range, only when $10^{-4}\, {\rm cm}^{-3}$ is considered, the expected number of events is very similar to that observed by Pierre Auger Observatory, otherwise we obtain an excessive luminosity. In addition, correlating the $\gamma$-ray and neutrino fluxes through pp interactions we calculate the number of high-energy neutrinos expected in IceCube. Our analysis indicates that neutrinos above 1 TeV cannot be produced in the lobes of Centaurus A, which is consistent with the results recently published by IceCube Collaboration.

Long gamma-ray bursts have been widely associated with collapsing massive stars in the framework of collapsar model. High-energy neutrinos and photons can be produced in the internal shocks of middle relativistic jets from core-collapse supernova. Although photons can hardly escape, high-energy neutrinos could be the only signature when the jets are hidden. We show that using suitable parameters, high-energy neutrinos in GeV - PeV range can be produced in the hidden jet inside the collapsar, thus demonstrating that these objects are candidates to produce neutrinos with energies between 1 - 10 PeV which were observed with IceCube. On the other hand, due to matter effects, high-energy neutrinos may oscillate resonantly from one flavor to another before leaving the star. Using two (solar, atmospheric and accelerator parameters) and three neutrino mixing, we study the possibility of resonant oscillation for these neutrinos created in internal shocks. Also we compute the probabilities of neutrino oscillations in the matter at different distances along the jet (before leaving the star) and after in vacuum, on their path to Earth. Finally, neutrino flavor ratios on Earth are estimated.

The High-Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory is designed to record air showers produced by cosmic rays and gamma rays between 100 GeV and 100 TeV. Because of its large field of view and high livetime, HAWC is well-suited to measure gamma rays from extended sources, diffuse emission, and transient sources. We describe the sensitivity of HAWC to emission from the extended Cygnus region as well as other types of galactic diffuse emission; searches for flares from gamma-ray bursts and active galactic nuclei; and the first measurement of the Crab Nebula with HAWC-30.

We describe measurements of GeV and TeV cosmic rays with the High-Altitude Water Cherenkov Gamma-Ray Observatory, or HAWC. The measurements include the observation of the shadow of the moon; the observation of small-scale and large-scale angular clustering of the TeV cosmic rays; the prospects for measurement of transient solar events with HAWC; and the observation of Forbush decreases with the HAWC engineering array and HAWC-30.

The High-Altitude Water Cherenkov Gamma Ray Observatory (HAWC) is under construction 4100 meters above sea level at Sierra Negra, Mexico. We describe the design and cabling of the detector, the characterization of the photomultipliers, and the timing calibration system. We also outline a next-generation detector based on the water Cherenkov technique.

The High-Altitude Water Cherenkov Gamma Ray Observatory (HAWC) is designed to perform a synoptic survey of the TeV sky. The high energy coverage of the experiment will enable studies of fundamental physics beyond the Standard Model, and the large field of view of the detector will enable detailed studies of cosmologically significant backgrounds and magnetic fields. We describe the sensitivity of the full HAWC array to these phenomena in five contributions shown at the 33rd International Cosmic Ray Conference in Rio de Janeiro, Brazil (July 2013).

The High Altitude Water Cherenkov (HAWC) observatory is an array of large water Cherenkov detectors sensitive to gamma rays and hadronic cosmic rays in the energy band between 100 GeV and 100 TeV. The observatory will be used to measure high-energy protons and cosmic rays via detection of the energetic secondary particles reaching the ground when one of these particles interacts in the atmosphere above the detector. HAWC is under construction at a site 4100 meters above sea level on the northern slope of the volcano Sierra Negra, which is located in central Mexico at 19 degrees N latitude. It is scheduled for completion in 2014. In this paper we estimate the sensitivity of the HAWC instrument to point-like and extended sources of gamma rays. The source fluxes are modeled using both unbroken power laws and power laws with exponential cutoffs. HAWC, in one year, is sensitive to point sources with integral power-law spectra as low as 5x10^-13 cm^-2 sec^-1 above 2 TeV (approximately 50 mCrab) over 5 sr of the sky. This is a conservative estimate based on simple event parameters and is expected to improve as the data analysis techniques are refined. We discuss known TeV sources and the scientific contributions that HAWC can make to our understanding of particle acceleration in these sources.

We present a leptonic model on the external shock context to describe the high-energy emission of GRB 940217, GRB 941017 and GRB 970217A. We argue that the emission consists of two components, one with a similar duration of the burst, and a second, longer-lasting GeV phase lasting hundred of seconds after the prompt phase. Both components can be described as synchrotron self-Compton emission from a reverse and forward shock respectively. For the reverse shock, we analyze the synchrotron self-Compton in the thick-shell case. The calculated fluxes and break energies are all consistent with the observed values.

Recent detections of Fanaroff-Riley Class I AGNs by HESS, MAGIC, and VERITAS suggest that very-high-energy gamma-rays (VHE, E > 100 GeV) may not have a leptonic origin. We present a hadronic model to describe the TeV photons as the neutral pion decay resulting from pgamma and pp interactions. For the pgamma interaction, we assume that the target photons are produced by leptonic processes and apparent at the second spectral peak. For the pp interaction we consider as targets the thermal particle densities in the lobes. We show that this model can describe the TeV spectra of the radio galaxies NCG 1275, M87 and Cen A

We present a leptonic model on the external shock framework to describe the long- and short- lasting GeV component of some GRBs. This model was already applied successfully to GRB 090926A, and we extend it to describe the high-energy emission of GRB 090902B and GRB 090510. We argue that the high-energy emission consists of two components, one at MeV energies with a duration of a few seconds during the prompt phase, and a second GeV component lasting hundred of seconds after the prompt phase. The short high-energy component can be described as SSC emission from a reverse shock and the longer component arises from SSC emission of the forward shock. The main assumption of our model is that the jet is magnetized and evolves in the thick-shell case. The calculated fluxes and break energies are all consistent with the observed values.

SSC emission from a reverse shock has been suggested as the origin for the high energy component lasting 2 s in the prompt phase of GRB98080923 (Fraija et al. 2012). The model describes spectral indices, fluxes and the duration of the high-energy component as well as a long keV tail present in the prompt phase of GRB980923. Here, we present an extension of this model to describe the high-energy emission of GRB090926A. We argue that the emission consist of two components, one with a duration less than 1s during the prompt phase, and a second, longer-lasting GeV phase lasting hundred of seconds after the prompt phase. The short high-energy phase can be described as SSC emission from a reverse shock similar to that observed in GRB980923, while the longer component arises from the forward shock. The main assumption is that the jet is magnetized and evolves in the thick-shell case, and the calculated fluxes and break energies are all consistent with the observed values. A comparison between the resulting parameters obtained for GRB980923 and GRB090926A suggests differences in burst tails that could be attributable to the circumburst medium, and this could account for previous analyses reported in the literature for other bursts. We find that the density of the surrounding medium inferred from the observed values associated to the forward shock agrees with standard values for host galaxies such as the one associated to GRB090926A.

The Pierre Auger Observatory has associated a few ultra high energy cosmic rays with the direction of Centaurus A. This source has been deeply studied in radio, infrared, X-ray and $\gamma$-rays (MeV-TeV) because it is the nearest radio-loud active galactic nuclei. Its spectral energy distribution or spectrum shows two main peaks, the low energy peak, at an energy of $10^{-2}$ eV, and the high energy peak, at about 150 keV. There is also a faint very high energy (E $\geq$ 100 GeV) $\gamma$-ray emission fully detected by the High Energy Stereoscopic System experiment. In this work we describe the entire spectrum, the two main peaks with a Synchrotron/Self-Synchrotron Compton model and, the Very High Energy emission with a hadronic model. We consider p$\gamma$ and $pp$ interactions. For the p$\gamma$ interaction, we assume that the target photons are those produced at 150 keV in the leptonic processes. On the other hand, for the pp interaction we consider as targets the thermal particle densities in the lobes. Requiring a satisfactory description of the spectra at very high energies with p$\gamma$ interaction we obtain an excessive luminosity in ultra high energy cosmic rays (even exceeding the Eddington luminosity). However, when considering pp interaction to describe the $\gamma$-spectrum, the obtained number of ultra high energy cosmic rays are in agreement with Pierre Auger observations. Moreover, we calculate the possible neutrino signal from pp interactions on a Km$^3 $ neutrino telescope using Monte Carlo simulations.

Centaurus A (Cen A) is the nearest radio-loud AGN and is detected from radio to very high energy gamma-rays. Its nuclear spectral energy distribution (SED) shows a double-peak feature, which is well explained by the leptonic synchrotron + synchrotron self-Compton model. This model however cannot account for the observed high energy photons in the TeV range, which display a distinct component. Here we show that ~ TeV photons can be well interpreted as the neutral pion decay products from p-gamma interactions of Fermi accelerated high energy protons in the jet with the seed photons around the second SED peak at ~170 keV. Extrapolating the inferred proton spectrum to high energies, we find that this same model is consistent with the detection of 2 ultra-high-energy cosmic ray events detected by Pierre Auger Observatory from the direction of Cen A. We also estimate the GeV neutrino flux from the same process, and find that it is too faint to be detected by current high-energy neutrino detectors.

GRB 980923 was one of the brightest bursts observed by the Burst and Transient Source Experiment (BATSE). Previous studies have detected two distinct components in addition to the main prompt episode, which is well described by a Band function. The first of these is a tail with a duration of approx. 400s, while the second is a high-energy component lasting approx. 2s. We summarize the observations, and argue for a unified model in which the tail can be understood as the early gamma-ray afterglow from forward shock synchrotron emission, while the high-energy component arises from synchrotron self-Comtpon (SSC) from the reverse shock. Consistency between the main assumption of thick shell emission and agreement between the observed and computed values for fluxes, break energies, starting times and spectral indices leads to a requirement that the ejecta be highly magnetized.

We present the sensitivity of HAWC to Gamma Ray Bursts (GRBs). HAWC is a very high-energy gamma-ray observatory currently under construction in Mexico at an altitude of 4100 m. It will observe atmospheric air showers via the water Cherenkov method. HAWC will consist of 300 large water tanks instrumented with 4 photomultipliers each. HAWC has two data acquisition (DAQ) systems. The main DAQ system reads out coincident signals in the tanks and reconstructs the direction and energy of individual atmospheric showers. The scaler DAQ counts the hits in each photomultiplier tube (PMT) in the detector and searches for a statistical excess over the noise of all PMTs. We show that HAWC has a realistic opportunity to observe the high-energy power law components of GRBs that extend at least up to 30 GeV, as it has been observed by Fermi LAT. The two DAQ systems have an energy threshold that is low enough to observe events similar to GRB 090510 and GRB 090902b with the characteristics observed by Fermi LAT. HAWC will provide information about the high-energy spectra of GRBs which in turn could help to understanding about e-pair attenuation in GRB jets, extragalactic background light absorption, as well as establishing the highest energy to which GRBs accelerate particles.

GRB 980923 was one of the brightest bursts observed by the Burst and Transient Source Experiment (BATSE). Previous studies have detected two distinct components in addition to the main prompt episode, which is well described by a Band function. The first of these is a tail with a duration of 400s, while the second is a high-energy component lasting 2s. After summarizing the observations, we present a model for this event and conclude that the tail can be understood as the early gamma-ray afterglow from forward shock synchrotron emission, while the high-energy component is described by the SSC emission from the reverse shock. The main assumption is that of a thick-shell case from highly magnetized ejecta. The calculated fluxes, break energies, starting times and spectral index are all consistent with the observed values.

Centaurus A is the nearest radio-loud AGN and is detected from radio to very high energy gamma-rays. Its nuclear spectral energy distribution shows two peaks, one in the far-infrared band and another at about 150 keV. By assuming the second peak is due to the electron synchrotron emission and the power index for the differential spectrum of the very high energy cosmic ray proton to be 2.7 we show that only pp interaction is responsible for the observed GeV-TeV emission from Centaurus A. We also found that indeed many very high energy cosmic ray protons from Centaurus A can arrive on Earth thus supporting the recent observation of two events by Pierre Auger Observatory.

Neutrinos of energy about 5-20 MeV are produced due to the stellar collapse or merger events that trigger the Gamma-Ray Burst. Also low energy MeV neutrinos are produced within the fireball due to electron positron annihilation and nucleonic bremsstrahlung. Many of these neutrinos will propagate through the dense and relativistic magnetized plasma of the fireball. We have studied the possibility of resonant oscillation of $\nu_e\leftrightarrow \nu_{\mu,\tau}$ by taking into account the neutrino oscillation parameters from SNO, SuperKamiokande and Liquid Scintillator Detector. Using the resonance condition we have calculated the resonance length for these neutrinos and also the fireball observables like lepton asymmetry and the baryon load are estimated based on the assumed fireball radius of 100 Km.