Showing 1–41 of 41 results for author: Bromberg, O

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

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

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

Submitted 28 May, 2024; originally announced May 2024.

Comments: 15 pages, 11 figures

Constraints on the accretion properties of quasi-periodic erupters from GRMHD simulations

Authors: Anna Chashkina, Omer Bromberg, Amir Levinson, Ehud Nakar

Abstract: Context: Some apparently quiescent supermassive black holes (BHs) at centers of galaxies show quasi-periodic eruptions (QPEs) in the X-ray band, the nature of which is still unknown. A possible origin for the eruptions is an accretion disk, however the properties of such disks are restricted by the timescales of reccurance and durations of the flares. Aims: In this work we test the possibility tha… ▽ More Context: Some apparently quiescent supermassive black holes (BHs) at centers of galaxies show quasi-periodic eruptions (QPEs) in the X-ray band, the nature of which is still unknown. A possible origin for the eruptions is an accretion disk, however the properties of such disks are restricted by the timescales of reccurance and durations of the flares. Aims: In this work we test the possibility that the known QPEs can be explained by accretion from a compact accretion disk with an outer radius $r_{\rm out}\sim 10-40 r_{\rm g}$, focusing on a particular object GSN 069. Methods: We run several 3D GRMHD simulations with the {\tt HARMPI} code of thin and thick disks and study how the initial disk parameters such as thickness, magnetic field configuration, magnetization and Kerr parameter affect the observational properties of QPEs. Results: We show that accretion onto a slowly rotating BH through a small, thick accretion disk with an initially low plasma $β$ can explain the observed flare duration, the time between outbursts and the lack of evidence for a variable jet emission. In order to form such a disk the accreting matter should have a low net angular momentum. A potential source for such low angular momentum matter with a quasi periodic feeding mechanism might be a tight binary of wind launching stars. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: Submitted to Astronomy and Astrophysics

Jets with a Twist: Emergence of FR0 Jets in 3D GRMHD Simulation of Zero Angular Momentum Black Hole Accretion

Authors: Aretaios Lalakos, Alexander Tchekhovskoy, Omer Bromberg, Ore Gottlieb, Jonatan Jacquemin-Ide, Matthew Liska, Haocheng Zhang

Abstract: Spinning supermassive black holes (BHs) in active galactic nuclei (AGN) magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within $3$ orders of magnitude in distance from the BH, well before reaching kpc-scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magneto… ▽ More Spinning supermassive black holes (BHs) in active galactic nuclei (AGN) magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within $3$ orders of magnitude in distance from the BH, well before reaching kpc-scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen much larger than the BH gravitational radius, $R_\text{B}=10^3R_\text{g}$. The BH develops dynamically-important large-scale magnetic fields, forms a magnetically-arrested disk (MAD), and launches relativistic jets that propagate well outside $R_\text{B}$ and suppress BH accretion to $1.5\%$ of the Bondi rate, $\dot{M}_\text{B}$. Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff-Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets $5-10\%$ of $\dot{M}_\text{B}$ reach the BH. Thereafter, the disk starts rocking back and forth by angles $90-180^\circ$: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to $\lesssim 2 \% ~ \dot{M}_\text{B}$. Because BAD and RAD states tangle up the jets and destroy them well inside $R_\text{B}$, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Relativistic Spherical Shocks in Expanding Media

Authors: Taya Govreen-Segal, Noam Youngerman, Ishika Palit, Ehud Nakar, Amir Levinson, Omer Bromberg

Abstract: We investigate the propagation of spherically symmetric shocks in relativistic homologously expanding media with density distributions following a power-law profile in their Lorentz factor. That is, $ρ_{ej} \propto t^{-3}γ_{e}(R,t)^{-α}$, where $ρ_{ej}$ is the medium proper density, $γ_{e}$ is its Lorentz factor, $α>0$ is constant and $t$, $R$ are the time and radius from the center. We find that… ▽ More We investigate the propagation of spherically symmetric shocks in relativistic homologously expanding media with density distributions following a power-law profile in their Lorentz factor. That is, $ρ_{ej} \propto t^{-3}γ_{e}(R,t)^{-α}$, where $ρ_{ej}$ is the medium proper density, $γ_{e}$ is its Lorentz factor, $α>0$ is constant and $t$, $R$ are the time and radius from the center. We find that the shocks behavior can be characterized by their proper velocity, $U'=Γ_s'β_s'$, where $Γ_s'$ is the shock Lorentz factor as measured in the immediate upstream frame and $β_s'$ is the corresponding 3-velocity. While generally, we do not expect the shock evolution to be self-similar, for every $α>0$ we find a critical value $U'_c$ for which a self-similar solution with constant $U'$ exists. We then use numerical simulations to investigate the behavior of general shocks. We find that shocks with $U'>U'_c$ have a monotonously growing $U'$, while those with $U'<U'_c$ have a decreasing $U'$ and will eventually die out. Finally, we present an analytic approximation, based on our numerical results, for the evolution of general shocks in the regime where $U'$ is ultra-relativistic. △ Less

Submitted 15 September, 2023; originally announced September 2023.

A kinetic study of black hole activation by local plasma injection into the inner magnetosphere

Authors: Idan Niv, Omer Bromberg, Amir Levinson, Benoit Cerutti, Benjamin Crinquand

Abstract: (Abridged) An issue of considerable interest in the theory of jet formation by the Blandford-Znajek mechanism, is how plasma is being continuously supplied to the magnetosphere to maintain it in a force-free state. Injection of electron-positron pairs via annihilation of MeV photons, emitted from a hot accretion flow, has been shown to be a viable possibility, but requires a high enough accretion… ▽ More (Abridged) An issue of considerable interest in the theory of jet formation by the Blandford-Znajek mechanism, is how plasma is being continuously supplied to the magnetosphere to maintain it in a force-free state. Injection of electron-positron pairs via annihilation of MeV photons, emitted from a hot accretion flow, has been shown to be a viable possibility, but requires a high enough accretion rate. At lower accretion rates, and in the absence of any other form of plasma supply, the magnetosphere becomes charge starved, forming intermittent spark gaps that can induce intense pair cascades via interactions with soft disk radiation, enabling outflow formation. It is often speculated that enough plasma can penetrate the inner magnetosphere from the accretion flow through some rearrangement of magnetic field lines (e.g., interchange instability). However, the question arises whether such episodes of plasma intrusion can prevent the formation of spark gaps. To address this question we conducted a suite of numerical experiments, by means of radiative, 2D axisymmetric general relativistic particle-in-cell simulations, in which plasma is injected into specified regions at a prescribed rate. We find that when pair production is switched off, nearly complete screening is achieved when the plasma is injected within the outer light cylinder at a high enough rate. Injection beyond the outer light cylinder results in either, the formation of large vacuum gaps, or coherent, large-amplitude oscillations of the magnetosphere, depending on the injection rate. Within the allowed dynamic range of our simulations, we see no evidence for the system to approach a steady state as the injection rate is increased. Switching on pair production results in nearly complete screening of the entire magnetosphere in all cases, with some fraction of the maximum Blandford-Znajek power emitted as TeV gamma-rays. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: 8 pages, 8 figures, for associated movie files, see https://www.youtube.com/playlist?list=PLsP9qR29TtwnAsBPx5xTPEUy6FZ_b-0PX

Interface instabilities in hydrodynamic relativistic jets

Authors: Pavel Abolmasov, Omer Bromberg

Abstract: Both the dynamics and the observational properties of relativistic jets are determined by their interaction with the ambient medium. A crucial role is played by the contact discontinuity at the jet boundary, which in the presence of jet collimation may become subject to Rayleigh-Taylor instability (RTI) and Richtmyer-Meshkov instability (RMI). Here, we study the evolution of these instabilities in… ▽ More Both the dynamics and the observational properties of relativistic jets are determined by their interaction with the ambient medium. A crucial role is played by the contact discontinuity at the jet boundary, which in the presence of jet collimation may become subject to Rayleigh-Taylor instability (RTI) and Richtmyer-Meshkov instability (RMI). Here, we study the evolution of these instabilities in non-magnetized relativistic jets using special relativistic three-dimensional hydrodynamic simulations. We show that the growth of initial perturbations is consistent with relativistic RTI operating in the jet collimation region. The contribution of RMI becomes important downstream from the collimation shock in agreement with theoretical expectations. Both instabilities reach non-linear scales above the shock convergence point and trigger strong turbulence, mixing jet with ambient matter. We devise an analytic solution for the mixing rate and show that it is sensitive to the external density gradients. Our results may be applied to different types of astrophysical objects. In particular, different contribution of interface instabilities is a natural explanation for the observed dichotomy between FR-I and FR-II radiogalaxies. The rapid slow-down in the jet of M87 is consistent with baryon entrainment from the circumnuclear matter with the observed density distribution. In microquasars, baryon loading triggered by interface instabilities is a probable reason for the low observed Lorentz factors. We show that the observed variability in gamma-ray bursts cannot come from mixing driven by interface instabilities and likely originates from the engine, suggesting the presence of magnetic fields in the jet. △ Less

