Searches for optical transients are usually performed with a cadence of days to weeks, optimised for supernova discovery. The optical fast transient sky is still largely unexplored, with only a few surveys to date having placed meaningful constraints on the detection of extragalactic transients evolving at sub-hour timescales. Here, we present the results of deep searches for dim, minute-timescale extragalactic fast transients using the Dark Energy Camera, a core facility of our all-wavelength and all-messenger Deeper, Wider, Faster programme. We used continuous 20s exposures to systematically probe timescales down to 1.17 minutes at magnitude limits $g > 23$ (AB), detecting hundreds of transient and variable sources. Nine candidates passed our strict criteria on duration and non-stellarity, all of which could be classified as flare stars based on deep multi-band imaging. Searches for fast radio burst and gamma-ray counterparts during simultaneous multi-facility observations yielded no counterparts to the optical transients. Also, no long-term variability was detected with pre-imaging and follow-up observations using the SkyMapper optical telescope. We place upper limits for minute-timescale fast optical transient rates for a range of depths and timescales. Finally, we demonstrate that optical $g$-band light curve behaviour alone cannot discriminate between confirmed extragalactic fast transients such as prompt GRB flashes and Galactic stellar flares.

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. A narrow feature at ~3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux ($>0.3$ keV) measured during the entire Swift mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3$\sigma$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$\sigma$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$\sigma$. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: a) a pileup in the EED due to stochastic acceleration; b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and c) an emission from an IC pair cascade.

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

The accretion-powered X-ray pulsar GX 301-2 was observed with the balloon-borne X-Calibur hard X-ray polarimeter during late December 2018, with contiguous observations by the NICER X-ray telescope, the Swift X-ray Telescope and Burst Alert Telescope, and the Fermi Gamma-ray Burst Monitor spanning several months. The observations detected the pulsar in a rare apastron flaring state coinciding with a significant spin-up of the pulsar discovered with the Fermi GBM. The X-Calibur, NICER, and Swift observations reveal a pulse profile strongly dominated by one main peak, and the NICER and Swift data show strong variation of the profile from pulse to pulse. The X-Calibur observations constrain for the first time the linear polarization of the 15-35 keV emission from a highly magnetized accreting neutron star, indicating a polarization degree of (27+38-27)% (90% confidence limit) averaged over all pulse phases. We discuss the spin-up and the X-ray spectral and polarimetric results in the context of theoretical predictions. We conclude with a discussion of the scientific potential of future observations of highly magnetized neutron stars with the more sensitive follow-up mission XL-Calibur.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

The Neil Gehrels Swift Observatory carried out prompt searches for gravitational wave (GW) events detected by the LIGO/Virgo Collaboration (LVC) during the second observing run ("O2"). Swift performed extensive tiling of eight LVC triggers, two of which had very low false-alarm rates (GW 170814 and the epochal GW 170817), indicating a high confidence of being astrophysical in origin; the latter was the first GW event to have an electromagnetic counterpart detected. In this paper we describe the follow-up performed during O2 and the results of our searches. No GW electromagnetic counterparts were detected; this result is expected, as GW 170817 remained the only astrophysical event containing at least one neutron star after LVC's later retraction of some events. A number of X-ray sources were detected, with the majority of identified sources being active galactic nuclei. We discuss the detection rate of transient X-ray sources and their implications in the O2 tiling searches. Finally, we describe the lessons learned during O2, and how these are being used to improve the \swift\ follow-up of GW events. In particular, we simulate a population of GRB afterglows to evaluate our source ranking system's ability to differentiate them from unrelated and uncatalogued X-ray sources. We find that $\approx$60-70% of afterglows whose jets are oriented towards Earth will be given high rank (i.e., "interesting" designation) by the completion of our second follow-up phase (assuming their location in the sky was observed), but that this fraction can be increased to nearly 100% by performing a third follow-up observation of sources exhibiting fading behavior.

GRB 160821B is a short duration gamma-ray burst (GRB) detected and localized by the Neil Gehrels Swift Observatory in the outskirts of a spiral galaxy at z=0.1613, at a projected physical offset of 16 kpc from the galaxy's center. We present X-ray, optical/nIR and radio observations of its counterpart and model them with two distinct components of emission: a standard afterglow, arising from the interaction of the relativistic jet with the surrounding medium, and a kilonova, powered by the radioactive decay of the sub-relativistic ejecta. Broadband modeling of the afterglow data reveals a weak reverse shock propagating backward into the jet, and a likely jet-break at 3.5 d. This is consistent with a structured jet seen slightly off-axis while expanding into a low-density medium. Analysis of the kilonova properties suggests a rapid evolution toward red colors, similar to AT2017gfo, and a low nIR luminosity, possibly due to the presence of a long-lived neutron star. The global properties of the environment, the inferred low mass (M_ej < 0.006 Msun) and velocities (v > 0.05 c) of lanthanide-rich ejecta are consistent with a binary neutron star merger progenitor.

