Gamma-ray Bursts Research Articles

A simulation study on the sub-threshold joint gravitational wave-electromagnetic wave observation on binary neutron star mergers

Abstract The coalescence of binary neutron stars (BNS) is a prolific source of gravitational waves (GWs) and electromagnetic (EM) radiation, offering a dual observational window into the Universe. Lowering the signal-to-noise ratio (S/N) threshold is a simple and cost-effective way to enhance the detection probability of GWs from BNS mergers. In this study, we introduce a metric of the purity of joint GW and EM detections Pjoint, which is in analogue to Pastro in GW only observations. By simulating BNS merger GWs jointly detected by the HLV network and EM counterparts (kilonovae and short Gamma-ray bursts, sGRBs) with an assumed merger rate density of BNS, we generate catalogs of GW events and EM counterparts. Through this simulation, we analyze joint detection pairs, both correct and misidentified. We find the following: 1. For kilonovae, requiring Pjoint > 95% instead of $P_{\rm astro}>95~{{\%}}$ reduces the S/N from 9.2 to 8.5-8.8, allowing 5-13 additional joint detections per year and increasing the GW detection volume by 9-17%; 2. For sGRBs, requiring Pjoint > 95% instead of Pastro reduces the S/N from 9.2 to 8.1-8.5; 3. Increasing kilonova or sGRB detection capability does not improve Pjoint due to a higher rate of misidentifications. We also show that sub-threshold GW and kilonova detections can reduce the uncertainty in measuring the Hubble constant to 89-92% of its original value, and sub-threshold GW and sGRB observations can enhance the precision of constraining the speed of GWs to 88% of previously established values.

DUNE: Dust depletion UNified method across cosmic time and Environments

We present a novel method to characterize dust depletion, namely, the depletion of metals into dust grains. We used observed correlations among relative abundances combining a total of 17 metals in diverse galactic environments, including the Milky Way (MW), Large Magellanic Cloud (LMC), Small Magellanic Cloud (SMC), and damped Lyman-alpha absorbers (DLAs) towards quasars and gamma-ray bursts (GRBs). We only considered the relative abundances of metals that qualify as tracers of dust and we used all available dust tracers. We find linear correlations among all studied dust tracers in a multidimensional space, where each dimension corresponds to an individual dust tracer. The fit to the linear correlations among the dust tracers describes the tendencies of different elements when depleting into dust grains. We determined the overall strength of dust depletion, Delta , along individual lines of sight, based on the correlations among different dust tracers. We avoided any preference for specific dust tracers or any other assumptions by including all available dust tracers in this multidimensional space. We also determined the dust depletion of Kr, C, O, Cl, P, Zn, Ge, Mg, Cu, Si, Fe, Ni, and Ti. Finally, we offer simple guidelines for the application of the method to the study of the observed patterns of abundances and relative abundances. This has allowed for a straightforward determination of the overall strength of depletion and the dust depletion of individual elements. We also obtained an estimate for the gas-phase metallicity and identified any additional deviations due to the nucleosynthesis of specific stellar populations. Thus, we have established a unified methodology for characterizing dust depletion across cosmic time and diverse galactic environments, offering a valuable new approach to the study of dust depletion in studies of the chemical evolution of galaxies.

Synchrotron self-Compton in a radiative-adiabatic fireball scenario: Modelling the multiwavelength observations in some Fermi/LAT bursts

Abstract Energetic GeV photons expected from the closest and the most energetic Gamma-ray bursts (GRBs) provide an unique opportunity to study the very-high-energy emission as well as the possible correlations with lower energy bands in realistic GRB afterglow models. In the standard GRB afterglow model, the relativistic homogeneous shock is usually considered to be fully adiabatic, however, it could be partially radiative. Based on the external forward-shock scenario in both stellar wind and constant-density medium. We present a radiative-adiabatic analytical model of the synchrotron self-Compton (SSC) and synchrotron processes considering an electron energy distribution with a power-law index of 1 < p < 2 and 2 ≤ p. We show that the SSC scenario plays a relevant role in the radiative parameter ε, leading to a prolonged evolution during the slow cooling regime. In a particular case, we derive the Fermi/LAT light curves together with the photons with energies ≥100 MeV in a sample of nine bursts from the second Fermi/LAT GRB catalog that exhibited temporal and spectral indices with ≳ 1.5 and ≈2, respectively. These events can hardly be described with closure relations of the standard synchrotron afterglow model, and also exhibit energetic photons above the synchrotron limit. We have modeled the multi-wavelength observations of our sample to constrain the microphysical parameters, the circumburst density, the bulk Lorentz factor and the mechanism responsible for explaining the energetic GeV photons.