Submitted 11 January, 2023; v1 submitted 20 October, 2022; originally announced October 2022.

Comments: accepted for publication to MNRAS; 18 pages, 10 figures

Black hole to photosphere: 3D GRMHD simulations of collapsars reveal wobbling and hybrid composition jets

Authors: Ore Gottlieb, Matthew Liska, Alexander Tchekhovskoy, Omer Bromberg, Aretaios Lalakos, Dimitrios Giannios, Philipp Mösta

Abstract: Long-duration $γ$-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth by a black-hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over 6 orders of magnitude in space and time. The collapsing stellar… ▽ More Long-duration $γ$-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth by a black-hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over 6 orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important and launches bipolar jets. The jets reach the photosphere at $\sim10^{12}$ cm with an opening angle $θ_j\sim6^\circ$ and a Lorentz factor $Γ_j\lesssim 30$, unbinding $\gtrsim90\%$ of the star. We find that (i) the disk-jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle $θ_t\sim12^\circ$, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection $θ_j+θ_t$ suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs may be rarer than currently thought and emerge in only $\sim 0.1\%$ of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where $\sim 10\%$ of the gas maintains magnetization $σ\gtrsim 0.1$. This indicates that both a photospheric component and reconnection may play a role in the prompt emission. △ Less

Submitted 20 June, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: For movies of the simulation, see https://oregottlieb.com/collapsar.html

Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients

Authors: Ore Gottlieb, Aretaios Lalakos, Omer Bromberg, Matthew Liska, Alexander Tchekhovskoy

Abstract: We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, $ l $, of the collapsing material and the magnetic field, $ B $, on the BH horizon: (i) Subrelativistic outflow (low $ l $ and high $ B $), (ii)… ▽ More We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, $ l $, of the collapsing material and the magnetic field, $ B $, on the BH horizon: (i) Subrelativistic outflow (low $ l $ and high $ B $), (ii) Stationary accretion shock instability (SASI; high $ l $ and low $ B $), (iii) Relativistic jets (high $ l $ and high $ B $). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices ($ \gtrsim 2 $), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the lightcurve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova+GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disk-jet flows, even in the absence of any tilt in the progenitor. △ Less

Submitted 26 December, 2021; v1 submitted 29 September, 2021; originally announced September 2021.

GRMHD simulations of BH activation by small scale magnetic loops: Formation of striped jets and active coronae

Authors: Anna Chashkina, Omer Bromberg, Amir Levinson

Abstract: We have performed a series of numerical experiments aimed at studying the activation of Kerr black holes (BHs) by advection of small scale magnetic fields. Such configurations may potentially give rise to the formation of quasi-striped Blandford-Znajek jets. It can also lead to enhanced dissipation and generation of plasmoids in current sheets formed in the vicinity of the BH horizon, which may co… ▽ More We have performed a series of numerical experiments aimed at studying the activation of Kerr black holes (BHs) by advection of small scale magnetic fields. Such configurations may potentially give rise to the formation of quasi-striped Blandford-Znajek jets. It can also lead to enhanced dissipation and generation of plasmoids in current sheets formed in the vicinity of the BH horizon, which may constitute a mechanism to power the hard X-ray emission seen in many accreting BH systems (a la lamppost models). Our analysis suggests that formation of quasi-striped jets with significant power may be possible provided loops with alternating polarity having sizes larger than $\sim 10 r_g$ or so can be maintained (either form sporadically or advected from outside) at a radius $\lesssim 10^2 r_g$. This conclusion is consistent with recent results of general relativistic force-free simulations. We also find that the accretion dynamics exhibits cyclic behaviour in MAD states, alternating between high accretion phases and quenched accretion phases during which the magnetosphere becomes force-free out to radii $\gtrsim 10r_g$. We suggest that such a behaviour should lead to notable variations of the observed luminosity and image of the inner disc (BH shadow image). Finally, we find that the transition between accreted loops on the BH gives rise to the formation of current sheets and energetic plasmoids on the jet boundary during intermittent periods when the jet becomes inactive, in addition to an equatorial current sheet that forms during peaks in the jet activity. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: 12 pages, 12 figures. Submitted to MNRAS

Modeling the Linear Polarization of GRB Afterglows Across the Electromagnetic Spectrum

Authors: Gal Birenbaum, Omer Bromberg

Abstract: Linear polarization measurements in the optical band show polarization degrees of a few percent at late times. Recently, polarization at sub-percent level was also detected in radio by ALMA, opening the window for multi-wavelength polarimetry and stressing the importance of properly modeling polarization in GRB afterglows across the EM spectrum. We introduce a numerical tool that can calculate the… ▽ More Linear polarization measurements in the optical band show polarization degrees of a few percent at late times. Recently, polarization at sub-percent level was also detected in radio by ALMA, opening the window for multi-wavelength polarimetry and stressing the importance of properly modeling polarization in GRB afterglows across the EM spectrum. We introduce a numerical tool that can calculate the polarization from relativistically moving surfaces by discretizing them to small patches of uniform magnetic field, calculating the polarized emission from each cell assuming synchrotron radiation and summing it to obtain the total degree of polarization. We apply this tool to afterglow shocks with random magnetic fields confined to the shock plane, considering electron radiative cooling. We analyze the observed polarization curves in several wavelengths above the cooling frequency and below the minimal synchrotron frequency and point to the characteristic differences between them. We present a method to constrain the jet opening angle and the viewing angle within the context of our model. Applying it to GRB 021004 we obtain angles of 10 and 8 degrees respectively and conclude that a non-negligible component of radial magnetic field is required to explain the 1% polarization level observed 3.5 days after the burst. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 14 pages, 10 figures

Intermittent mildly magnetized jets as the source of GRBs

Authors: Ore Gottlieb, Omer Bromberg, Amir Levinson, Ehud Nakar

Abstract: Gamma-ray bursts (GRBs) are powered by relativistic jets that exhibit intermittency over a broad range of timescales - from $ \sim $ ms to seconds. Previous numerical studies have shown that hydrodynamic (i.e., unmagnetized) jets that are expelled from a variable engine are subject to strong mixing of jet and cocoon material, which strongly inhibits the GRB emission. In this paper we conduct 3D RM… ▽ More Gamma-ray bursts (GRBs) are powered by relativistic jets that exhibit intermittency over a broad range of timescales - from $ \sim $ ms to seconds. Previous numerical studies have shown that hydrodynamic (i.e., unmagnetized) jets that are expelled from a variable engine are subject to strong mixing of jet and cocoon material, which strongly inhibits the GRB emission. In this paper we conduct 3D RMHD simulations of mildly magnetized jets with power modulation over durations of 0.1 s and 1 s, and a steady magnetic field at injection. We find that when the jet magnetization at the launching site is $σ\sim 0.1$, the initial magnetization is amplified by shocks formed in the flow to the point where it strongly suppresses baryon loading. We estimate that a significant contamination can be avoided if the magnetic energy at injection constitutes at least a few percent of the jet energy. The variability timescales of the jet after it breaks out of the star are then governed by the injection cycles rather than by the mixing process, suggesting that in practice jet injection should fluctuate on timescales as short as $ \sim 10 $ ms in order to account for the observed light curves. Better stability is found for jets with shorter modulations. We conclude that for sufficiently hot jets, the Lorentz factor near the photosphere can be high enough to allow efficient photospheric emission. Our results imply that jets with $ 10^{-2} < σ< 1 $ injected by a variable engine with $ \sim 10 $ ms duty cycle are plausible sources of long GRBs. △ Less

Submitted 14 April, 2021; v1 submitted 29 January, 2021; originally announced February 2021.

The structure of weakly-magnetized $ γ$-ray burst jets

Authors: Ore Gottlieb, Omer Bromberg, Chandra B. Singh, Ehud Nakar