The first simultaneous detection of a short gamma-ray burst (SGRB) with a gravitational-wave (GW) signal ushered in a new era of multi-messenger astronomy. In order to increase the number of SGRB-GW simultaneous detections, we need full sky coverage in the gamma-ray regime. BurstCube, a CubeSat for Gravitational Wave Counterparts, aims to expand sky coverage in order to detect and localize gamma-ray bursts (GRBs). BurstCube will be comprised of 4 Cesium Iodide scintillators coupled to arrays of Silicon photo-multipliers on a 6U CubeSat bus (a single U corresponds to cubic unit $\sim$10 cm $\times$ 10 cm $\times$ 10 cm) and will be sensitive to gamma-rays between 50 keV and 1 MeV, the ideal energy range for GRB prompt emission. BurstCube will assist current observatories, such as $Swift$ and $Fermi$, in the detection of GRBs as well as provide astronomical context to gravitational wave events detected by Advanced LIGO, Advanced Virgo, and KAGRA. BurstCube is currently in its development and testing phase to prepare for launch readiness in the fall of 2021. We present the mission concept, preliminary performance, and status.

We report on a search for high-energy counterparts to fast radio bursts (FRBs) with the Fermi Gamma-ray Burst Monitor (GBM), Fermi Large Area Telescope (LAT), and the Neil Gehrels Swift Observatory Burst Alert Telescope (BAT). We find no significant associations for any of the 23 FRBs in our sample, but report upper limits to the high-energy fluence for each on timescales of 0.1, 1, 10, and 100 s. We report lower limits on the ratio of the radio to high-energy fluence, $\frac{f_{r}}{f_{\gamma}}$, for timescales of 0.1 and 100 s. We discuss the implications of our non-detections on various proposed progenitor models for FRBs, including analogs of giant pulses from the Crab pulsar and hyperflares from magnetars. This work demonstrates the utility of analyses of high-energy data for FRBs in tracking down the nature of these elusive sources.

A Milky-Way Type Ia Supernova (SNIa) could be unidentified or even initially unnoticed, being dim in radio, X-rays, and neutrinos, and suffering large optical/IR extinction in the Galactic plane. But SNIa emit nuclear gamma-ray lines from $^{56}{\rm Ni}\to ^{56}{\rm Co}\to ^{56}{\rm Fe}$ radioactive decays. These lines fall within the Fermi/GBM energy range, and the $^{56}{\rm Ni}$ 158 keV line is detectable by Swift/BAT. Both instruments frequently monitor the Galactic plane, which is transparent to gamma rays. Thus GBM and BAT are ideal Galactic SNIa early warning systems. We simulate SNIa MeV light curves and spectra to show that GBM and BAT could confirm a Galactic SNIa explosion, followed by Swift localization and observation in X-rays and UVOIR band. The time of detection depends sensitively on the $^{56}{\rm Ni}$ distribution, and can be as early as a few days if $\gtrsim$ 10% of the $^{56}{\rm Ni}$ is present in the surface as suggested by SN2014J gamma data.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Though they play a key role in cultivating the cosmological environment and/or enabling our studies of it, there is still much we do not know about AGNs and GRBs, particularly the avenue in which and through which they supply radiation and energetic particles, namely their jets. This White Paper is the second of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is active galactic nuclei and gamma-ray bursts.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Energetic particles in our own galaxy, galactic cosmic rays (GCRs), engage in a complex interplay with the interstellar medium and magnetic fields in the galaxy, giving rise to many of its key characteristics. This White Paper is the first of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is galactic cosmic rays, supernova remnants, protostellar jets and superbubbles, and colliding wind binaries.

The bright transient AT2018cow has been unlike any other known type of transient. Its high brightness, rapid rise and decay and initially nearly featureless spectrum are unprecedented and difficult to explain using models for similar burst sources. We present evidence for faint gamma-ray emission continuing for at least 8 days, and featureless spectra in the ultraviolet bands -- both unusual for eruptive sources. The X-ray variability of the source has a burst-like character. The UV-optical spectrum does not show any CNO line but is well described by a blackbody. We demonstrate that a model invoking the tidal disruption of a 0.1 - 0.4 Msun Helium White Dwarf (WD) by a 100,000 to one million solar mass Black Hole (BH) located in the outskirts of galaxy Z~137-068 could provide an explanation for most of the characteristics shown in the multi-wavelength observations. A blackbody-like emission is emitted from an opaque photosphere, formed by the debris of the WD disruption. Broad features showing up in the optical/infrared spectra in the early stage are probably velocity broadened lines produced in a transient high-velocity outward moving cocoon. The asymmetric optical/infrared lines that appeared at a later stage are emission from an atmospheric layer when it detached from thermal equilibrium with the photosphere, which undergoes more rapid cooling. The photosphere shrinks when its temperature drops, and the subsequent infall of the atmosphere produced asymmetric line profiles. Additionally, a non-thermal jet might be present, emitting X-rays in the 10-150 keV band.

