Search

The Einstein Probe (EP) achieved its first detection and localization of a bright X-ray flare, EP240219a, on February 19, 2024, during its commissioning phase. Subsequent targeted searches triggered by the EP240219a alert identified a faint, untriggered gamma-ray burst (GRB) in the archived data of Fermi/GBM, Swift/BAT, Insight-HXMT/HE and INTEGRAL/SPI-ACS. The EP/WXT light curve reveals a long duration of approximately 160 seconds with a slow decay, whereas the Fermi/GBM light curve shows a total duration of approximately 70 seconds. The peak in the Fermi/GBM light curve occurs slightly later with respect to the peak seen in the EP/WXT light curve. Our spectral analysis shows that a single cutoff power-law model effectively describes the joint EP/WXT-Fermi/GBM spectra in general, indicating coherent broad emission typical of GRBs. The model yielded a photon index of $\sim -1.70 \pm 0.05$ and a peak energy of $\sim 257 \pm 134$ keV. After detection of GRB 240219A, long-term observations identified several candidates in optical and radio wavelengths, none of which was confirmed as the afterglow counterpart during subsequent optical and near-infrared follow-ups. The analysis of GRB 240219A classifies it as an X-ray rich GRB with a high peak energy, presenting both challenges and opportunities for studying the physical origins of X-ray flashes (XRFs), X-ray rich GRBs (XRRs), and classical GRBs (C-GRBs). Furthermore, linking the cutoff power-law component to non-thermal synchrotron radiation suggests that the burst is driven by a Poynting flux-dominated outflow., Comment: 14 pages, 8 figures, 3 tables

We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wave Transient Catalogs (GWTC-3). Targeted searches were carried out on the entire GW sample using the maximum--likelihood NITRATES pipeline on the BAT data made available via the GUANO infrastructure. We do not detect any significant electromagnetic emission that is temporally and spatially coincident with any of the GW candidates. We report flux upper limits in the 15-350 keV band as a function of sky position for all the catalog candidates. For GW candidates where the Swift-BAT false alarm rate is less than 10$^{-3}$ Hz, we compute the GW--BAT joint false alarm rate. Finally, the derived Swift-BAT upper limits are used to infer constraints on the putative electromagnetic emission associated with binary black hole mergers., Comment: 50 pages, 10 figures, 4 tables

GW230529 is the first compact binary coalescence detected by the LIGO-Virgo-KAGRA collaboration with at least one component mass confidently in the lower mass-gap, corresponding to the range 3-5$M_{\odot}$. If interpreted as a neutron star-black hole merger, this event has the most symmetric mass ratio detected so far and therefore has a relatively high probability of producing electromagnetic (EM) emission. However, no EM counterpart has been reported. At the merger time $t_0$, Swift-BAT and Fermi-GBM together covered 100$\%$ of the sky. Performing a targeted search in a time window $[t_0-20 \text{s},t_0+20 \text{s}]$, we report no detection by the Swift-BAT and the Fermi-GBM instruments. Combining the position-dependent $\gamma-$ray flux upper limits and the gravitational-wave posterior distribution of luminosity distance, sky localization and inclination angle of the binary, we derive constraints on the characteristic luminosity and structure of the jet possibly launched during the merger. Assuming a top-hat jet structure, we exclude at 90$\%$ credibility the presence of a jet which has at the same time an on-axis isotropic luminosity $\gtrsim 10^{48}$ erg s$^{-1}$, in the bolometric band 1 keV-10 MeV, and a jet opening angle $\gtrsim 15$ deg. Similar constraints are derived testing other assumptions about the jet structure profile. Excluding GRB 170817A, the luminosity upper limits derived here are below the luminosity of any GRB observed so far., Comment: 18 pages, 1 table, 11 figures

The Advanced X-ray Imaging Satellite (AXIS) promises revolutionary science in the X-ray and multi-messenger time domain. AXIS will leverage excellent spatial resolution (<1.5 arcsec), sensitivity (80x that of Swift), and a large collecting area (5-10x that of Chandra) across a 24-arcmin diameter field of view to discover and characterize a wide range of X-ray transients from supernova-shock breakouts to tidal disruption events to highly variable supermassive black holes. The observatory's ability to localize and monitor faint X-ray sources opens up new opportunities to hunt for counterparts to distant binary neutron star mergers, fast radio bursts, and exotic phenomena like fast X-ray transients. AXIS will offer a response time of <2 hours to community alerts, enabling studies of gravitational wave sources, high-energy neutrino emitters, X-ray binaries, magnetars, and other targets of opportunity. This white paper highlights some of the discovery science that will be driven by AXIS in this burgeoning field of time domain and multi-messenger astrophysics., Comment: This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at http://axis.astro.umd.edu

The arrival of gravitational wave astronomy and a growing number of time-domain focused observatories are set to lead to a increasing number of detections of short gamma-ray bursts (GRBs) launched with a moderate inclination to Earth. Being nearby events, these are also prime candidates for very long-term follow-up campaigns and very-long-baseline interferometry (VLBI), which has implications for multi-messenger modelling, data analysis, and statistical inference methods applied to these sources. Here we present a comprehensive modelling update that directly incorporates into afterglowpy astrometric observations of the GRB position, Poissonian statistics for faint sources, and modelling of a trans-relativistic population of electrons. We use the revolutionary event GW170817 to demonstrate the impact of these extensions both for the best-fit physics parameters and model selection methods that assess the statistical significance of additional late-time emission components. By including in our analysis the latest Chandra X-ray observations of GRB 170817A, we find only weak evidence (less than two sigma) for a new emission component at late times, which makes for a slightly more natural fit to the centroid evolution and prediction for the external medium density., Comment: 23 pages, 12 figures, submitted to ApJ