Polarization from a Radially Stratified GRB Outflow

While the dominant radiation mechanism of gamma-ray bursts (GRBs) remains a question of debate, synchrotron emission is one of the foremost candidates to describe the multi-wavelength afterglow observations. As such, it is expected that GRBs should present some degree of polarization across their evolution—presenting a feasible means of probing these bursts’ energetic and angular properties. Although obtaining polarization data is difficult due to the inherent complexities regarding GRB observations, advances are being made, and theoretical modeling of synchrotron polarization is now more relevant than ever. In this manuscript, we present the polarization for a fiduciary model, where the synchrotron FS emission evolving in the radiative–adiabatic regime is described by a radially stratified off-axis outflow. This is parameterized with a power-law velocity distribution and decelerated in a constant-density and wind-like external environment. We apply this theoretical polarization model for two select GRBs, presenting upper limits in their polarization—GRB 170817A, a known off-axis GRB with radio polarization upper limits, and GRB 190014C, an on-axis case, where the burst was seen from within the half-opening angle of the jet, with observed optical polarization—in an attempt to constrain their magnetic field geometry in the emitting region.

Diversity in Fermi/GBM Gamma-Ray Bursts: New Insights from Machine Learning

Classification of gamma-ray bursts (GRBs) has been a long-standing puzzle in high-energy astrophysics. Recent observations challenge the traditional short versus long viewpoint, where long GRBs are thought to originate from the collapse of massive stars and short GRBs from compact binary mergers. Machine learning (ML) algorithms have been instrumental in addressing this problem, revealing five distinct GRB groups within the Swift Burst Alert Telescope (BAT) light-curve data, two of which are associated with kilonovae (KNe). In this work, we extend our analysis to the Fermi Gamma-ray Burst Monitor catalog and identify five clusters using unsupervised ML techniques, consistent with the Swift/BAT results. These five clusters are well separated in the fluence-duration plane, hinting at a potential link between fluence, duration, and complexities (or structures) in the light curves of GRBs. Further, we confirm two distinct classes of KN-associated GRBs. The presence of GRB 170817A in one of the two KN-associated clusters lends evidence to the hypothesis that this class of GRBs could potentially be produced by binary neutron star mergers. The second KN-associated GRB cluster could potentially originate from neutron star–black hole mergers. Future multimessenger observations of compact binaries in gravitational waves and electromagnetic waves can be paramount in understanding these clusters better.

Varying linear polarisation in the dust-free gamma-ray burst 210610B

Context. Long gamma-ray bursts (GRBs) are produced by the collapse of some very massive stars, that emit ultra-relativistic jets. When the jets collide with the interstellar medium they decelerate and generate the so-called afterglow emission, which has been observed to be polarised. Aims. We study the polarimetric evolution of the GRB 210610B afterglow, at z = 1.1341. This allows us to evaluate the role of geometric and/or magnetic mechanisms in the GRB afterglow polarisation. Methods. We observed GRB 210610B using imaging polarimetry with CAFOS on the 2.2 m Calar Alto Telescope and FORS2 on the 4 × 8.1 m Very Large Telescope. Complementary optical spectroscopy was obtained with OSIRIS on the 10.4 m Gran Telescopio Canarias. We studied the GRB light-curve from X-rays to the optical bands and the Spectral Energy Distribution (SED). This allowed us to strongly constrain the line-of-sight extinction. Finally, we studied the GRB host galaxy using optical to NIR data to fit the SED and derive its integrated properties. Results. GRB 210610B had a bright afterglow with a negligible line-of-sight extinction. Polarimetry was obtained at three epochs: during an early plateau phase, at the time when the light curve breaks, and after the light curve steepened. We observe an initial polarisation of ∼4% that goes to zero at the time of the break, and it then again increases to ∼2%, with a change in the position angle of 54 ± 9 deg. The spectrum shows features with very low equivalent widths. This indicate a small amount of material in the line of sight within the host. Conclusions. The lack of dust and the low amount of material in the line of sight to GRB 210610B allowed us to study the intrinsic polarisation of the GRB optical afterglow. The GRB polarisation signals are consistent with ordered magnetic fields in refreshed shock or/and hydrodynamics-scale turbulent fields in the forward shock.

