Search

This study presents a semi-analytical model that models the effects that a static magnetized medium with a tangled field produces in relativistic collimated jets. The model is a first approximation that addresses the magnetic field present in the medium and is based on pressure equilibrium principles between the jet, cocoon, and external medium. A fraction of the ambient medium field is allowed to be entrained in the cocoon. We find that the jet and cocoon properties are affected at high magnetic fields and that they are sensitive to the mixing in the cocoon. For low-mixing a slower-broader jet with a broader and more energetic cocoon would be produced. On the other hand, high-mixing would produce a faster-narrower jet with a narrow and less-energetic cocoon. Two-dimensional hydrodynamical simulations are used to validate the model and to constrain the mixing parameter. We found good qualitative agreement between the model and the simulations. For high magnetization, the results were found to be consistent with the low mixing case of our semi-analytical model., Comment: 9 pages, 6 figures, 1 table, sent to MNRAS

The majority of massive stars are members of binary systems, a fraction of which is expected to remain bound after the first star goes off as a supernova. When the second star also explodes, the SN ejecta are bound to interact with the compact object companion. We explore the consequences of this interaction in the case of a black hole (BH) companion. We show that accretion of the SN ejecta by the BH generally occurs with highly super-Eddington rates, which can be strongly modulated if the ejecta are clumpy, as typically observed in supernova remnants. This late accretion produces transient flaring activity with a time delay from the SN explosion of days to months, depending on the orbital separation and the velocity of the ejecta. The flares are expected to have a non-thermal, broad band spectrum, but their high-frequency emission (UV and X-rays) would be absorbed within the remnant. Observed flares should therefore be in the optical and near-infrared. We propose this model as an explanation to the late-time flaring activity observed in the fast blue optical transient AT2022tsd., Comment: Submitted to ApJL, comments most welcome

Long-duration gamma-ray bursts (LGRBs), thought to be produced during core-collapse supernovae, may have a prominent neutron component in the outflow material. If present, neutrons can change how photons scatter in the outflow by reducing its opacity, thereby allowing the photons to decouple sooner than if there were no neutrons present. Understanding the details of this process could therefore allow us to probe the central engine of LGRBs, which is otherwise hidden. Here, we present results of the photospheric emission from an LGRB jet, using a combination of relativistic hydrodynamic simulations and radiative transfer post-processing using the Monte Carlo Radiation Transfer (MCRaT) code. We control the size of the neutron component in the jet material by varying the equilibrium electron fraction $Y_{e}$, and we find that the presence of neutrons in the GRB fireball affects the Band parameters $\alpha$ and $E_{0}$, while the picture with the $\beta$ parameter is less clear. In particular, the break energy $E_{0}$ is shifted to higher energies. Additionally, we find that increasing the size of the neutron component also increases the total radiated energy of the outflow across multiple viewing angles. Our results not only shed light on LGRBs, but are also relevant to short-duration gamma-ray bursts associated with binary neutron star mergers, due to the likelihood of a prominent neutron component in such systems.

Inverse reinforcement learning (IRL) aims to recover the reward function of an expert agent from demonstrations of behavior. It is well-known that the IRL problem is fundamentally ill-posed, i.e., many reward functions can explain the demonstrations. For this reason, IRL has been recently reframed in terms of estimating the feasible reward set (Metelli et al., 2021), thus, postponing the selection of a single reward. However, so far, the available formulations and algorithmic solutions have been proposed and analyzed mainly for the online setting, where the learner can interact with the environment and query the expert at will. This is clearly unrealistic in most practical applications, where the availability of an offline dataset is a much more common scenario. In this paper, we introduce a novel notion of feasible reward set capturing the opportunities and limitations of the offline setting and we analyze the complexity of its estimation. This requires the introduction an original learning framework that copes with the intrinsic difficulty of the setting, for which the data coverage is not under control. Then, we propose two computationally and statistically efficient algorithms, IRLO and PIRLO, for addressing the problem. In particular, the latter adopts a specific form of pessimism to enforce the novel desirable property of inclusion monotonicity of the delivered feasible set. With this work, we aim to provide a panorama of the challenges of the offline IRL problem and how they can be fruitfully addressed., Comment: International Conference on Machine Learning 41 (ICML 2024)

In online Inverse Reinforcement Learning (IRL), the learner can collect samples about the dynamics of the environment to improve its estimate of the reward function. Since IRL suffers from identifiability issues, many theoretical works on online IRL focus on estimating the entire set of rewards that explain the demonstrations, named the feasible reward set. However, none of the algorithms available in the literature can scale to problems with large state spaces. In this paper, we focus on the online IRL problem in Linear Markov Decision Processes (MDPs). We show that the structure offered by Linear MDPs is not sufficient for efficiently estimating the feasible set when the state space is large. As a consequence, we introduce the novel framework of rewards compatibility, which generalizes the notion of feasible set, and we develop CATY-IRL, a sample efficient algorithm whose complexity is independent of the cardinality of the state space in Linear MDPs. When restricted to the tabular setting, we demonstrate that CATY-IRL is minimax optimal up to logarithmic factors. As a by-product, we show that Reward-Free Exploration (RFE) enjoys the same worst-case rate, improving over the state-of-the-art lower bound. Finally, we devise a unifying framework for IRL and RFE that may be of independent interest.