We present analytical and numerical models of the bright long GRB 210822A at $z=1.736$. The intrinsic extreme brightness exhibited in the optical, which is very similar to other bright GRBs (e.g., GRBs 080319B, 130427A, 160625A 190114C, and 221009A), makes GRB 210822A an ideal case for studying the evolution of this particular kind of GRB. We use optical data from the RATIR instrument starting at $T+315.9$ s, with publicly available optical data from other ground-based observatories, as well as Swift/UVOT, and X-ray data from the Swift/XRT instrument. The temporal profiles and spectral properties during the late stages align consistently with the conventional forward shock model, complemented by a reverse shock element that dominates optical emissions during the initial phases ($T<300$ s). Furthermore, we observe a break at $T=80000$s that we interpreted as evidence of a jet break, which constrains the opening angle to be about $\theta_\mathrm{j}=(3-5)$ degrees. Finally, we apply a machine-learning technique to model the multi-wavelength light curve of GRB 210822A using the AFTERGLOWPY library. We estimate the angle of sight $\theta_{obs}=(6.4 \pm 0.1) \times 10^{-1}$ degrees, the energy $E_0=(7.9 \pm 1.6)\times 10^{53}$ ergs, the electron index $p=2.54 \pm 0.10$, the thermal energy fraction in electrons $\epsilon_\mathrm{e}=(4.63 \pm 0.91) \times 10^{-5}$ and in the magnetic field $\epsilon_\mathrm{B}= (8.66 \pm 1.01) \times 10^{-6}$, the efficiency $\chi = 0.89 \pm 0.01$, and the density of the surrounding medium $n_\mathrm{0} = 0.85 \pm 0.01 cm^{-3}$., Comment: Resubmitted to MNRAS after moderate revision, 12 pages, 6 figures

Kilonovae are a rare class of astrophysical transients powered by the radioactive decay of nuclei heavier than iron, synthesized in the merger of two compact objects. Over the first few days, the kilonova evolution is dominated by a large number of radioactive isotopes contributing to the heating rate. On timescales of weeks to months, its behavior is predicted to differ depending on the ejecta composition and merger remnant. However, late-time observations of known kilonovae are either missing or limited. Here we report observations of a luminous red transient with a quasi-thermal spectrum, following an unusual gamma-ray burst of long duration. We classify this thermal emission as a kilonova and track its evolution up to two months after the burst. At these late times, the recession of the photospheric radius and the rapidly-decaying bolometric luminosity ($L_{\rm bol}\propto t^{-2.7\pm 0.4}$) support the recombination of lanthanide-rich ejecta as they cool., Comment: 47 pages, 14 figures, 9 tables; submitted; a minor typo fixed

Kilonovae are a novel class of astrophysical transients, and the only observationally-confirmed site of rapid neutron capture nucleosynthesis (the r-process) in the Universe. To date, only a handful of kilonovae have been detected, with just a single spectroscopically-observed event (AT 2017gfo). Spectra of AT 2017gfo provided evidence for the formation of elements heavier than iron; however, these spectra were collected during the first ~ 10 days, when emission from light r-process elements dominates the observations. Heavier elements, if synthesised, are expected to shape the late-time evolution of the kilonova, beyond the phases for which we have spectral observations. Here we present spectroscopic observations of a rapidly-reddening thermal transient, following the gamma-ray burst, GRB 230307A. Early (2.4 day) optical spectroscopy identifies the presence of a hot (T ~ 6700 K) thermal continuum. By 29 days, this component has expanded and cooled significantly (T ~ 640 K), yet it remains optically thick, indicating the presence of high-opacity ejecta. We show that these properties can only be explained by the merger of compact objects, and further, leads us to infer the production of the heavy lanthanide elements. We identify several spectral features (in both absorption and emission), whose cause can be explained by newly-synthesised heavy elements. This event marks only the second recorded spectroscopic evidence for the synthesis of r-process elements, and the first to be observed at such late times., Comment: 41 pages, 11 figures, 4 tables. Submitted - comments welcome!

Swift has now completed 18 years of mission, during which it discovered thousands of gamma-ray bursts (GRBs) as well as new classes of high-energy transient phenomena. Its first breakthrough result was the localization of short duration GRBs, which enabled for redshift measurements and kilonova searches. Swift, in synergy with the Hubble Space Telescope and a wide array of ground-based telescopes, provided the first tantalizing evidence of a kilonova in the aftermath of a short GRB. In 2017, Swift observations of the gravitational wave event GW170817 captured the early UV photons from the kilonova AT2017gfo, opening a new window into the physics of kilonovae. Since then, Swift has continued to expand the sample of known kilonovae, leading to the surprising discovery of a kilonova in a long duration GRB. This article will discuss recent advances in the study of kilonovae driven by the fundamental contribution of Swift., Comment: 26 pages, 7 figures, 2 tables. Published as part of the Special Issue "18 Years of Science with the Neil Gehrels Swift Observatory's Ultra-Violet/Optical Telescope"

