Showing 1–50 of 86 results for author: Beniamini, P

Magnetospheric origin of a fast radio burst constrained using scintillation

Authors: Kenzie Nimmo, Ziggy Pleunis, Paz Beniamini, Pawan Kumar, Adam E. Lanman, D. Z. Li, Robert Main, Mawson W. Sammons, Shion Andrew, Mohit Bhardwaj, Shami Chatterjee, Alice P. Curtin, Emmanuel Fonseca, B. M. Gaensler, Ronniy C. Joseph, Zarif Kader, Victoria M. Kaspi, Mattias Lazda, Calvin Leung, Kiyoshi W. Masui, Ryan Mckinven, Daniele Michilli, Ayush Pandhi, Aaron B. Pearlman, Masoud Rafiei-Ravandi , et al. (4 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are micro-to-millisecond duration radio transients that originate mostly from extragalactic distances. The emission mechanism responsible for these high luminosity, short duration transients remains debated. The models are broadly grouped into two classes: physical processes that occur within close proximity to a central engine; and central engines that release energy whic… ▽ More Fast radio bursts (FRBs) are micro-to-millisecond duration radio transients that originate mostly from extragalactic distances. The emission mechanism responsible for these high luminosity, short duration transients remains debated. The models are broadly grouped into two classes: physical processes that occur within close proximity to a central engine; and central engines that release energy which moves to large radial distances and subsequently interacts with surrounding media producing radio waves. The expected emission region sizes are notably different between these two types of models. FRB emission size constraints can therefore be used to distinguish between these competing models and inform on the physics responsible. Here we present the measurement of two mutually coherent scintillation scales in the frequency spectrum of FRB 20221022A: one originating from a scattering screen located within the Milky Way, and the second originating from a scattering screen located within its host galaxy or local environment. We use the scattering media as an astrophysical lens to constrain the size of the lateral emission region, $R_{\star\mathrm{obs}} \lesssim 3\times10^{4}$ km. We find that this is inconsistent with the expected emission sizes for the large radial distance models, and is more naturally explained with an emission process that operates within or just beyond the magnetosphere of a central compact object. Recently, FRB 20221022A was found to exhibit an S-shaped polarisation angle swing, supporting a magnetospheric emission process. The scintillation results presented in this work independently support this conclusion, while highlighting scintillation as a useful tool in our understanding of FRB emission physics and progenitors. △ Less

Submitted 16 June, 2024; originally announced June 2024.

Comments: Submitted, comments welcome

X-ray Afterglow limits on the viewing angles of short gamma-ray bursts

Authors: Brendan O'Connor, Paz Beniamini, Ramandeep Gill

Abstract: The observed behavior of a short gamma-ray burst (sGRB) afterglow lightcurve can reveal the angular structure of the relativistic jet and constrain the observer's viewing angle $θ_\textrm{obs}$. Regardless of viewing angle, the afterglow emission is produced by the interaction between the relativistic jet and its surrounding environment. However, the observed deceleration time of the jet, and, the… ▽ More The observed behavior of a short gamma-ray burst (sGRB) afterglow lightcurve can reveal the angular structure of the relativistic jet and constrain the observer's viewing angle $θ_\textrm{obs}$. Regardless of viewing angle, the afterglow emission is produced by the interaction between the relativistic jet and its surrounding environment. However, the observed deceleration time of the jet, and, therefore, the time of the afterglow peak, depends on the observer's viewing angle. A larger viewing angle leads to a later peak of the afterglow and a lower flux at peak. We use the the earliest X-ray afterglow detections of 58 cosmological sGRBs detected with the Neil Gehrels Swift Observatory X-ray Telescope to set an upper limit on the ratio of the viewing angle $θ_\textrm{obs}$ to the jet's half-opening angle $θ_\textrm{c}$. We adopt a power-law angular jet structure in both energy $E(θ)\proptoθ^{-a}$ and Lorentz factor $Γ(θ)\proptoθ^{-b}$ beyond the core. For this structured jet scenario we find that either sGRBs are viewed within $θ_\textrm{obs}/θ_\textrm{c}<1$ or the initial Lorentz factor of material in their jet's core is extremely high ($Γ_0>500$). If we consider a tophat jet structure, we constrain 90% of our sample to be viewed within $θ_\textrm{obs}/θ_\textrm{c}<1.06$ and 1.15 for our canonical and conservative afterglow scenarios. For a subset of events with measurements of the jet break, under the assumption that they on-axis we can constrain $Γ_0θ_\textrm{c}\gtrsim 30$. This confirmation that cosmological sGRBs are viewed either on-axis or very close to their jet's core has significant implications for the nature of the prompt gamma-ray production mechanism and for the rate of future sGRB detections coincident with gravitational waves (GWs), implying that they are extremely rare. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS

Afterglow Linear Polarization Signatures from Shallow GRB Jets: Implications for Energetic GRBs

Authors: Gal Birenbaum, Ramandeep Gill, Omer Bromberg, Paz Beniamini, Jonathan Granot

Abstract: Gamma-ray bursts (GRBs) are powered by ultra-relativistic jets. The launching sites of these jets are surrounded by dense media, which the jets must cross before they can accelerate and release the high energy emission. Interaction with the medium leads to the formation of a mildly relativistic sheath around the jet resulting in an angular structures in the jet's asymptotic Lorentz factor and ener… ▽ More Gamma-ray bursts (GRBs) are powered by ultra-relativistic jets. The launching sites of these jets are surrounded by dense media, which the jets must cross before they can accelerate and release the high energy emission. Interaction with the medium leads to the formation of a mildly relativistic sheath around the jet resulting in an angular structures in the jet's asymptotic Lorentz factor and energy per solid angle, which modifies the afterglow emission. We build a semi-analytical tool to analyze the afterglow light curve and polarization signatures of jets observed from a wide range of viewing angles, and focus on ones with slowly declining energy profiles known as shallow jets. We find overall lower polarization compared to the classical top-hat jet model. We provide an analytical expression for the peak polarization degree as a function of the energy profile power-law index, magnetic field configuration and viewing angle, and show that it occurs near the light curve break time for all viewers. When applying our tool to GRB 221009A, suspected to originate from a shallow jet, we find that the suggested jet structures for this event agree with the upper limits placed on the afterglow polarization in the optical and X-ray bands. We also find that at early times the polarization levels may be significantly higher, allowing for a potential distinction between different jet structure models and possibly constraining the magnetization in both forward and reverse shocks at that stage. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 15 pages, 11 figures

Joint gravitational wave-short GRB detection of Binary Neutron Star mergers with existing and future facilities

Authors: Soumyadeep Bhattacharjee, Smaranika Banerjee, Varun Bhalerao, Paz Beniamini, Sukanta Bose, Kenta Hotokezaka, Archana Pai, Muhammed Saleem, Gaurav Waratkar

Abstract: We explore the joint detection prospects of short gamma-ray bursts (sGRBs) and their gravitational wave (GW) counterparts by the current and upcoming high-energy GRB and GW facilities from binary neutron star (BNS) mergers. We consider two GW detector networks: (1) A four-detector network comprising LIGO Hanford, Livingston, Virgo, and Kagra, (IGWN4) and (2) a future five-detector network includin… ▽ More We explore the joint detection prospects of short gamma-ray bursts (sGRBs) and their gravitational wave (GW) counterparts by the current and upcoming high-energy GRB and GW facilities from binary neutron star (BNS) mergers. We consider two GW detector networks: (1) A four-detector network comprising LIGO Hanford, Livingston, Virgo, and Kagra, (IGWN4) and (2) a future five-detector network including the same four detectors and LIGO India (IGWN5). For the sGRB detection, we consider existing satellites Fermi and Swift and the proposed all-sky satellite Daksha. Most of the events for the joint detection will be off-axis, hence, we consider a broad range of sGRB jet models predicting the off-axis emission. Also, to test the effect of the assumed sGRB luminosity function, we consider two different functions for one of the emission models. We find that for the different jet models, the joint sGRB and GW detection rates for Fermi and Swift with IGWN4 (IGWN5) lie within 0.07-0.62$\mathrm{\ yr^{-1}}$ (0.8-4.0$\mathrm{\ yr^{-1}}$) and 0.02-0.14$\mathrm{\ yr^{-1}}$ (0.15-1.0$\mathrm{\ yr^{-1}}$), respectively, when the BNS merger rate is taken to be 320$\mathrm{\ Gpc^{-3}~yr^{-1}}$. With Daksha, the rates increase to 0.2-1.3$\mathrm{\ yr^{-1}}$ (1.3-8.3$\mathrm{\ yr^{-1}}$), which is 2-9 times higher than the existing satellites. We show that such a mission with higher sensitivity will be ideal for detecting a higher number of fainter events observed off-axis or at a larger distance. Thus, Daksha will boost the joint detections of sGRB and GW, especially for the off-axis events. Finally, we find that our detection rates with optimal SNRs are conservative, and noise in GW detectors can increase the rates further. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: 9 pages, 6 figures, 2 tables; accepted for publication in MNRAS. The definitive version will be available on the journal page

Ultrafast compact binary mergers

Authors: Paz Beniamini, Tsvi Piran

Abstract: The duration of orbital decay induced by gravitational waves (GWs) is often the bottleneck of the evolutionary phases going from star formation to a merger. We show here that kicks imparted to the newly born compact object during the second collapse generically result in a GW merger time distribution behaving like $dN/d\log t \propto t^{2/7}$ at short durations, leading to ultrafast mergers. Namel… ▽ More The duration of orbital decay induced by gravitational waves (GWs) is often the bottleneck of the evolutionary phases going from star formation to a merger. We show here that kicks imparted to the newly born compact object during the second collapse generically result in a GW merger time distribution behaving like $dN/d\log t \propto t^{2/7}$ at short durations, leading to ultrafast mergers. Namely, a nonnegligible fraction of neutron star binaries, formed in this way, will merge on a time scale as short as 10 Myr, and a small fraction will merge even on a time scale less than 10 kyr. The results can be applied to different types of compact binaries. We discuss here the implications for binary neutron star mergers. These include unique short gamma-ray bursts (GRBs), eccentric and misaligned mergers, $r$-process enrichment in the very early Universe and in highly star-forming regions, and possible radio precursors. Interestingly, we conclude that among the few hundred short GRBs detected so far, a few must have formed via this ultrafast channel. △ Less

Submitted 2 April, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 19 pages, 8 figures. Accepted for publication in ApJ

Internal Shocks Hydrodynamics: the Collision of Two Cold Shells in Detail

Authors: Sk. Minhajur Rahaman, Jonathan Granot, Paz Beniamini

Abstract: Emission in many astrophysical transients originates from a shocked fluid. A central engine typically produces an outflow with varying speeds, leading to internal collisions within the outflow at finite distances from the source. Each such collision produces a pair of forward and reverse shocks with the two shocked regions separated by a contact discontinuity (CD). As a useful approximation, we co… ▽ More Emission in many astrophysical transients originates from a shocked fluid. A central engine typically produces an outflow with varying speeds, leading to internal collisions within the outflow at finite distances from the source. Each such collision produces a pair of forward and reverse shocks with the two shocked regions separated by a contact discontinuity (CD). As a useful approximation, we consider the head-on collision between two cold and uniform shells (a slower leading shell and a faster trailing shell) of finite radial width, and study the dynamics of shock propagation in planar geometry. We find significant differences between the forward and reverse shocks, in terms of their strength, internal energy production efficiency, and the time it takes for the shocks to sweep through the respective shells. We consider the subsequent propagation of rarefaction waves in the shocked regions and explore the cases where these waves can catch up with the shock fronts and thereby limit the internal energy dissipation. We demonstrate the importance of energy transfer from the trailing to leading shell through $pdV$ work across the CD. We outline the parameter space regions relevant for models of different transients,e.g., Gamma-ray burst (GRB) internal shock model, fast radio burst (FRB) blastwave model, Giant flare due to magnetars, and superluminous supernovae (SLSN) ejecta. We find that the reverse shock likely dominates the internal energy production for many astrophysical transients. △ Less

Submitted 22 December, 2023; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS

Ultra High Energy Cosmic Rays from Tidal Disruption Events

Authors: Tsvi Piran, Paz Beniamini

Abstract: The tidal disruption event AT2018hyz was a regular optically detected one with no special prompt features. However, it suddenly displayed a fast-rising radio flare almost three years after the disruption. The flare is most naturally interpreted as arising from an off-axis relativistic jet. We didn't see the jet at early times as its emission was relativistically beamed away from us. However, we co… ▽ More The tidal disruption event AT2018hyz was a regular optically detected one with no special prompt features. However, it suddenly displayed a fast-rising radio flare almost three years after the disruption. The flare is most naturally interpreted as arising from an off-axis relativistic jet. We didn't see the jet at early times as its emission was relativistically beamed away from us. However, we could see the radiation once the jet has slowed down due to interaction with the surrounding matter. Analysis of the radio data enabled estimates of the jet's kinetic energy and opening angle, as well as the conditions (size and magnetic field) within the radio-emitting region. We show here that such a jet satisfies the Hillas condition for the acceleration of UHECRs to the highest energies. We also show that the rate and total power of this event are consistent with the observed luminosity density of UHECRs. These results strongly support earlier suggestions that TDEs are the sources of UHECRs. △ Less

Submitted 4 December, 2023; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: Accepted for publication in JCAP. The revised version includes a rebuttle to Cendes et al., 2023 arXiv:2308.13595 and additional references to earlier work on the relation between UHECRs and TDEs

A Numerical Study of Relativistic Oblique Shock Reflection

Authors: Prasanta Bera, Jonathan Granot, Michael Rabinovich, Paz Beniamini

Abstract: Shocks are ubiquitous in astrophysical sources, many of which involve relativistic bulk motions, leading to the formation of relativistic shocks. Such relativistic shocks have so far been studied mainly in one dimension, for simplicity, but the complex nature of the relevant astrophysical flows often requires higher dimensional studies. Here we study the two-dimensional problem of the reflection o… ▽ More Shocks are ubiquitous in astrophysical sources, many of which involve relativistic bulk motions, leading to the formation of relativistic shocks. Such relativistic shocks have so far been studied mainly in one dimension, for simplicity, but the complex nature of the relevant astrophysical flows often requires higher dimensional studies. Here we study the two-dimensional problem of the reflection of a planer shock off of a wall for a general incidence angle and a cold unshocked medium. We use primarily relativistic hydrodynamic numerical simulations, and elaborately compare the results to an analytic treatment. The simulations are performed both in the rest frame S of the unshocked fluid, where the dimensionless proper speed of the singly shocked fluid is $u_1=Γ_1β_1$ and the shock incidence angle is $α_1$, and in the rest frame S$^\prime$ of the point P of intersection of the incident shock and the wall for regular reflection (RR). Good agreement is obtained between the simulations in these two frames and with the analytic solution. The establishment of a steady flow in frame S$^\prime$ is explored, along with the transition between the strong and weak shock RR solutions. The transition line between RR and Mach reflection (MR) is studied numerically in the $u_1$-$α_1$ plane and found to coincide with the analytic detachment/sonic line. The flow properties along the sonic line are investigated in detail focusing on how they vary between the Newtonian and relativistic limits. △ Less

Submitted 22 September, 2023; v1 submitted 20 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS

Constraining the long-lived supramassive neutron stars by magnetar boosted kilonovae

Authors: Hao Wang, Paz Beniamini, Dimitrios Giannios

Abstract: Kilonovae are optical transients following the merger of neutron star binaries, which are powered by the r-process heating of merger ejecta. However, if a merger remnant is a long-lived supramassive neutron star supported by its uniform rotation, it will inject energy into the ejecta through spindown power. The energy injection can boost the peak luminosity of a kilonova by many orders of magnitud… ▽ More Kilonovae are optical transients following the merger of neutron star binaries, which are powered by the r-process heating of merger ejecta. However, if a merger remnant is a long-lived supramassive neutron star supported by its uniform rotation, it will inject energy into the ejecta through spindown power. The energy injection can boost the peak luminosity of a kilonova by many orders of magnitudes, thus significantly increasing the detectable volume. Therefore, even if such events are only a small fraction of the kilonovae population, they could dominate the detection rates. However, after many years of optical sky surveys, no such event has been confirmed. In this work, we build a boosted kilonova model with rich physical details, including the description of the evolution and stability of a proto neutron star, and the energy absorption through X-ray photoionization. We simulate the observation prospects and find the only way to match the absence of detection is to limit the energy injection by the newly born magnetar to only a small fraction of the neutron star rotational energy, thus they should collapse soon after the merger. Our result indicates that most supramassive neutron stars resulting from binary neutron star mergers are short lived and they are likely to be rare in the Universe. △ Less

Submitted 24 November, 2023; v1 submitted 17 August, 2023; originally announced August 2023.

Comments: 17 pages, 10 figures. Accepted by MNRAS

An off-axis relativistic jet seen in the long lasting delayed radio flare of the TDE AT 2018hyz

Authors: Itai Sfaradi, Paz Beniamini, Assaf Horesh, Tsvi Piran, Joe Bright, Lauren Rhodes, David R. A. Willians, Rob Fender, James K. Leung, Tara Murphy, Dave A. Green

Abstract: The Tidal Disruption Event (TDE) AT 2018hyz exhibited a delayed radio flare almost three years after the stellar disruption. Here we report new radio observations of the TDE AT 2018hyz with the AMI-LA and ATCA spanning from a month to more than four years after the optical discovery and 200 days since the last reported radio observation. We detected no radio detection from 30-220 days after the op… ▽ More The Tidal Disruption Event (TDE) AT 2018hyz exhibited a delayed radio flare almost three years after the stellar disruption. Here we report new radio observations of the TDE AT 2018hyz with the AMI-LA and ATCA spanning from a month to more than four years after the optical discovery and 200 days since the last reported radio observation. We detected no radio detection from 30-220 days after the optical discovery in our observations at 15.5 GHz down to a $3σ$ level of < 0.14 mJy. The fast-rising, delayed, radio flare is observed in our radio data set and continues to rise almost ~1580 days after the optical discovery. We find that the delayed radio emission, first detected $972$ days after optical discovery, evolves as $t^{4.2 \pm 0.9}$, at 15.5 GHz. Here, we present an off-axis jet model that can explain the full set of radio observations. In the context of this model, we require a powerful narrow jet with an isotropic equivalent kinetic energy $E_{\rm k,iso} \sim 10^{55}$ erg, an opening angle of $ \rm \sim 7^{\circ}$, and a relatively large viewing angle of $ \rm \sim 42^{\circ}$, launched at the time of the stellar disruption. Within our framework, we find that the minimal collimated energy possible for an off-axis jet from AT 2018hyz is $E_k \geq 3 \times 10^{52}$ erg. Finally, we provide predictions based on our model for the light curve turnover time, and for the proper motion of the radio emitting source. △ Less

Submitted 3 August, 2023; originally announced August 2023.

Comments: 8 pages, 6 figures, submitted to MNRAS

Prompt Gamma-Ray Burst Emission from Internal Shocks -- New Insights

Authors: Sk. Minhajur Rahaman, Jonathan Granot, Paz Beniamini

Abstract: Internal shocks are a leading candidate for the dissipation mechanism that powers the prompt $γ$-ray emission in gamma-ray bursts (GRBs). In this scenario a compact central source produces an ultra-relativistic outflow with varying speeds, causing faster parts or shells to collide with slower ones. Each collision produces a pair of shocks -- a forward shock (FS) propagating into the slower leading… ▽ More Internal shocks are a leading candidate for the dissipation mechanism that powers the prompt $γ$-ray emission in gamma-ray bursts (GRBs). In this scenario a compact central source produces an ultra-relativistic outflow with varying speeds, causing faster parts or shells to collide with slower ones. Each collision produces a pair of shocks -- a forward shock (FS) propagating into the slower leading shell and a reverse shock (RS) propagating into the faster trailing shell. The RS's lab-frame speed is always smaller, while the RS is typically stronger than the FS, leading to different conditions in the two shocked regions that both contribute to the observed emission. We show that optically-thin synchrotron emission from both (weaker FS + stronger RS) can naturally explain key features of prompt GRB emission such as the pulse shapes, time-evolution of the $νF_ν$ peak flux and photon-energy, and the spectrum. Particularly, it can account for two features commonly observed in GRB spectra: (i) a sub-dominant low-energy spectral component (often interpreted as ``photospheric''-like), or (ii) a doubly-broken power-law spectrum with the low-energy spectral slope approaching the slow cooling limit. Both features can be obtained while maintaining high overall radiative efficiency without any fine-tuning of the physical conditions. △ Less

Submitted 22 December, 2023; v1 submitted 1 August, 2023; originally announced August 2023.

Comments: Published in MNRAS Letters

Constraining the FRB mechanism from scintillation in the host galaxy

Authors: Pawan Kumar, Paz Beniamini, Om Gupta, James M. Cordes

Abstract: Most FRB models can be divided into two groups based on the distance of the radio emission region from the central engine. The first group of models, the so-called `nearby' or magnetospheric models, invoke FRB emission at distances of 10$^9$ cm or less from the central engine, while the second `far-away' models involve emission from distances of 10$^{11}$ cm or greater. The lateral size for the em… ▽ More Most FRB models can be divided into two groups based on the distance of the radio emission region from the central engine. The first group of models, the so-called `nearby' or magnetospheric models, invoke FRB emission at distances of 10$^9$ cm or less from the central engine, while the second `far-away' models involve emission from distances of 10$^{11}$ cm or greater. The lateral size for the emission region for the former class of models ($\lesssim$ 10$^7$ cm) is much smaller than the second class of models ($\gtrsim 10^9$ cm). We propose that an interstellar scattering screen in the host galaxy is well-suited to differentiate between the two classes of models, particularly based on the level of modulations in the observed intensity with frequency, in the regime of strong diffractive scintillation. This is because the diffractive length scale for the host galaxy's ISM scattering screen is expected to lie between the transverse emission-region sizes for the `nearby' and the `far-away' class of models. Determining the strength of flux modulation caused by scintillation (scintillation modulation index) across the scintillation bandwidth ($\sim 1/2πδt_s$) would provide a strong constraint on the FRB radiation mechanism when the scatter broadening ($δt_s$) is shown to be from the FRB host galaxy. The scaling of the scintillation bandwidth as $\sim ν^{4.4}$ may make it easier to determine the modulation index at $\gtrsim$ 1 GHz. △ Less

Submitted 29 September, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: Paper accepted for publication in MNRAS

The Swift Deep Galactic Plane Survey (DGPS) Phase-I Catalog

Authors: B. O'Connor, C. Kouveliotou, P. A. Evans, N. Gorgone, A. J. van Kooten, S. Gagnon, H. Yang, M. G. Baring, E. Bellm, P. Beniamini, J. Brink, D. A. H. Buckley, S. B. Cenko, O. D. Egbo, E. Gogus, J. Granot, C. Hailey, J. Hare, F. Harrison, D. Hartmann, A. J. van der Horst, D. Huppenkothen, L. Kaper, O. Kargaltsev, J. A. Kennea , et al. (8 additional authors not shown)

Abstract: The \textit{Swift} Deep Galactic Plane Survey is a \textit{Swift} Key Project consisting of 380 tiled pointings covering 40 deg$^{2}$ of the Galactic Plane between longitude $10$\,$<$\,$|l|$\,$<$\,$30$ deg and latitude $|b|$\,$<$\,$0.5$ deg. Each pointing has a $5$ ks exposure, yielding a total of 1.9 Ms spread across the entire survey footprint. Phase-I observations were carried out between March… ▽ More The \textit{Swift} Deep Galactic Plane Survey is a \textit{Swift} Key Project consisting of 380 tiled pointings covering 40 deg$^{2}$ of the Galactic Plane between longitude $10$\,$<$\,$|l|$\,$<$\,$30$ deg and latitude $|b|$\,$<$\,$0.5$ deg. Each pointing has a $5$ ks exposure, yielding a total of 1.9 Ms spread across the entire survey footprint. Phase-I observations were carried out between March 2017 and May 2021. The Survey is complete to depth $L_X$\,$>$\,$10^{34}$ erg s$^{-1}$ to the edge of the Galaxy. The main Survey goal is to produce a rich sample of new X-ray sources and transients, while also covering a broad discovery space. Here, we introduce the Survey strategy and present a catalog of sources detected during Phase-I observations. In total, we identify 928 X-ray sources, of which 348 are unique to our X-ray catalog. We report on the characteristics of sources in our catalog and highlight sources newly classified and published by the DGPS team. △ Less

Submitted 21 September, 2023; v1 submitted 25 June, 2023; originally announced June 2023.

Comments: Accepted in ApJS. This is the final version

The Signature of Refreshed Shocks in the of Afterglow of GRB030329

Authors: Michael J. Moss, Robert Mochkovitch, Frédéric Daigne, Paz Beniamini, Sylvain Guiriec

Abstract: GRB030329 displays one clear and, possibly, multiple less intense fast-rising ($Δt / t \sim 0.3$) jumps in its optical afterglow light curve. The decay rate of the optical light curve remains the same before and after the photon flux jumps. This may be the signature of energy injection into the forward and reverse shocked material at the front of the jet. In this study, we model the Gamma-Ray Burs… ▽ More GRB030329 displays one clear and, possibly, multiple less intense fast-rising ($Δt / t \sim 0.3$) jumps in its optical afterglow light curve. The decay rate of the optical light curve remains the same before and after the photon flux jumps. This may be the signature of energy injection into the forward and reverse shocked material at the front of the jet. In this study, we model the Gamma-Ray Burst (GRB) ejecta as a series of shells of material. We follow the dynamical evolution of the ejecta as it interacts with itself (i.e., internal shocks) and with the circumburst medium (i.e., external forward and reverse shocks), and we calculate the emission from each shock event assuming synchrotron emission. We confirm the viability of the model proposed by \citet{2003Natur.426..138G} in which the jumps in the optical afterglow light curve of GRB030329 are produced via refreshed shocks. The refreshed shocks may be the signatures of the collisions between earlier ejected material with an average Lorentz factor $\barΓ\gtrsim 100$ and later ejected material with $\barΓ \sim 10$ once the early material has decelerated due to interaction with the circumburst medium. We show that even if the late material is ejected with a spread of Lorentz factors, internal shocks naturally produce a narrow distribution of Lorentz factors ($ΔΓ/Γ\lesssim0.1$), which is a necessary condition to produce the observed quick rise times of the jumps. These results imply a phase of internal shocks at some point in the dynamical evolution of the ejecta, which requires a low magnetization in the outflow. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Comments: 11 pages, 6 figures

Report number: MN-23-1927-MJ.R1

Journal ref: Mon. Not. Roy. Astron. Soc. stad2594 (2023)

Swift J1644+57 as an off-axis Jet

Authors: Paz Beniamini, Tsvi Piran, Tatsuya Matsumoto

Abstract: One of the intriguing puzzles concerning Swift J1644+57, the first jetted tidal disruption event (TDE) discovered, is the constant increase in its jet energy, as implied by radio observations. During the first two hundred days the jet energy has increased by an order of magnitude. We suggest that the jet was viewed slightly off-axis. In this case, the apparent energy increase arises due to the slo… ▽ More One of the intriguing puzzles concerning Swift J1644+57, the first jetted tidal disruption event (TDE) discovered, is the constant increase in its jet energy, as implied by radio observations. During the first two hundred days the jet energy has increased by an order of magnitude. We suggest that the jet was viewed slightly off-axis. In this case, the apparent energy increase arises due to the slowing down of the jet and the corresponding broadening of its beaming cone. Using equipartition analysis, we infer an increasing jet energy as a larger region of the jet is observed. A simple off-axis model accounts nicely for the multi-wavelength radio observations, resolving this long-standing puzzle. The model allows us to self-consistently evolve the synchrotron signature from an off-axis jet as a function of time. It also allows us to estimate, for the first time, the beaming angle of the jet, $θ_0 \approx 21^{\circ}$. Considering existing limits on the black hole mass, $\lesssim 10^{7}M_{\odot}$, this angle implies that the prompt phase beaming corrected luminosity of Swift J1644+57, $\sim 10^{47}$ergs/sec, was super Eddington. We also present a closure relation between the spectral and temporal flux for off-axis jets, which can be used to test whether a given radio transient is off-axis or not. △ Less

Submitted 12 July, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: 10 pages, 8 figures. Published in MNRAS

A structured jet explains the extreme GRB 221009A

Authors: B. O'Connor, E. Troja, G. Ryan, P. Beniamini, H. van Eerten, J. Granot, S. Dichiara, R. Ricci, V. Lipunov, J. H. Gillanders, R. Gill, M. Moss, S. Anand, I. Andreoni, R. L. Becerra, D. A. H. Buckley, N. R. Butler, S. B. Cenko, A. Chasovnikov, J. Durbak, C. Francile, E. Hammerstein, A. J. van der Horst, M. Kasliwal, C. Kouveliotou , et al. (7 additional authors not shown)

Abstract: Long duration gamma-ray bursts (GRBs) are powerful cosmic explosions, signaling the death of massive stars. Among them, GRB 221009A is by far the brightest burst ever observed. Due to its enormous energy ($E_\textrm{iso}\!\approx$10$^{55}$ erg) and proximity ($z\!\approx$0.15), GRB 221009A is an exceptionally rare event that pushes the limits of our theories. We present multi-wavelength observatio… ▽ More Long duration gamma-ray bursts (GRBs) are powerful cosmic explosions, signaling the death of massive stars. Among them, GRB 221009A is by far the brightest burst ever observed. Due to its enormous energy ($E_\textrm{iso}\!\approx$10$^{55}$ erg) and proximity ($z\!\approx$0.15), GRB 221009A is an exceptionally rare event that pushes the limits of our theories. We present multi-wavelength observations covering the first three months of its afterglow evolution. The X-ray brightness decays as a power-law with slope $\approx\!t^{-1.66}$, which is not consistent with standard predictions for jetted emission. We attribute this behavior to a shallow energy profile of the relativistic jet. A similar trend is observed in other energetic GRBs, suggesting that the most extreme explosions may be powered by structured jets launched by a common central engine. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: Submitted version. 53 pages, 9 figures, 6 tables

Science with the Daksha High Energy Transients Mission

Authors: Varun Bhalerao, Disha Sawant, Archana Pai, Shriharsh Tendulkar, Santosh Vadawale, Dipankar Bhattacharya, Vikram Rana, Hitesh Kumar L. Adalja, G C Anupama, Suman Bala, Smaranika Banerjee, Judhajeet Basu, Hrishikesh Belatikar, Paz Beniamini, Mahesh Bhaganagare, Ankush Bhaskar, Soumyadeep Bhattacharjee, Sukanta Bose, Brad Cenko, Mehul Vijay Chanda, Gulab Dewangan, Vishal Dixit, Anirban Dutta, Priyanka Gawade, Abhijeet Ghodgaonkar , et al. (50 additional authors not shown)

Abstract: We present the science case for the proposed Daksha high energy transients mission. Daksha will comprise of two satellites covering the entire sky from 1~keV to $>1$~MeV. The primary objectives of the mission are to discover and characterize electromagnetic counterparts to gravitational wave source; and to study Gamma Ray Bursts (GRBs). Daksha is a versatile all-sky monitor that can address a wide… ▽ More We present the science case for the proposed Daksha high energy transients mission. Daksha will comprise of two satellites covering the entire sky from 1~keV to $>1$~MeV. The primary objectives of the mission are to discover and characterize electromagnetic counterparts to gravitational wave source; and to study Gamma Ray Bursts (GRBs). Daksha is a versatile all-sky monitor that can address a wide variety of science cases. With its broadband spectral response, high sensitivity, and continuous all-sky coverage, it will discover fainter and rarer sources than any other existing or proposed mission. Daksha can make key strides in GRB research with polarization studies, prompt soft spectroscopy, and fine time-resolved spectral studies. Daksha will provide continuous monitoring of X-ray pulsars. It will detect magnetar outbursts and high energy counterparts to Fast Radio Bursts. Using Earth occultation to measure source fluxes, the two satellites together will obtain daily flux measurements of bright hard X-ray sources including active galactic nuclei, X-ray binaries, and slow transients like Novae. Correlation studies between the two satellites can be used to probe primordial black holes through lensing. Daksha will have a set of detectors continuously pointing towards the Sun, providing excellent hard X-ray monitoring data. Closer to home, the high sensitivity and time resolution of Daksha can be leveraged for the characterization of Terrestrial Gamma-ray Flashes. △ Less

Submitted 27 January, 2024; v1 submitted 22 November, 2022; originally announced November 2022.

Comments: 19 pages, 7 figures. Accepted in Experimental Astronomy. More details about the mission at https://www.dakshasat.in/

Hybrid Pulsar-Magnetar Model for FRB 20191221A

Authors: Paz Beniamini, Pawan Kumar

Abstract: We show that the 216.8$\pm$0.1 ms periodicity reported for the fast radio burst (FRB) 20191221A is very constraining for burst models. The high accuracy of burst periodicity (better than one part in 10$^3$), and the 2\% duty cycle (ratio of burst-duration and inter-burst interval), suggest a pulsar-like rotating beam model for the observed activity; the radio waves are produced along open field li… ▽ More We show that the 216.8$\pm$0.1 ms periodicity reported for the fast radio burst (FRB) 20191221A is very constraining for burst models. The high accuracy of burst periodicity (better than one part in 10$^3$), and the 2\% duty cycle (ratio of burst-duration and inter-burst interval), suggest a pulsar-like rotating beam model for the observed activity; the radio waves are produced along open field lines within $\sim 10^7$ cm of the neutron star surface, and the beam periodically sweeps across the observer as the star spins. According to this picture, FRB 20191221A is a factor $\sim 10^{12}$ scaled up version of galactic pulsars with one major difference: whereas pulsars convert rotational kinetic energy to EM waves, the outbursts of 20191221A require conversion of magnetic energy to radiation. △ Less

Submitted 5 January, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 7 pages, 3 figures. Accepted for publication in MNRAS

Constraints on Electron Acceleration in Gamma-Ray Bursts Afterglows from Radio Peaks

Authors: Ruby A. Duncan, Alexander J. van der Horst, Paz Beniamini

Abstract: Studies of gamma-ray bursts (GRBs) and their multi-wavelength afterglows have led to insights in electron acceleration and emission properties from relativistic, high-energy astrophysical sources. Broadband modeling across the electromagnetic spectrum has been the primary means of investigating the physics behind these sources, although independent diagnostic tools have been developed to inform an… ▽ More Studies of gamma-ray bursts (GRBs) and their multi-wavelength afterglows have led to insights in electron acceleration and emission properties from relativistic, high-energy astrophysical sources. Broadband modeling across the electromagnetic spectrum has been the primary means of investigating the physics behind these sources, although independent diagnostic tools have been developed to inform and corroborate assumptions made in particle acceleration simulations and broadband studies. We present a methodology to constrain three physical parameters related to electron acceleration in GRB blast waves: the fraction of shock energy in electrons, $ε_e$; the fraction of electrons that gets accelerated into a power-law distribution of energies, $ξ_e$; and the minimum Lorentz factor of the accelerated electrons, $γ_m$. These parameters are constrained by observations of the peaks in radio afterglow light curves and spectral energy distributions. From a sample of 49 radio afterglows, we are able to find narrow distributions for these parameters, hinting at possible universality of the blast wave microphysics, although observational bias could play a role in this. Using radio peaks and considerations related to the prompt gamma-ray emission efficiency, we constrain the allowed parameter ranges for both $ε_e$ and $ξ_e$ to within about one order of magnitude, $0.01\lesssimε_e\lesssim0.2$ and $0.1\lesssimξ_e\lesssim1$. Such stringent constraints are inaccessible for $ξ_e$ from broadband studies due to model degeneracies. △ Less

Submitted 1 November, 2022; originally announced November 2022.

Comments: 10 pages, 6 figures, 2 tables; accepted for publication in MNRAS

Evidence for an abundant old population of Galactic ultra long period magnetars and implications for fast radio bursts

Authors: P. Beniamini, Z. Wadiasingh, J. Hare, K. Rajwade, G. Younes, A. J. van der Horst

Abstract: Two recent discoveries, namely PSR J0901-4046 and GLEAM-X J162759.5-523504.3 (hereafter GLEAM-X J1627), have corroborated an extant population of radio-loud periodic sources with long periods (76 s and 1091 s respectively) whose emission can hardly be explained by rotation losses. We argue that GLEAM-X J1627 is a highly-magnetized object consistent with a magnetar (an ultra long period magnetar -… ▽ More Two recent discoveries, namely PSR J0901-4046 and GLEAM-X J162759.5-523504.3 (hereafter GLEAM-X J1627), have corroborated an extant population of radio-loud periodic sources with long periods (76 s and 1091 s respectively) whose emission can hardly be explained by rotation losses. We argue that GLEAM-X J1627 is a highly-magnetized object consistent with a magnetar (an ultra long period magnetar - ULPM), and demonstrate it is unlikely to be either a magnetically or a rotationally-powered white dwarf. By studying these sources together with previously detected objects, we find there are at least a handful of promising candidates for Galactic ULPMs. The detections of these objects imply a substantial number, $N \gtrsim 13000$ and $N \gtrsim 500$ for PSR J0901--4046 like and GLEAM-X J1627 like objects, respectively, within our Galaxy. These source densities, as well as cooling age limits from non-detection of thermal X-rays, Galactic offsets, timing stability and dipole spindown limits, all imply the ULPM candidates are substantially older than confirmed Galactic magnetars and that their formation channel is a common one. Their existence implies widespread survival of magnetar-like fields for several Myr, distinct from the inferred behaviour in confirmed Galactic magnetars. ULPMs may also constitute a second class of FRB progenitors which could naturally exhibit very long periodic activity windows. Finally, we show that existing radio campaigns are biased against detecting objects like these and discuss strategies for future radio and X-ray surveys to identify more such objects. We estimate that ${\cal O}(100)$ more such objects should be detected with SKA-MID and DSA-2000. △ Less

Submitted 16 February, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 22 pages, 10 figures. Published in MNRAS

Gravitational lensing in the presence of plasma scattering with application to Fast Radio Bursts

Authors: Pawan Kumar, Paz Beniamini

Abstract: We describe how gravitational lensing of fast radio bursts (FRBs) is affected by a plasma screen in the vicinity of the lens or somewhere between the source and the observer. Wave passage through a turbulent medium affects gravitational image magnification, lensing probability (particularly for strong magnification events), and the time delay between images. The magnification is suppressed because… ▽ More We describe how gravitational lensing of fast radio bursts (FRBs) is affected by a plasma screen in the vicinity of the lens or somewhere between the source and the observer. Wave passage through a turbulent medium affects gravitational image magnification, lensing probability (particularly for strong magnification events), and the time delay between images. The magnification is suppressed because of the broadening of the angular size of the source due to scattering by the plasma. The time delay between images is modified as the result of different dispersion measure (DM) along photon trajectories for different images. Each of the image lightcurve is also broadened due to wave scattering so that the images could have distinct temporal profiles. The first two effects are most severe for stellar and sub-stellar mass lens, and the last one (scatter broadening) for lenses and plasma screens at cosmological distances from the source/observer. This could limit the use of FRBs to measure their cosmic abundance. On the other hand, when the time delay between images is large, such that the lightcurve of a transient source has two or more well separated peaks, the different DMs along the wave paths of different images can probe density fluctuations in the IGM on scales $\lesssim 10^{-6}$ rad and explore the patchy reionization history of the universe using lensed FRBs at high redshifts. Different rotation measure (RM) along two image paths can convert linearly polarized radiation from a source to partial circular polarization. △ Less

Submitted 17 August, 2022; v1 submitted 5 August, 2022; originally announced August 2022.

Comments: 11 pages, 5 figures

A deep survey of short GRB host galaxies over $z\sim0-2$: implications for offsets, redshifts, and environments

Authors: B. O'Connor, E. Troja, S. Dichiara, P. Beniamini, S. B. Cenko, C. Kouveliotou, J. Becerra Gonzalez, J. Durbak, P. Gatkine, A. Kutyrev, T. Sakamoto, R. Sanchez-Ramirez, S. Veilleux

Abstract: A significant fraction (30\%) of well-localized short gamma-ray bursts (sGRBs) lack a coincident host galaxy. This leads to two main scenarios: \textit{i}) that the progenitor system merged outside of the visible light of its host, or \textit{ii}) that the sGRB resided within a faint and distant galaxy that was not detected by follow-up observations. Discriminating between these scenarios has impo… ▽ More A significant fraction (30\%) of well-localized short gamma-ray bursts (sGRBs) lack a coincident host galaxy. This leads to two main scenarios: \textit{i}) that the progenitor system merged outside of the visible light of its host, or \textit{ii}) that the sGRB resided within a faint and distant galaxy that was not detected by follow-up observations. Discriminating between these scenarios has important implications for constraining the formation channels of neutron star mergers, the rate and environments of gravitational wave sources, and the production of heavy elements in the Universe. In this work, we present the results of our observing campaign targeted at 31 sGRBs that lack a putative host galaxy. Our study effectively doubles the sample of well-studied sGRB host galaxies, now totaling 72 events of which $28\%$ lack a coincident host to deep limits ($r$\,$\gtrsim$\,$26$ or $F110W$\,$\gtrsim$\,$27$ AB mag), and represents the largest homogeneously selected catalog of sGRB offsets to date. We find that 70\% of sub-arcsecond localized sGRBs occur within 10 kpc of their host's nucleus, with a median projected physical offset of $5.6$ kpc. Using this larger population, we discover an apparent redshift evolution in their locations: bursts at low-$z$ occur at $2\times$ larger offsets compared to those at $z$\,$>$\,$0.5$. This evolution could be due to a physical evolution of the host galaxies themselves or a bias against faint high-$z$ galaxies. Furthermore, we discover a sample of hostless sGRBs at $z$\,$\gtrsim$\,$1$ that are indicative of a larger high-$z$ population, constraining the redshift distribution and disfavoring log-normal delay time models. △ Less