Abstract: The interaction of gamma-ray burst (GRB) jets with the dense media into which they are launched promote the growth of local hydrodynamic instabilities along the jet boundary. In a companion paper we study the evolution of hydrodynamic (unmagnetized) jets, finding that mixing of jet-cocoon material gives rise to an interface layer, termed jet-cocoon interface (JCI), which contains a significant fra… ▽ More The interaction of gamma-ray burst (GRB) jets with the dense media into which they are launched promote the growth of local hydrodynamic instabilities along the jet boundary. In a companion paper we study the evolution of hydrodynamic (unmagnetized) jets, finding that mixing of jet-cocoon material gives rise to an interface layer, termed jet-cocoon interface (JCI), which contains a significant fraction of the system energy. We find that the angular structure of the jet + JCI, when they reach the homologous phase, can be approximated by a flat core (the jet) + a power-law function (the JCI) with indices that depend on the degree of mixing. In this paper we examine the effect of subdominant toroidal magnetic fields on the jet evolution and morphology. We find that weak fields can stabilize the jet against local instabilities. The suppression of the mixing diminishes the JCI and thus reshapes the jet's post-breakout structure. Nevertheless, the overall shape of the outflow can still be approximated by a flat core + a power-law function, although the JCI power-law decay is steeper. The effect of weak fields is more prominent in long GRB jets, where the mixing in hydrodynamic jets is stronger. In short GRB jets there is small mixing in both weakly magnetized and unmagnetized jets. This result influences the expected jet emission which is governed by the jet's morphology. Therefore, prompt and afterglow observations in long GRBs may be used as probes for the magnetic nature at the base of the jets. △ Less

Submitted 19 August, 2020; v1 submitted 22 July, 2020; originally announced July 2020.

The structure of hydrodynamic $ γ$-ray burst jets

Authors: Ore Gottlieb, Ehud Nakar, Omer Bromberg

Abstract: After being launched, GRB jets propagate through dense media prior to their breakout. The jet-medium interaction results in the formation of a complex structured outflow, often referred to as a "structured jet". The underlying physics of the jet-medium interaction that sets the post-breakout jet morphology has never been explored systematically. Here we use a suite of 3D simulations to follow the… ▽ More After being launched, GRB jets propagate through dense media prior to their breakout. The jet-medium interaction results in the formation of a complex structured outflow, often referred to as a "structured jet". The underlying physics of the jet-medium interaction that sets the post-breakout jet morphology has never been explored systematically. Here we use a suite of 3D simulations to follow the evolution of hydrodynamic long and short gamma-ray bursts (GRBs) jets after breakout to study the post-breakout structure induced by the interaction. Our simulations feature Rayleigh-Taylor fingers that grow from the cocoon into the jet, mix cocoon with jet material and destabilize the jet. The mixing gives rise to a previously unidentified region sheathing the jet from the cocoon, which we denote the jet-cocoon interface (JCI). long GRBs undergo strong mixing, resulting in most of the jet energy to drift into the JCI, while in short GRBs weaker mixing is possible, leading to a comparable amount of energy in the two components. Remarkably, the jet structure (jet-core plus JCI) can be characterized by simple universal angular power-law distributions, with power-law indices that depend solely on the mixing level. This result supports the commonly used power-law angular distribution, and disfavors Gaussian jets. At larger angles, where the cocoon dominates, the structure is more complex. The mixing shapes the prompt emission light curve and implies that typical long GRB afterglows are different from those of short GRBs. Our predictions can be used to infer jet characteristics from prompt and afterglow observations. △ Less

Submitted 1 November, 2020; v1 submitted 3 June, 2020; originally announced June 2020.

Particle acceleration in kink-unstable jets

Authors: Jordy Davelaar, Alexander A. Philippov, Omer Bromberg, Chandra B. Singh

Abstract: Magnetized jets in GRBs and AGNs are thought to be efficient accelerators of particles, however, the process responsible for the acceleration is still a matter of active debate. In this work, we study the kink-instability in non-rotating force-free jets using first-principle particle-in-cell simulations. We obtain similar overall evolution of the instability as found in MHD simulations. The instab… ▽ More Magnetized jets in GRBs and AGNs are thought to be efficient accelerators of particles, however, the process responsible for the acceleration is still a matter of active debate. In this work, we study the kink-instability in non-rotating force-free jets using first-principle particle-in-cell simulations. We obtain similar overall evolution of the instability as found in MHD simulations. The instability first generates large scale current sheets, which at later times break up into small-scale turbulence. Reconnection in these sheets proceeds in the strong guide field regime, which results in a formation of steep power laws in the particle spectra. Later evolution shows heating of the plasma, which is driven by small-amplitude turbulence induced by the kink instability. These two processes energize particles due to a combination of ideal and non-ideal electric fields. △ Less

Submitted 26 May, 2020; v1 submitted 29 October, 2019; originally announced October 2019.

Comments: 9 pages, 5 figures, accepted for publication in ApJL

Kink instability: evolution and energy dissipation in Relativistic Force-Free Non-Rotating Jets

Authors: Omer Bromberg, Chandra B. Singh, Jordy Davelaar, Alexander A. Philippov

Abstract: We study the evolution of kink instability in a force-free, non-rotating plasma column of high magnetization. The main dissipation mechanism is identified as reconnection of magnetic field-lines with various intersection angles, driven by the compression of the growing kink lobes. We measure dissipation rates ${\rm d} U_{Bφ}/{{\rm d}t} \approx -0.1 U_{Bφ}/τ$, where $τ$ is the linear growth time of… ▽ More We study the evolution of kink instability in a force-free, non-rotating plasma column of high magnetization. The main dissipation mechanism is identified as reconnection of magnetic field-lines with various intersection angles, driven by the compression of the growing kink lobes. We measure dissipation rates ${\rm d} U_{Bφ}/{{\rm d}t} \approx -0.1 U_{Bφ}/τ$, where $τ$ is the linear growth time of the kink instability. This value is consistent with the expansion velocity of the kink mode, which drives the reconnection. The relaxed state is close to a force-free Taylor state. We constraint the energy of that state using considerations from linear stability analysis. Our results are important for understanding magnetic field dissipation in various extreme astrophysical objects, most notably in relativistic jets. We outline the evolution of the kink instability in such jets and derive constrains on the conditions that allow for the kink instability to grow in these systems. △ Less

Submitted 27 August, 2019; v1 submitted 22 August, 2019; originally announced August 2019.

Comments: 18 pages, 13 figures

Eccentric binary black hole mergers in globular clusters hosting intermediate-mass black holes

Authors: Giacomo Fragione, Omer Bromberg

Abstract: Globular clusters (GCs) may harbour intermediate-mass black holes (IMBHs) at their centres. In these dynamically active environments stellar-mass black holes (SBHs) sink to the center soon after formation, due to dynamical friction and start interacting among themselves and with the central IMBH. Likely, some of the SBHs will form bound systems with the IMBH. A fraction of those will be triple sys… ▽ More Globular clusters (GCs) may harbour intermediate-mass black holes (IMBHs) at their centres. In these dynamically active environments stellar-mass black holes (SBHs) sink to the center soon after formation, due to dynamical friction and start interacting among themselves and with the central IMBH. Likely, some of the SBHs will form bound systems with the IMBH. A fraction of those will be triple systems composed of binary SBHs and the IMBH acting as a third distant perturber. If the SBH binary orbit is sufficiently inclined it can develop Lidov-Kozai (LK) oscillations, which can drive the system to high eccentricities and eventually to a merger due to gravitational wave (GW) emission on short timescales. In this work, we focus on the dynamics of the IMBH-SBH-SBH triples and illustrate that these systems can be possible sources of GWs. A distinctive signature of this scenario is that a considerable fraction of these mergers are highly eccentric when entering the LIGO band ($10$ Hz). Assuming that $\sim 20\%$ of GCs host IMBHs and a GC density in the range $n_{_{\rm GC}}=0.32$-$2.31\,\mathrm{Mpc}^{-3}$, we have estimated a rate $Γ=0.06$-$0.46\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ of these events. This suggests that dynamically-driven binary SBH mergers in this scenario could contribute to the merger events observed by LIGO/VIRGO. Full $N$-body simulations of GCs harbouring IMBHs are highly desirable to give a more precise constrain on this scenario. △ Less

Submitted 20 July, 2019; v1 submitted 22 March, 2019; originally announced March 2019.

Comments: 9 pages, 6 figures, 1 table, accepted by MNRAS

Energetic eruptions leading to a peculiar hydrogen-rich explosion of a massive star