A simple model of dust extinction in gamma-ray burst host galaxies

Context. Gamma-ray burst (GRB) afterglows are powerful probes for studying the different properties of their host galaxies (e.g., the interstellar dust) at all redshifts. By fitting their spectral energy distribution (SED) over a large range of wavelengths, we can gain direct insights into the properties of the interstellar dust by studying the extinction curves. Unlike the dust extinction templates, such as those of the average Milky Way (MW) or the Small and Large Magellanic Cloud (SMC and LMC), the extinction curves of galaxies outside the Local Group exhibit deviation from these laws. Altogether, X-ray and gamma-ray satellites as well as ground-based telescopes, such as Neil Gehrels Swift Observatory (Swift) and Gamma-Ray Optical and Near-Infrared Detector (GROND), provide measurements of the afterglows from the X-ray to the NIR, which can be used to extract information on dust extinction curves along their lines of sight (LoS). The study presented in this paper undertakes such a photometric study, comprising a preparatory work for the SVOM mission and its ground-based follow-up telescope COLIBRI. Aims. We propose a simple approach to parameterize the dust extinction curve of GRB host galaxies. The model used in this analysis is based on a power law form with the addition of a Loretzian-like Drude profile with two parameters: the extinction slope, γ, and the 2175 Å bump amplitude, Eb. Methods. Using the g′r′i′z′JHKs GROND filter bands, we tested our dust extinction model and explored the parameter space in extinction and redshift by fitting SEDs of simplified simulations of GRB afterglow spectra based on different extinction curve templates. From a final sample of 10 real Swift/GROND extinguished GRBs, we determined the quantities of the dust extinction in their host and measured their extinction curves. Results. We find that our derived extinction curves are in agreement with the spectroscopic measurements reported for four GRBs in the literature. We compared four other GRBs to the results of photometric studies where fixed laws were used to fit their data. We additionally derived two new GRB extinction curves. The measured average extinction curve is given by a slope of γ = 1.051 ± 0.129 and Eb = 0.070 ± 0.036, which is equivalent to a quasi-featureless in-between SMC-LMC template. This is consistent with previous studies aimed at deriving the dust host galaxy extinction where we expect that small dust grains dominate in GRB environment, yielding a steeper curve than the mean MW extinction curve.

Investigating the off-axis GRB afterglow scenario for extragalactic fast X-ray transients