During its first observing run, in late 2015, the advanced LIGO facility announced 3 gravitational wave (GW) triggers to electromagnetic follow-up partners. Two of these have since been confirmed as being of astrophysical origin: both are binary black hole mergers at ~500 Mpc; the other trigger was later found not to be astrophysical. In this paper we report on the Swift follow up observations of the second and third triggers, including details of 21 X-ray sources detected; none of which can be associated with the GW event. We also consider the challenges that the next GW observing run will bring as the sensitivity and hence typical distance of GW events will increase. We discuss how to effectively use galaxy catalogues to prioritise areas for follow up, especially in the presence of distance estimates from the GW data. We also consider two galaxy catalogues and suggest that the high completeness at larger distances of the 2MASS Photometric Redshift Catalogue (2MPZ) makes it very well suited to optimise Swift follow-up observations.

Gamma-ray bursts (GRBs) occurring in the local Universe constitute an interesting sub-class of the GRB family, since their luminosity is on average lower than that of their cosmological analogs. We aim to contribute to the study of local bursts by reporting the case of GRB 171205A. This source was discovered by Swift Burst Alert Telescope (BAT) on 2017, December 5 and soon associated with a low redshift host galaxy (z=0.037), and an emerging SN (SN 2017iuk). We analyzed the full Swift, dataset, comprising the UV-Optical Telescope (UVOT), X-ray Telescope (XRT) and BAT data. In addition, we employed the Konus-Wind high energy data as a valuable extension at gamma-ray energies. The photometric SN signature is clearly visible in the UVOT u, b and v filters. The maximum emission is reached at ~ 13 (rest frame) days, and the whole bump resembles that of SN 2006aj, but lower in magnitude and with a shift in time of +2 d. A prebump in the v-band is also clearly visible, and this is the first time that such a feature is not observed achromatically in GRB-SNe. Its physical origin cannot be easily explained. The X-ray spectrum shows an intrinsic Hydrogen column density N_H,int = 7.4(+4.1 -3.6) X 10^20 / cm^2$, which is at the low end of the N_H,int, even considering just low redshift GRBs. The spectrum also features a thermal component, which is quite common in GRBs associated with SNe, but whose origin is still a matter of debate. Finally, the isotropic energy in the gamma-ray band, E_iso = 2.18(+0.63 -0.50) X 10^49 erg, is lower than those of cosmological GRBs. Combining this value with the peak energy in the same band, E_p=125(+141 -37) keV, implies that GRB 171205A is an outlier of the Amati relation, as are some other low redshift GRBs, and its emission mechanism should be different from that of canonical, farther away GRBs.

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $\gamma$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

We have observed a large glitch in the Crab pulsar (PSR B0531+21). The glitch occurred around MJD 58064 (2017 November 8) when the pulsar underwent an increase in the rotation rate of $\Delta \nu = 1.530 \times 10^{-5}$ Hz, corresponding to a fractional increase of $\Delta \nu / \nu = 0.516 \times 10^{-6}$ making this event the largest glitch ever observed in this source. Due to our high-cadence and long-dwell time observations of the Crab pulsar we are able to partially resolve a fraction of the total spin-up of the star. This delayed spin-up occurred over a timescale of $\sim$1.7 days and is similar to the behaviour seen in the 1989 and 1996 large Crab pulsar glitches. The spin-down rate also increased at the glitch epoch by $\Delta \dot{\nu} / \dot{\nu} = 7 \times 10^{-3}$. In addition to being the largest such event observed in the Crab, the glitch occurred after the longest period of glitch inactivity since at least 1984 and we discuss a possible relationship between glitch size and waiting time. No changes to the shape of the pulse profile were observed near the glitch epoch at 610 MHz or 1520 MHz, nor did we identify any changes in the X-ray flux from the pulsar. The long-term recovery from the glitch continues to progress as $\dot{\nu}$ slowly rises towards pre-glitch values. In line with other large Crab glitches, we expect there to be a persistent change to $\dot{\nu}$. We continue to monitor the long-term recovery with frequent, high quality observations.

We present a catalog of hard X-ray sources detected in the first 105 months of observations with the Burst Alert Telescope (BAT) coded mask imager on board the Swift observatory. The 105 month Swift-BAT survey is a uniform hard X-ray all-sky survey with a sensitivity of $8.40\times 10^{-12}\ {\rm erg\ s^{-1}\ cm^{-2}}$ over 90% of the sky and $7.24\times 10^{-12}\ {\rm erg\ s^{-1}\ cm^{-2}}$ over 50% of the sky in the 14-195 keV band. The Swift-BAT 105 month catalog provides 1632 (422 new detections) hard X-ray sources in the 14-195 keV band above the 4.8{\sigma} significance level. Adding to the previously known hard X-ray sources, 34% (144/422) of the new detections are identified as Seyfert AGN in nearby galaxies (z<0.2). The majority of the remaining identified sources are X-ray binaries (7%, 31) and blazars/BL Lac objects (10%, 43). As part of this new edition of the Swift-BAT catalog, we release eight-channel spectra and monthly sampled light curves for each object in the online journal and at the Swift-BAT 105 month Web site.

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.