Authors: Iair Arcavi, D. Andrew Howell, Daniel Kasen, Lars Bildsten, Griffin Hosseinzadeh, Curtis McCully, Zheng Chuen Wong, Sarah Rebekah Katz, Avishay Gal-Yam, Jesper Sollerman, Francesco Taddia, Giorgos Leloudas, Christoffer Fremling, Peter E. Nugent, Assaf Horesh, Kunal Mooley, Clare Rumsey, S. Bradley Cenko, Melissa L. Graham, Daniel A. Perley, Ehud Nakar, Nir J. Shaviv, Omer Bromberg, Ken J. Shen, Eran O. Ofek , et al. (28 additional authors not shown)

Abstract: Every supernova hitherto observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curv… ▽ More Every supernova hitherto observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curve peak, or a plateau in luminosity, lasting approximately 100 days before declining. Here we report observations of iPTF14hls, an event that has spectra identical to a hydrogen-rich core-collapse supernova, but characteristics that differ extensively from those of known supernovae. The light curve has at least five peaks and remains bright for more than 600 days; the absorption lines show little to no decrease in velocity; and the radius of the line-forming region is more than an order of magnitude bigger than the radius of the photosphere derived from the continuum emission. These characteristics are consistent with a shell of several tens of solar masses ejected by the star at supernova-level energies a few hundred days before a terminal explosion. Another possible eruption was recorded at the same position in 1954. Multiple energetic pre-supernova eruptions are expected to occur in stars of 95-130 solar masses, which experience the pulsational pair instability. That model, however, does not account for the continued presence of hydrogen, or the energetics observed here. Another mechanism for the violent ejection of mass in massive stars may be required. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: Published in Nature

The gamma-rays that accompanied GW170817 and the observational signature of a magnetic jet breaking out of NS merger ejecta

Authors: Omer Bromberg, Alexander Tchekhovskoy, Ore Gottlieb, Ehud Nakar, Tsvi Piran

Abstract: We present the first relativistic MHD numerical simulation of a magnetic jet that propagates through and emerges from the dynamical ejecta of a binary neutron star merger. Generated by the magnetized rotation of the merger remnant, the jet propagates through the ejecta and produces an energetic cocoon that expands at mildly relativistic velocities and breaks out of the ejecta. We show that if the… ▽ More We present the first relativistic MHD numerical simulation of a magnetic jet that propagates through and emerges from the dynamical ejecta of a binary neutron star merger. Generated by the magnetized rotation of the merger remnant, the jet propagates through the ejecta and produces an energetic cocoon that expands at mildly relativistic velocities and breaks out of the ejecta. We show that if the ejecta has a low-mass ($\sim10^{-7} M_\odot$) high-velocity ($v\sim0.85$ c) tail, the cocoon shock breakout will generate $γ$-ray emission that is comparable to the observed short GRB170817A that accompanied the recent gravitational wave event GW170817. Thus, we propose that this GRB, which is quite different from all other short GRBs observed before, was produced by a different mechanism. We expect, however, that such events are numerous and many will be detected in coming LIGO-Virgo runs. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: 7 pages 5 figures

Why Are We Not Detecting Any Type Ib Supernova Associated with Long Gamma Ray Bursts?

Authors: E. Sobacchi, J. Granot, O. Bromberg

Abstract: Core-collapse supernovae (SNe) of Types Ib and Ic arise from hydrogen-stripped stars, while the latter are also stripped of their helium. Both SN types have a similar temporal evolution, suggesting broadly similar progenitors. However, while some Ic SNe are associated with relativistic jets of gamma-ray bursts (GRBs), no GRB has yet been found in association with Type Ib SNe. Here we find that, ev… ▽ More Core-collapse supernovae (SNe) of Types Ib and Ic arise from hydrogen-stripped stars, while the latter are also stripped of their helium. Both SN types have a similar temporal evolution, suggesting broadly similar progenitors. However, while some Ic SNe are associated with relativistic jets of gamma-ray bursts (GRBs), no GRB has yet been found in association with Type Ib SNe. Here we find that, even if GRB-like central engines operate in both SNe Ib and Ic, different properties of their envelopes may accommodate this potential tension. In particular, we focus on the case of a low-mass, extended envelope surrounding the progenitor star (as produced, for example, by strong mass losses prior to explosion). If the envelopes of Type Ib SNe are sufficiently massive ($M_{\rm ext}\sim(0.3-1)M_\odot$) and extended ($R_{\rm ext}\sim 10^{13}\text{ cm}$), we show that (i) GRB-like jets cannot break out of the star; (ii) the SN light curve is compatible with current observations. Different envelope properties of Type Ib SNe with respect to (at least some) Type Ic SNe may be connected to the presence of a helium layer surrounding the star. △ Less

Submitted 5 May, 2017; v1 submitted 2 May, 2017; originally announced May 2017.

Comments: We realised that eq. 4 and 5 are not appropriate if a large amount of energy is deposited directly into the envelope. We are trying to fix this issue, but this is not immediate and can affect the conclusions. Hence, we prefer to withdraw the paper

A Common Central Engine for Long Gamma Ray Bursts and Type Ib/c Supernovae?

Authors: E. Sobacchi, J. Granot, O. Bromberg, M. C. Sormani

Abstract: Long-duration, spectrally-soft Gamma-Ray Bursts (GRBs) are associated with Type Ic Core Collapse (CC) Supernovae (SNe), and thus arise from the death of massive stars. In the collapsar model, the jet launched by the central engine must bore its way out of the progenitor star before it can produce a GRB. Most of these jets do not break out, and are instead "choked" inside the star, as the central-e… ▽ More Long-duration, spectrally-soft Gamma-Ray Bursts (GRBs) are associated with Type Ic Core Collapse (CC) Supernovae (SNe), and thus arise from the death of massive stars. In the collapsar model, the jet launched by the central engine must bore its way out of the progenitor star before it can produce a GRB. Most of these jets do not break out, and are instead "choked" inside the star, as the central-engine activity time, $t_{\rm e}$, is not long enough. Modelling the long-soft GRB duration distribution assuming a power-law distribution for their central-engine activity times, $\propto t_{\rm e}^{-α}$ for $t_{\rm e}>t_{\rm b}$, we find a steep distribution ($α\sim4$) and a typical GRB jet breakout time of $t_{\rm b}\sim 60\text{ s}$ in the star's frame. The latter suggests the presence of a low-density, extended envelope surrounding the progenitor star, similar to that previously inferred for low-luminosity GRBs. Extrapolating the range of validity of this power law below what is directly observable, to $t_{\rm e}<t_{\rm b}$, by only a factor of $\sim$4-5 produces enough events to account for all Type Ib/c SNe. Such extrapolation is necessary to avoid fine-tuning the distribution of central engine activity times with the breakout time, which are presumably unrelated. We speculate that central engines launching relativistic jets may operate in all Type Ib/c SNe. In this case, the existence of a common central engine would imply that (i) the jet may significantly contribute to the energy of the SN; (ii) various observational signatures, like the asphericity of the explosion, could be directly related to jet's interaction with the star. △ Less

Submitted 10 August, 2017; v1 submitted 30 April, 2017; originally announced May 2017.

Comments: 12 pages, 2 figures, accepted by MNRAS

Three-dimensional Relativistic MHD Simulations of Active Galactic Nuclei Jets: Magnetic Kink Instability and Fanaroff-Riley Dichotomy

Authors: Alexander Tchekhovskoy, Omer Bromberg

Abstract: Energy deposition by active galactic nuclei jets into the ambient medium can affect galaxy formation and evolution, the cooling of gas flows at the centres of galaxy clusters, and the growth of the supermassive black holes. However, the processes that couple jet power to the ambient medium and determine jet morphology are poorly understood. For instance, there is no agreement on the cause of the w… ▽ More Energy deposition by active galactic nuclei jets into the ambient medium can affect galaxy formation and evolution, the cooling of gas flows at the centres of galaxy clusters, and the growth of the supermassive black holes. However, the processes that couple jet power to the ambient medium and determine jet morphology are poorly understood. For instance, there is no agreement on the cause of the well-known Fanaroff-Riley (FR) morphological dichotomy of jets, with FRI jets being shorter and less stable than FRII jets. We carry out global 3D magnetohydrodynamic simulations of relativistic jets propagating through the ambient medium. We show that the flat density profiles of galactic cores slow down and collimate the jets, making them susceptible to the 3D magnetic kink instability. We obtain a critical power, which depends on the galaxy core mass and radius, below which jets become kink-unstable within the core, stall, and inflate cavities filled with relativistically-hot plasma. Jets above the critical power stably escape the core and form powerful backflows. Thus, the kink instability controls the jet morphology and can lead to the FR dichotomy. The model-predicted dependence of the critical power on the galaxy optical luminosity agrees well with observations. △ Less