Context. Extragalactic fast X-ray transients (FXTs) are short-duration (∼ks) X-ray flashes of unknown origin, potentially arising from binary neutron star (BNS) mergers, tidal disruption events, or supernova shock breakouts. Aims. In the context of the BNS scenario, we investigate the possible link between FXTs and the afterglows of off-axis merger-induced gamma-ray bursts (GRBs). Methods. By modelling well-sampled broadband afterglows of 13 merger-induced GRBs, we make predictions for their X-ray light curve behaviour had they been observed off-axis, considering both a uniform jet with core angle θC and a Gaussian-structured jet whose edge lies at an angle θW = 2θC. We compare their peak X-ray luminosity, duration, and temporal indices α (where F ∝ tα) with those of the currently known extragalactic FXTs. Results. Our analysis reveals that a slightly off-axis observing angle of θobs ≈ (2.2 − 3)θC and a structured jet are required to explain the shallow (|α|≲0.3) temporal indices of the FXT light curves, which cannot be reproduced in the uniform-jet case at any viewing angle. In the case of a structured jet with truncation angle θW = 2θC, the distributions of the duration of the FXTs are consistent with those of the off-axis afterglows for the same range of observing angles, θobs ≈ (2.2 − 3)θC. While the distributions of the off-axis peak X-ray luminosity are consistent only for θobs = 2.2θC, focussing on individual events with different intrinsic luminosities reveals that the match of all three properties (peak X-ray luminosity, duration and temporal indices) of the FXTs at the same viewing angle is possible in the range θobs ∼ (2.2 − 2.6)θC. Despite the small sample of GRBs analysed, these results show that there is a region of the parameter space – although quite limited – where the observational properties of off-axis GRB afterglow can be consistent with those of the newly discovered FXTs. Future observations of FXTs discovered by the recently launched Einstein Probe mission and GRB population studies combined with more complex afterglow models will shed light on this possible GRB-FXT connection, and eventually unveil the progenitors of some FXTs.

Photospheric Prompt Emission from Long Gamma Ray Burst Simulations. III. X-Ray Spectropolarimetry

While gamma-ray bursts (GRBs) have the potential to shed light on the astrophysics of jets, compact objects, and cosmology, a major set back in their use as probes of these phenomena stems from our incomplete knowledge surrounding their prompt emission. There are numerous models that can account for various observations of GRBs in the gamma-ray and X-ray energy ranges, due to the flexibility in the number of parameters that can be tuned to increase agreement with data. Furthermore, these models lack predictive power that can test future spectropolarimetric observations of GRBs across the electromagnetic spectrum. In this work, we use the MCRaT radiative transfer code to calculate the X-ray spectropolarimetric signatures expected from the photospheric model for two unique hydrodynamic simulations of long GRBs. We make time-resolved and time-integrated comparisons between the X-ray and gamma-ray mock observations, shedding light on the information that can be obtained from X-ray prompt emission signatures. Our results show that the T 90 derived from the X-ray light curve is the best diagnostic for the time that the central engine is active. We also find that our simulations reproduce the observed characteristics of the Einstein Probe–detected GRB 240315C. Based on our simulations, we are also able to make predictions for future X-ray spectropolarimetric measurements. Our results show the importance of conducting global radiative transfer calculations of GRB jets to better contextualize the prompt emission observations and constrain the mechanisms that produce the prompt emission.

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

Gamma-ray bursts (GRBs) are powered by ultrarelativistic jets. The launching sites of these jets are surrounded by dense media, which the jets must cross before they can accelerate and release high-energy emission. Interaction with the medium leads to the formation of a mildly relativistic sheath around the jet, resulting in an angular structure to 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.

Towards a new model-independent calibration of Gamma-Ray Bursts

Current data on baryon acoustic oscillations and Supernovae of Type Ia (SNIa) cover up to z∼2.5. These low-redshift observations play a very important role in the determination of cosmological parameters and have been widely used to constrain the ΛCDM and models beyond the standard, such as the ones with open curvature. To extend this investigation to higher redshifts, Gamma-Ray Bursts (GRBs) stand out as one of the most promising observables. In spite of being transient, they are extremely energetic and can be used to probe the universe up to z∼9.4. They exhibit characteristics that suggest they are potentially standardizable candles and this allows their use to extend the distance ladder beyond SNIa. The use of GRB correlations is still a challenge due to the spread in their intrinsic properties. One of the correlations that can be employed for the standardization is the fundamental plane relation between the peak prompt luminosity, the rest-frame end time of the plateau phase, and its corresponding luminosity, also known as the three-dimensional Dainotti correlation. In this work, we propose an innovative method of calibration of the Dainotti relation which is independent of cosmology. We employ state-of-the-art data on Cosmic Chronometers (CCH) at z≲2 and use the Gaussian Processes Bayesian reconstruction tool. To match the CCH redshift range, we select 20 long GRBs in the range 0.553≤z≤1.96 from the Platinum sample, which consists of well-defined GRB plateau properties that obey the fundamental plane relation. To ensure the generality of our method, we verify that the choice of priors on the parameters of the Dainotti relation and the modelling of CCH uncertainties and covariance have negligible impact on our results. Moreover, we consider the case in which the redshift evolution of the physical features of the plane is accounted for. We find that the use of CCH allows us to identify a sub-sample of GRBs that adhere even more closely to the fundamental plane relation, with an intrinsic scatter of σint=0.20−0.05+0.03 obtained in this analysis when evolutionary effects are considered. In an epoch in which we strive to reduce uncertainties on the variables of the GRB correlations in order to tighten constraints on cosmological parameters, we have found a novel model-independent approach to pinpoint a sub-sample that can thus represent a valuable set of standardizable candles. This allows us to extend the cosmic distance ladder presenting a new catalogue of calibrated luminosity distances up to z=5.