In 2022, researchers from the University of Melbourne undertook a project exploring how Australia-based migrants (diasporas) help in times of humanitarian crises overseas. The project involved community researchers from eight diaspora communities (Afghanistan, Indonesia, Lebanon, Myanmar, Nepal, Pacific Islands, South Sudan and Syria) to co-design methodology, engage communities and collect data to find out how, why and what was done by communities in Australia in response to specific crises or events. The project aims to better understand the strengths and challenges faced by Australia-based migrants responding to different kinds of crises (disasters, conflicts and complex crises), and to identify potential tools that can support diaspora communities in their responses in the future. More information about this project can be found at https://diasporahumanitarians. com/.

In 2022, researchers from the University of Melbourne and Australian National University undertook a project exploring how Australia-based migrants (diasporas) help in times of humanitarian crises overseas. The project involved community researchers from eight diaspora communities (Afghanistan, Indonesia, Lebanon, Myanmar, Nepal, Pacific Islands, South Sudan and Syria) to co-design methodology, engage communities and collect data to find out how, why and what was done by communities in Australia in response to specific crises or events. The project aims to better understand the strengths and challenges faced by Australia-based migrants responding to different kinds of crises (disasters, conflicts and complex crises), and to identify potential tools that can support diaspora communities in their responses in the future. More information about this project can be found at https://diasporahumanitarians.com/.

In 2022, researchers from the University of Melbourne undertook a project exploring how Australia-based migrants (diasporas) help in times of humanitarian crises overseas. The project involved community researchers from eight diaspora communities (Afghanistan, Indonesia, Lebanon, Myanmar, Nepal, Pacific Islands, South Sudan and Syria) to co-design methodology, engage communities and collect data to find out how, why and what was done by communities in Australia in response to specific crises or events. The project aims to better understand the strengths and challenges faced by Australia-based migrants responding to different kinds of crises (disasters, conflicts and complex crises), and to identify potential tools that can support diaspora communities in their responses in the future. More information about this project can be found at https://diasporahumanitarians. com/.

Contains fulltext : 294450.pdf (Publisher’s version ) (Open Access)

Long and short gamma-ray bursts (GRBs), canonically separated at around 2 seconds 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_{90}\sim64$ s). Despite the relatively small redshift z=0.248 and HST followup 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 duration t=1.1 s and of energy $E_{\rm{iso}}=10^{51}$ erg seen moderately off-axis, exploding in a medium of density $10^7-10^8$ cm$^{-3}$. 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., Comment: 6 pages, 2 figures. Accepted for publication in ApJL

Inverse reinforcement learning (IRL) denotes a powerful family of algorithms for recovering a reward function justifying the behavior demonstrated by an expert agent. A well-known limitation of IRL is the ambiguity in the choice of the reward function, due to the existence of multiple rewards that explain the observed behavior. This limitation has been recently circumvented by formulating IRL as the problem of estimating the feasible reward set, i.e., the region of the rewards compatible with the expert's behavior. In this paper, we make a step towards closing the theory gap of IRL in the case of finite-horizon problems with a generative model. We start by formally introducing the problem of estimating the feasible reward set, the corresponding PAC requirement, and discussing the properties of particular classes of rewards. Then, we provide the first minimax lower bound on the sample complexity for the problem of estimating the feasible reward set of order ${\Omega}\Bigl( \frac{H^3SA}{\epsilon^2} \bigl( \log \bigl(\frac{1}{\delta}\bigl) + S \bigl)\Bigl)$, being $S$ and $A$ the number of states and actions respectively, $H$ the horizon, $\epsilon$ the desired accuracy, and $\delta$ the confidence. We analyze the sample complexity of a uniform sampling strategy (US-IRL), proving a matching upper bound up to logarithmic factors. Finally, we outline several open questions in IRL and propose future research directions., Comment: Submitted to: ICML23

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.

Both long and short gamma-ray bursts (GRBs) are expected to occur in the dense environments of active galactic nuclei (AGN) accretion disks. As these bursts propagate through the disks 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 $10^5$ and $5 \times 10^7$ solar masses in the case of long GRBs, and between $10^4$ and $10^7$ solar masses in the case of short GRBs. Our findings can lead to the identification of bursts exploding in AGN disk 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 disk 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, and the launching of a collimated and relativistic jet, which will eventually produce a short gamma-ray burst (SGRB). 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.5 dimensional 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 from 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., Comment: 7 pages, 5 figures, 1 table. Submitted to MNRAS

Long and short gamma-ray bursts are traditionally associated with galactic environments, where circumburst densities are small or moderate (few to hundreds of protons per cubic cm). 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 semi-analytical 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 multi-pulse light curves are not impossible, they would require the central engine to go dormant for a long time before re-igniting. In addition, short GRB engines would produce bursts with prompt duration that would exceed the canonical 2 s separation threshold and would 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., Comment: 9 pages, 5 figures, under review for AAS Journals