Submitted 3 June, 2016; v1 submitted 14 December, 2015; originally announced December 2015.

Comments: MNRAS Letters, in press. 5 pages, 4 figures. This version matches the proofs. Movies available at https://goo.gl/o9FTwc

Relativistic MHD simulations of core-collapse GRB jets: 3D instabilities and magnetic dissipation

Authors: Omer Bromberg, Alexander Tchekhovskoy

Abstract: Relativistic jets naturally occur in astrophysical systems that involve accretion onto compact objects, such as core collapse of massive stars in gamma-ray bursts (GRBs) and accretion onto supermassive black holes in active galactic nuclei (AGN). It is generally accepted that these jets are powered electromagnetically, by the magnetised rotation of a central compact object. However, how they produ… ▽ More Relativistic jets naturally occur in astrophysical systems that involve accretion onto compact objects, such as core collapse of massive stars in gamma-ray bursts (GRBs) and accretion onto supermassive black holes in active galactic nuclei (AGN). It is generally accepted that these jets are powered electromagnetically, by the magnetised rotation of a central compact object. However, how they produce the observed emission and survive the propagation for many orders of magnitude in distance without being disrupted by current-driven non-axisymmetric instabilities is the subject of active debate. We carry out time-dependent 3D relativistic magnetohydrodynamic simulations of relativistic, Poynting flux dominated jets. The jets are launched self-consistently by the rotation of a strongly magnetised central compact object. This determines the natural degree of azimuthal magnetic field winding, a crucial factor that controls jet stability. We find that the jets are susceptible to two types of instability: (i) a global, external kink mode that grows on long time scales and causes the jets to bodily bend sideways. Whereas this mode does not cause jet disruption over the simulated distances, it substantially reduces jet propagation speed. We show, via an analytic model, that the growth of the external kink mode depends on the slope of the ambient medium density profile. In flat density distributions characteristic of galactic cores, an AGN jet may stall, whereas in stellar envelopes the external kink weakens as the jet propagates outward; (ii) a local, internal kink mode that grows over short time scales and causes small-angle magnetic reconnection and conversion of about half of jet electromagnetic energy flux into heat. Based on the robustness and energetics of the internal kink mode, we suggest that this instability is the main dissipation mechanism responsible for powering GRB prompt emission. △ Less

Submitted 11 August, 2015; originally announced August 2015.

Comments: 21 pages, 17 figures, submitted to MNRAS

Gamma-Ray Bursts as Sources of Strong Magnetic Fields

Authors: Jonathan Granot, Tsvi Piran, Omer Bromberg, Judith L. Racusin, Frédéric Daigne

Abstract: Gamma-Ray Bursts (GRBs) are the strongest explosions in the Universe, which due to their extreme character likely involve some of the strongest magnetic fields in nature. This review discusses the possible roles of magnetic fields in GRBs, from their central engines, through the launching, acceleration and collimation of their ultra-relativistic jets, to the dissipation and particle acceleration t… ▽ More Gamma-Ray Bursts (GRBs) are the strongest explosions in the Universe, which due to their extreme character likely involve some of the strongest magnetic fields in nature. This review discusses the possible roles of magnetic fields in GRBs, from their central engines, through the launching, acceleration and collimation of their ultra-relativistic jets, to the dissipation and particle acceleration that power their $γ$-ray emission, and the powerful blast wave they drive into the surrounding medium that generates their long-lived afterglow emission. An emphasis is put on particular areas in which there have been interesting developments in recent years. △ Less

Submitted 30 July, 2015; originally announced July 2015.

Comments: 51 pages; 12 figures; invited topical review, to be published in "The Strongest Magnetic Fields in the Universe" (Space Sciences Series of ISSI, Springer), Space Science Reviews

On the Composition of GRBs' Collapsar Jets

Authors: Omer Bromberg, Jonathan Granot, Tsvi Piran

Abstract: The duration distribution of long Gamma Ray Bursts reveals a plateau at durations shorter than ~20 s (in the observer frame) and a power-law decline at longer durations (Bromberg et al., 2012). Such a plateau arises naturally in the Collapsar model. In this model the engine has to operate long enough to push the jet out of the stellar envelope and the observed duration of the burst is the differen… ▽ More The duration distribution of long Gamma Ray Bursts reveals a plateau at durations shorter than ~20 s (in the observer frame) and a power-law decline at longer durations (Bromberg et al., 2012). Such a plateau arises naturally in the Collapsar model. In this model the engine has to operate long enough to push the jet out of the stellar envelope and the observed duration of the burst is the difference between the engine's operation time and the jet breakout time. We compare the jet breakout time inferred from the duration distribution (~10 s in the burst's frame) to the breakout time of a hydrodynamic jet (~10 s for typical parameters) and of a Poynting flux dominated jet with the same overall energy (<~1 s). As only the former is compatible with the duration of the plateau in the GRB duration distribution, we conclude that the jet is hydrodynamic during most of the time that its head is within the envelope of the progenitor star and around the time when it emerges from the star. This would naturally arise if the jet forms as a hydrodynamic jet in the first place or if it forms Poynting flux dominated but dissipates most of its magnetic energy early on within the progenitor star and emerges as a hydrodynamic jet. △ Less

Submitted 1 July, 2014; originally announced July 2014.

Comments: 17 pages, 3 figures

The dynamics of a highly magnetized jet propagating inside a star

Authors: Omer Bromberg, Jonathan Granot, Yuri Lyubarsky, Tsvi Piran

Abstract: The collapsar model explains the association of long duration gamma-Ray Bursts (GRBs) with stellar collapse. It involves a relativistic jet that forms at the core of a collapsing massive star. The jet penetrates the stellar envelope and the prompt GRB emission is produced once the jet is well outside the star. Most current models for generation of relativistic jets involve Poynting flux dominated… ▽ More The collapsar model explains the association of long duration gamma-Ray Bursts (GRBs) with stellar collapse. It involves a relativistic jet that forms at the core of a collapsing massive star. The jet penetrates the stellar envelope and the prompt GRB emission is produced once the jet is well outside the star. Most current models for generation of relativistic jets involve Poynting flux dominated outflows. We explore here the propagation of such a jet through a stellar envelope. The jet forms a bow shock around it. Energy dissipation at the head of this shock supplies energy to a cocoon that surrounds the jet. This cocoon exerts pressure on the jet and collimates it. While this description resembles the propagation of a hydrodynamic jets there are significant qualitative differences. Two Strong shocks, the reverse shock that slows down the hydrodynamic jet and the collimation shock that collimates it, cannot form within the Poynting flux dominated jet. As a result this jet moves much faster and dissipates much less energy while it crosses the stellar envelope. We construct here a simple analytic model that explores, self consistently, the jet-cocoon interaction and dynamics. Using this model we determine the properties of the jet, including its velocity, propagation time and shape. △ Less

Submitted 17 February, 2014; originally announced February 2014.

Comments: 31 pages, 6 figures

Mass entrainment and turbulence-driven acceleration of ultra-high energy cosmic rays in Centaurus A

Authors: Sarka Wykes, Judith H. Croston, Martin J. Hardcastle, Jean A. Eilek, Peter L. Biermann, Abraham Achterberg, Justin D. Bray, Alex Lazarian, Marijke Haverkorn, Ray J. Protheroe, Omer Bromberg

Abstract: Observations of the FR I radio galaxy Centaurus A in radio, X-ray and gamma-ray bands provide evidence for lepton acceleration up to several TeV and clues about hadron acceleration to tens of EeV. Synthesising the available observational constraints on the physical conditions and particle content in the jets, inner lobes and giant lobes of Centaurus A, we aim to evaluate its feasibility as an ultr… ▽ More Observations of the FR I radio galaxy Centaurus A in radio, X-ray and gamma-ray bands provide evidence for lepton acceleration up to several TeV and clues about hadron acceleration to tens of EeV. Synthesising the available observational constraints on the physical conditions and particle content in the jets, inner lobes and giant lobes of Centaurus A, we aim to evaluate its feasibility as an ultra-high-energy cosmic-ray source. We apply several methods of determining jet power and affirm the consistency of various power estimates of ~ 1 x 10^43 erg s^-1. Employing scaling relations based on previous results for 3C 31, we estimate particle number densities in the jets, encompassing available radio through X-ray observations. Our model is compatible with the jets ingesting ~ 3 x 10^21 g s^-1 of matter via external entrainment from hot gas and ~ 7 x 10^22 g s^-1 via internal entrainment from jet-contained stars. This leads to an imbalance between the internal lobe pressure available from radiating particles and magnetic field, and our derived external pressure. Based on knowledge of the external environments of other FR I sources, we estimate the thermal pressure in the giant lobes as 1.5 x 10^-12 dyn cm^-2, from which we deduce a lower limit to the temperature of ~ 1.6 x 10^8 K. Using dynamical and buoyancy arguments, we infer ~ 440-645 Myr and ~ 560 Myr as the sound-crossing and buoyancy ages of the giant lobes respectively, inconsistent with their spectral ages. We re-investigate the feasibility of particle acceleration via stochastic processes in the lobes, placing new constraints on the energetics and on turbulent input to the lobes. The same 'very hot' temperatures that allow self-consistency between the entrainment calculations and the missing pressure also allow stochastic UHECR acceleration models to work. △ Less