High energy neutrino production in gamma-ray bursts: dependence of the neutrino signal on the jet composition

Heavy nuclei can be synthetized or entrained in gamma-ray bursts (GRBs) with implications on the high-energy neutrino emission. By means of a Monte-Carlo algorithm, we model nuclear cascades and investigate their impact on the neutrino production considering kinetic dominated jets (in the internal shock model, including a dissipative photosphere) as well as Poynting flux dominated jets (for a jet model invoking internal-collision-induced magnetic reconnection and turbulence, ICMART). We find that the ICMART model allows for efficient nuclear cascades leading to an overall larger neutrino fluence than in the other two jet models. The survival of nuclei and inefficient nuclear cascades lead to an overall reduction of the neutrino fluence up to one order of magnitude. However, if nuclei are disintegrated, the neutrino fluence may be comparable to the one emitted from a jet loaded with protons. Exploring the parameter space of jet properties, we conclude that the composition and the bulk Lorentz factor have significant impact on the efficiency of nuclear cascades as well as the spectral shape of the expected neutrino fluence. On the other hand, the neutrino spectral distribution is less sensitive to the power-law index of the accelerated population of protons or heavier nuclei. For what concerns the diffuse emission of neutrinos from GRBs, we find that the uncertainty due to the jet composition can be at most comparable to the one related to the GRB cosmological rate.

High-energy Neutrinos from Late-time Jets of Gamma-Ray Bursts Seeded with Cocoon Photons

In gamma-ray bursts (GRBs), ∼100–1000 s after the prompt emission, afterglow observations have consistently shown X-ray excesses detected in the form of flares (in long GRBs) or extended emission (in short GRBs). These observations are interpreted as emissions from jets launched by late central engine activity. However, the characteristics of these late-time jets, particularly the dissipation radius (r diss), Lorentz factor (Γ), and cosmic-ray loading factor (ξ p ), remain unknown despite their importance. Here, in order to understand the properties of the late-time jets with future multimessenger observations, we estimate the detectability of neutrinos associated with late-time emissions for a wide range of r diss and Γ, assuming ξ p = 10. We take into account external seed photons from the cocoon around the jets, which can enhance the neutrino production through photohadronic interaction in the jet dissipation region. Our results are still consistent with the upper limit obtained by IceCube. Our calculations indicate a promising prospect for neutrino detection with IceCube-Gen2 through the stacking of ∼1000–2000 events, for a wide range of r diss and Γ. We found that setting an optimal energy threshold of 10 TeV can significantly reduce noise without negatively affecting neutrino detection. Furthermore, even in the case of nondetection, we show that meaningful constraints on the characteristics of the late-time jets can be obtained.

Multipolar Electromagnetic Emission of Newborn Magnetars

It is generally recognized that the electromagnetic multipolar emission from magnetars can be used to explain radiation from soft gamma repeaters or anomalous X-ray pulsars, but they have little impact on the spin-down of magnetars. We here present an analytical solution for the neutron star multipolar electromagnetic fields and their associated expected luminosities. We find that for newborn millisecond magnetars, the spin-down luminosity from higher multipolar components can match or even exceed that from the dipole component. Such high-intensity radiation will undoubtedly affect related astrophysical phenomena at the birth of a magnetar. We show that the spin-down luminosity from multipoles can well explain the majority of gamma-ray burst (GRB) afterglows, from the plateau starting at several hundred seconds until the normal decay phase lasting for many years. The fitted magnetar parameters for GRB afterglows are all typical values, with spins in the millisecond range and magnetic field strengths on the order of 1014–1015 G. Our results, in turn, provide support for the hypothesis that GRBs originate from the birth of magnetars with a period of a few milliseconds, thus deepening our understanding of the complex magnetic field structure and the equation of state of magnetars.

