Search for: All records

The origin of short gamma-ray bursts is associated with outflows powered by the remnant of a binary neutron star merger. This remnant can be either a black hole or a highly magnetized, fast-spinning neutron star, also known as a magnetar. Here we present the results of two relativistic magnetohydrodynamical simulations aimed at investigating the large-scale dynamics and propagation of magnetar collimated outflows through the medium surrounding the remnant. The first simulation evolves a realistic jet by injecting external simulation data, while the second evolves an analytical model jet with similar properties for comparison. We find that both outflows remain collimated and successfully emerge through the static medium surrounding the remnant. However, they fail to attain relativistic velocities and only reach a mean maximum speed of ∼0.7cfor the realistic jet and ∼0.6cfor the analytical jet. We also find that the realistic jet has a much more complex structure. The lack of highly relativistic speeds, which makes these jets unsuitable as short gamma-ray burst sources, is due to numerical limitations and is not general to all possible magnetar outflows. A jet like the one we study, however, could give rise to or augment a blue kilonova component. In addition, it would make the propagation of a relativistic jet easier, should one be launched after the neutron star collapses into a black hole.

Long and short gamma-ray bursts (GRBs), canonically separated at around 2 s duration, are associated with different progenitors: the collapse of a massive star and the merger of two compact objects, respectively. GRB 191019A was a long GRB (T₉₀∼ 64 s). Despite the relatively small redshiftz= 0.248 and Hubble Space Telescope follow-up observations, an accompanying supernova was not detected. In addition, the host galaxy did not have significant star formation activity. Here we propose that GRB 191019A was produced by a binary compact merger, whose prompt emission was stretched in time by the interaction with a dense external medium. This would be expected if the burst progenitor was located in the disk of an active galactic nucleus, as supported by the burst localization close to the center of its host galaxy. We show that the light curve of GRB 191019A can be well modeled by a burst of intrinsic durationt_eng= 1.1 s and of energyE_iso= 10⁵¹erg seen moderately off axis, exploding in a medium of density ∼10⁷–10⁸cm⁻³. The double-peaked light curve carries the telltale features predicted for GRBs in high-density media, where the first peak is produced by the photosphere and the second by the overlap of reverse shocks that take place before the internal shocks could happen. This would make GRB 191019A the first confirmed stellar explosion from within an accretion disk, with important implications for the formation and evolution of stars in accretion flows and for gravitational-waves source populations.

Both long and short gamma-ray bursts (GRBs) are expected to occur in the dense environments of active galactic nucleus (AGN) accretion discs. As these bursts propagate through the discs they live in, they photoionize the medium causing time-dependent opacity that results in transients with unique spectral evolution. In this paper, we use a line-of-sight radiation transfer code coupling metal and dust evolution to simulate the time-dependent absorption that occurs in the case of both long and short GRBs. Through these simulations, we investigate the parameter space in which dense environments leave a potentially observable imprint on the bursts. Our numerical investigation reveals that time-dependent spectral evolution is expected for central supermassive black hole masses between 105 and 5 × 107 solar masses in the case of long GRBs, and between 104 and 107 solar masses in the case of short GRBs. Our findings can lead to the identification of bursts exploding in AGN disc environments through their unique spectral evolution coupled with a central location. In addition, the study of the time-dependent evolution would allow for studying the disc structure, once the identification with an AGN has been established. Finally, our findings lead to insight into whether GRBs contribute to the AGN emission, and which kind, thus helping to answer the question of whether GRBs can be the cause of some of the as-of-yet unexplained AGN time variability.