Submitted 5 August, 2013; v1 submitted 13 May, 2013; originally announced May 2013.

Comments: 13 pages, 2 figures. Minor updates to match accepted version. Accepted by A&A

Journal ref: Astronomy & Astrophysics, Volume 558, id. A19, 13 pp. (2013)

Short vs Long and Collapsars vs. non-Collapsar: a quantitative classification of GRBs

Authors: Omer Bromberg, Ehud Nakar, Tsvi Piran, Re'em Sari

Abstract: Gamma-Ray Bursts (GRBs) are traditionally divided to long and short according to their durations (>/< 2 sec). It was generally believed that this reflects a different physical origin: Collapsars (long) and non-Collapsars (short). We have recently shown that the duration distribution of Collapsars is flat, namely independent of the duration, at short durations. Using this model for the distribution… ▽ More Gamma-Ray Bursts (GRBs) are traditionally divided to long and short according to their durations (>/< 2 sec). It was generally believed that this reflects a different physical origin: Collapsars (long) and non-Collapsars (short). We have recently shown that the duration distribution of Collapsars is flat, namely independent of the duration, at short durations. Using this model for the distribution of Collapsars we determine the duration distribution of non-Collapsars and estimate the probability that a burst with a given duration (and hardness) is a Collapsar or not. We find that this probability depends strongly on the spectral window of the observing detector. While the commonly used limit of 2 sec is conservative and suitable for BATSE bursts, 40% of swift's bursts shorter than 2 sec are Collapsars and division >/<0.8 sec is more suitable for swift. We find that the duration overlap of the two populations is very large. On the one hand there is a non-negligible fraction of non-Collapsars longer than 10 sec, while on the other hand even bursts shorter than 0.5 sec in the swift sample have a non-negligible probability to be Collapsars. Our results enable the construction of non-Collapsar samples while controlling the Collapsar contamination. They also highlight that no firm conclusions can be drawn based on a single burst and they have numerous implications concerning previous studies of non-Collapsar properties that were based on the current significantly contaminated swift samples of localized short GRBs. Specifically: (i) all known short bursts with z>1 are most likely Collapsars, (ii) the only short burst with a clear jet break is most likely a Collapsar, indicating our lack of knowledge concerning non-Collapsar beaming (iii) the existence of non-Collapsars with durations up to 10 sec impose new challenges to non-Collapsar models. △ Less

Submitted 28 September, 2012; originally announced October 2012.

Comments: 39 pages, 9 figures

The Swift short gamma-ray burst rate density: prospects for detecting binary neutron star mergers by aLIGO

Authors: David Coward, Eric Howell, Tsvi Piran, Giulia Stratta, Marica Branchesi, Omer Bromberg, Bruce Gendre, Ronald Burman, Dafne Guetta

Abstract: Presently only 30% of short gamma ray bursts (SGRBs) have accurate redshifts, and this sample is highly biased by the limited sensitivity of {\it Swift} to detect SGRBs. We account for the dominant biases to calculate a realistic SGRB rate density out to $z = 0.5$ using the {\it Swift} sample of peak fluxes, redshifts, and those SGRBs with a beaming angle constraint from X-ray/optical observations… ▽ More Presently only 30% of short gamma ray bursts (SGRBs) have accurate redshifts, and this sample is highly biased by the limited sensitivity of {\it Swift} to detect SGRBs. We account for the dominant biases to calculate a realistic SGRB rate density out to $z = 0.5$ using the {\it Swift} sample of peak fluxes, redshifts, and those SGRBs with a beaming angle constraint from X-ray/optical observations. Assuming a significant fraction of binary neutron star mergers produce SGRBs, we calculate lower and upper detection rate limits of (1-180) per Yr by an advanced LIGO and Virgo coincidence search. Our detection rate is compatible with extrapolations using Galactic pulsar observations and population synthesis. △ Less

Submitted 22 June, 2012; originally announced June 2012.

Comments: 6 pages, 1 figure, to appear in: Proceedings of Science, Gamma-Ray Bursts 2012 Conference - GRB2012, May 07-11, 2012, Munich Germany

The Long, the Short and the Weak - the origin of GRBs

Authors: Tsvi Piran, Omer Bromberg, Ehud Nakar, Re'em Sari

Abstract: The origin of Gamma-Ray Bursts is one of the most interesting puzzles in recent astronomy. During the last decade a consensus formed that long GRBs (LGRBs) arise from the collapse of massive stars and that short GRBs (SGRBs) have a different origin, most likely neutron star mergers. A key ingredient of the Collapsar model that explains how the collapse of massive stars produces a GRB is the emerge… ▽ More The origin of Gamma-Ray Bursts is one of the most interesting puzzles in recent astronomy. During the last decade a consensus formed that long GRBs (LGRBs) arise from the collapse of massive stars and that short GRBs (SGRBs) have a different origin, most likely neutron star mergers. A key ingredient of the Collapsar model that explains how the collapse of massive stars produces a GRB is the emergence of a relativistic jet that penetrates the stellar envelope. The condition that the emerging jet penetrates the envelope poses strong constraints on the system. Using these constraints we show that: (i) Low luminosity GRBs (llGRBs), a sub population of GRBs with a very low luminosities (and other peculiar properties: single peaked, smooth and soft) cannot be formed by Collapsars. llGRBs must have a different origin (most likely a shock breakout). (ii) On the other hand regular LGRBs must be formed by Collapsars. (iii) While for BATSE the dividing duration between Collapsars and non-Collapsar is indeed at $\sim 2$ sec, the dividing duration is different for other GRBs detectors. In particular most Swift bursts longer than 0.8 sec are of a Collapsar origin. This last results requires a revision of many conclusions concerning the origin of Swift SGRBs which were based on the commonly used 2 sec limit. △ Less

Submitted 27 August, 2012; v1 submitted 4 June, 2012; originally announced June 2012.

Comments: Revised version accepted for publication in the Philosophical transactions - the proceedings of "New windows on transients across the Universe" a discussion meeting of the Royal Society, London, England and in the proceedings of "Xth International Conference on Gravitation, Astrophysics and Cosmology", Quy-Nhon, Vietnam. (Talks given by Tsvi Piran)

The Swift short gamma-ray burst rate density: implications for binary neutron star merger rates

Authors: David Coward, Eric Howell, Tsvi Piran, Giulia Stratta, Marica Branchesi, Omer Bromberg, Bruce Gendre, Ronald Burman, Dafne Guetta

Abstract: Short gamma-ray bursts (SGRBs) observed by {\it Swift} are potentially revealing the first insight into cataclysmic compact object mergers. To ultimately acquire a fundamental understanding of these events requires pan-spectral observations and knowledge of their spatial distribution to differentiate between proposed progenitor populations. Up to April 2012 there are only some 30% of SGRBs with re… ▽ More Short gamma-ray bursts (SGRBs) observed by {\it Swift} are potentially revealing the first insight into cataclysmic compact object mergers. To ultimately acquire a fundamental understanding of these events requires pan-spectral observations and knowledge of their spatial distribution to differentiate between proposed progenitor populations. Up to April 2012 there are only some 30% of SGRBs with reasonably firm redshifts, and this sample is highly biased by the limited sensitivity of {\it Swift} to detect SGRBs. We account for the dominant biases to calculate a realistic SGRB rate density out to $z\approx0.5$ using the {\it Swift} sample of peak fluxes, redshifts, and those SGRBs with a beaming angle constraint from X-ray/optical observations. We find an SGRB lower rate density of $8^{+5}_{-3}$ $\mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$ (assuming isotropic emission), and a beaming corrected upper limit of $1100^{+700}_{-470} $ $\mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$. Assuming a significant fraction of binary neutron star mergers produce SGRBs, we calculate lower and upper detection rate limits of $(1-180)$ yr$^{-1}$ by an aLIGO and Virgo coincidence search. Our detection rate is similar to the lower and realistic rates inferred from extrapolations using Galactic pulsar observations and population synthesis. △ Less

