Source author record

Omer Bromberg

Omer Bromberg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2020arXiv

Particle acceleration in kink-unstable jets

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.

preprint2020arXiv

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

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.

preprint2016arXiv

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

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.

preprint2015arXiv

Gamma-Ray Bursts as Sources of Strong Magnetic Fields

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.

preprint2015arXiv

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

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.

preprint2014arXiv

On the Composition of GRBs' Collapsar Jets

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.

preprint2014arXiv

The dynamics of a highly magnetized jet propagating inside a star

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.

preprint2013arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2011arXiv

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

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.

preprint2011arXiv

Are low luminosity GRBs generated by relativistic jets?

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.

preprint2011arXiv

New insights on the duration distribution of long GRBs from Collapsars

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.

preprint2011arXiv

Sub-Photospheric Emission from Relativistic Radiation Mediated Shocks in GRBs

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.

preprint2011arXiv

The propagation of relativistic jets in external media

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.

preprint2007arXiv

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

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.