Short gamma-ray bursts within 200 Mpc

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

All-particle cosmic ray energy spectrum measured by the HAWC experiment from 10 to 500 TeV

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

Very high energy particle acceleration powered by the jets of the microquasar SS 433

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.Comment: Preprint version of Nature paper. Contacts: S. BenZvi, B. Dingus, K. Fang, C.D. Rho , H. Zhang, H. Zho

A Reverse Shock and Unusual Radio Properties in GRB 160625B

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

Measurement of the Crab Nebula Spectrum Past 100 TeV with HAWC

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.Comment: published in Ap

Constraining the pˉ/p\bar{p}/p Ratio in TeV Cosmic Rays with Observations of the Moon Shadow by HAWC

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.Comment: 10 pages, 5 figures. Accepted by Physical Review

GRB Fermi-LAT Afterglows: Explaining Flares, Breaks, and Energetic Photons

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 (≈10² 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 an 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_B ≃ 30) and it decelerates in a uniform-density medium with a very low density (n = 4.554_(-1.121)^(+1.128) x 10⁻⁴ cm⁻³