Submitted 22 June, 2012; v1 submitted 9 February, 2012; originally announced February 2012.

Comments: version 2: 7 pages, 1 figure, updated table & references, submitted to MNRAS

Journal ref: 2012, MNRAS, 425, 1365

New insights on the duration distribution of long GRBs from Collapsars

Authors: Omer Bromberg, Ehud Nakar, Tsvi Piran, Re'em Sari

Abstract: According the Collapsar model long gamma-ray bursts (LGRBs) involve relativistic jets that puncture the envelope of a collapsing star, and produced the γ-rays after they break out. This model provides a theoretical framework for the well known association between LGRBs and massive stars. However although this association is supported by a wealth of observations, to this date there is no direct obs… ▽ More According the Collapsar model long gamma-ray bursts (LGRBs) involve relativistic jets that puncture the envelope of a collapsing star, and produced the γ-rays after they break out. This model provides a theoretical framework for the well known association between LGRBs and massive stars. However although this association is supported by a wealth of observations, to this date there is no direct observational evidence for the emergence of the jet from the star. In other words there is no direct evidence for the Collapsar model. Here we show that a distinct signature of the Collapsar model is the appearance of a plateau in the duration distribution of the prompt GRB emission at times much shorter than the typical breakout time of the jet. This plateau is evident in the data of all major GRB satellites, and provides a direct evidence supporting the Collapsar model. It also enables us to place limits on the sizes and masses of LGRB progenitors; suggests the existence of a large population of choked (failed) GRBs; and indicates that the 2 s duration commonly used to separate Collapsars and non-Collapsars holds for BATSE and possibly Fermi GBM GRBs, but it is inconsistent with the duration distributions of Swift GRBs. △ Less

Submitted 27 December, 2011; originally announced December 2011.

Comments: 7 pages, 2 figures. To appear in: proceedings of the INTEGRAL conference "The Extreme and Variable High Energy Sky", held on September 19-23, 2011 in Chia Laguna, Cagliari, Italy; Proceedings of Science, SISSA, in press

An observational imprint of the Collapsar model of long Gamma Ray Bursts

Authors: Omer Bromberg, Ehud Nakar, Tsvi Piran, Re'em Sari

Abstract: The Collapsar model provides a theoretical framework for the well known association between long gamma-ray bursts (GRBs) and collapsing massive stars. A bipolar relativistic jet, launched at the core of a collapsing star, drills its way through the stellar envelope and breaks out of the surface before producing the observed gamma-rays. While a wealth of observations associate GRBs with the death o… ▽ More The Collapsar model provides a theoretical framework for the well known association between long gamma-ray bursts (GRBs) and collapsing massive stars. A bipolar relativistic jet, launched at the core of a collapsing star, drills its way through the stellar envelope and breaks out of the surface before producing the observed gamma-rays. While a wealth of observations associate GRBs with the death of massive stars, as yet there is no direct evidence for the Collapsar model itself. Here we show that a distinct signature of the Collapsar model is the appearance of a plateau in the duration distribution of the prompt GRB emission at times much shorter than the typical breakout time of the jet. This plateau is evident in the data of all three major satellites. These findings provide an evidence that directly supports the Collapsar model. Additionally, it suggests the existence of a large population of choked (failed) GRBs and that the 2 s duration commonly used to separate Collapsars and non-Collapasars is inconsistent with the duration distributions of {\it Swift} and Fermi GRBs and only holds for BATSE GRBs. △ Less

Submitted 13 November, 2011; originally announced November 2011.

Comments: 12 pages, 1 figure

Are low luminosity GRBs generated by relativistic jets?

Authors: Omer Bromberg, Ehud Nakar, Tsvi Piran

Abstract: Low luminosity gamma-ray bursts (ll-GRBs) constitute a sub-class of gamma-ray bursts (GRBs) that plays a central role in the GRB-supernova connection. While ll-GRBs differ from typical long GRBs (LGRBs) in many aspects, they also share some common features. Therefore, the question whether the gamma-ray emission of ll-GRBs and LGRBs has a common origin is of great interest. Here we address this que… ▽ More Low luminosity gamma-ray bursts (ll-GRBs) constitute a sub-class of gamma-ray bursts (GRBs) that plays a central role in the GRB-supernova connection. While ll-GRBs differ from typical long GRBs (LGRBs) in many aspects, they also share some common features. Therefore, the question whether the gamma-ray emission of ll-GRBs and LGRBs has a common origin is of great interest. Here we address this question by testing whether ll-GRBs, like LGRBs according to the Collapsar model, can be generated by relativistic jets that punch holes in the envelopes of their progenitor stars. The collapsar model predicts that the durations of most observed bursts will be comparable to, or longer than, the time it takes the jets to breakout of the star. We calculate the jet breakout times of ll-GRBs and compare them to the observed durations. We find that there is a significant access of ll-GRBs with durations that are much shorter than the jet breakout time and that these are inconsistent with the Collapsar model. We conclude that the processes that dominate the gamma-ray emission of ll-GRBs and of LGRBs are most likely fundamentally different. △ Less

Submitted 7 July, 2011; originally announced July 2011.

Comments: Submitted to ApJ letters. 12 pages, 1 figure

The propagation of relativistic jets in external media

Authors: Omer Bromberg, Ehud Nakar, Tsvi Piran, Re'em Sari

Abstract: Relativistic jets are ubiquitous in astrophysical systems that contain compact objects. They transport large amounts of energy to large distances from the source, and their interaction with the ambient medium has a crucial effect on the evolution of the system. The propagation of the jet is characterized by the formation of a shocked "head" at the front of the jet which dissipates the jet's energy… ▽ More Relativistic jets are ubiquitous in astrophysical systems that contain compact objects. They transport large amounts of energy to large distances from the source, and their interaction with the ambient medium has a crucial effect on the evolution of the system. The propagation of the jet is characterized by the formation of a shocked "head" at the front of the jet which dissipates the jet's energy and a cocoon that surrounds the jet and potentially collimates it. We present here a self consistent, analytic model that follows the evolution of the jet and its cocoon, and describes their interaction. We show that the critical parameter that determines the properties of the jet-cocoon system is the dimensionless ratio between the jet's energy density and the rest-mass energy density of the ambient medium. This parameter, together with the jet's injection angle, also determines whether the jet is collimated by the cocoon or not. The model is applicable to relativistic, unmagnetized, jets on all scales and may be used to determine the conditions in AGNs jets as well as in GRBs or microquasars. It shows that AGN and microquasar jets are hydrodynamically collimated due to the interaction with the ambient medium, while GRB jets can be collimated only inside a star and become uncollimated once they breakout. △ Less

Submitted 7 July, 2011; originally announced July 2011.

Comments: 35 pages, 4 figures

Sub-Photospheric Emission from Relativistic Radiation Mediated Shocks in GRBs

Authors: Omer Bromberg, Ziv Mikolitzky, Amir Levinson

Abstract: It is proposed that the prompt emission observed in bursts that exhibit a thermal component originates from relativistic radiation mediated shocks that form below the photosphere of the GRB outflow. It is argue that such shocks are expected to form in luminous bursts via collisions of shells that propagate with moderate Lorentz factors $Γ\lesssim 500$. Faster shells will collide above the photosph… ▽ More It is proposed that the prompt emission observed in bursts that exhibit a thermal component originates from relativistic radiation mediated shocks that form below the photosphere of the GRB outflow. It is argue that such shocks are expected to form in luminous bursts via collisions of shells that propagate with moderate Lorentz factors $Γ\lesssim 500$. Faster shells will collide above the photosphere to form collisionless shocks. We demonstrate that in events like GRB 090902B a substantial fraction of the explosion energy is dissipated below the photosphere, in a region of moderate optical depth $τ\lesssim300$, whereas in GRB 080916C the major fraction of the energy dissipates above the photosphere. We show that under conditions anticipated in many GRBs, such relativistic radiation mediated shocks convect enough radiation upstream to render photon production in the shock transition negligible, unlike the case of shock breakout in supernovae. The resulting spectrum, as measured in the shock frame, has a relatively low thermal peak, followed by a broad, nonthermal component extending up to the KN limit. △ Less