We report observations of the optical counterpart of the long gamma-ray burst (LGRB) GRB 221009A. Due to the extreme rarity of being both nearby ($z = 0.151$) and highly energetic ($E_{\gamma,\mathrm{iso}} \geq 10^{54}$ erg), GRB 221009A offers a unique opportunity to probe the connection between massive star core collapse and relativistic jet formation across a very broad range of $\gamma$-ray properties. Adopting a phenomenological power-law model for the afterglow and host galaxy estimates from high-resolution Hubble Space Telescope imaging, we use Bayesian model comparison techniques to determine the likelihood of an associated SN contributing excess flux to the optical light curve. Though not conclusive, we find moderate evidence ($K_{\rm{Bayes}}=10^{1.2}$) for the presence of an additional component arising from an associated supernova, SN 2022xiw, and find that it must be substantially fainter ($<$ 67% as bright at the 99% confidence interval) than SN 1998bw. Given the large and uncertain line-of-sight extinction, we attempt to constrain the supernova parameters ($M_{\mathrm{Ni}}$, $M_{\mathrm{ej}}$, and $E_{\mathrm{KE}}$) under several different assumptions with respect to the host galaxy's extinction. We find properties that are broadly consistent with previous GRB-associated SNe: $M_{\rm{Ni}}=0.05$ - $0.25 \, \rm{M_\odot}$, $M_{\rm{ej}}=3.5$ - $11.1 \, \rm{M_\odot}$, and $E_{\rm{KE}} = (1.6$ - $5.2) \times 10^{52} \, \rm{erg}$. We note that these properties are weakly constrained due to the faintness of the supernova with respect to the afterglow and host emission, but we do find a robust upper limit on the $M_{\rm{Ni}}$ of $M_{\rm{Ni}}<0.36\, \rm{M_\odot}$. Given the tremendous range in isotropic gamma-ray energy release exhibited by GRBs (7 orders of magnitude), the SN emission appears to be decoupled from the central engine in these systems., Comment: 18 pages, accepted to ApJL, 4 tables, 5 figures. Updated abstract in Preview

On August 17, 2017, Advanced LIGO and Virgo observed GW170817, the first gravitational-wave (GW) signal from a binary neutron star merger. It was followed by a short-duration gamma-ray burst, GRB 170817A, and by a non-thermal afterglow emission. In this work, a combined simultaneous fit of the electromagnetic (EM, specifically, afterglow) and GW domains is implemented, both using the posterior distribution of a GW standalone analysis as prior distribution to separately process the EM data, and fitting the EM and GW domains simultaneously. These approaches coincide mathematically, as long as the actual posterior of the GW analysis, and not an approximation, is used as prior for the EM analysis. We treat the viewing angle, $\theta_v$, as shared parameter across the two domains. In the afterglow modelling with a Gaussian structured jet this parameter and the jet core angle, $\theta_c$, are correlated, leading to high uncertainties on their values. The joint EM+GW analysis relaxes this degeneracy, reducing the uncertainty compared to an EM-only fit. We also apply our methodology to hypothetical GW170817-like events occurring in the next GW observing run at $\sim$140 and 70 Mpc. At 70 Mpc the existing EM degeneracy is broken, thanks to the inclusion of the GW domain in the analysis. At 140 Mpc, the EM-only fit cannot constrain $\theta_v$ nor $\theta_c$ because of the lack of detections in the afterglow rising phase. Folding the GW data into the analysis leads to tighter constraints on $\theta_v$, still leaving $\theta_c$ unconstrained, requiring instruments with higher sensitivities, such as Athena., Comment: 15 pages, 10 figures. Accepted for publication in MNRAS

FRB 180916.J0158+65 is one of the nearest, periodically repeating, and actively bursting fast radio burst (FRB) which has been localized to the outskirts of a spiral galaxy. In this work we study the FRB with the hard X-ray $14-195$ keV data from the Burst Alert Telescope (BAT) on board The Neil Gehrels Swift Observatory. BAT uses coded mask technology giving a localization of $\lesssim 3$ arc-minute in the hard X-ray band, along with an accurate background estimation. BAT has been observing the source location in survey mode since February 2020. The survey mode observations involves background subtracted spectra, integrated over a time span ranging $300-2000$ seconds, at the source location (from Feb 2020-Jan 2022). We analyzed all the $\sim 230$ survey mode observations from BAT and checked for any signal in any of the observations. We did not detect any signal at $>5\sigma$ confidence level in any of the observations. We could estimate a $5\sigma$ upper limit on the $14-195$ keV flux, which ranged between $4.5\times 10^{-10} - 7.6\times 10^{-9}\, \rm erg\, cm^{-2}\, s^{-1}$. At the source distance this relates to a $5\sigma$ upper limit on luminosity of $5.08\times 10^{44}- 8.5\times 10^{45} \rm erg\, s^{-1}$. With this estimate, we could rule out any persistent X-ray emission, at the source location for these snapshots of BAT observations., Comment: Accepted for publication in ApJ