The merger of two neutron stars (NSs) produces the emission of gravitational waves, the formation of a compact object surrounded by a dense and magnetized environment. If the binary undergoes delayed collapse a collimated and relativistic jet, which will eventually produce a short gamma-ray burst (SGRB), may be launched. The interaction of the jet with the environment has been shown to play a major role in shaping the structure of the outflow that eventually powers the gamma-ray emission. In this paper, we present a set of 2.5D RMHD simulations that follow the evolution of a relativistic non-magnetized jet through a medium with different magnetization levels, as produced after the merger of two NSs. We find that the predominant consequence of a magnetized ambient medium is that of suppressing instabilities within the jet and preventing the formation of a series of collimation shocks. One implication of this is that internal shocks lose efficiency, causing bursts with low-luminosity prompt emission. On the other hand, the jet-head velocity and the induced magnetization within the jet are fairly independent of the magnetization of the ambient medium. Future numerical studies with a larger domain are necessary to obtain light curves and spectra in order to better understand the role of magnetized media.

Long and shortγ-ray bursts (GRBs) are traditionally associated with galactic environments, where circumburst densities are small or moderate (few to hundreds of protons per cubic centimeter). However, both are also expected to occur in the disks of active galactic nuclei, where the ambient medium density can be much larger. In this work we study, via semianalytical methods, the propagation of the GRB outflow, its interaction with the external material, and the ensuing prompt radiation. In particular, we focus on the case in which the external shock develops early in the evolution at a radius that is smaller than the internal shock one. We find that bursts in such high-density environments are likely characterized by a single, long emission episode that is due to the superposition of individual pulses, with a characteristic hard-to-soft evolution irrespective of the light-curve luminosity. While multipulse light curves are not impossible, they would require the central engine to go dormant for a long time before reigniting. In addition, short GRB engines would produce bursts with prompt duration that would exceed the canonical 2 s separation threshold and likely be incorrectly classified as long events, even though they would not be accompanied by a simultaneous supernova. Finally, these events have a large dynamical efficiency, which would produce a bright prompt emission followed by a somewhat dim afterglow.

The discs of active galactic nuclei (AGNs) have emerged as rich environments for the production and capture of stars and the compact objects that they leave behind. These stars produce long gamma-ray bursts (GRBs) at their deaths, while frequent interactions among compact objects form binary neutron stars and neutron star–black hole binaries, leading to short GRBs upon their merger. Predicting the properties of these transients as they emerge from the dense environments of AGN discs is key to their proper identification and to better constrain the star and compact object population in AGN discs. Some of these transients would appear unusual because they take place in much higher densities than the interstellar medium. Others, which are the subject of this paper, would additionally be modified by radiation diffusion, since they are generated within optically thick regions of the accretion discs. Here, we compute the GRB afterglow light curves for diffused GRB sources for a representative variety of central black hole masses and disc locations. We find that the radiation from radio to ultraviolet and soft X-rays can be strongly suppressed by synchrotron self-absorption in the dense medium of the AGN disc. In addition, photon diffusion can significantly delay the emergence of the emission peak, turning a beamed, fast transient into a slow, isotropic, and dimmer one. These would appear as broad-band correlated AGN variability with a dominance at the higher frequencies. Their properties can constrain both the stellar populations within AGN discs and the disc structure.

The association of GRB170817A with a binary neutron star (BNS) merger has revealed that BNSs produce at least a fraction of short gamma-ray bursts (SGRBs). As gravitational wave (GW) detectors push their horizons, it is important to assess coupled electromagnetic (EM)/GW probabilities and maximize observational prospects. Here, we perform BNS population synthesis calculations with the code mobse, seeding the binaries in galaxies at three representative redshifts, $z$ = 0.01, 0.1, and 1 of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the microphysics parameters of GRB170817A, we numerically compute the broad-band light curves of jets from BNS mergers, with the afterglow brightness dependent on the local medium density at the merger site. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find a gamma-ray detection probability of $\sim\!2{{\rm per\ cent}},10{{\rm per\ cent}},\mathrm{and}\ 40{{\rm per\ cent}}$ for BNSs at $z$ = 1, 0.1, and 0.01, respectively, for the random case, rising to $\sim\!75{{\rm per\ cent}}$ for the $z$ = 0.01, GW-triggered aligned case. Afterglow detection probabilities of GW-triggered BNS mergers vary in the range of $\sim \! 0.3 \!-\! 0.5{{\rm per\ cent}}$, with higher values for aligned jets, and are comparable across the high- and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses.