Formation of GW230529 from Isolated Binary Evolution

In this paper, we explore the formation of the mass-gap black hole-neutron star (mgBHNS) merger detected in gravitational wave (GW) event, i.e., GW230529, from the isolated binary evolution channel, and study potential signatures of its electromagnetic counterparts. By adopting the “delayed” supernova prescription and reasonable model realizations, our population synthesis simulation results can simultaneously match the rate densities of mgBHNS and total BHNS mergers inferred from the population analyses, along with the population distribution of the BH mass in BHNS mergers reported by the LIGO–Virgo–KAGRA Collaboration. Because GW230529 contributes significantly to the inferred mgBHNS rate densities, we suggest that GW230529 can be explained through the isolated binary evolution channel. Considering the AP4 (DD2) equation of state, the probability that GW230529 can make tidal disruption is 12.8% (63.2%). If GW230529 is a disrupted event, its kilonova peak apparent magnitude is predicted ∼23–24 mag, and hence, can be detected by the present survey projects and Large Synoptic Survey Telescope. Since GW230529 could be an off-axis event inferred from the GW observation, its associated gamma-ray burst (GRB) might be too dim to be observed by γ-ray detectors, interpreting the lack of GRB observations. Our study suggests the existence of mgBHNS mergers formed through the isolated binary evolution channel due to the discovery of GW230529, indicating that BHNS mergers are still likely to be multimessenger sources that emit GWs, GRBs, and kilonovae. Although mgBHNS mergers account for ∼50% of the cosmological BHNS population, we find that ≳90% of disrupted BHNS mergers are expected to originate from mgBHNS mergers.

Testing the standardizability of, and deriving cosmological constraints from, a new Amati-correlated gamma-ray burst data compilation

By using gamma-ray burst (GRB) data to simultaneously constrain Amati correlation parameters and cosmological parameters in six spatially flat and nonflat dark energy cosmological models, we show that an updated 220 GRB version of the Jia et al. [1] GRB data compilation are standardizable through the Amati correlation and so can be used for cosmological analyses. However, the resulting GRB data constraints on the current value of the nonrelativistic matter density parameter, Ω m 0, are in > 2σ tension with those from a joint analysis of better-established Hubble parameter [H(z)] and baryon acoustic oscillation (BAO) data for most of the cosmological models we consider, indicating that these GRB data cannot be jointly used with better-established H(z) + BAO data to constrain cosmological parameters.

Studies of GAGG:Ce scintillators for space missions dedicated to Terrestrial Gamma-Ray Flashes and Gamma-Ray Bursts observation

Fast Gamma ray Spectrometer (FGS) is an on-going detector development for the observation of transient gamma ray events such as Terrestrial Gamma ray Flashes (TGFs) produced in Earth’s thunderstorms and Gamma Ray Bursts (GRBs) produced in the far Universe. In this work, we compare different Cerium-doped Gadolinium Aluminium Gallium Garnet (GAGG:Ce) scintillators, associated with Silicon Photo-Multipliers (SiPMs), quantifying their performances at atmospheric pressure and under vacuum, with different surface finitions. We also study their afterglow following proton irradiation, similar to what will be observed in space environment. Our measurements with GAGG:Ce scintillators after proton irradiation confirm that GAGG:Ce can be used in space.

Disk mass after a binary neutron star merger as a constraining parameter for short gamma-ray bursts