Submitted 21 January, 2011; originally announced January 2011.

Comments: 27 pages, 5 figures

Recollimation and Radiative Focusing of Relativistic Jets: Applications to Blazars and M87

Authors: Omer Bromberg, Amir Levinson

Abstract: Recent observations of M87 and some blazars reveal violent activity in small regions located at relatively large distances from the central engine. Motivated by these considerations, we study the hydrodynamic collimation of a relativistic cooling outflow using a semi-analytical model developed earlier. We first demonstrate that radiative cooling of the shocked outflow layer can lead to a focusin… ▽ More Recent observations of M87 and some blazars reveal violent activity in small regions located at relatively large distances from the central engine. Motivated by these considerations, we study the hydrodynamic collimation of a relativistic cooling outflow using a semi-analytical model developed earlier. We first demonstrate that radiative cooling of the shocked outflow layer can lead to a focusing of the outflow and its reconfinement in a region having a very small cross-sectional radius. Such a configuration can produce rapid variability at large distances from the central engine via reflections of the converging recollimation shock. Possible applications of this model to TeV blazars are discussed. We then apply our model to M87. The low radiative efficiency of the M87 jet renders focusing unlikely. However, the shallow profile of the ambient medium pressure inferred from observations, results in extremely good collimation that can explain the reported variability of the X-ray flux emitted from the HST-1 knot. △ Less

Submitted 3 May, 2009; v1 submitted 3 October, 2008; originally announced October 2008.

Comments: 17 pages, 4 figures. Accepted by ApJ

Journal ref: Astrophys.J.699:1274-1280,2009

Collimation and Radiative Deceleration of Jets in TeV AGNs

Authors: Amir Levinson, Omer Bromberg

Abstract: We consider some implications of the rapid X-ray and TeV variability observed in M87 and the TeV blazars. We outline a model for jet focusing and demonstrate that modest radiative cooling can lead to recollimation of a relativistic jet in a nozzle having a very small cross-sectional radius. Such a configuration can produce rapid variability at large distances from the central engine and may expl… ▽ More We consider some implications of the rapid X-ray and TeV variability observed in M87 and the TeV blazars. We outline a model for jet focusing and demonstrate that modest radiative cooling can lead to recollimation of a relativistic jet in a nozzle having a very small cross-sectional radius. Such a configuration can produce rapid variability at large distances from the central engine and may explain recent observations of the HST-1 knot in M87. Possible applications of this model to TeV blazars are discussed. We also discuss a scenario for the very rapid TeV flares observed with HESS and MAGIC in some blazars, that accommodates the relatively small Doppler factors inferred from radio observations. △ Less

Submitted 17 December, 2007; originally announced December 2007.

Comments: 5 pages. Proceedings of "High Energy Phenomena in Relativistic Outflows", held in Dublin, Ireland, September 24-28, 2007

Report number: Int.J.Mod.Phys.D17,2008,1603-1610

Journal ref: Int.J.Mod.Phys.D17:1603-1610,2008

Relativistic Photon Mediated Shocks

Authors: Amir Levinson, Omer Bromberg

Abstract: A system of equations governing the structure of a steady, relativistic radiation dominated shock is derived, starting from the general form of the transfer equation obeyed by the photon distribution function. Closure is obtained by truncating the system of moment equations at some order. The anisotropy of the photon distribution function inside the shock is shown to increase with increasing sho… ▽ More A system of equations governing the structure of a steady, relativistic radiation dominated shock is derived, starting from the general form of the transfer equation obeyed by the photon distribution function. Closure is obtained by truncating the system of moment equations at some order. The anisotropy of the photon distribution function inside the shock is shown to increase with increasing shock velocity, approaching nearly perfect beaming at upstream Lorentz factors $Γ_{-}>>1$. Solutions of the shock equations are presented for some range of upstream conditions. These solutions are shown to converge as the truncation order is increased. △ Less

Submitted 24 March, 2008; v1 submitted 21 November, 2007; originally announced November 2007.

Comments: 5 pages, a shorter version will appear in PRL

Journal ref: Phys.Rev.Lett.100:131101,2008

Hydrodynamic Collimation of Relativistic Outflows: Semianalytic Solutions and Application to Gamma-Ray Bursts

Authors: Omer Bromberg, Amir Levinson

Abstract: A model is developed for the confinement and collimation of a baryon poor outflow by its surrounding medium. Both, confinement by kinetic pressure of a static corona, and confinement by the ram pressure of a supersonic wind emanating from a disk surrounding the inner source are considered. Solutions are presented for the structure of the shocked layers of a deflected baryon poor jet (BPJ) and ex… ▽ More A model is developed for the confinement and collimation of a baryon poor outflow by its surrounding medium. Both, confinement by kinetic pressure of a static corona, and confinement by the ram pressure of a supersonic wind emanating from a disk surrounding the inner source are considered. Solutions are presented for the structure of the shocked layers of a deflected baryon poor jet (BPJ) and exterior wind. The dependence of the opening angle of the BPJ on the parameters of the confining medium are carefully examined. It is found that the BPJ shock may either converge to the symmetry axis or diverge away from it, depending on the opening angle of the BPJ injection cone. In the latter case the inner flow exhibits a non-uniform structure, consisting of an ultra-relativistic core containing the unshocked BPJ enveloped by the slower, shocked BPJ layer. The implications of our results to the prompt GRB emission are briefly discussed. △ Less

Submitted 11 November, 2007; v1 submitted 15 May, 2007; originally announced May 2007.

Comments: 22 pages, 7 figures. Added a section on the model limitation, and possible applications of this model in systems other then GRBs. Accepted by Apj

Journal ref: Astrophys.J.671:678-688,2007

An Observational Limit on the Earliest GRBs

Authors: Smadar Naoz, Omer Bromberg

Abstract: We predict the redshift of the first observable (i.e., in our past light cone) Gamma Ray Burst (GRB) and calculate the GRB-rate redshift distribution of the Population III stars at very early times (z=20-60). Using the last 2 years of data from Swift we place an upper limit on the efficiency (η_{GRB}) of GRB production per solar mass from the first generation of stars. We find that the first obs… ▽ More We predict the redshift of the first observable (i.e., in our past light cone) Gamma Ray Burst (GRB) and calculate the GRB-rate redshift distribution of the Population III stars at very early times (z=20-60). Using the last 2 years of data from Swift we place an upper limit on the efficiency (η_{GRB}) of GRB production per solar mass from the first generation of stars. We find that the first observable GRB is most likely to have formed at redshift 60. The observed rate of extremely high redshift GRBs (XRGs) is a subset of a group of 15 long GRBs per year, with no associated redshift and no optical afterglow counterparts, detected by Swift. Taking this maximal rate we get that η_{GRB}<1.1~10^{-4} GRBs per solar mass in stars. A more realistic evaluation, e.g., taking a subgroup of 5% of the total sample of Swift gives an upper limit of η_{GRB}<3.2~10^{-5} GRBs per solar mass. △ Less

Submitted 13 February, 2007; originally announced February 2007.

Comments: 6 Pages, 3 figures, submitted to MNRAS

The gravitational-wave spectrum of a non-axisymmetric torus around a rapidly spinning black hole

Authors: Omer Bromberg, Amir Levinson, Maurice van Putten

Abstract: The gravitational-wave spectrum emitted by a non-axisymmetric torus rotating at an angular velocity $Ω_T$, is derived in terms of a structure function representing a combination of sausage-tilt modes in the torus in the limit of an incompressible fluid. The analysis of the gravitational-wave spectrum is then applied to a model proposed recently, in which a highly magnetized torus interacts with… ▽ More The gravitational-wave spectrum emitted by a non-axisymmetric torus rotating at an angular velocity $Ω_T$, is derived in terms of a structure function representing a combination of sausage-tilt modes in the torus in the limit of an incompressible fluid. The analysis of the gravitational-wave spectrum is then applied to a model proposed recently, in which a highly magnetized torus interacts with a stellar mass, Kerr black hole via poloidal field lines that connect the torus and the horizon. The deformation of the torus results from global magnetic instabilities when the magnetic field strength inside the torus exceeds a few times $10^{15}$ Gauss. The dynamics of the system is illustrated using a non-MHD toy model. It is found that, quite generally, most of the gravitational-wave energy emitted is in the frequency range of sensitivity of LIGO and Virgo. △ Less

Submitted 4 July, 2005; originally announced July 2005.

Comments: 20 pages, 5 figures

Journal ref: New Astron.11:619-627,2006