Although gamma ray bursts (GRBs) have been detected for many decades, the lack of knowledge regarding the radiation mechanism that produces the energetic flash of radiation, or prompt emission, from these events has prevented the full use of GRBs as probes of high-energy astrophysical processes. While there are multiple models that attempt to describe the prompt emission, each model can be tuned to account for observed GRB characteristics in the gamma and X-ray energy bands. One energy range that has not been fully explored for the purpose of prompt emission model comparison is that of the optical band, especially with regard to polarization. Here, we use an improved Monte Carlo radiation transfer code to calculate the expected photospheric optical and gamma-ray polarization signatures (Π_optand Π_γ, respectively) from a set of two relativistic hydrodynamic long GRB simulations, which emulate a constant and variable jet. We find that time-resolved Π_optcan be large (∼75%) while time-integrated Π_optcan be smaller due to integration over the asymmetries in the GRB jet where optical photons originate; Π_γfollows a similar evolution as Π_optwith smaller polarization degrees. We also show that Π_optand Π_γagree well with observations in each energy range. Additionally, we make predictions for the expected polarization of GRBs based on their location within the Yonetoku relationship. While improvements can be made to our analyses and predictions, they exhibit the insight that global radiative transfer simulations of GRB jets can provide with respect to current and future observations.

A complete understanding of gamma-ray bursts (GRBs) has been difficult to achieve, due to our incomplete knowledge of the radiation mechanism that is responsible for producing the prompt emission. This emission, which is detected in the first tens of seconds of the GRB, is typically dominated by hard X-ray and gamma-ray photons, although there have also been a few dozen prompt optical detections. These optical detections have the potential to discriminate between plausible prompt emission models, such as the photospheric and synchrotron shock models. In this work, we use an improved MCRaT code, which includes cyclo-synchrotron emission and absorption, to conduct radiative transfer calculations from optical to gamma-ray energies under the photospheric model. The calculations are conducted using a set of two-dimensional relativistic hydrodynamic long GRB jet simulations, consisting of a constant and a variable jet. We predict the correlations between the optical and gamma-ray light curves as functions of observer angle and jet variability, and find that there should be extremely dim optical prompt precursors for large viewing angles. Additionally, the detected optical emission originates from dense regions of the outflow, such as shock interfaces and the jet-cocoon interface. Our results also show that the photospheric model is unable to account for the current set of optical prompt detections that have been made and therefore additional radiative mechanisms are needed to explain these prompt optical observations. These findings show the importance of conducting global radiative transfer simulations using hydrodynamically calculated jet structures.

The association of GRB170817A with GW170817 has confirmed the long-standing hypothesis that binary neutron star (BNS) mergers are the progenitors of at least some short gamma-ray bursts (SGRBs). This connection has ushered in an era in which broadband observations of SGRBs, together with measurements of the time delay between the gravitational waves and the electromagnetic radiation, allow for probing the properties of the emitting outflow and its engine to an unprecedented detail. Because the structure of the radiating outflow is molded by the interaction of a relativistic jet with the binary ejecta, it is of paramount importance to study the system in a realistic setting. Here we present a three-dimensional hydrodynamic simulation of a relativistic jet propagating in the ejecta of a BNS merger, which were computed with a general relativistic magnetohydrodynamic simulation. We find that the jet’s centroid oscillates around the axis of the system, due to inhomogeneities encountered in the propagation. These oscillations allow the jet to find the path of least resistance and travel faster than an identical jet in smooth ejecta. In our setup the breakout time is ∼0.6 s, which is comparable to the expected central engine duration in SGRBs and possibly a non-negligible fraction of the total delay between the gravitational and gamma-ray signals. Our simulation also shows that energy is carried in roughly equal amounts by the jet and by the cocoon, and that about 20% of the injected energy is transferred to the ejecta via mechanical work.