FRB180301 is one of the most actively repeating fast radio bursts (FRBs) which has shown polarization angle changes in its radio burst emission, an indication for their likely origin in the magnetosphere of a highly-magnetized neutron star. We carried out a multi-wavelength campaign with the FAST radio telescope and NICER X-ray observatory to investigate any possible X-ray emission temporally coincident with the bright radio bursts. The observations took place on 2021 March 4, 9 and 19. We detected five bright radio bursts with FAST, four of which were strictly simultaneous with the NICER observations. The peak flux-density of the radio bursts ranged between $28-105$ mJy, the burst fluence between $27-170$ mJy-ms, and the burst durations between $1.7-12.3$ ms. The radio bursts from FRB~180301 exhibited complex time domain structure, and sub-pulses were detected in individual bursts, with no significant circular polarisation. The linear degree of polarisation in L-band reduced significantly compared to the 2019 observations. We do not detect any X-ray emission in excess of the background during the 5ms, 10ms, 100ms, 1sec and 100sec time intervals at/around the radio-burst barycenter-corrected arrival times, at a $>5\sigma$ confidence level. The $5\sigma$ upper limits on the X-ray a) persistent flux is $<7.64\times 10^{-12}\, \rm erg\, cm^{-2}\, s^{-1}$ , equivalent to $L_{\rm X}<2.50 \times 10^{45} \rm erg\, s^{-1}$ and b) 5 ms fluence is $<2\times 10^{-11} \rm erg\, cm^{-2}$, at the radio burst regions. Using the $5$ ms X-ray fluence upper limit, we can estimate the radio efficiency $\eta_{R/X} \equiv L_{\rm Radio}/L_{\rm X-ray} \gtrsim 10^{-8}$. The derived upper limit on $\eta_{R/X}$ is consistent with both magnetospheric models and synchrotron maser models involving relativistic shocks., Comment: Resubmitted to ApJ following referee's minor comments

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources., Comment: accepted in ApJS; follow-up paper to 1902.10045

The workshop Dawn VI: Next Generation Observatories took place online over three days, 5-7 October, 2021. More than 200 physicists and astronomers attended to contribute to, and learn from, a discussion of next-generation ground-based gravitational-wave detectors. The program was centered on the next generation of ground-based gravitational-wave observatories and their synergy with the greater landscape of scientific observatories of the 2030s. Cosmic Explorer (CE), a concept developed with US National Science Foundation support, was a particular focus; Einstein Telescope (ET), the European next generation concept, is an important complement and partner in forming a network. The concluding summary of the meeting expressed the sentiment that the observational science accessible to CE and ET, also in combination with data from other non-GW observatories, will stimulate a very broad community of analysts and yield insights which are exciting given the access to GWs from the entire universe. The need, and desire, for closer collaboration between ET and CE was expressed; a three-detector network is optimal for delivering much of the science. The science opportunities afforded by CE and ET are broad and compelling, impacting a wide range of disciplines in physics and high energy astrophysics. There was a consensus that CE is a concept that can deliver the promised science. A strong endorsement of Cosmic Explorer, as described in the CE Horizon Study, is a primary outcome of DAWN VI., Comment: Proceedings of the Dawn VI workshop

Kilonovae are ultraviolet, optical, and infrared transients powered by the radioactive decay of heavy elements following a neutron star merger. Joint observations of kilonovae and gravitational waves can offer key constraints on the source of Galactic $r$-process enrichment, among other astrophysical topics. However, robust constraints on heavy element production requires rapid kilonova detection (within $\sim 1$ day of merger) as well as multi-wavelength observations across multiple epochs. In this study, we quantify the ability of 13 wide field-of-view instruments to detect kilonovae, leveraging a large grid of over 900 radiative transfer simulations with 54 viewing angles per simulation. We consider both current and upcoming instruments, collectively spanning the full kilonova spectrum. The Roman Space Telescope has the highest redshift reach of any instrument in the study, observing kilonovae out to $z \sim 1$ within the first day post-merger. We demonstrate that BlackGEM, DECam, GOTO, the Vera C. Rubin Observatory's LSST, ULTRASAT, and VISTA can observe some kilonovae out to $z \sim 0.1$ ($\sim$475 Mpc), while DDOTI, MeerLICHT, PRIME, $Swift$/UVOT, and ZTF are confined to more nearby observations. Furthermore, we provide a framework to infer kilonova ejecta properties following non-detections and explore variation in detectability with these ejecta parameters., Comment: 22 pages, 8 figures, 1 appendix (with 9 pages, 10 figures), submitted to ApJ

Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the universe. The duration and hardness distribution of GRBs has two clusters, now understood to reflect (at least) two different progenitors. Short-hard GRBs (SGRBs; T90 <2 s) arise from compact binary mergers, while long-soft GRBs (LGRBs; T90 >2 s) have been attributed to the collapse of peculiar massive stars (collapsars). The discovery of SN 1998bw/GRB 980425 marked the first association of a LGRB with a collapsar and AT 2017gfo/GRB 170817A/GW170817 marked the first association of a SGRB with a binary neutron star merger, producing also gravitational wave (GW). Here, we present the discovery of ZTF20abwysqy (AT2020scz), a fast-fading optical transient in the Fermi Satellite and the InterPlanetary Network (IPN) localization regions of GRB 200826A; X-ray and radio emission further confirm that this is the afterglow. Follow-up imaging (at rest-frame 16.5 days) reveals excess emission above the afterglow that cannot be explained as an underlying kilonova (KN), but is consistent with being the supernova (SN). Despite the GRB duration being short (rest-frame T90 of 0.65 s), our panchromatic follow-up data confirms a collapsar origin. GRB 200826A is the shortest LGRB found with an associated collapsar; it appears to sit on the brink between a successful and a failed collapsar. Our discovery is consistent with the hypothesis that most collapsars fail to produce ultra-relativistic jets., Comment: submitted to Nature Astronomy