Submitted 11 July, 2022; v1 submitted 19 April, 2022; originally announced April 2022.

Comments: Accepted by MNRAS

Robust Features of Off-Axis Gamma-Ray Burst Afterglow Lightcurves

Authors: Paz Beniamini, Ramandeep Gill, Jonathan Granot

Abstract: The ultra-relativistic outflows powering gamma-ray bursts (GRBs) acquire angular structure through their interaction with external material. They are often characterized by a compact, nearly uniform narrow core (with half-opening angle $θ_{c,\{ε,Γ\}}$) surrounded by material with energy per unit solid angle ($ε=ε_cΘ_ε^{-a}$, where $Θ_{\{ε,Γ\}}=[1+θ^2/θ_{c,\{ε,Γ\}}^2]^{1/2}$) and initial specific k… ▽ More The ultra-relativistic outflows powering gamma-ray bursts (GRBs) acquire angular structure through their interaction with external material. They are often characterized by a compact, nearly uniform narrow core (with half-opening angle $θ_{c,\{ε,Γ\}}$) surrounded by material with energy per unit solid angle ($ε=ε_cΘ_ε^{-a}$, where $Θ_{\{ε,Γ\}}=[1+θ^2/θ_{c,\{ε,Γ\}}^2]^{1/2}$) and initial specific kinetic energy ($Γ_0-1=[Γ_c-1]Θ_Γ^{-b}$) declining as power laws. Multi-wavelength afterglow lightcurves of off-axis jets (with viewing angle $θ_{\rm obs} > θ_c$) offer robust ways to constrain $a$, $b$ and the external density radial profile ($ρ\propto R^{-k}$), even while other burst parameters may remain highly degenerate. We extend our previous work on such afterglows to include more realistic angular structure profiles derived from three-dimensional hydrodynamic simulations of both long and short GRBs (addressing also jets with shallow angular energy profiles, whose emission exhibits unique evolution). We present afterglow lightcurves based on our parameterized power-law jet angular profiles for different viewing angles $θ_{\rm obs}$ and $k=\{0,1,2\}$. We identify a unique evolutionary power-law phase of the characteristic synchrotron frequencies ($ν_m$ and $ν_c$) that manifests when the lightcurve is dominated by emission sensitive to the angular structure of the outflow. We calculate the criterion for obtaining single or double peaked light-curves in the general case when $θ_{c,Γ}\neqθ_{c,ε}$. We emphasize how the shape of the lightcurve and the temporal evolution of $ν_m$ and $ν_c$ can be used to constrain the outflow structure and potentially distinguish between magnetic and hydrodynamic jets. △ Less