A long-standing paradigm in astrophysics is that collisions- or mergers- of two neutron stars (NSs) form highly relativistic and collimated outflows (jets) powering gamma-ray bursts (GRBs) of short (< 2 s) duration. However, the observational support for this model is only indirect. A hitherto outstanding prediction is that gravitational wave (GW) events from such mergers should be associated with GRBs, and that a majority of these GRBs should be off-axis, that is, they should point away from the Earth. Here we report the discovery of the X-ray counterpart associated with the GW event GW170817. While the electromagnetic counterpart at optical and infrared frequencies is dominated by the radioactive glow from freshly synthesized r-process material in the merger ejecta, known as kilonova, observations at X-ray and, later, radio frequencies exhibit the behavior of a short GRB viewed off-axis. Our detection of X-ray emission at a location coincident with the kilonova transient provides the missing observational link between short GRBs and GWs from NS mergers, and gives independent confirmation of the collimated nature of the GRB emission.

With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading ultraviolet emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultra-relativistic, highly collimated ejecta (a gamma-ray burst afterglow).

We present Hubble Space Telescope and Chandra imaging, combined with Very Large Telescope MUSE integral field spectroscopy of the counterpart and host galaxy of the first binary neutron star merger detected via gravitational wave emission by LIGO & Virgo, GW170817. The host galaxy, NGC 4993, is an S0 galaxy at z=0.009783. There is evidence for large, face-on spiral shells in continuum imaging, and edge-on spiral features visible in nebular emission lines. This suggests that NGC 4993 has undergone a relatively recent (<1 Gyr) ``dry'' merger. This merger may provide the fuel for a weak active nucleus seen in Chandra imaging. At the location of the counterpart, HST imaging implies there is no globular or young stellar cluster, with a limit of a few thousand solar masses for any young system. The population in the vicinity is predominantly old with <1% of any light arising from a population with ages <500 Myr. Both the host galaxy properties and those of the transient location are consistent with the distributions seen for short-duration gamma-ray bursts, although the source position lies well within the effective radius (r_e ~ 3 kpc), providing an r_e-normalized offset that is closer than ~90% of short GRBs. For the long delay time implied by the stellar population, this suggests that the kick velocity was significantly less than the galaxy escape velocity. We do not see any narrow host galaxy interstellar medium features within the counterpart spectrum, implying low extinction, and that the binary may lie in front of the bulk of the host galaxy.