We use a sample of 27 GRBs at redshift $z=2-6$ to probe the outflows in their respective host galaxies ($\mathrm{log(M_*/M_{\odot})}~\sim~9-11$) and search for possible relations between the outflow properties and those of the host galaxies such as $\mathrm{M_*}$, SFR, and specific SFR. First, we consider three outflow properties $-$ outflow column density ($\mathrm{N_{out}}$), maximum outflow velocity ($\mathrm{V_{max}}$), and normalized maximum velocity ($\mathrm{V_{norm}}$ = $\mathrm{V_{max}/V_{circ, halo}}$, where $\mathrm{V_{circ,halo}}$ is the halo circular velocity). We observe clear trends of $\mathrm{N_{out}}$ and $\mathrm{V_{max}}$ with increasing SFR in high-ion-traced outflows, with a stronger ($>~3\sigma$) $\mathrm{V_{max}}-$SFR correlation. We find that the estimated mass outflow rate and momentum flux of the high-ion outflows scale with SFR and can be supported by the momentum imparted by star formation (supernovae and stellar winds). The kinematic correlations of high-ion-traced outflows with SFR are similar to those observed for star-forming galaxies at low redshifts. The correlations with SFR are weaker in low-ions. This, along with the lower detection fraction in low-ions, indicates that the outflow is primarily high-ion dominated. We also observe a strong ($>~3\sigma$) trend of normalized velocity ($\mathrm{V_{norm}}$) decreasing with halo mass and increasing with sSFR, suggesting that outflows from low-mass halos and high-sSFR galaxies are most likely to escape and enrich the outer CGM and IGM with metals. By comparing the CGM-GRB stacks with those of starbursts at $z\sim2$ and $z\sim0.1$, we find that over a broad redshift range, the outflow strength strongly depends on the main-sequence offset at the respective redshifts rather than simply the SFR., Comment: Published in Astrophysical Journal, 28 pages, 20 figures. Comments are welcome

The transition from prompt to the afterglow emission is one of the most exciting and least understood phases in gamma-ray bursts (GRBs). Correlations among optical, X-ray and gamma-ray emission in GRBs have been explored, to attempt to answer whether the earliest optical emission comes from internal and/or external shocks. We present optical photometric observations of GRB 180325A collected with the TAROT and RATIR ground-based telescopes. These observations show two strong optical flashes with separate peaks at $\sim50\;$s and $\sim120\;$s, followed by a temporally extended optical emission. We also present X-rays and gamma-ray observations of GRB 180325A, detected by the Burst Alert Telescope (BAT) and X-ray Telescope (XRT), on the Neil Gehrels Swift observatory, which both observed a narrow flash at $\sim80\;$s. We show that the prompt gamma- and X-ray early emission shares similar temporal and spectral features consistent with internal dissipation within the relativistic outflow (e.g. by internal shocks or magnetic reconnection), while the early optical flashes are likely generated by the reverse shock that decelerates the ejecta as it sweeps up the external medium., Comment: Accepted for publication in ApJ

Neutron star mergers produce a substantial amount of fast-moving ejecta, expanding outwardly for years after the merger. The interaction of these ejecta with the surrounding medium may produce a weak isotropic radio remnant, detectable in relatively nearby events. We use late-time radio observations of short duration gamma-ray bursts (sGRBs) to constrain this model. Two samples of events were studied: four sGRBs that are possibly in the local (<200 Mpc) universe were selected to constrain the remnant non-thermal emission from the sub-relativistic ejecta, whereas 17 sGRBs at cosmological distances were used to constrain the presence of a proto-magnetar central engine, possibly re-energezing the merger ejecta. We consider the case of GRB~170817A/GW170817, and find that in this case the early radio emission may be quenched by the jet blast-wave. In all cases, for ejecta mass range of M_ej \lesssim 10^{-2} (5 * 10^{-2}) M_sun, we can rule out very energetic merger ejecta E_ej \gtrsim 5 * 10^{52}(10^{53}) erg, thus excluding the presence of a powerful magnetar as a merger remnant., Comment: 13 pages, 8 figures, 3 tables. Accepted for pubblication in MNRAS

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

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., Comment: 13 pages, 2 figures, 2 tables, submitted to ApJ

Gamma-ray bursts (GRBs) associated with gravitational wave events are, and will likely continue to be, viewed at a larger inclination than GRBs without gravitational wave detections. As demonstrated by the afterglow of GW170817, this requires an extension of the common GRB afterglow models which typically assume emission from an on-axis top hat jet. We present a characterization of the afterglows arising from structured jets, providing a framework covering both successful and choked jets. We compute new closure relations for structured jets and compare with the established relations for energy injection and refreshed shock models. The temporal slope before the jet break is found to be a simple function of the ratio between the viewing angle and effective opening angle of the jet. A numerical model to calculate synthetic light curves and spectra is publicly available as the open source Python package afterglowpy., Comment: 29 pages, 13 figures, submitted to ApJ

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