The origin of short gamma-ray bursts (sGRBs) 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, fastly spinning neutron star, also known as a magnetar. Here, we present the results of two relativistic magnetohydrodynamical (RMHD) 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.7c for the realistic jet, and ~0.6c for the analytical jet. We also find that the realistic jet has a much more complex structure. The lack of highly relativistic speeds, that makes these jets unsuitable as short GRB sources, is due to numerical limitations and 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., Comment: 11 pages, 7 figures, accepted to ApJ

The disks 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 (LGRBs) at their deaths, while frequent interactions among compact objects form binary neutron stars and neutron star-black hole binaries, leading to short GRBs (SGRBs) upon their merger. Predicting the properties of these transients as they emerge from the dense environments of AGN disks is key to their proper identification and to better constrain the star and compact object population in AGN disks. 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 disks. Here we compute the GRB afterglow light curves for diffused GRB sources for a representative variety of central black-hole masses and disk locations. We find that the radiation from radio to UV and soft X-rays can be strongly suppressed by synchrotron self-absorption in the dense medium of the AGN disk. 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 broadband-correlated AGN variability with dominance at the higher frequencies. Their properties can constrain both the stellar populations within AGN disks as well as the disk structure., Comment: Accepted for publication in MNRAS

GW170817 is the first binary neutron star (NS) merger detected in gravitational waves (GWs) and photons, and so far remains the only GW event of its class with a definitive electromagnetic (EM) counterpart. Radio emission from the structured jet associated with GW170817 has faded below the sensitivity achievable via deep radio observations with the most sensitive radio arrays currently in operation. Hence, we now have the opportunity to probe the radio re-brightening that some models predict, should emerge at late times from the interaction of the dynamically-stripped merger ejecta with the interstellar medium. Here we present the latest results from our deep radio observations of the GW170817 field with the Karl G. Jansky Very Large Array (VLA), 4.5 years after the merger. Our new data at $3\,$GHz do not show any compelling evidence for emission in excess to the tail of the jet afterglow ($<3.3\,\mu$Jy), confirming our previous results. We thus set new constraints on the dynamical ejecta afterglow models. These constraints favor single-speed ejecta with energy $\lesssim 10^{50}\,$erg (for an ejecta speed of $\beta_0=0.5$), or steeper energy-speed distributions of the kilonova ejecta. Our results also suggest larger values of the cold, non-rotating maximum NS mass in equal mass scenarios. However, without a detection of the dynamical ejecta afterglow, obtaining precise constraints on the NS equation of state remains challenging., Comment: 6 pages, 5 figures. Published in ApJ

Background: There is a paucity of data comparing 30-day morbidity and mortality of sleeve gastrectomy (SG), Roux-en-Y gastric bypass (RYGB), and one anastomosis gastric bypass (OAGB). This study aimed to compare the 30-day safety of SG, RYGB, and OAGB in propensity score-matched cohorts. Materials and methods: This analysis utilised data collected from the GENEVA study which was a multicentre observational cohort study of bariatric and metabolic surgery (BMS) in 185 centres across 42 countries between 01/05/2022 and 31/10/2020 during the Coronavirus Disease-2019 (COVID-19) pandemic. 30-day complications were categorised according to the Clavien–Dindo classification. Patients receiving SG, RYGB, or OAGB were propensity-matched according to baseline characteristics and 30-day complications were compared between groups. Results: In total, 6770 patients (SG 3983; OAGB 702; RYGB 2085) were included in this analysis. Prior to matching, RYGB was associated with highest 30-day complication rate (SG 5.8%; OAGB 7.5%; RYGB 8.0% (p = 0.006)). On multivariate regression modelling, Insulin-dependent type 2 diabetes mellitus and hypercholesterolaemia were associated with increased 30-day complications. Being a non-smoker was associated with reduced complication rates. When compared to SG as a reference category, RYGB, but not OAGB, was associated with an increased rate of 30-day complications. A total of 702 pairs of SG and OAGB were propensity score-matched. The complication rate in the SG group was 7.3% (n = 51) as compared to 7.5% (n = 53) in the OAGB group (p = 0.68). Similarly, 2085 pairs of SG and RYGB were propensity score-matched. The complication rate in the SG group was 6.1% (n = 127) as compared to 7.9% (n = 166) in the RYGB group (p = 0.09). And, 702 pairs of OAGB and RYGB were matched. The complication rate in both groups was the same at 7.5 % (n = 53; p = 0.07). Conclusions: This global study found no significant difference in the 30-day morbidity and mortality of SG

The detection of electromagnetic radiation (EM) accompanying the gravitational wave (GW) signal from the binary neutron star (BNS) merger GW170817 has revealed that these systems constitute at least a fraction of the progenitors of short gamma-ray bursts (SGRBs). As gravitational wave detectors keep pushing their detection horizons, it is important to assess coupled GW/EM probabilities, and how to maximize observational prospects. Here we perform population synthesis calculations of BNS evolution with the code MOBSE, and seed the binaries in galaxies at three representative redshifts (z=0.01,0.1,1) of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the astrophysical parameters of GRB170817A as a prototype, we numerically compute the broadband lightcurves of jets from BNS mergers, with the afterglow brightness depending on the local medium density at the merger sites. 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 that ~70-80% of BNSs from z=0.01 should be detectable in gamma-rays. The afterglow detection probabilities of GW-triggered BNS mergers vary between ~0.3-0.7%, with higher values for jets aligned with the BNS angular momentum, 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., Comment: 15 pages, 12 figures, submitted to MNRAS