Submitted 28 June, 2022; v1 submitted 12 April, 2022; originally announced April 2022.

Comments: 15 pages, 13 figures. Accepted for publication in MNRAS

Magnetic field reversal in the turbulent environment around a repeating fast radio burst

Authors: Reshma Anna-Thomas, Liam Connor, Shi Dai, Yi Feng, Sarah Burke-Spolaor, Paz Beniamini, Yuan-Pei Yang, Yongkun Zhang, Kshitij Aggarwal, Casey J. Law, Di Li, Chenhui Niu, Shami Chatterjee, Marilyn Cruces, Ran Duan, Miroslav D. Filipovi, George Hobbs, Ryan S. Lynch, Chenchen Miao, Jiarui Niu, Stella K. Ocker, Chao-Wei Tsai, Pei Wang, Mengyao Xue, Jumei Yao , et al. (5 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are brief, intense flashes of radio waves from unidentified extragalactic sources. Polarized FRBs originate in highly magnetized environments. We report observations of the repeating FRB 20190520B spanning seventeen months , which show its amount of Faraday rotation is highly variable and twice changes its sign. The FRB also depolarizes below radio frequencies around 1 to… ▽ More Fast radio bursts (FRBs) are brief, intense flashes of radio waves from unidentified extragalactic sources. Polarized FRBs originate in highly magnetized environments. We report observations of the repeating FRB 20190520B spanning seventeen months , which show its amount of Faraday rotation is highly variable and twice changes its sign. The FRB also depolarizes below radio frequencies around 1 to 3 GHz. We interpret these properties as due to change in the parallel component of the integrated magnetic field along the line-of-sight, including reversals. This could result from propagation through a turbulent, magnetized screen of plasma located between $10^{-5}$ to 100 parsecs of the FRB source. This is consistent with the bursts passing through the stellar wind of a binary companion of the FRB source. △ Less

Submitted 12 May, 2023; v1 submitted 22 February, 2022; originally announced February 2022.

Comments: arXiv:2203.08151 has been merged into this

Journal ref: Science 380,599-603(2023)

Plug Disintegration in GRB Jet Eruption

Authors: Almog Yalinewich, Paz Beniamini

Abstract: In this work we consider the eruption of a tenuous relativistic hydrodynamic jet from a dense baryonic envelope. As the jet moves out and away, it carries along and continues to accelerate a layer of baryonic material which we refer to as the plug. We solve the relativistic equations of motion for the trajectory of the plug, and verify it using a relativistic hydrodynamic simulation. We show that… ▽ More In this work we consider the eruption of a tenuous relativistic hydrodynamic jet from a dense baryonic envelope. As the jet moves out and away, it carries along and continues to accelerate a layer of baryonic material which we refer to as the plug. We solve the relativistic equations of motion for the trajectory of the plug, and verify it using a relativistic hydrodynamic simulation. We show that under these conditions, the plug breaks up at a radius larger by a factor of a few from the radius of the envelope, due to the onset of the Rayleigh Taylor instability. After breakup the jet continues to accelerate to higher Lorentz factors while the plug fragments maintain a moderate Lorentz factor. The presence of slower moving ejecta can explain late time features of GRBs such as X ray flares without recourse to a long lived engine. △ Less

Submitted 26 January, 2022; originally announced January 2022.

Flares in gamma-ray burst X-ray afterglows as prompt emission from slightly misaligned structured jets

Authors: Raphaël Duque, Paz Beniamini, Frédéric Daigne, Robert Mochkovitch

Abstract: We develop a model to explain the flaring activity in gamma-ray burst X-ray afterglows within the framework of slightly misaligned observers to structured jets. We suggest that flares could be the manifestation of prompt dissipation within the core of the jet, appearing to a misaligned observer in the X-ray band because of less favorable Doppler boosting. These flares appear during the afterglow p… ▽ More We develop a model to explain the flaring activity in gamma-ray burst X-ray afterglows within the framework of slightly misaligned observers to structured jets. We suggest that flares could be the manifestation of prompt dissipation within the core of the jet, appearing to a misaligned observer in the X-ray band because of less favorable Doppler boosting. These flares appear during the afterglow phase because of core--observer light travel delays. In this picture, the prompt emission recorded by this observer comes from material along their line of sight, in the lateral structure of the jet, outside the jet's core. We start by laying down the basic analytical framework to determine the flares characteristics as a function of those of the gamma-ray pulse an aligned observer would see. We show that there is viable parameter space to explain flares with typical observing times and luminosities. We then analytically explore this model, showing that it naturally produces flares with small aspect ratios, as observed. We perform fits of our model to two $Swift$/XRT flares representing two different types of morphology, to show that our model can capture both. The ejection time of the core jet material responsible of the flare is a critical parameter. While it always remains small compared to the observed time of the flare, confirming that our model does not require very late central engine activity, late ejection times are strongly favored, sometimes larger than the observed duration of the parent gamma-ray burst's prompt emission as measured by $T_{90}$. △ Less

Submitted 1 April, 2022; v1 submitted 6 December, 2021; originally announced December 2021.

Comments: Main text 11 pages, 4 figures and 2 tables. Final version to appear in MNRAS

A Toy Model for the Time-Frequency Structure of Fast Radio Bursts: Implications for the CHIME Burst Dichotomy

Authors: Brian D. Metzger, Navin Sridhar, Ben Margalit, Paz Beniamini, Lorenzo Sironi

Abstract: We introduce a toy model for the time-frequency structure of fast radio bursts (FRB), in which the observed emission is produced as a narrowly-peaked intrinsic spectral energy distribution sweeps down in frequency across the instrumental bandpass as a power-law in time. Though originally motivated by emission models which invoke a relativistic shock, the model could in principle apply to a wider r… ▽ More We introduce a toy model for the time-frequency structure of fast radio bursts (FRB), in which the observed emission is produced as a narrowly-peaked intrinsic spectral energy distribution sweeps down in frequency across the instrumental bandpass as a power-law in time. Though originally motivated by emission models which invoke a relativistic shock, the model could in principle apply to a wider range of emission scenarios. We quantify the burst's detectability using the frequency bandwidth over which most of its signal-to-noise ratio (SNR) is accumulated. We demonstrate that by varying just a single parameter of the toy model-the power-law index βof the frequency drift rate-one can transform a long (and hence preferentially time-resolved) burst with a narrow time-integrated spectrum into a shorter burst with a broad power-law time-integrated spectrum. We suggest that burst-to-burst diversity in the value of βcould generate the dichotomy between burst duration and frequency-width recently found by CHIME. In shock models, the value of βis related to the radial density profile of external medium, which in light of the preferentially longer duration of bursts from repeating sources may point to diversity in the external environments surrounding repeating versus one-off FRB sources. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 11 pages, 7 figures, to be submitted to ApJ, comments appreciated

Faraday depolarization and induced circular polarization by multi-path propagation with application to FRBs

Authors: Paz Beniamini, Pawan Kumar, Ramesh Narayan

Abstract: We describe how the observed polarization properties of an astronomical object are related to its intrinsic polarization properties and the finite temporal and spectral resolutions of the observing device. Moreover, we discuss the effect that a scattering screen, with non-zero magnetic field, between the source and observer has on the observed polarization properties. We show that the polarization… ▽ More We describe how the observed polarization properties of an astronomical object are related to its intrinsic polarization properties and the finite temporal and spectral resolutions of the observing device. Moreover, we discuss the effect that a scattering screen, with non-zero magnetic field, between the source and observer has on the observed polarization properties. We show that the polarization properties are determined by the ratio of observing bandwidth and coherence bandwidth of the scattering screen and the ratio of temporal resolution of the instrument and the variability time of screen, as long as the length over which the Faraday rotation induced by the screen changes by $\simπ$ is smaller than the size of the screen visible to the observer. We describe the conditions under which a source that is 100\% linearly polarized intrinsically might be observed as partially depolarized, and how the source's temporal variability can be distinguished from the temporal variability induced by the scattering screen. In general, linearly polarized waves passing through a magnetized scattering screen can develop a significant circular polarization. We apply the work to the observed polarization properties of a few fast radio bursts (FRBs), and outline potential applications to pulsars. △ Less

Submitted 28 December, 2021; v1 submitted 30 September, 2021; originally announced October 2021.

Comments: 15 pages, 9 figures. Accepted for publication in MNRAS

Implications of a rapidly varying FRB in a globular cluster of M81

Authors: Wenbin Lu, Paz Beniamini, Pawan Kumar

Abstract: A repeating source of fast radio bursts (FRBs) is recently discovered from a globular cluster of M81. Association to a globular cluster (or other old stellar systems) suggests that strongly magnetized neutron stars, which are the most likely objects responsible for FRBs, are born not only when young massive stars undergo core-collapse, but also by mergers of old white dwarfs. We find that the frac… ▽ More A repeating source of fast radio bursts (FRBs) is recently discovered from a globular cluster of M81. Association to a globular cluster (or other old stellar systems) suggests that strongly magnetized neutron stars, which are the most likely objects responsible for FRBs, are born not only when young massive stars undergo core-collapse, but also by mergers of old white dwarfs. We find that the fractional contribution to the total FRB rate by old stellar populations is at least a few percent, and the precise fraction can be constrained by FRB searches in the directions of nearby galaxies, both star-forming and elliptical ones. Using very general arguments, we show that the activity time of the M81-FRB object is between 10^4 and 10^6 years, and more likely of order 10^5 years. The energetics of radio outbursts puts a lower limit on the magnetic field strength of 10^{13} G, and the spin period > 0.2 sec, thereby ruling out the source being a milli-second pulsar. The upper limit on the persistent X-ray luminosity (provided by Chandra), together with the high FRB luminosity and frequent repetitions, severely constrains (or rules out) the possibility that the M81-FRB is a scaled-up version of giant pulses from Galactic pulsars. Finally, the 50 ns variability time of the FRB lightcurve suggests that the emission is produced in a compact region inside the neutron star magnetosphere, as it cannot be accounted for when the emission is at distances > 10^{10} cm. △ Less

Submitted 3 December, 2021; v1 submitted 8 July, 2021; originally announced July 2021.

Comments: MNRAS accepted after minor corrections

Survival times of supramassive neutron stars resulting from binary neutron star mergers

Authors: Paz Beniamini, Wenbin Lu

Abstract: A binary neutron star (BNS) merger can lead to various outcomes, from indefinitely stable neutron stars, through supramassive (SMNS) or hypermassive (HMNS) neutron stars supported only temporarily against gravity, to black holes formed promptly after the merger. Up-to-date constraints on the BNS total mass and the neutron star equation of state suggest that a long-lived SMNS may form in… ▽ More A binary neutron star (BNS) merger can lead to various outcomes, from indefinitely stable neutron stars, through supramassive (SMNS) or hypermassive (HMNS) neutron stars supported only temporarily against gravity, to black holes formed promptly after the merger. Up-to-date constraints on the BNS total mass and the neutron star equation of state suggest that a long-lived SMNS may form in $\sim 0.45-0.9$ of BNS mergers. A maximally rotating SMNS needs to lose $\sim 3-6\times 10^{52}$ erg of it's rotational energy before it collapses, on a fraction of the spin-down timescale. A SMNS formation imprints on the electromagnetic counterparts to the BNS merger. However, a comparison with observations reveals tensions. First, the distribution of collapse times is too wide and that of released energies too narrow (and the energy itself too large) to explain the observed distributions of internal X-ray plateaus, invoked as evidence for SMNS-powered energy injection. Secondly, the immense energy injection into the blastwave should lead to extremely bright radio transients which previous studies found to be inconsistent with deep radio observations of short gamma-ray bursts. Furthermore, we show that upcoming all-sky radio surveys will constrain the extracted energy distribution, independently of a GRB jet formation. Our results can be self-consistently understood, provided that most BNS merger remnants collapse shortly after formation (even if their masses are low enough to allow for SMNS formation). This naturally occurs if the remnant retains half or less of its initial energy by the time it enters solid body rotation. △ Less

Submitted 20 July, 2021; v1 submitted 2 April, 2021; originally announced April 2021.

Comments: Accepted for publication in ApJ

Evidence of extended emission in GRB 181123B and other high-redshift short GRBs

Authors: S. Dichiara, E. Troja, P. Beniamini, B. O'Connor, M. Moss, A. Y. Lien, R. Ricci, L. Amati, G. Ryan, T. Sakamoto

Abstract: We study the high-energy properties of GRB 181123B, a short gamma-ray burst (sGRB) at redshift $z\approx$1.75. We show that, despite its nominal short duration with $T_{90}<$2 s, this burst displays evidence of a temporally extended emission (EE) at high energies and that the same trend is observed in the majority of sGRBs at $z\gtrsim$1. We discuss the impact of instrumental selection effects on… ▽ More We study the high-energy properties of GRB 181123B, a short gamma-ray burst (sGRB) at redshift $z\approx$1.75. We show that, despite its nominal short duration with $T_{90}<$2 s, this burst displays evidence of a temporally extended emission (EE) at high energies and that the same trend is observed in the majority of sGRBs at $z\gtrsim$1. We discuss the impact of instrumental selection effects on the GRB classification, stressing that the measured $T_{90}$ is not an unambiguous indicator of the burst physical origin. By examining their environment (e.g. stellar mass, star formation, offset distribution), we find that these high-$z$ sGRBs share many properties of long GRBs at a similar distance and are consistent with a short-lived progenitor system. If produced by compact binary mergers, these sGRBs with EE may be easier to localize at large distances and herald a larger population of sGRBs in the early universe. △ Less

Submitted 4 July, 2022; v1 submitted 3 March, 2021; originally announced March 2021.

Comments: 10 pages, 7 figures, 2 tables. ApJ Letters, in press

Periodic Fast Radio Bursts from Luminous X-ray Binaries

Authors: Navin Sridhar, Brian D. Metzger, Paz Beniamini, Ben Margalit, Mathieu Renzo, Lorenzo Sironi, Konstantinos Kovlakas

Abstract: The discovery of periodicity in the arrival times of the fast radio bursts (FRBs) poses a challenge to the oft-studied magnetar scenarios. However, models that postulate that FRBs result from magnetized shocks or magnetic reconnection in a relativistic outflow are not specific to magnetar engines; instead, they require only the impulsive injection of relativistic energy into a dense magnetized med… ▽ More The discovery of periodicity in the arrival times of the fast radio bursts (FRBs) poses a challenge to the oft-studied magnetar scenarios. However, models that postulate that FRBs result from magnetized shocks or magnetic reconnection in a relativistic outflow are not specific to magnetar engines; instead, they require only the impulsive injection of relativistic energy into a dense magnetized medium. Motivated thus, we outline a new scenario in which FRBs are powered by short-lived relativistic outflows (``flares'') from accreting black holes or neutron stars, which propagate into the cavity of the pre-existing (``quiescent'') jet. In order to reproduce FRB luminosities and rates, we are driven to consider binaries of stellar-mass compact objects undergoing super-Eddington mass-transfer, similar to ultraluminous X-ray (ULX) sources. Indeed, the host galaxies of FRBs, and their spatial offsets within their hosts, show broad similarities with ULXs. Periodicity on timescales of days to years could be attributed to precession (e.g., Lens-Thirring) of the polar accretion funnel, along which the FRB emission is geometrically and relativistically beamed, which sweeps across the observer line of sight. Accounting for the most luminous FRBs via accretion power may require a population of binaries undergoing brief-lived phases of unstable (dynamical-timescale) mass-transfer. This will lead to secular evolution in the properties of some repeating FRBs on timescales of months to years, followed by a transient optical/IR counterpart akin to a luminous red nova, or a more luminous accretion-powered optical/X-ray transient. We encourage targeted FRB searches of known ULX sources. △ Less

Submitted 1 July, 2021; v1 submitted 11 February, 2021; originally announced February 2021.

Comments: Accepted for publication in The Astrophysical Journal; 25 pages, 8 figures

Evidence for r-process delay in very metal-poor stars

Authors: Yuta Tarumi, Kenta Hotokezaka, Paz Beniamini

Abstract: The abundances of r-process elements of very metal-poor stars capture the history of the r-process enrichment in the early stage of star formation in a galaxy. Currently, various types of astrophysical sites including neutron star mergers, magneto-rotational supernovae, and collapsars, are suggested as the origin of r-process elements. The time delay between the star formation and the production o… ▽ More The abundances of r-process elements of very metal-poor stars capture the history of the r-process enrichment in the early stage of star formation in a galaxy. Currently, various types of astrophysical sites including neutron star mergers, magneto-rotational supernovae, and collapsars, are suggested as the origin of r-process elements. The time delay between the star formation and the production of r-process elements is the key to distinguish these scenarios with the caveat that the diffusion of r-process elements in the interstellar medium may induce the delay in r-process enrichment because r-process events are rare. Here we study the observed Ba abundance data of very metal-poor stars as the tracer of the early enrichment history of r-process elements. We find that the gradual increase of [Ba/Mg] with [Fe/H], which is remarkably similar among the Milky Way and classical dwarfs, requires a significant time delay (100 Myr -- 1 Gyr) of r-process events from star formation rather than the diffusion-induced delay. We stress that this conclusion is robust to the assumption regarding s-process contamination in the Ba abundances because the sources with no delay would overproduce Ba at very low metallicities even without the contribution from the s-process. Therefore we conclude that sources with a delay, possibly neutron star mergers, are the origins of r-process elements. △ Less

Submitted 18 May, 2021; v1 submitted 5 February, 2021; originally announced February 2021.

Comments: 10 pages, 6 figures, ApJ Letters in press

No velocity-kicks are required to explain large-distance offsets of Ca-rich supernovae and short-GRBs

Authors: Hagai B. Perets, Paz Beniamini

Abstract: The environments of explosive transients link their progenitors to the underlying stellar population, providing critical clues for their origins. However, some Ca-rich supernovae (SNe) and short gamma ray burst (sGRBs) appear to be located at remote locations, far from the stellar population of their host galaxy, challenging our understanding of their origin and/or physical evolution. These findin… ▽ More The environments of explosive transients link their progenitors to the underlying stellar population, providing critical clues for their origins. However, some Ca-rich supernovae (SNe) and short gamma ray burst (sGRBs) appear to be located at remote locations, far from the stellar population of their host galaxy, challenging our understanding of their origin and/or physical evolution. These findings instigated models suggesting that either large velocity kicks were imparted to their progenitors, allowing them to propagate to large distances and attain their remote locations; or that they formed in dense globular clusters residing in the halos. Here we show that instead, the large spatial-offsets of these transients are naturally explained by the observations of highly extended underlying stellar populations in (mostly early type) galaxy halos, typically missed since they can only be identified through ultra-deep/stacked images. Consequently, no large velocity kicks, nor halo globular cluster environments are required in order to explain the origin of these transients. These findings support thermonuclear explosions on white-dwarfs, for the origins of Ca-rich SNe progenitors, and no/small-natal-kicks given to sGRB progenitors. Since early-type galaxies contain older stellar populations, transient arising from older stellar populations would have larger fractions of early-type galaxy hosts, and consequently larger fractions of transients at large offsets. This is verified by our results for sGRBs and Ca-rich SNe showing different offset distributions in early vs. late-type galaxies. Furthermore, once divided between early and late type galaxies, the offsets' distributions of the different transients are consistent with each other. Finally, we point out that studies of other transients' offset distribution (e.g. Ia-SNe or FRBs) should similarly consider the host galaxy-type. △ Less

Submitted 27 January, 2021; originally announced January 2021.

Comments: Comments are welcome

Exploring the epoch of hydrogen reionization using FRBs

Authors: Paz Beniamini, Pawan Kumar, Xiangcheng Ma, Eliot Quataert

Abstract: We describe three different methods for exploring the hydrogen reionization epoch using fast radio bursts (FRBs) and provide arguments for the existence of FRBs at high redshift (z). The simplest way, observationally, is to determine the maximum dispersion measure (DM$_{\rm max}$) of FRBs for an ensemble that includes bursts during the reionization. The DM$_{\rm max}$ provides information regardin… ▽ More We describe three different methods for exploring the hydrogen reionization epoch using fast radio bursts (FRBs) and provide arguments for the existence of FRBs at high redshift (z). The simplest way, observationally, is to determine the maximum dispersion measure (DM$_{\rm max}$) of FRBs for an ensemble that includes bursts during the reionization. The DM$_{\rm max}$ provides information regarding reionization much like the optical depth of the CMB to Thomson scattering does, and it has the potential to be more accurate than constraints from Planck, if DM$_{\rm max}$ can be measured to a precision better than 500 $\mbox{pc cm}^{-3}$. Another method is to measure redshifts of about 40 FRBs between z of 6-10 with$\sim10\%$ accuracy to obtain the average electron density in 4 different z-bins with $\sim4\%$ accuracy. These two methods don't require knowledge of the FRB luminosity function and its possible redshift evolution. Finally, we show that the reionization history is reflected in the number of FRBs per unit DM, given a fluence limited survey of FRBs that includes bursts during the reionization epoch; we show using FIRE simulations that the contributions to DM from the FRB host galaxy $\&$ CGM during the reionization era is a small fraction of the observed DM. This third method requires no redshift information but does require knowledge of the FRB luminosity function. △ Less

Submitted 1 February, 2021; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: 13 pages, 10 figures. Accepted for publication in MNRAS

On the formation of GW190814

Authors: Wenbin Lu, Paz Beniamini, Clément Bonnerot

Abstract: The LIGO-Virgo collaboration recently reported a puzzling event, GW190814, with component masses of 23 and 2.6 solar masses. Motivated by the relatively small rate of such a coalescence and the fact that the mass of the secondary is close to the total mass of known binary neutron star (bNS) systems, we propose that GW190814 was a second-generation merger from a hierarchical triple system, i.e., th… ▽ More The LIGO-Virgo collaboration recently reported a puzzling event, GW190814, with component masses of 23 and 2.6 solar masses. Motivated by the relatively small rate of such a coalescence and the fact that the mass of the secondary is close to the total mass of known binary neutron star (bNS) systems, we propose that GW190814 was a second-generation merger from a hierarchical triple system, i.e., the remnant from the bNS coalescence was able to merge again with the 23 Msun black hole (BH) tertiary. We show that this occurs at a sufficiently high probability provided that the semimajor axis of the outer orbit is less than a few AU at the time of bNS coalescence. It remains to be explored whether the conditions for the formation of such tight triple systems are commonly realized in the Universe, especially in low metallicity (less than 0.1 solar) environments. Our model provides a number of predictions. (1) The spin of the secondary in GW190814-like systems is 0.6 to 0.7. (2) The component mass distribution from a large sample of LIGO sources should have a narrow peak between 2.5 and ~3.5 Msun, whereas the range between ~3.5 and ~5 Msun stays empty (provided that stellar evolution does not generate such BHs in the "mass gap"). (3) About 90% (10%) of GW190814-like events have an eccentricity of e > 2x10^{-3} (> 0.1) near gravitational wave frequency of 10 mHz. (4) A significant fraction (> 10%) of bNS mergers should have signatures of a massive tertiary at a distance of a few AU in the gravitational waveform. (5) There are 10^5 undetected radio-quiet bNS systems with a massive BH tertiary in the Milky Way. △ Less

Submitted 26 October, 2020; v1 submitted 21 September, 2020; originally announced September 2020.

Comments: MNRAS accepted

GRB Spectrum from Gradual Dissipation in a Magnetized Outflow

Authors: Ramandeep Gill, Jonathan Granot, Paz Beniamini

Abstract: Modeling of gamma-ray burst (GRB) prompt emission spectra sometimes requires a (quasi-) thermal spectral component in addition to the Band function. In photospheric emission models, a prominent thermal component broadened by sub-photospheric dissipation is expected to be released at the photospheric radius, $r_{\rm ph}\sim10^{12}\,$cm. We consider an ultra-relativistic strongly magnetized outflow… ▽ More Modeling of gamma-ray burst (GRB) prompt emission spectra sometimes requires a (quasi-) thermal spectral component in addition to the Band function. In photospheric emission models, a prominent thermal component broadened by sub-photospheric dissipation is expected to be released at the photospheric radius, $r_{\rm ph}\sim10^{12}\,$cm. We consider an ultra-relativistic strongly magnetized outflow with a striped-wind magnetic-field structure undergoing gradual and continuous magnetic energy dissipation at $r<r_s$ that heats and accelerates the flow, leading to a bulk Lorentz factor $Γ(r)=Γ_\infty\min[1,(r/r_s)^{1/3}]$, where typically $r_{\rm ph}<r_s$. Similar dynamics and energy dissipation rates are also expected in highly-variable magnetized outflows without stripes/field-reversals. Two modes of particle energy injection are considered: (a) power-law electrons, e.g. accelerated by magnetic reconnection, and (b) continuous distributed heating of all electrons (and $e^\pm$-pairs), e.g. due to MHD instabilities. Time-resolved energy spectra are obtained using a numerical code that evolves coupled kinetic equations for a photon-electron-positron plasma. We find that (i) the thermal component peaks at $(1+z)E_{\rm pk}\sim0.2-1\,$MeV, for a source at redshift $z$, and becomes subdominant if the total injected energy density exceeds the thermal one, (ii) power-law electrons cool mainly by synchrotron emission whereas mildly relativistic and almost monoenergetic electrons in the distributed heating scenario cool by Comptonization on thermal peak photons, (iii) both scenarios can yield a low-energy break at $E_{\rm br}\approx E_{\rm th}$, and (iv) the $0.5(1+z)^{-1}\,$keV X-ray emission is suppressed in the power-law injection case, but it is expected for the distributed heating scenario. Energy-dependent linear polarization can differentiate between the two energy injection cases. △ Less

Submitted 16 September, 2020; v1 submitted 24 August, 2020; originally announced August 2020.

Comments: 18 pages, 8 figures, Accepted for publication in MNRAS

What does FRB light-curve variability tell us about the emission mechanism?

Authors: Paz Beniamini, Pawan Kumar

Abstract: A few fast radio bursts' (FRBs) light-curves have exhibited large intrinsic modulations of their flux on extremely short ($t_{\rm r}\sim 10μ$s) time scales, compared to pulse durations ($t_{\rm FRB}\sim1$ms). Light-curve variability timescales, the small ratio of rise time of the flux to pulse duration, and the spectro-temporal correlations in the data constrain the compactness of the source and t… ▽ More A few fast radio bursts' (FRBs) light-curves have exhibited large intrinsic modulations of their flux on extremely short ($t_{\rm r}\sim 10μ$s) time scales, compared to pulse durations ($t_{\rm FRB}\sim1$ms). Light-curve variability timescales, the small ratio of rise time of the flux to pulse duration, and the spectro-temporal correlations in the data constrain the compactness of the source and the mechanism responsible for the powerful radio emission. The constraints are strongest when radiation is produced far ($\gtrsim 10^{10}$cm) from the compact object. We describe different physical set-ups that can account for the observed $t_{\rm r}/t_{\rm FRB}\ll 1$ despite having large emission radii. The result is either a significant reduction in the radio production efficiency or distinct light-curves features that could be searched for in observed data. For the same class of models, we also show that due to high-latitude emission, if a flux $f_1(ν_1)$ is observed at $t_1$ then at a lower frequency $ν_2<ν_1$ the flux should be at least $(ν_2/ν_1)^2f_1$ at a slightly later time ($t_2=t_1ν_1/ν_2$) independent of the duration and spectrum of the emission in the comoving frame. These features can be tested, once light-curve modulations due to scintillation are accounted for. We provide the timescales and coherence bandwidths of the latter for a range of possibilities regarding the physical screens and the scintillation regime. Finally, if future highly resolved FRB light-curves are shown to have intrinsic variability extending down to $\sim μ$s timescales, this will provide strong evidence in favor of magnetospheric models. △ Less

Submitted 14 August, 2020; v1 submitted 14 July, 2020; originally announced July 2020.

Comments: 14 pages, 9 figures. Accepted for publication in MNRAS

Modeling Synchrotron Self-Compton and Klein-Nishina Effects in Gamma-Ray Burst Afterglows

Authors: Taylor Jacovich, Paz Beniamini, Alexander van der Horst

Abstract: We present a self-consistent way of modeling synchrotron self-Compton (SSC) effects in gamma-ray burst afterglows, with and without approximated Klein-Nishina suppressed scattering. We provide an analytic approximation of our results, so that it can be incorporated into the afterglow modeling code \texttt{boxfit}, which is currently based on pure synchrotron emission. We discuss the changes in spe… ▽ More We present a self-consistent way of modeling synchrotron self-Compton (SSC) effects in gamma-ray burst afterglows, with and without approximated Klein-Nishina suppressed scattering. We provide an analytic approximation of our results, so that it can be incorporated into the afterglow modeling code \texttt{boxfit}, which is currently based on pure synchrotron emission. We discuss the changes in spectral shape and evolution due to SSC effects, and comment on how these changes affect physical parameters derived from broadband modeling. We show that SSC effects can have a profound impact on the shape of the X-ray light curve using simulations including these effects. This leads to data that cannot be simultaneously fit well in both the X-ray and radio bands when considering synchrotron-only fits, and an inability to recover the correct physical parameters, with some fitted parameters deviating orders of magnitude from the simulated input parameters. This may have a significant impact on the physical parameter distributions based on previous broadband modeling efforts. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: 15 pages, 11 figures, 6 tables. Submitted to MNRAS

ZTF20aajnksq (AT2020blt): A Fast Optical Transient at $z \approx 2.9$ With No Detected Gamma-Ray Burst Counterpart

Authors: Anna Y. Q. Ho, Daniel A. Perley, Paz Beniamini, S. Bradley Cenko, S. R. Kulkarni, Igor Andreoni, Leo P. Singer, Kishalay De, Mansi M. Kasliwal, Christoffer Fremling, Eric C. Bellm, Richard Dekany, Alexandre Delacroix, Dmitry A. Duev, Daniel A. Goldstein, V. Zach Golkhou, Ariel Goobar, Matthew Graham, David Hale, Thomas Kupfer, Russ R. Laher, Frank J. Masci, A. A. Miller, James D. Neill, Reed Riddle , et al. (5 additional authors not shown)

Abstract: We present ZTF20aajnksq (AT2020blt), a fast-fading ($Δr=2.4$ mag in $Δt=1.3$ days) red ($g-r\approx0.6$ mag) and luminous ($M_{1626}=-25.9$) optical transient at $z=2.9$ discovered by the Zwicky Transient Facility (ZTF). AT2020blt shares several features in common with afterglows to long-duration gamma-ray bursts (GRBs): (1) an optical light curve well-described by a broken power-law with a break… ▽ More We present ZTF20aajnksq (AT2020blt), a fast-fading ($Δr=2.4$ mag in $Δt=1.3$ days) red ($g-r\approx0.6$ mag) and luminous ($M_{1626}=-25.9$) optical transient at $z=2.9$ discovered by the Zwicky Transient Facility (ZTF). AT2020blt shares several features in common with afterglows to long-duration gamma-ray bursts (GRBs): (1) an optical light curve well-described by a broken power-law with a break at $t_\mathrm{j}=1$ day (observer-frame); (2) a luminous $(L_X = 10^{46}$ $\mathrm{erg}$ $\mathrm{s}^{-1})$ X-ray counterpart; and (3) luminous ($L_ν= 4 \times 10^{31}$ $\mathrm{erg}$ $\mathrm{sec}^{-1}$ $\mathrm{Hz}^{-1}$ at 10 GHz) radio emission. However, no GRB was detected in the 0.74d between the last ZTF non-detection ($r > 20.64$) and the first ZTF detection ($r = 19.57$), with an upper limit on the isotropic-equivalent gamma-ray energy release of $E_{γ,\mathrm{iso}} < 7 \times 10^{52}$ erg. AT2020blt is thus the third afterglow-like transient discovered without a detected GRB counterpart (after PTF11agg and ZTF19abvizsw) and the second (after ZTF19abvizsw) with a redshift measurement. We conclude that the properties of AT2020blt are consistent with a classical (initial Lorentz factor $Γ_0 \gtrsim 100$) on-axis GRB that was missed by high-energy satellites. Furthermore, by estimating the rate of transients with light curves similar to that of AT2020blt in ZTF high-cadence data, we agree with previous results that there is no evidence for an afterglow-like phenomenon that is significantly more common than classical GRBs. We conclude by discussing the status and future of fast-transient searches in wide-field high-cadence optical surveys. △ Less

Submitted 19 October, 2020; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 19 pages, 4 figures, 1 table. Version accepted by ApJ on 17 Oct 2020

GRB Fermi-LAT afterglows: explaining flares, breaks, and energetic photons

Authors: N. Fraija, T. Laskar, S. Dichiara, P. Beniamini, R. Barniol Duran, M. G. Dainotti, R. L. Becerra

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

Submitted 21 November, 2020; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 19 pages, 5 figures. Accepted for publication in ApJ. New affiliations to the author list are added

Deciphering the properties of the central engine in GRB collapsars

Authors: M. Petropoulou, P. Beniamini, G. Vasilopoulos, D. Giannios, R. Barniol Duran

Abstract: The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic pr… ▽ More The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic properties of central engines that are consistent with prompt observations of long GRBs detected by the Burst Alert Telescope (BAT) on board the Neil Gehrels Swift Observatory. Adopting a generic model for the central engine, in which the engine power and activity timescale are independent of each other, we perform Monte Carlo simulations of long GRBs produced by jets that successfully breakout from the star. Our simulations consider the dependence of the jet breakout timescale on the engine luminosity and the effects of the detector's flux threshold. The two-dimensional (2D) distribution of simulated detectable bursts in the gamma-ray luminosity versus gamma-ray duration plane is consistent with the observed one for a range of parameter values describing the central engine. The intrinsic 2D distribution of simulated collapsar GRBs peaks at lower gamma-ray luminosities and longer durations than the observed one, a prediction that can be tested in the future with more sensitive detectors. Black-hole accretors, whose power and activity time are set by the large-scale magnetic flux through the progenitor star and stellar structure, respectively, are compatible with the properties of the central engine inferred by our model. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Comments: 13 pages, 6 figures, accepted for publication in MNRAS

Deceleration of relativistic jets with lateral expansion

Authors: Wenbin Lu, Paz Beniamini, Austin McDowell

Abstract: We present a model for the hydrodynamics of a relativistic jet interacting with the circum-stellar medium (CSM). The shocked CSM and the jet material are assumed to be in an infinitely thin surface, so the original 2D problem is effectively reduced to 1D. From general conservation laws, we derive the equation of motion for each fluid element along this surface, taking into account the deceleration… ▽ More We present a model for the hydrodynamics of a relativistic jet interacting with the circum-stellar medium (CSM). The shocked CSM and the jet material are assumed to be in an infinitely thin surface, so the original 2D problem is effectively reduced to 1D. From general conservation laws, we derive the equation of motion for each fluid element along this surface, taking into account the deceleration along the surface normal due to newly swept-up mass and lateral expansion due to pressure gradient in the tangential direction. The pressure and energy density of the shocked CSM are given by the jump conditions at the forward shock. The method is implemented with a finite-differencing numerical scheme, along with calculation of synchrotron emission and absorption from shock-accelerated electrons, in a new code $\mathtt{Jedi}$ (for "jet dynamics"). We present a number of test cases, including top-hat jet, power-law structured jet, "boosted fireball" profile, and CSM with density jump at the wind termination shock. Based on the agreement with other analytical and numerical calculations, we conclude that our simplified method provides a good approximation for the hydrodynamics and afterglow emission for a wide variety of jet structures and CSM density profiles. Efficient modeling of the afterglow from e.g., neutron star mergers, will provide important information on the jet energetics, CSM properties, and the viewing angle. △ Less

Submitted 20 May, 2020; originally announced May 2020.

Comments: 12 pages, 12 figures, plus appendices. Comments welcome!

Implications of a "Fast Radio Burst" from a Galactic Magnetar

Authors: Ben Margalit, Paz Beniamini, Navin Sridhar, Brian D. Metzger

Abstract: A luminous radio burst was recently detected in temporal coincidence with a hard X-ray flare from the Galactic magnetar SGR 1935+2154 with a time and frequency structure consistent with cosmological fast radio bursts (FRB) and a fluence within a factor of $\lesssim 10$ of the least energetic extragalactic FRB previously detected. Although active magnetars are commonly invoked FRB sources, several… ▽ More A luminous radio burst was recently detected in temporal coincidence with a hard X-ray flare from the Galactic magnetar SGR 1935+2154 with a time and frequency structure consistent with cosmological fast radio bursts (FRB) and a fluence within a factor of $\lesssim 10$ of the least energetic extragalactic FRB previously detected. Although active magnetars are commonly invoked FRB sources, several distinct mechanisms have been proposed for generating the radio emission which make different predictions for the accompanying higher frequency radiation. We show that the properties of the coincident radio and X-ray flares from SGR 1935+2154, including their approximate simultaneity and relative fluence $E_{\rm radio}/E_{\rm X} \sim 10^{-5}$, as well as the duration and spectrum of the X-ray emission, are consistent with extant predictions for the synchrotron maser shock model. Rather than arising from the inner magnetosphere, the X-rays are generated by (incoherent) synchrotron radiation from thermal electrons heated at the same shocks which produce the coherent maser emission. Although the rate of SGR 1935+2154-like bursts in the local universe is not sufficient to contribute appreciably to the extragalactic FRB rate, the inclusion of an additional population of more active magnetars with stronger magnetic fields than the Galactic population can explain both the FRB rate as well as the repeating fraction, however only if the population of active magnetars are born at a rate that is at least two-orders of magnitude lower than that of SGR 1935+2154-like magnetars. This may imply that the more active magnetar sources are not younger magnetars formed in a similar way to the Milky Way population (e.g. via ordinary supernovae), but instead through more exotic channels such as superluminous supernovae, accretion-induced collapse or neutron star mergers. △ Less

Submitted 17 May, 2020; v1 submitted 11 May, 2020; originally announced May 2020.

Comments: 21 pages, 9 figures; submitted to ApJL; comments welcome!

Constraints on the circumburst environments of short gamma-ray bursts

Authors: Brendan O'Connor, Paz Beniamini, Chryssa Kouveliotou

Abstract: Observational follow-up of well localized short gamma-ray bursts (SGRBs) has left $20-30\%$ of the population without a coincident host galaxy association to deep optical and NIR limits ($\gtrsim 26$ mag). These SGRBs have been classified as observationally hostless due to their lack of strong host associations. It has been argued that these hostless SGRBs could be an indication of the large dista… ▽ More Observational follow-up of well localized short gamma-ray bursts (SGRBs) has left $20-30\%$ of the population without a coincident host galaxy association to deep optical and NIR limits ($\gtrsim 26$ mag). These SGRBs have been classified as observationally hostless due to their lack of strong host associations. It has been argued that these hostless SGRBs could be an indication of the large distances traversed by the binary neutron star system (due to natal kicks) between its formation and its merger (leading to a SGRB). The distances of GRBs from their host galaxies can be indirectly probed by the surrounding circumburst densities. We show that a lower limit on those densities can be obtained from early afterglow lightcurves. We find that $\lesssim16\%$ of short GRBs in our sample took place at densities $\lesssim10^{-4}$ cm$^{-3}$. These densities represent the expected range of values at distances greater than the host galaxy's virial radii. We find that out of the five SGRBs in our sample that have been found to be observationally hostless, none are consistent with having occurred beyond the virial radius of their birth galaxies. This implies one of two scenarios. Either these observationally hostless SGRBs occurred outside of the half-light radius of their host galaxy, but well within the galactic halo, or in host galaxies at moderate to high redshifts ($z\gtrsim 2$) that were missed by follow-up observations. △ Less

Submitted 20 May, 2020; v1 submitted 31 March, 2020; originally announced April 2020.

Comments: 20 pages, 9 figures, accepted for publication in MNRAS

Periodicity in recurrent fast radio bursts and the origin of ultra long period magnetars

Authors: Paz Beniamini, Zorawar Wadiasingh, Brian D. Metzger

Abstract: The recurrent fast radio burst FRB 180916 was recently shown to exhibit a 16 day period (with possible aliasing) in its bursting activity. Given magnetars as widely considered FRB sources, this period has been attributed to precession of the magnetar spin axis or the orbit of a binary companion. Here, we make the simpler connection to a {\it rotational period}, an idea observationally motivated by… ▽ More The recurrent fast radio burst FRB 180916 was recently shown to exhibit a 16 day period (with possible aliasing) in its bursting activity. Given magnetars as widely considered FRB sources, this period has been attributed to precession of the magnetar spin axis or the orbit of a binary companion. Here, we make the simpler connection to a {\it rotational period}, an idea observationally motivated by the 6.7 hour period of the Galactic magnetar candidate, 1E 161348--5055. We explore three physical mechanisms that could lead to the creation of ultra long period magnetars: (i) enhanced spin-down due to episodic mass-loaded charged particle winds (e.g. as may accompany giant flares), (ii) angular momentum kicks from giant flares and (iii) fallback leading to long lasting accretion disks. We show that particle winds and fallback accretion can potentially lead to a sub-set of the magnetar population with ultra long periods, sufficiently long to accommodate FRB 180916 or 1E 161348--5055. If confirmed, such periods implicate magnetars in relatively mature states (ages $1-10$ kyr) and which possessed large internal magnetic fields at birth $B_{\rm int}\gtrsim 10^{16}$ G. In the low-twist magnetar model for FRBs, such long period magnetars may dominate FRB production for repeaters at lower isotropic-equivalent energies and broaden the energy distribution beyond that expected for a canonical population of magnetars which terminate their magnetic activity at shorter periods $P \lesssim 10$ s. △ Less

Submitted 19 June, 2020; v1 submitted 27 March, 2020; originally announced March 2020.

Comments: 12 pages, 5 figures. Accepted for publication in MNRAS

On the origin of the multi-GeV photons from the closest burst with intermediate luminosity: GRB 190829A

Authors: N. Fraija, P. Veres, P. Beniamini, A. Galvan-Gamez, B. D. Metzger, R. Barniol Duran, R. L. Becerra

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

Submitted 25 March, 2020; originally announced March 2020.

Comments: 26 pages, 7 figures, Submitted to ApJ

Turbulent mixing of r-process elements in the Milky Way

Authors: Paz Beniamini, Kenta Hotokezaka

Abstract: We study turbulent gas diffusion affects on $r$-process abundances in Milky Way stars, by a combination of an analytical approach and a Monte Carlo simulation. Higher $r$-process event rates and faster diffusion, lead to more efficient mixing corresponding to a reduced scatter of $r$-process abundances and causing $r$-process enriched stars to start appearing at lower metallicities. We use three i… ▽ More We study turbulent gas diffusion affects on $r$-process abundances in Milky Way stars, by a combination of an analytical approach and a Monte Carlo simulation. Higher $r$-process event rates and faster diffusion, lead to more efficient mixing corresponding to a reduced scatter of $r$-process abundances and causing $r$-process enriched stars to start appearing at lower metallicities. We use three independent observations to constrain the model parameters: (i) the scatter of radioactively stable $r$-process element abundances, (ii) the largest $r$-process enrichment values observed in any solar neighborhood stars and (iii) the isotope abundance ratios of different radioactive $r$-process elements ($^{244}$Pu/$^{238}$U and $^{247}$Cm/$^{238}$U) at the early solar system as compared to their formation. Our results indicate that the Galactic $r$-process rate and the diffusion coefficient are respectively $r<4\times 10^{-5}\mbox{ yr}^{-1}, D>0.1 \mbox{ kpc}^2\mbox{Gyr}^{-1}$ ($r<4\times 10^{-6}\mbox{ yr}^{-1}, D>0.5 \mbox{ kpc}^2\mbox{Gyr}^{-1}$ for collapsars or similarly prolific $r$-process sources) with allowed values satisfying an approximate anti-correlation such that $D\approx r^{-2/3}$, implying that the time between two $r$-process events that enrich the same location in the Galaxy, is $τ_{\rm mix}\approx 100-200\mbox{ Myr}$. This suggests that a fraction of $\sim 0.8$ ($\sim 0.5$) of the observed $^{247}$Cm ($^{244}$Pu) abundance is dominated by one $r$-process event in the early solar system. Radioactively stable element abundances are dominated by contributions from $\sim 10$ different events in the early solar system. For metal poor stars (with [Fe/H]$\lesssim -2$), their $r$-process abundances are dominated by either a single or several events, depending on the star formation history. △ Less

Submitted 16 June, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: 11 pages, 9 figures, accepted for publication in MNRAS

Afterglow Lightcurves from Misaligned Structured Jets

Authors: Paz Beniamini, Jonathan Granot, Ramandeep Gill

Abstract: GRB 170817A/GW 170817 is the first GRB clearly viewed far from the GRB jet's symmetry axis. Its afterglow was densely monitored over a wide range of frequencies and times. It has been modeled extensively, primarily numerically, and although this endeavour was very fruitful, many of the underlying model parameters remain undetermined. We provide analytic modelling of GRB afterglows observed off-axi… ▽ More GRB 170817A/GW 170817 is the first GRB clearly viewed far from the GRB jet's symmetry axis. Its afterglow was densely monitored over a wide range of frequencies and times. It has been modeled extensively, primarily numerically, and although this endeavour was very fruitful, many of the underlying model parameters remain undetermined. We provide analytic modelling of GRB afterglows observed off-axis, considering jets with a narrow core (of half-opening angle $θ_{\rm c}$) and power-law wings in energy per unit solid angle ($ε=ε_cΘ^{-a}$ where $Θ=[1+(θ/θ_{\rm c})^2]^{1/2}$) and initial specific kinetic energy ($Γ_0-1=[Γ_{\rm c,0}-1]Θ^{-b}$), as well as briefly discuss Gaussian jets. Our study reveals qualitatively different types of lightcurves that can be viewed in future off-axis GRBs, with either single or double peaks, depending on the jet structure and the viewing angle. Considering the lightcurve shape rather than the absolute normalizations of times and / or fluxes, removes the dependence of the lightcurve on many of the highly degenerate burst parameters. This study can be easily used to determine the underlying jet structure, significantly reduce the effective parameter space for numerical fitting attempts and provide physical insights. As an illustration, we show that for GRB 170817A, there is a strong correlation between the allowed values of $Γ_{\rm c,0}$ and $b$, leading to a narrow strip of allowed solutions in the $Γ_{\rm c,0}$-$b$ plane above some minimal values $Γ_{\rm c,0}\gtrsim 40, b\gtrsim1.2$. Furthermore, the Lorentz factor of the material dominating the early lightcurve can be constrained by three independent techniques to be $Γ_{0}(θ_{\rm min,0})\approx5-7$. △ Less

Submitted 4 March, 2020; v1 submitted 7 January, 2020; originally announced January 2020.

Comments: 14 pages, 17 figures. Published in MNRAS

Ready, set, launch: time interval between BNS merger and short GRB jet formation

Authors: Paz Beniamini, Rodolfo Barniol Duran, Maria Petropoulou, Dimitrios Giannios

Abstract: The joint detection of GW~170817/GRB 170817 confirmed the long-standing theory that binary neutron star mergers produce short gamma-ray burst (sGRB) jets that can successfully break out of the surrounding ejecta. At the same time, the association with a kilonova provided unprecedented information regarding the physical properties (such as masses and velocities) of the different ejecta constituents… ▽ More The joint detection of GW~170817/GRB 170817 confirmed the long-standing theory that binary neutron star mergers produce short gamma-ray burst (sGRB) jets that can successfully break out of the surrounding ejecta. At the same time, the association with a kilonova provided unprecedented information regarding the physical properties (such as masses and velocities) of the different ejecta constituents. Combining this knowledge with the observed luminosities and durations of cosmological sGRBs detected by the Burst Alert Telescope (BAT) onboard the Neil Gehrels Swift Observatory, we revisit the breakout conditions of sGRB jets. Assuming self-collimation of sGRB jets does not play a critical role, we find that the time interval between the binary merger and the launching of a typical sGRB jet is $\lesssim0.1$~s. We also show that for a fraction of at least $\sim 30\%$ of sGRBs, the usually adopted assumption of static ejecta is inconsistent with observations, even if the polar ejecta mass is an order of magnitude smaller than the one in GRB 170817. Our results disfavour magnetar central engines for powering cosmological sGRBs, limit the amount of energy deposited in the cocoon prior to breakout, and suggest that the observed delay of $\sim 1.$7~s in GW 170817 /GRB 170817 between the gravitational wave and $γ$-ray signals is likely dominated by the propagation time of the jet to the $γ$-ray production site. △ Less

Submitted 18 May, 2020; v1 submitted 3 January, 2020; originally announced January 2020.

Comments: 15 pages, 8 figures. Accepted for publication in ApJ Letters