We present 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\%., Comment: Accepted por publication in ApJ

Using deep and high-cadence gamma-ray burst (GRB) afterglow data from RATIR, we observe a sharp and achromatic light curve break 12.6 days after the GRB, accompanied by an approximately achromatic bump. Fitting of the optical, NIR, and X-ray data suggest a very narrow (2 degree) jet which remains collimated at late-time. We argue that the sharp light curve bump suggests an edge brightened jet, perhaps emitting only during a brief period of lateral jet expansion. The lightcurve also exhibits a gradual spectral evolution lasting $>10$ days. The evolution of the flux can be modeled as $\textrm{Flux} \sim \big(\frac{t}{[20 \textrm{days}]}\big)^\alpha \big(\frac{\lambda}{[800 \textrm{nm}]}\big)^\beta$, with a temporal slope $\alpha=-0.956 \pm 0.003$ and a gradually time-varying spectral slope $\beta = (0.60 \pm 0.07)+(0.26 \pm 0.06) \textrm{log}\big(\frac{t}{20 \rm{days}}\big)$., Comment: 10 pages, 3 figures, accepted to ApJ

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.

Wide-field optical surveys have begun to uncover large samples of fast (t_rise < 5d), luminous (M_peak < -18), blue transients. While commonly attributed to the breakout of a supernova shock into a dense wind, the great distances to the transients of this class found so far have hampered detailed investigation of their properties. We present photometry and spectroscopy from a comprehensive worldwide campaign to observe AT2018cow (ATLAS18qqn), the first fast-luminous optical transient to be found in real time at low redshift. Our first spectra (<2 days after discovery) are entirely featureless. A very broad absorption feature suggestive of near-relativistic velocities develops between 3-8 days, then disappears. Broad emission features of H and He develop after >10 days. The spectrum remains extremely hot throughout its evolution, and the photospheric radius contracts with time (receding below R<10^14 cm after 1 month). This behaviour does not match that of any known supernova, although a relativistic jet within a fallback supernova could explain some of the observed features. Alternatively, the transient could originate from the disruption of a star by an intermediate-mass black hole, although this would require long-lasting emission of highly super-Eddington thermal radiation. In either case, AT2018cow suggests that the population of fast luminous transients represents a new class of astrophysical event. Intensive follow-up of this event in its late phases, and of any future events found at comparable distance, will be essential to better constrain their origins., Comment: Corrected Figure 8 / Table 4 to use final fits. Includes machine-readable photometry table (hopefully for real this time)

IC443 is a supernova remnant located in a quite complex environment since it interacts with nearby clouds. Indications for the presence of overionized plasma have been found though the possible physical causes of overionization are still debated. Moreover, because of its peculiar position and proper motion, it is not clear if the pulsar wind nebula (PWN) within the remnant is the relic of the IC443 progenitor star or just a rambling one seen in projection on the remnant. Here we address the study of IC443 plasma in order to clarify the relationship PWN-remnant, the presence of overionization and the origin of the latter. We analyzed two \emph{XMM-Newton} observations producing background-subtracted, vignetting-corrected and mosaicked images in two different energy bands and we performed a spatially resolved spectral analysis of the X-ray emission. We identified an elongated (jet-like) structure with Mg-rich plasma in overionization. The head of the jet is interacting with a molecular cloud and the jet is aligned with the position of the PWN at the instant of the supernova explosion. Interestingly, the direction of the jet of ejecta is somehow consistent with the direction of the PWN jet. Our discovery of a jet of ejecta in IC443 enlarge the sample of core-collapse SNRs with collimated ejecta structures. IC443's jet is the first one which shows overionized plasma, possibly associated with the adiabatic expansion of ejecta. The match between the jet's direction and the original position of the PWN strongly supports the association between the neutron star and IC443., Comment: Accepted for publication in A&A

Recently, the optical counterpart of a gravitational wave source GW170817 has been identified in NGC 4993 galaxy. Together with evidence from observations in electromagnetic waves, the event has been suggested as a result of a merger of two neutron stars. We analyze the multi-wavelength data to characterize the host galaxy property and its distance to examine if the properties of NGC 4993 are consistent with this picture. Our analysis shows that NGC 4993 is a bulge-dominated galaxy with reff ~ 2-3 kpc and the Sersic index of n = 3-4 for the bulge component. The spectral energy distribution from 0.15 to 24 micron indicates that this galaxy has no significant ongoing star formation, the mean stellar mass of (0.3 - 1.2) times 10^11 Msun,the mean stellar age greater than ~3 Gyr, and the metallicity of about 20% to 100% of solar abundance. Optical images reveal dust lanes and extended features that suggest a past merging activity. Overall, NGC 4993 has characteristics of normal, but slightly disturbed elliptical galaxies. Furthermore, we derive the distance to NGC 4993 with the fundamental plane relation using 17 parameter sets of 7 different filters and the central stellar velocity dispersion from literature, finding an angular diameter distance of 37.7 +- 8.7 Mpc. NGC 4993 is similar to some host galaxies of short gamma-ray bursts but much different from those of long gamma-ray bursts, supporting the picture of GW170817 as a result of a merger of two NSs., Comment: 7 pages, 2 figures, 3 tables, accepted for publication in ApJL