Context. The coincident detection of GW170817 and gamma-ray burst GRB170817A marked a milestone for the connection between binary neutron star (BNS) mergers and short gamma-ray bursts (sGRBs). These mergers can lead to the formation of a black hole that is surrounded by a disk and to the generation of a powerful jet. It spends energy to break free from the merger ejecta, and then a portion of it is dissipated to produce observable emissions. Aims. Our primary goal is to enhance our comprehension of BNS mergers by constraining the disk mass for a selection of sGRBs. To do this, we used the isotropic gamma-ray luminosity and corresponding emission times as key indicators. Methods. We leveraged data from GW170817 to estimate the disk mass surrounding the BNS merger remnant, and we subsequently inferred the efficiency of the accretion onto the jet. We then statistically examined other sGRB observations to estimate whether they might have been induced by BNS mergers Results. Our findings suggest that when similar physical parameters are employed as in the only observed BNS-powered GRB event, GRB170817A, a substantial fraction of sGRBs would need an unrealistically massive disk remnant. Conclusions. This observation raises the possibility that either a different mechanism powered those events or that the post-collapse disk efficiency varies significantly in different BNS merger scenarios.

Host Galaxy Properties of Gamma-Ray Bursts Involving Neutron Star Binary Mergers and Their Impact on Kilonovae Rates

In the upcoming gravitational-wave (GW) observing runs, identifying host galaxies is crucial as it provides essential redshift information and enables the use of GW events as standard sirens. However, pinpointing host galaxies remains challenging due to the large localization uncertainties and the rapidly fading nature of their optical counterparts. Analyzing the host galaxies of short gamma-ray bursts (sGRBs) offers an alternative approach to deepen our understanding of the environments where binary neutron stars (BNS) primarily merge. This study compiles archival photometric data for the host galaxies of 76 sGRBs and four hybrid GRBs that are long GRBs accompanied by kilonova-like signals. We use this data to evaluate their physical properties through spectral energy distribution fitting. To assess the characteristics of the host galaxies, we utilized a volume-limited sample (z < 0.5) from the COSMOS field as a control group. Contrary to expectations that the BNS merger rate is proportional to host stellar mass, the short and hybrid GRB population appears less massive than the mass-weighted distribution of the control sample. Instead, we propose a formulation for the expected BNS merger rate from a galaxy as log(nBNSGyr-1)=0.86×log(M*/M⊙) + 0.44×log(sSFRyr-1)+0.857 , which optimally explains the deviation between the stellar mass distributions of the GRB host galaxies and the control sample. These insights provide a strategic framework for targeted GW follow-ups and enhance our ability to identify potential host galaxies for future GW events.

Late Jets, Early Sparks: Illuminating the Premaximum Bumps in Superluminous Supernovae

Superluminous supernovae (SLSNe) radiate ≳10–100 times more energy than ordinary stellar explosions, implicating a novel power source behind these enigmatic events. One frequently discussed source, particularly for hydrogen-poor (Type I) SLSNe, is a central engine such as a millisecond magnetar or accreting black hole. Both black hole and magnetar engines are expected to channel a fraction of their luminosity into a collimated relativistic jet. Using 3D relativistic hydrodynamical simulations, we explore the interaction of a relativistic jet, endowed with a luminosity L j ≈ 1045.5 erg s−1 and duration t eng ≈ 10 days compatible with those needed to power SLSNe, launched into the envelope of the exploding star. The jet successfully breaks through the expanding ejecta, and its shocked cocoon powers ultraviolet/optical emission lasting several days after the explosion and reaching a peak luminosity ≳1044 erg s−1, corresponding to a sizable fraction of L j . This high radiative efficiency is the result of the modest adiabatic losses the cocoon experiences owing to the low optical depths of the enlarged ejecta at these late times, e.g., compared to the more compact stars in gamma-ray bursts. The luminosity and temperature of the cocoon emission match those of the “bumps” in SLSN light curves observed weeks prior to the optical maximum in many SLSNe. Confirmation of jet breakout signatures by future observations (e.g., days-long to weeks-long internal X-ray emission from the jet for on-axis observers, spectroscopy confirming large photosphere velocities v/c ≳ 0.1, or detection of a radio afterglow) would offer strong evidence for central engines powering SLSNe.