The study of Gamma Ray Bursts (GRBs) has the potential to improve our understanding of high energy astrophysical phenomena. In order to reliably use GRBs to this end, we first need to have a well-developed grasp of the mechanism that produces the radiation within GRB jets and how that relates to their structure. One model for the emission mechanism of GRBs invokes radiation produced deep in the jet which eventually escapes the jet at its photosphere. While this model has been able to explain a number of observed GRB characteristics, it is currently lacking in predictive power and in ability to fully reproduce GRB spectra. In order to address these shortcomings of the model, we have expanded the capabilities of the MCRaT code, a state of the art radiative transfer code that can now simulate optical to gamma ray radiation propagating in a hydrodynamically simulated GRB jet. Using the MCRaT code, we have constructed mock observed light curves, spectra, and polarization from optical to gamma ray energies for the simulated GRBs. Using these mock observables, we have compared our simulations of photospheric emission to observations and found much agreement between the two. Furthermore, the MCRaT calculations combined with the hydrodynamical simulations allow us to connect the mock observables to the structure of the simulated GRB jet in a way that was not previously possible. While there are a number of improvements that can be made to the analyses, the steps taken here begin to pave the way for us to fully understand the connection between the structure of a given GRB jet and the radiation that would be expected from it., Comment: 10 pages, 3 figures, accepted as a conference proceeding to the Sixteenth Marcel Grossmann Meeting

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 regards to polarization. Here, we use an improved MCRaT code to calculate the expected photospheric optical and gamma-ray polarization signatures ($\Pi_\mathrm{opt}$ and $\Pi_\gamma$, respectively) from a set of two relativistic hydrodynamic long GRB simulations, which emulate a constant and variable jet. We find that time resolved $\Pi_\mathrm{opt}$ can be large ($\sim 75\%$) while time-integrated $\Pi_\mathrm{opt}$ can be smaller due to integration over the asymmetries in the GRB jet where optical photons originate; $\Pi_\gamma$ follows a similar evolution as $\Pi_\mathrm{opt}$ with smaller polarization degrees. We also show that $\Pi_\mathrm{opt}$ and $\Pi_\gamma$ 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., Comment: 16 pages; 7 figures including 2 animations (animations can be obtained from: https://github.com/parsotat/SpectropolarimetryAnimations); accepted to ApJ; the MCRaT and ProcessMCRaT codes are open source and available at: github.com/lazzati-astro/MCRaT and at: github.com/parsotat/ProcessMCRaT

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 to probe the properties of the emitting outflow and its engine to an unprecedented detail. Since 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 sec, 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 per cent of the injected energy is transferred to the ejecta via mechanical work., Comment: Under review for publication in ApJL

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 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 not able to account for the current set of optical prompt detections that have been made and 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., Comment: 22 pages, with 5 pages for the appendix; 11 figures; accepted to ApJ; the MCRaT code is open source and available at: https://github.com/lazzati-astro/MCRaT and the ProcessMCRaT code is also open source and available at: https://github.com/parsotat/ProcessMCRaT

The detection of electromagnetic radiation (EM) accompanying the gravitational wave (GW) signal from the binary neutron star (BNS) merger GW170817 has revealed that these systems constitute at least a fraction of the progenitors of short gamma-ray bursts (SGRBs). As gravitational wave detectors keep pushing their detection horizons, it is important to assess coupled GW/EM probabilities, and how to maximize observational prospects. Here we perform population synthesis calculations of BNS evolution with the code MOBSE, and seed the binaries in galaxies at three representative redshifts (z=0.01,0.1,1) of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the astrophysical parameters of GRB170817A as a prototype, we numerically compute the broadband lightcurves of jets from BNS mergers, with the afterglow brightness depending on the local medium density at the merger sites. 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 that ~70-80% of BNSs from z=0.01 should be detectable in gamma-rays. The afterglow detection probabilities of GW-triggered BNS mergers vary between ~0.3-0.7%, with higher values for jets aligned with the BNS angular momentum, 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., Comment: 15 pages, 12 figures, submitted to MNRAS

The study of Gamma Ray Bursts (GRBs) has the potential to improve our understanding of high energy astrophysical phenomena. In order to reliably use GRBs to this end, we first need to have a well-developed grasp of the mechanism that produces the radiation within GRB jets and how that relates to their structure. One model for the emission mechanism of GRBs invokes radiation produced deep in the jet which eventually escapes the jet at its photosphere. While this model has been able to explain a number of observed GRB characteristics, it is currently lacking in predictive power and in ability to fully reproduce GRB spectra. In order to address these shortcomings of the model, we have expanded the capabilities of the MCRaT code, a state of the art radiative transfer code that can now simulate optical to gamma ray radiation propagating in a hydrodynamically simulated GRB jet. Using the MCRaT code, we have constructed mock observed light curves, spectra, and polarization from optical to gamma ray energies for the simulated GRBs. Using these mock observables, we have compared our simulations of photospheric emission to observations and found much agreement between the two. Furthermore, the MCRaT calculations combined with the hydrodynamical simulations allow us to connect the mock observables to the structure of the simulated GRB jet in a way that was not previously possible. While there are a number of improvements that can be made to the analyses, the steps taken here begin to pave the way for us to fully understand the connection between the structure of a given GRB jet and the radiation that would be expected from it., Comment: 10 pages, 3 figures, accepted as a conference proceeding to the Sixteenth Marcel Grossmann Meeting

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 regards to polarization. Here, we use an improved MCRaT code to calculate the expected photospheric optical and gamma-ray polarization signatures ($\Pi_\mathrm{opt}$ and $\Pi_\gamma$, respectively) from a set of two relativistic hydrodynamic long GRB simulations, which emulate a constant and variable jet. We find that time resolved $\Pi_\mathrm{opt}$ can be large ($\sim 75\%$) while time-integrated $\Pi_\mathrm{opt}$ can be smaller due to integration over the asymmetries in the GRB jet where optical photons originate; $\Pi_\gamma$ follows a similar evolution as $\Pi_\mathrm{opt}$ with smaller polarization degrees. We also show that $\Pi_\mathrm{opt}$ and $\Pi_\gamma$ 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., Comment: 16 pages; 7 figures including 2 animations (animations can be obtained from: https://github.com/parsotat/SpectropolarimetryAnimations); accepted to ApJ; the MCRaT and ProcessMCRaT codes are open source and available at: github.com/lazzati-astro/MCRaT and at: github.com/parsotat/ProcessMCRaT

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 to probe the properties of the emitting outflow and its engine to an unprecedented detail. Since 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 sec, 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 per cent of the injected energy is transferred to the ejecta via mechanical work., Comment: Under review for publication in ApJL

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 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 not able to account for the current set of optical prompt detections that have been made and 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., Comment: 22 pages, with 5 pages for the appendix; 11 figures; accepted to ApJ; the MCRaT code is open source and available at: https://github.com/lazzati-astro/MCRaT and the ProcessMCRaT code is also open source and available at: https://github.com/parsotat/ProcessMCRaT

We present new radio observations of the binary neutron star merger GW170817 carried out with the Karl G. Jansky Very large Array (VLA) more than 3\,yrs after the merger. Our combined dataset is derived by co-adding more than $\approx32$\,hours of VLA time on-source, and as such provides the deepest combined observation (rms sensitivity $\approx 0.99\,\mu$Jy) of the GW170817 field obtained to date at 3\,GHz. We find no evidence for a late-time radio re-brightening at a mean epoch of $t\approx 1200$\,d since merger, in contrast to a $\approx 2.1\,\sigma$ excess observed at X-ray wavelengths at the same mean epoch. Our measurements agree with expectations from the post-peak decay of the radio afterglow of the GW170817 structured jet. Using these results, we constrain the parameter space of models that predict a late-time radio re-brightening possibly arising from the high-velocity tail of the GW170817 kilonova ejecta, which would dominate the radio and X-ray emission years after the merger (once the structured jet afterglow fades below detection level). Our results point to a steep energy-speed distribution of the kilonova ejecta (with energy-velocity power law index $\alpha \gtrsim 5$). We suggest possible implications of our radio analysis, when combined with the recent tentative evidence for a late-time re-brightening in the X-rays, and highlight the need for continued radio-to-X-ray monitoring to test different scenarios., Comment: 6 pages, 2 figures. Version published in ApJL (includes minor updates to figures, data and refs)

GRB 190114C was a bright burst that occurred in the local Universe (z=0.425). It was the first gamma-ray burst (GRB) ever detected at TeV energies, thanks to MAGIC. We characterize the ambient medium properties of the host galaxy through the study of the absorbing X-ray column density. Joining Swift, XMM-Newton, and NuSTAR observations, we find that the GRB X-ray spectrum is characterized by a high column density that is well in excess of the expected Milky Way value and decreases, by a factor of ~2, around ~$10^5$ s. Such a variability is not common in GRBs. The most straightforward interpretation of the variability in terms of photoionization of the ambient medium is not able to account for the decrease at such late times, when the source flux is less intense. Instead, we interpret the decrease as due to a clumped absorber, denser along the line of sight and surrounded by lower-density gas. After the detection at TeV energies of GRB 190114C, two other GRBs were promptly detected. They share a high value of the intrinsic column density and there are hints for a decrease of the column density, too. We speculate that a high local column density might be a common ingredient for TeV-detected GRBs., Comment: 6 pages, 4 figures. Accepted for publication in A&A

The binary neutron star (BNS) merger GW170817 was the first astrophysical source detected in gravitational waves and multiwavelength electromagnetic radiation. The almost simultaneous observation of a pulse of gamma rays proved that BNS mergers are associated with at least some short gamma-ray bursts (GRBs). However, the gamma-ray pulse was faint, casting doubt on the association of BNS mergers with the luminous, highly relativistic outflows of canonical short GRBs. Here we show that structured jets with a relativistic, energetic core surrounded by slower and less energetic wings produce afterglow emission that brightens characteristically with time, as recently seen in the afterglow of GW170817. Initially, we only see the relatively slow material moving towards us. As time passes, larger and larger sections of the outflow become visible, increasing the luminosity of the afterglow. The late appearance and increasing brightness of the multiwavelength afterglow of GW170817 allow us to constrain the geometry of its ejecta and thus reveal the presence of an off-axis jet pointing about 30° away from Earth. Our results confirm a single origin for BNS mergers and short GRBs: GW170817 produced a structured outflow with a highly relativistic core and a canonical short GRB. We did not see the bright burst because it was beamed away from Earth. However, approximately one in 20 mergers detected in gravitational waves will be accompanied by a bright, canonical short GRB.

The binary neutron star (BNS) merger GW170817 was the first astrophysical source detected in gravitational waves and multiwavelength electromagnetic radiation. The almost simultaneous observation of a pulse of gamma rays proved that BNS mergers are associated with at least some short gamma-ray bursts (GRBs). However, the gamma-ray pulse was faint, casting doubt on the association of BNS mergers with the luminous, highly relativistic outflows of canonical short GRBs. Here we show that structured jets with a relativistic, energetic core surrounded by slower and less energetic wings produce afterglow emission that brightens characteristically with time, as recently seen in the afterglow of GW170817. Initially, we only see the relatively slow material moving towards us. As time passes, larger and larger sections of the outflow become visible, increasing the luminosity of the afterglow. The late appearance and increasing brightness of the multiwavelength afterglow of GW170817 allow us to constrain the geometry of its ejecta and thus reveal the presence of an off-axis jet pointing about 30° away from Earth. Our results confirm a single origin for BNS mergers and short GRBs: GW170817 produced a structured outflow with a highly relativistic core and a canonical short GRB. We did not see the bright burst because it was beamed away from Earth. However, approximately one in 20 mergers detected in gravitational waves will be accompanied by a bright, canonical short GRB.

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have evaluated the impact of differences in data processing or flux determination methods used by different groups, and attempted to mitigate these differences to provide a more uniform dataset. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.92\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.584\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.168\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al. (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle., Comment: Accepted by ApJ. Revised X-ray and up-to-date afterglow measurements are available in ASCII format on GitHub (https://github.com/kmooley/GW170817/) and mirrored on tauceti.caltech.edu (http://www.tauceti.caltech.edu/kunal/gw170817/). The online dataset will be continuously updated as new measurements become available

GRB 170817A was a subluminous short gamma-ray burst detected about 1.74 s after the gravitational wave signal GW170817 from a binary neutron star (BNS) merger. It is now understood as an off-axis event powered by the cocoon of a relativistic jet pointing 15 to 30 degrees away from the direction of observation. The cocoon was energized by the interaction of the incipient jet with the non-relativistic baryon wind from the merger remnant, resulting in a structured outflow with a narrow core and broad wings. In this paper, we couple the observational constraints on the structured outflow with a model for the jet-wind interaction to constrain the intrinsic properties with which the jet was launched by the central engine, including its time delay from the merger event. Using wind prescriptions inspired by magnetized BNS merger simulations, we find that the jet was launched within about 0.4 s from the merger, implying that the 1.74 s observed delay was dominated by the fireball propagation up to the photospheric radius. We also constrain, for the first time for any gamma-ray burst, the jet opening angle at injection and set a lower limit to its asymptotic Lorentz factor. These findings suggest an initially Poynting-flux dominated jet, launched via electromagnetic processes. If the jet was powered by an accreting black hole, they also provide a significant constraint on the survival time of the metastable neutron star remnant., Comment: Accepted for publication in ApJ

A comprehensive understanding of Gamma Ray Bursts (GRBs) has been elusive due to the variety of questions surrounding the radiation mechanism at play in these events. Polarization measurements of GRBs can heavily constrain the relevant radiation mechanisms and the structure of the GRB jet; however, there is a limited number of theoretical predictions that observed GRB polarizations can be compared against. Here, we conduct radiative transfer calculations of a set of two dimensional relativistic hydrodynamic long GRB (LGRB) jet simulations, of a constant and a variable jet, using the Monte Carlo Radiation Transport (MCRaT) code. MCRaT has been enhanced by the inclusion of polarization; it has been first verified by reproducing a variety of results in the literature and then used to obtain the time integrated and resolved polarization degrees and angles of the synthetic LGRBs. While the obtained time-integrated polarization degrees {($\lesssim 1$\%)} are consistent with the constraints from the POLAR experiment, they are lower than other theoretical studies due to the lack of strong gradients in the model jet profiles that we use. The time resolved results suggests that GRBs with wide jets observed on axis will have small polarization degrees ($\lesssim 2\%$) and {constant polarization angles}, during the brightest portion of the light curve. GRBs observed off axis will have larger polarization degrees and polarization angles that change with the temporal structure of radiating shells in the outflow. We then place our results in the context of GRB prompt emission models and future LEAP and POLAR-2 GRB polarimetry detections., Comment: 14 pages, 7 figures, final version after it has been accepted to ApJ

We study the ensemble of linear polarization measurement in the optical afterglows of long-duration gamma-ray bursts. We assume a non sideways-expanding top-hat jet geometry and use the relatively large number of measurements under the assumption that they represent a statistically unbiased sample. This allows us to constrain the ratio between the maximum predicted polarization and the measured one, which is an indicator of the geometry of the magnetic field in the downstream region of the external shock. We find that the measured polarization is substantially suppressed with respect to the maximum possible for either a completely ordered magnetic field parallel to the shock normal or to a field that is entirely contained in the shock plane. The measured polarization is limited, on average, to between 25 and 30% of the maximum theoretically possible value. This reduction requires the perpendicular component of the magnetic field to be dominant in energy with respect to the component parallel to the shock front, as expected for a shock generated and/or shock compressed field. We find, however, that the data only marginally support the assumption of a simple top-hat jet, pointing towards a more complex geometry for the outflow., Comment: 6 pages, submitted to AAS journals

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a t^(0.86±0.04) rise, a t^(−1.90±0.12) decline, and a peak occurring at 155±4 days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index (−0.569±0.002) across all epochs. This gives a precise and updated estimate of the electron power-law index, p=2.138±0.004. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, <0.01 cm⁻³, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10^(10.4) G and 10¹⁶ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

The disks of Active Galactic Nuclei (AGNs), traditionally studied as the feeders of the supermassive black holes (SMBHs) at their centers, have recently triggered a lot of interest also as hosts to massive stars and hence their neutron star (NS) and black hole (BH) remnants. Migration traps and gas torques in these disks favor binary formation and enhance the rate of compact object mergers. In these environments both long gamma-ray bursts (GRBs) from the death of massive stars and short GRBs from NS-NS and NS-BH mergers are expected. However, their properties in the unique environments of AGN disks have never been studied. Here we show that GRBs in AGNs can display unique features, owing to the unusual relative position of the shocks that characterize the burst evolution and the Thomson photosphere of the AGN disk. In dense environments, for example, the external shock develops before the internal shocks, leading to prompt emission powered by a relativistic reverse shock. The transient's time evolution is also compressed, yielding afterglow emission that is much brighter and peaks much earlier than for GRBs in the interstellar medium. Additionally, in regions of the disk that are sufficiently dense and extended, the light curves are dominated by diffusion, since the fireball is trapped inside the disc photosphere. These diffusion-dominated transients emerge on timescales of days in disks around SMBHs of $\sim 10^6 M_\odot$ to years for SMBHs of $\sim 10^8 M_\odot$. Finally, a large fraction of events, especially in AGNs with SMBHs $\lesssim 10^7 M_\odot$, display time-variable absorption in the X-ray band., Comment: Submitted to ApJ Letters

The detection of GW170817, it's extensive multi-wavelength follow-up campaign, and the large amount of theoretical development and interpretation that followed, have resulted in a significant step forward in the understanding of the binary neutron star merger phenomenon as a whole. One of its aspects is seeing the merger as a progenitor of short gamma-ray bursts (SGRB), which will be the subject of this review. On the one hand, GW170817 observations have confirmed some theoretical expectations, exemplified by the confirmation that binary neutron star mergers are the progenitors of SGRBs. In addition, the multimessenger nature of GW170817 has allowed for gathering of unprecedented data, such as the trigger time of the merger, the delay with which the gamma-ray photons were detected, and the brightening afterglow of an off-axis event. All together, the incomparable richness of the data from GW170817 has allowed us to paint a fairly detailed picture of at least one SGRB. I will detail what we learned, what new questions have arisen, and the perspectives for answering them when a sample of GW170817-comparable events have been studied., Comment: Mini-review for Frontiers. Final published version

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have evaluated the impact of differences in data processing or flux determination methods used by different groups, and attempted to mitigate these differences to provide a more uniform dataset. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.92\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.584\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.168\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al. (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle., Comment: Accepted by ApJ. Revised X-ray and up-to-date afterglow measurements are available in ASCII format on GitHub (https://github.com/kmooley/GW170817/) and mirrored on tauceti.caltech.edu (http://www.tauceti.caltech.edu/kunal/gw170817/). The online dataset will be continuously updated as new measurements become available

GRB 170817A was a subluminous short gamma-ray burst detected about 1.74 s after the gravitational wave signal GW170817 from a binary neutron star (BNS) merger. It is now understood as an off-axis event powered by the cocoon of a relativistic jet pointing 15 to 30 degrees away from the direction of observation. The cocoon was energized by the interaction of the incipient jet with the non-relativistic baryon wind from the merger remnant, resulting in a structured outflow with a narrow core and broad wings. In this paper, we couple the observational constraints on the structured outflow with a model for the jet-wind interaction to constrain the intrinsic properties with which the jet was launched by the central engine, including its time delay from the merger event. Using wind prescriptions inspired by magnetized BNS merger simulations, we find that the jet was launched within about 0.4 s from the merger, implying that the 1.74 s observed delay was dominated by the fireball propagation up to the photospheric radius. We also constrain, for the first time for any gamma-ray burst, the jet opening angle at injection and set a lower limit to its asymptotic Lorentz factor. These findings suggest an initially Poynting-flux dominated jet, launched via electromagnetic processes. If the jet was powered by an accreting black hole, they also provide a significant constraint on the survival time of the metastable neutron star remnant., Comment: Accepted for publication in ApJ

A comprehensive understanding of Gamma Ray Bursts (GRBs) has been elusive due to the variety of questions surrounding the radiation mechanism at play in these events. Polarization measurements of GRBs can heavily constrain the relevant radiation mechanisms and the structure of the GRB jet; however, there is a limited number of theoretical predictions that observed GRB polarizations can be compared against. Here, we conduct radiative transfer calculations of a set of two dimensional relativistic hydrodynamic long GRB (LGRB) jet simulations, of a constant and a variable jet, using the Monte Carlo Radiation Transport (MCRaT) code. MCRaT has been enhanced by the inclusion of polarization; it has been first verified by reproducing a variety of results in the literature and then used to obtain the time integrated and resolved polarization degrees and angles of the synthetic LGRBs. While the obtained time-integrated polarization degrees {($\lesssim 1$\%)} are consistent with the constraints from the POLAR experiment, they are lower than other theoretical studies due to the lack of strong gradients in the model jet profiles that we use. The time resolved results suggests that GRBs with wide jets observed on axis will have small polarization degrees ($\lesssim 2\%$) and {constant polarization angles}, during the brightest portion of the light curve. GRBs observed off axis will have larger polarization degrees and polarization angles that change with the temporal structure of radiating shells in the outflow. We then place our results in the context of GRB prompt emission models and future LEAP and POLAR-2 GRB polarimetry detections., Comment: 14 pages, 7 figures, final version after it has been accepted to ApJ

We study the ensemble of linear polarization measurement in the optical afterglows of long-duration gamma-ray bursts. We assume a non sideways-expanding top-hat jet geometry and use the relatively large number of measurements under the assumption that they represent a statistically unbiased sample. This allows us to constrain the ratio between the maximum predicted polarization and the measured one, which is an indicator of the geometry of the magnetic field in the downstream region of the external shock. We find that the measured polarization is substantially suppressed with respect to the maximum possible for either a completely ordered magnetic field parallel to the shock normal or to a field that is entirely contained in the shock plane. The measured polarization is limited, on average, to between 25 and 30% of the maximum theoretically possible value. This reduction requires the perpendicular component of the magnetic field to be dominant in energy with respect to the component parallel to the shock front, as expected for a shock generated and/or shock compressed field. We find, however, that the data only marginally support the assumption of a simple top-hat jet, pointing towards a more complex geometry for the outflow., Comment: 6 pages, submitted to AAS journals

It has been proposed that blazar jets are structured, with a fast spine surrounded by a slower sheath or layer. This structured jet model explains some properties of their emission and morphology. Because of their relative motion, the radiation produced by one component is seen amplified by the other, thus enhancing the inverse Compton emission of both. Radiation is emitted anisotropically in the comoving frames and causes the emitting plasma to recoil. As seen in the observer frame, this corresponds to a deceleration of the fastest component (the spine) and an acceleration of the slower one (the layer). While the deceleration of the spine has already been investigated, here we study for the first time the acceleration of the sheath and find self-consistent velocity profile solutions for both the spine and the sheath while accounting for radiative cooling. We find that the sheath can be accelerated to the velocities required by the observations if its leptons remain energetic in the acceleration region, assumed to be of the order of ~100 Schwarzschild radii, demanding continuous injection of energetic particles in that region.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations., Comment: submitted to Astro2020

The discovery of GW170817, the merger of a binary neutron star (NS) triggered by a gravitational wave detection by LIGO and Virgo, has opened a new window of exploration in the physics of NSs and their cosmological role. Among the important quantities to measure are the mass and velocity of the ejecta produced by the tidally disrupted NSs and the delay - if any - between the merger and the launching of a relativistic jet. These encode information on the equation of state of the NS, the nature of the merger remnant, and the jet launching mechanism, as well as yielding an estimate of the mass available for r-process nucleosynthesis. Here we derive analytic estimates for the structure of jets expanding in environments with different density, velocity, and radial extent. We compute the jet-cocoon structure and the properties of the broadband afterglow emission as a function of the ejecta mass, velocity, and time delay between merger and launch of the jet. We show that modeling of the afterglow light curve can constrain the ejecta properties and, in turn, the physics of neutron density matter. Our results increase the interpretative power of electromagnetic observations by allowing for a direct connection with the merger physics., Comment: ApJ in press. Updated references and discussion of BNS ejecta

GHz radio astronomy has played a fundamental role in the recent dazzling discovery of GW170817, a neutron star (NS)-NS merger observed in both gravitational waves (GWs) and light at all wavelengths. Here we show how the expected progress in sensitivity of ground-based GW detectors over the next decade calls for U.S.-based GHz radio arrays to be improved beyond current levels. We discuss specifically how several new scientific opportunities would emerge in multi-messenger time-domain astrophysics if a next generation GHz radio facility with sensitivity and resolution $10\times$ better than the current Jansky Very Large Array (VLA) were to work in tandem with ground-based GW detectors. These opportunities include probing the properties, structure, and size of relativistic jets and wide-angle ejecta from NS-NS mergers, as well as unraveling the physics of their progenitors via host galaxy studies., Comment: 8 pages, 5 figures, Submitted to the Astro2020 Decadal Survey

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science., Comment: 50 pages, Probe class mission concept study report submitted to NASA for Astro2020 Decadal Survey