Recently we have witnessed the first multi-messenger detection of colliding neutron stars through Gravitational Waves (GWs) and Electromagnetic (EM) waves (GW170817), thanks to the joint efforts of LIGO/Virgo and Space/Ground-based telescopes. In this paper, we report on the RATIR followup observation strategies and show the results for the trigger G194575. This trigger is not of astrophysical interest; however, is of great interests to the robust design of a followup engine to explore large sky error regions. We discuss the development of an image-subtraction pipeline for the 6-color, optical/NIR imaging camera RATIR. Considering a two band ($i$ and $r$) campaign in the Fall of 2015, we find that the requirement of simultaneous detection in both bands leads to a factor $\sim$10 reduction in false alarm rate, which can be further reduced using additional bands. We also show that the performance of our proposed algorithm is robust to fluctuating observing conditions, maintaining a low false alarm rate with a modest decrease in system efficiency that can be overcome utilizing repeat visits. Expanding our pipeline to search for either optical or NIR detections (3 or more bands), considering separately the optical $riZ$ and NIR $YJH$ bands, should result in a false alarm rate $\approx 1\%$ and an efficiency $\approx 90\%$. RATIR's simultaneous optical/NIR observations are expected to yield about one candidate transient in the vast 100 $\mathrm{deg^2}$ LIGO error region for prioritized followup with larger aperture telescopes., Comment: 10 pages, 8 figures, 1 table, Accepted to AJ

We report on the first open-use based Atacama Large Millimeter/submm Array (ALMA) 345-GHz observation for the late afterglow phase of GRB131030A. The ALMA observation constrained a deep limit at 17.1 d for the afterglow and host galaxy. We also identified a faint submillimeter source (ALMAJ2300-0522) near the GRB131030A position. The deep limit at 345 GHz and multifrequency observations obtained using {\it Swift} and RATIR yielded forward shock modeling with a two-dimensional relativistic hydrodynamic jet simulation and described X-ray excess in the afterglow. The excess was inconsistent with the synchrotron self-inverse Compton radiation from the forward shock. The host galaxy of GRB131030A and optical counterpart of ALMAJ2300-0522 were also identified in the SUBARU image. Based on the deep ALMA limit for the host galaxy, the 3-$\sigma$ upper limits of IR luminosity and the star formation rate (SFR) is estimated as $L_{IR}<1.11\times10^{11} L_{\odot}$ and SFR$<18.7$ ($M_{\odot}$~yr$^{-1}$), respectively. Although the separation angle from the burst location (3.5 arcsec) was rather large, ALMAJ2300-0522 may be one component of the GRB131030A host galaxy, according to previous host galaxy cases., Comment: 16 pages, 7 figures, 1 table, PASJ in press

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 ), Comment: 110 pages, 35 figures, accepted for publication on ApJS

DDOTI will be a wide-field robotic imager consisting of six 28-cm telescopes with prime focus CCDs mounted on a common equatorial mount. Each telescope will have a field of view of 12 square degrees, will have 2 arcsec pixels, and will reach a 10-sigma limiting magnitude in 60 seconds of r = 18.7 in dark time and r = 18.0 in bright time. The set of six will provide an instantaneous field of view of about 72 square degrees. DDOTI uses commercial components almost entirely. The first DDOTI will be installed at the Observatorio Astron\'omico Nacional in Sierra San Pedro Mart\'ir, Baja California, M\'exico in early 2017. The main science goals of DDOTI are the localization of the optical transients associated with GRBs detected by the GBM instrument on the Fermi satellite and with gravitational-wave transients. DDOTI will also be used for studies of AGN and YSO variability and to determine the occurrence of hot Jupiters. The principal advantage of DDOTI compared to other similar projects is cost: a single DDOTI installation costs only about US$500,000. This makes it possible to contemplate a global network of DDOTI installations. Such geographic diversity would give earlier access and a higher localization rate. We are actively exploring this option., Comment: To appear in the proceedings of SPIE conference 9910 "Observatory Operations: Strategies, Processes, and Systems VI". 12 pages

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., Comment: Submitted for publication in the Astrophysical Journal

To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays., Comment: Resubmitted to ApJS

We present a sample of 27 GRBs with detailed Swift light curves supplemented by late time Chandra observations. To answer the missing jet-break problem in general, we develop a numerical simulation based model which can be directly fit to the data using Monte Carlo methods. Our numerical model takes into account all the factors that can shape a jet-break: (i) lateral expansion (ii) edge effects and (iii) off-axis effects. Our results provide improved fits to the light curves and constraints on physical parameters. More importantly, our results suggest that off-axis effects are important and must be included in interpretations of GRB jet breaks., Comment: 11 pages, 4 figures, 1 table. ApJ in press. On-line figures (33 MB) can be found at https://psu.box.com/s/8o7vwenwkp7sidzebuywss9vvaqx69ha