S5 0716+714 is a well known BL Lac object, one of the brightest and most active blazars. The discovery in the Very High Energy band (VHE, E > 100 GeV) by MAGIC happened in 2008. In January 2015 the source went through the brightest optical state ever observed, triggering MAGIC follow-up and a VHE detection with 13{\sigma} significance (ATel 6999). Rich multiwavelength coverage of the flare allowed us to construct the broad-band spectral energy distribution of S5 0716+714 during its brightest outburst. In this work we will present the preliminary analysis of MAGIC and Fermi-LAT data of the flaring activity in January and February 2015 for the HE (0.1 < HE < 300 GeV) and VHE band, together with radio (Mets\"ahovi, OVRO, VLBA, Effelsberg), sub-millimeter (SMA), optical (Tuorla, Perkins, Steward, AZT-8+ST7, LX-200, Kanata), X-ray and UV (Swift-XRT and UVOT), in the same time-window and discuss the time variability of the multiwavelength light curves during this impressive outburst.

We report the identification of a bright hard X-ray source dominating the M31 bulge above 25 keV from a simultaneous NuSTAR-Swift observation. We find that this source is the counterpart to Swift J0042.6+4112, which was previously detected in the Swift BAT All-sky Hard X-ray Survey. This Swift BAT source had been suggested to be the combined emission from a number of point sources; our new observations have identified a single X-ray source from 0.5 to 50 keV as the counterpart for the first time. In the 0.5-10 keV band, the source had been classified as an X-ray binary candidate in various Chandra and XMM studies; however, since it was not clearly associated with Swift J0042.6+4112, the previous E < 10 keV observations did not generate much attention. This source has a spectrum with a soft X-ray excess (kT~ 0.2 keV) plus a hard spectrum with a power law of Gamma ~ 1 and a cutoff around 15-20 keV, typical of the spectral characteristics of accreting pulsars. Unfortunately, any potential pulsation was undetected in the NuSTAR data, possibly due to insufficient photon statistics. The existing deep HST images exclude high-mass (>3 Msun) donors at the location of this source. The best interpretation for the nature of this source is an X-ray pulsar with an intermediate-mass (<3 Msun) companion or a symbiotic X-ray binary. We discuss other possibilities in more detail.

We report the results of our observing campaign on GRB140903A, a nearby (z=0.351) short duration (T90~0.3 s) gamma-ray burst discovered by Swift. We monitored the X-ray afterglow with Chandra up to 21 days after the burst, and detected a steeper decay of the X-ray flux after approximately 1 day. Continued monitoring at optical and radio wavelengths showed a similar decay in flux at nearly the same time, and we interpret it as evidence of a narrowly collimated jet. By using the standard fireball model to describe the afterglow evolution, we derive a jet opening angle of 5 deg and a collimation-corrected total energy release of 2E50 erg. We further discuss the nature of the GRB progenitor system. Three main lines disfavor a massive star progenitor: the properties of the prompt gamma-ray emission, the age and low star-formation rate of the host galaxy, and the lack of a bright supernova. We conclude that this event was likely originated by a compact binary merger.

To date, the Burst Alert Telescope (BAT) onboard Swift has detected ~ 1000 gamma-ray bursts (GRBs), of which ~ 360 GRBs have redshift measurements, ranging from z = 0.03 to z = 9.38. We present the analyses of the BAT-detected GRBs for the past ~ 11 years up through GRB151027B. We report summaries of both the temporal and spectral analyses of the GRB characteristics using event data (i.e., data for each photon within approximately 250 s before and 950 s after the BAT trigger time), and discuss the instrumental sensitivity and selection effects of GRB detections. We also explore the GRB properties with redshift when possible. The result summaries and data products are available at http://swift.gsfc.nasa.gov/results/batgrbcat/index.html . In addition, we perform searches for GRB emissions before or after the event data using the BAT survey data. We estimate the false detection rate to be only one false detection in this sample. There are 15 ultra-long GRBs (~ 2% of the BAT GRBs) in this search with confirmed emission beyond ~ 1000 s of event data, and only two GRBs (GRB100316D and GRB101024A) with detections in the survey data prior to the starting of event data. (Some figures shown here are in lower resolution due to the size limit on arXiv. The full resolution version can be found at http://swift.gsfc.nasa.gov/results/batgrbcat/3rdBATcatalog.pdf )

The gamma-ray burst (GRB) rate is essential for revealing the connection between GRBs, supernovae and stellar evolution. Additionally, the GRB rate at high redshift provides a strong probe of star formation history in the early universe. While hundreds of GRBs are observed by Swift, it remains difficult to determine the intrinsic GRB rate due to the complex trigger algorithm of Swift. Current studies of the GRB rate usually approximate the Swift trigger algorithm by a single detection threshold. However, unlike the previously flown GRB instruments, Swift has over 500 trigger criteria based on photon count rate and additional image threshold for localization. To investigate possible systematic biases and explore the intrinsic GRB properties, we develop a program that is capable of simulating all the rate trigger criteria and mimicking the image threshold. Our simulations show that adopting the complex trigger algorithm of Swift increases the detection rate of dim bursts. As a result, our simulations suggest bursts need to be dimmer than previously expected to avoid over-producing the number of detections and to match with Swift observations. Moreover, our results indicate that these dim bursts are more likely to be high redshift events than low-luminosity GRBs. This would imply an even higher cosmic GRB rate at large redshifts than previous expectations based on star-formation rate measurements, unless other factors, such as the luminosity evolution, are taken into account. The GRB rate from our best result gives a total number of 4571^{+829}_{-1584} GRBs per year that are beamed toward us in the whole universe. SPECIAL NOTE (2015.05.16): This new version incorporates an erratum. All the GRB rate normalizations ($R_{\rm GRB}(z=0)$) should be a factor of 2 smaller than previously reported. Please refer to the Appendix for more details. We sincerely apologize for the mistake.

To draw inferences about gamma-ray burst (GRB) source populations based on Swift observations, it is essential to understand the detection efficiency of the Swift burst alert telescope (BAT). This study considers the problem of modeling the Swift/BAT triggering algorithm for long GRBs, a computationally expensive procedure, and models it using machine learning algorithms. A large sample of simulated GRBs from Lien 2014 is used to train various models: random forests, boosted decision trees (with AdaBoost), support vector machines, and artificial neural networks. The best models have accuracies of $\gtrsim97\%$ ($\lesssim 3\%$ error), which is a significant improvement on a cut in GRB flux which has an accuracy of $89.6\%$ ($10.4\%$ error). These models are then used to measure the detection efficiency of Swift as a function of redshift $z$, which is used to perform Bayesian parameter estimation on the GRB rate distribution. We find a local GRB rate density of $n_0 \sim 0.48^{+0.41}_{-0.23} \ {\rm Gpc}^{-3} {\rm yr}^{-1}$ with power-law indices of $n_1 \sim 1.7^{+0.6}_{-0.5}$ and $n_2 \sim -5.9^{+5.7}_{-0.1}$ for GRBs above and below a break point of $z_1 \sim 6.8^{+2.8}_{-3.2}$. This methodology is able to improve upon earlier studies by more accurately modeling Swift detection and using this for fully Bayesian model fitting. The code used in this is analysis is publicly available online (https://github.com/PBGraff/SwiftGRB_PEanalysis).

Compact binary system mergers are expected to generate gravitational radiation detectable by ground-based interferometers. A subset of these, the merger of a neutron star with another neutron star or a black hole, are also the most popular model for the production of short gamma-ray bursts (GRBs). The Swift Burst Alert Telescope (BAT) and the Fermi Gamma-ray Burst Monitor (GBM) trigger on short GRBs (SGRBs) at rates that reflect their relative sky exposures, with the BAT detecting 10 per year compared to about 45 for GBM. We examine the SGRB populations detected by Swift BAT and Fermi GBM. We find that the Swift BAT triggers on weaker SGRBs than Fermi GBM, providing they occur close to the center of the BAT field-of-view, and that the Fermi GBM SGRB detection threshold remains flatter across its field-of-view. Overall, these effects combine to give the instruments the same average sensitivity, and account for the SGRBs that trigger one instrument but not the other. We do not find any evidence that the BAT and GBM are detecting significantly different populations of SGRBs. Both instruments can detect untriggered SGRBs using ground searches seeded with time and position. The detection of SGRBs below the on-board triggering sensitivities of Swift BAT and Fermi GBM increases the possibility of detecting and localizing the electromagnetic counterparts of gravitational wave events seen by the new generation of gravitational wave detectors.

Isolated Swift gamma-ray burst (GRB) pulses, like their higher-energy BATSE counterparts, emit the bulk of their pulsed emission as a hard-to-soft component that can be fitted by the Norris et al. (2005) empirical pulse model. This signal is overlaid by a fainter, three-peaked signal that can be modeled by an empirical wave-like function (Hakkila and Preece, 2014): the two fits combine to reproduce GRB pulses with distinctive three-peaked shapes. The precursor peak appears on or before the pulse rise and is often the hardest component, the central peak is the brightest, and the decay peak converts exponentially decaying emission into a long, soft, power-law tail. Accounting for systematic instrumental differences, the general characteristics of the fitted pulses are remarkably similar. Isolated GRB pulses are dominated by hard-to-soft evolution; this is more pronounced for asymmetric pulses than for symmetric ones. Isolated GRB pulses can also exhibit intensity tracking behaviors that, when observed, are tied to the timing of the three peaks: pulses with the largest maximum hardnesses are hardest during the precursor, those with smaller maximum hardnesses are hardest during the central peak, and all pulses can re-harden during the central peak and/or during the decay peak. Since these behaviors are essentially seen in all isolated pulses, the distinction between "hard-to-soft" and "intensity-tracking" pulses really no longer applies. Additionally, the triple-peaked nature of isolated GRB pulses seems to indicate that energy is injected on three separate occasions during the pulse duration: theoretical pulse models need to account for this.

Gamma-ray bursts (GRBs) are the most luminous explosions in the universe, yet the nature and physical properties of their energy sources are far from understood. Very important clues, however, can be inferred by studying the afterglows of these events. We present optical and X-ray observations of GRB 130831A obtained by Swift, Chandra, Skynet, RATIR, Maidanak, ISON, NOT, LT and GTC. This burst shows a steep drop in the X-ray light-curve at $\simeq 10^5$ s after the trigger, with a power-law decay index of $\alpha \sim 6$. Such a rare behaviour cannot be explained by the standard forward shock (FS) model and indicates that the emission, up to the fast decay at $10^5$ s, must be of "internal origin", produced by a dissipation process within an ultrarelativistic outflow. We propose that the source of such an outflow, which must produce the X-ray flux for $\simeq 1$ day in the cosmological rest frame, is a newly born magnetar or black hole. After the drop, the faint X-ray afterglow continues with a much shallower decay. The optical emission, on the other hand, shows no break across the X-ray steep decrease, and the late-time decays of both the X-ray and optical are consistent. Using both the X-ray and optical data, we show that the emission after $\simeq 10^5$ s can be explained well by the FS model. We model our data to derive the kinetic energy of the ejecta and thus measure the efficiency of the central engine of a GRB with emission of internal origin visible for a long time. Furthermore, we break down the energy budget of this GRB into the prompt emission, the late internal dissipation, the kinetic energy of the relativistic ejecta, and compare it with the energy of the associated supernova, SN 2013fu.

We report here the discovery by the Intermediate Palomar Transient Factory (iPTF) of iPTF14yb, a luminous ($M_{r}\approx-27.8$ mag), cosmological (redshift 1.9733), rapidly fading optical transient. We demonstrate, based on probabilistic arguments and a comparison with the broader population, that iPTF14yb is the optical afterglow of the long-duration gamma-ray burst GRB 140226A. This marks the first unambiguous discovery of a GRB afterglow prior to (and thus entirely independent of) an associated high-energy trigger. We estimate the rate of iPTF14yb-like sources (i.e., cosmologically distant relativistic explosions) based on iPTF observations, inferring an all-sky value of $\Re_{\mathrm{rel}}=610$ yr$^{-1}$ (68% confidence interval of 110-2000 yr$^{-1}$). Our derived rate is consistent (within the large uncertainty) with the all-sky rate of on-axis GRBs derived by the Swift satellite. Finally, we briefly discuss the implications of the nondetection to date of bona fide "orphan" afterglows (i.e., those lacking detectable high-energy emission) on GRB beaming and the degree of baryon loading in these relativistic jets.

GRB 130925A was an unusual GRB, consisting of 3 distinct episodes of high-energy emission spanning $\sim$20 ks, making it a member of the proposed category of `ultra-long' bursts. It was also unusual in that its late-time X-ray emission observed by Swift was very soft, and showed a strong hard-to-soft spectral evolution with time. This evolution, rarely seen in GRB afterglows, can be well modelled as the dust-scattered echo of the prompt emission, with stringent limits on the contribution from the normal afterglow (i.e. external shock) emission. We consider and reject the possibility that GRB 130925A was some form of tidal disruption event, and instead show that if the circumburst density around GRB 130925A is low, the long duration of the burst and faint external shock emission are naturally explained. Indeed, we suggest that the ultra-long GRBs as a class can be explained as those with low circumburst densities, such that the deceleration time (at which point the material ejected from the nascent black hole is decelerated by the circumburst medium) is $\sim$20 ks, as opposed to a few hundred seconds for the normal long GRBs. The increased deceleration radius means that more of the ejected shells can interact before reaching the external shock, naturally explaining both the increased duration of GRB 130925A, the duration of its prompt pulses, and the fainter-than-normal afterglow.

The Swift/Burst Alert Telescope (BAT) hard X-ray transient monitor provides near real-time coverage of the X-ray sky in the energy range 15-50 keV. The BAT observes 88% of the sky each day with a detection sensitivity of 5.3 mCrab for a full-day observation and a time resolution as fine as 64 seconds. The three main purposes of the monitor are (1) the discovery of new transient X-ray sources, (2) the detection of outbursts or other changes in the flux of known X-ray sources, and (3) the generation of light curves of more than 900 sources spanning over eight years. The primary interface for the BAT transient monitor is a public web page. Between 2005 February 12 and 2013 April 30, 245 sources have been detected in the monitor, 146 of them persistent and 99 detected only in outburst. Among these sources, 17 were previously unknown and were discovered in the transient monitor. In this paper, we discuss the methodology and the data processing and filtering for the BAT transient monitor and review its sensitivity and exposure. We provide a summary of the source detections and classify them according to the variability of their light curves. Finally, we review all new BAT monitor discoveries; for the new sources that are previously unpublished, we present basic data analysis and interpretations.

The long gamma-ray burst (GRB) rate is essential for revealing the connection between GRBs, supernovae and stellar evolution. Additionally, the GRB rate at high redshift provides a strong probe of star formation history in the early universe. While hundreds of GRBs are observed by Swift, it remains difficult to determine the intrinsic GRB rate due to the complex trigger algorithm of Swift. Current studies usually approximate the Swift trigger algorithm by a single detection threshold. However, unlike the previously flown GRB instruments, Swift has over 500 trigger criteria based on photon count rate and additional image threshold for localization. To investigate possible systematic biases and explore the intrinsic GRB properties, we developed a program that is capable of simulating all the rate trigger criteria and mimicking the image trigger threshold. We use this program to search for the intrinsic GRB rate. Our simulations show that adopting the complex trigger algorithm of Swift increases the detection rate of dim bursts. As a result, we find that either the GRB rate is much higher than previously expected at large redshift, or the luminosity evolution is non-negligible. We will discuss the best results of the GRB rate in our search, and their impact on the star-formation history.

The origin of the diffuse extragalactic gamma-ray background (EGB) has been intensively studied but remains unsettled. Current popular source candidates include unresolved star-forming galaxies, starburst galaxies, and blazars. In this paper we calculate the EGB contribution from the interactions of cosmic rays accelerated by Type Ia supernovae (SNe), extending earlier work which only included core-collapse SNe. We consider Type Ia events in star-forming galaxies, but also in quiescent galaxies that lack star formation. For star-forming galaxies, consistently including Type Ia events makes little change to the star-forming EGB prediction, so long as both SN types have the same cosmic-ray acceleration efficiencies in star-forming galaxies. Thus, our updated EGB estimate continues to show that star-forming galaxies can represent a substantial portion of the signal measured by Fermi. For quiescent galaxies, conversely, we find a wide range of possibilities for the EGB contribution. The dominant uncertainty we investigated comes from the mass in hot gas, which provides targets for cosmic rays; total gas masses are as yet poorly known, particularly at larger radii. Additionally, the EGB estimation is very sensitive to the cosmic-ray acceleration efficiency and confinement, especially in quiescent galaxies. In the most optimistic allowed scenarios, quiescent galaxies can be an important source of the EGB. In this case, star-forming galaxies and quiescent galaxies together will dominate the EGB and leave little room for other contributions. If other sources, such as blazars, are found to have important contributions to the EGB, then either the gas mass or cosmic-ray content of quiescent galaxies must be significantly lower than in their star-forming counterparts. In any case, improved Fermi EGB measurements will provide important constraints on hot gas and cosmic rays in quiescent galaxies.

Radio properties of supernova outbursts remain poorly understood despite longstanding campaigns following events discovered at other wavelengths. After ~ 30 years of observations, only ~ 50 supernovae have been detected at radio wavelengths, none of which are Type Ia. Even the most radio-loud events are ~ 10^4 fainter in the radio than in the optical; to date, such intrinsically dim objects have only been visible in the very local universe. The detection and study of radio supernovae (RSNe) will be fundamentally altered and dramatically improved as the next generation of radio telescopes comes online, including EVLA, ASKAP, and MeerKAT, and culminating in the Square Kilometer Array (SKA); the latter should be > 50 times more sensitive than present facilities. SKA can repeatedly scan large (> 1 deg^2) areas of the sky, and thus will discover RSNe and other transient sources in a new, automatic, untargeted, and unbiased way. We estimate SKA will be able to detect core-collapse RSNe out to redshift z ~ 5, with an all-redshift rate ~ 620 events yr^-1 deg^-2, assuming a survey sensitivity of 50 nJy and radio lightcurves like those of SN 1993J. Hence SKA should provide a complete core-collapse RSN sample that is sufficient for statistical studies of radio properties of core-collapse supernovae. EVLA should find ~ 160 events yr^-1 deg^-2 out to redshift z ~ 3, and other SKA precursors should have similar detection rates. We also provided recommendations of the survey strategy to maximize the RSN detections of SKA. This new radio core-collapse supernovae sample will complement the detections from the optical searches, such as the LSST, and together provide crucial information on massive star evolution, supernova physics, and the circumstellar medium, out to high redshift. Additionally, SKA may yield the first radio Type Ia detection via follow-up of nearby events discovered at other wavelengths.

The cumulative (anti)neutrino production from all core-collapse supernovae within our cosmic horizon gives rise to the diffuse supernova neutrino background (DSNB), which is on the verge of detectability. The observed flux depends on supernova physics, but also on the cosmic history of supernova explosions; currently, the cosmic supernova rate introduces a substantial (+/-40%) uncertainty, largely through its absolute normalization. However, a new class of wide-field, repeated-scan (synoptic) optical sky surveys is coming online, and will map the sky in the time domain with unprecedented depth, completeness, and dynamic range. We show that these surveys will obtain the cosmic supernova rate by direct counting, in an unbiased way and with high statistics, and thus will allow for precise predictions of the DSNB. Upcoming sky surveys will substantially reduce the uncertainties in the DSNB source history to an anticipated +/-5% that is dominated by systematics, so that the observed high-energy flux thus will test supernova neutrino physics. The portion of the universe (z < 1) accessible to upcoming sky surveys includes the progenitors of a large fraction (~ 87%) of the expected 10-26 MeV DSNB event rate. We show that precision determination of the (optically detected) cosmic supernova history will also make the DSNB into a strong probe of an extra flux of neutrinos from optically invisible supernovae, which may be unseen either due to unexpected large dust obscuration in host galaxies, or because some core-collapse events proceed directly to black hole formation and fail to give an optical outburst.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

Large synoptic (repeated scan) imaging sky surveys are poised to observe enormous numbers of core-collapse supernovae. We quantify the discovery potential of upcoming projects, including DES, Pan-STARRS, and LSST. These surveys will map out the cosmic core-collapse supernova redshift distribution via direct counting, with very small statistical uncertainties out to a redshift depth which is a strong function of the survey limiting magnitude. This supernova redshift history encodes rich information about cosmology, star formation, and supernova astrophysics and phenomenology; the large statistics of the supernova sample will be crucial to disentangle possible degeneracies among these issues. For example, the cosmic supernova rate can be measured to high precision out to z ~ 0.5 for all core-collapse types, and out to redshift z ~ 1 for Type IIn events. Because of the tight link between supernovae and star formation, synoptic sky surveys will provide precision measurements of the normalization and z < 1 history of cosmic star-formation rate in a manner independent of and complementary to current methods. Furthermore, Type II supernovae can serve as distance indicators and would independently cross-check Type Ia distances measured in the same surveys. Arguably the largest and least-controlled uncertainty in these efforts comes from dust obscuration of supernovae in their host galaxies; we outline a strategy to determine empirically the obscuration properties by leveraging the large supernova samples over a broad range of redshift. [abstract abridged]