GRB 130925A is a peculiar event characterized by an extremely long gamma-ray duration ($\approx$7 ks), as well as dramatic flaring in the X-rays for $\approx$20 ks. After this period, its X-ray afterglow shows an atypical soft spectrum with photon index $\Gamma$$\sim$4, as observed by Swift and Chandra, until $\approx 10^7$ s, when XMM-Newton observations uncover a harder spectral shape with $\Gamma$$\sim$2.5, commonly observed in GRB afterglows. We find that two distinct emission components are needed to explain the X-ray observations: a thermal component, which dominates the X-ray emission for several weeks, and a non-thermal component, consistent with a typical afterglow. A forward shock model well describes the broadband (from radio to X-rays) afterglow spectrum at various epochs. It requires an ambient medium with a very low density wind profile, consistent with that expected from a low-metallicity blue supergiant (BSG). The thermal component has a remarkably constant size and a total energy consistent with those expected by a hot cocoon surrounding the relativistic jet. We argue that the features observed in this GRB (its ultralong duration, the thermal cocoon, and the low density wind environment) are associated with a low metallicity BSG progenitor and, thus, should characterize the class of ultralong GRBs., Comment: 6 pgs, 3 figs, fig1 revised, ApJL in press

We report the serendipitous discovery of a bright point source flare in the Abell cluster 1795 with archival EUVE and Chandra observations. Assuming the EUVE emission is associated with the Chandra source, the X-ray 0.5-7 keV flux declined by a factor of ~2300 over a time span of 6 years, following a power-law decay with index ~2.44+-0.40. The Chandra data alone vary by a factor of ~20. The spectrum is well fit by a blackbody with a constant temperature of kT~0.09 keV (~10^6 K). The flare is spatially coincident with the nuclear region of a faint, inactive galaxy with a photometric redshift consistent at the one sigma level with the cluster (z=0.062476). We argue that these properties are indicative of a tidal disruption of a star by a black hole with log(M_BH/M_sun)~5.5+-0.5. If so, such a discovery indicates that tidal disruption flares may be used to probe black holes in the intermediate mass range, which are very difficult to study by other means., Comment: 51 pages (single column), 4 figures; Accepted for publication in ApJ

In recent years, detailed observations and accurate numerical simulations have provided support to the idea that mergers of compact binaries containing either two neutron stars (NSs) or an NS and a black hole (BH) may constitute the central engine of short gamma-ray bursts (SGRBs). The merger of such compact binaries is expected to lead to the production of a spinning BH surrounded by an accreting torus. Several mechanisms can extract energy from this system and power the SGRBs. Here we connect observations and numerical simulations of compact binary mergers, and use the current sample of SGRBs with measured energies to constrain the mass of their powering tori. By comparing the masses of the tori with the results of fully general-relativistic simulations, we are able to infer the properties of the binary progenitors which yield SGRBs. By assuming a constant efficiency in converting torus mass into jet energy, epsilon_{jet}=10%, we find that most of the tori have masses smaller than 0.01M_{sun}, favoring "high-mass" binary NSs mergers, i.e., binaries with total masses >~1.5 the maximum mass of an isolated NS. This has important consequences for the gravitational-wave signals that may be detected in association with SGRBs, since "high-mass" systems do not form a long-lived hypermassive NS after the merger. While NS-BH systems cannot be excluded to be the engine of at least some of the SGRBs, the BH would need to have an initial spin of ~0.9, or higher., Comment: 6 pages, 3 figures. Extended the discussion on jet efficiency. Accepted for publication on The Astrophysical Journal Letters

We present broad-band observations of the afterglow and environment of the short GRB 111020A. An extensive X-ray light curve from Swift/XRT, XMM-Newton and Chandra, spanning ~100 seconds to 10 days after the burst, reveals a significant break at t~2 days with pre- and post-break decline rates of alphaX,1 ~ -0.78 and alphaX,2<-1.7, respectively. Interpreted as a jet break, we infer a collimated outflow with an opening angle of ~3-8 degrees. The resulting beaming-corrected gamma-ray (10-1000 keV band) and blastwave kinetic energies are (2-3)e48 erg and (0.3-2)e49 erg, respectively, with the range depending on the unknown redshift of the burst. We report a radio afterglow limit of <39 microJy (3-sigma) from EVLA observations which, along with our finding that vc24.4 mag at 18 hours after the burst, and reveal a potential host galaxy with i~24.3 mag. The sub-arcsecond localization from Chandra provides a precise offset of 0.80"+/-0.11" (1-sigma) from this galaxy corresponding to an offset of 5-7 kpc for z=0.5-1.5. We find a high excess neutral Hydrogen column density of (7.5+/-2.0)e21 cm^(-2) (z=0). Our observations demonstrate that a growing fraction of short GRBs are collimated which may lead to a true event rate of >100-1000 Gpc^(-3) yr^(-1), in good agreement with the NS-NS merger rate of ~200-3000 Gpc^(-3) yr^(-1). This consistency is promising for coincident short GRB-gravitational wave searches in the forthcoming era of Advanced LIGO/VIRGO., Comment: Submitted to ApJ; emulateapj style; 20 pages, 5 tables, 6 figures

We carried out a systematic search of precursors on the sample of short GRBs observed by Swift. We found that ~8-10% of short GRBs display such early episode of emission. One burst (GRB 090510) shows two precursor events, the former ~13 s and the latter ~0.5 s before the GRB. We did not find any substantial difference between the precursor and the main GRB emission, and between short GRBs with and without precursors. We discuss possible mechanisms to reproduce the observed precursor emission within the scenario of compact object mergers. The implications of our results on quantum gravity constraints are also discussed., Comment: 8 pages, 2 figures, accepted for publication in The Astrophysical Journal