A collapsar is a theoretical astrophysical model describing the core collapse of a rapidly rotating massive star, typically a Wolf–Rayet star with a mass exceeding 10 solar masses, into a black hole, where the failure of a successful supernova shock leads to the formation of a hyper-accreting disk that powers highly collimated relativistic jets responsible for long-duration gamma-ray bursts (GRBs).[1][2]The model was first proposed by Stan Woosley in 1993 as a mechanism linking gamma-ray bursts to the accretion of stellar material onto a newborn black hole, emphasizing the role of rapid rotation in enabling disk formation and beamed energy release.[1] In this scenario, the star's iron core collapses under gravity without producing an outgoing shock wave, resulting in a "failed supernova" where the outer envelope falls inward rather than being expelled.[3] The ensuing accretion disk, fed at rates of approximately 0.01–0.1 solar masses per second, generates intense neutrino emissions; annihilation of these neutrinos in the polar regions deposits energy of order 10^{51} erg, launching bipolar jets with Lorentz factors exceeding 100 that pierce the stellar envelope in seconds.[2] These jets, confined by the star's density stratification, produce the characteristic prompt emission of GRBs through internal shocks and subsequent afterglow via external shocks with the circumstellar medium.[4]Collapsars are closely associated with broad-lined Type Ic supernovae, often termed "hypernovae" due to their extreme energies, as evidenced by events like SN 1998bw linked to GRB 980425, where the jet energy also contributes to the supernova's luminosity.[3] While the core mechanism relies on specific progenitor conditions—low-metallicity environments to help preserve angular momentum, with models favoring metallicities below ~0.3 solar, though recent binary pathways enable formation near solarmetallicity—simulations show variability in outcomes, from energetic GRBs viewed on-axis to off-axis "orphan" afterglows or X-ray flashes.[5][6] Recent three-dimensional general relativistic magnetohydrodynamic simulations confirm the jet-launching process from the black hole-disk system, highlighting magnetic fields' role in collimation and the potential for r-process nucleosynthesis of heavy elements in the accretion flows.[7] As of 2025, binaryprogenitor models and updated simulations further extend the scenarios to solar-metallicity environments and enhanced heavy element production. This model unifies GRBs, supernovae, and black hole formation, providing a cornerstone for understanding the endpoints of massive star evolution.[8]
Definition and Characteristics
Core Concept
A collapsar refers to the theoretical scenario in which the core of a rapidly rotating massive star undergoes direct collapse into a black hole, bypassing neutron star formation and instead producing a hyper-accreting black hole surrounded by a centrifugally supported accretion disk that launches bipolar relativistic jets.[3] This process requires progenitor stars with zero-age main sequence masses typically in the range of 25 to 140 solar masses, where the iron core exceeds the Tolman-Oppenheimer-Volkoff limit for neutron star stability, leading to immediate black hole formation.[9] The rapid rotation of the progenitor is essential, providing the angular momentum necessary to form the disk and collimate the outflows along the rotational axis.[3]In the basic schematic of a collapsar, the central black hole—initially around 2 to 3 solar masses—is enveloped by a toroidal structure of infalling stellar material, where the disk's inner regions radiate efficiently via neutrinos.[3] Hyper-accretion onto the black hole occurs at rates of approximately 0.1 to 1 solar mass per second, converting a significant fraction of the gravitational binding energy into relativistic jets that propagate through the stellar envelope.[3] These jets, powered by mechanisms such as neutrino annihilation or magnetohydrodynamic processes, emerge from an evacuated polar channel created by the initial uninhibited infall along the rotation axis.[3]The collapsar model fundamentally differs from standard core-collapse supernovae, in which a neutrino-heated shock wave successfully explodes the star and leaves a neutron star remnant with isotropic ejecta.[3] Instead, collapsars constitute "failed supernovae," where the lack of a viable explosion directs the binding energy—potentially up to several times 10^51 ergs—predominantly into the beamed jets rather than a spherical outburst.[3] This model is closely linked to the production of long-duration gamma-ray bursts, observed in association with some such events.[3]
Progenitor Requirements
Collapsars arise from the core collapse of massive stars that develop iron cores sufficiently heavy to form black holes directly, bypassing the formation of a stable neutron star. The required initial mass range for these progenitors is typically 25 to 140 solar masses (M⊙), as this allows the stars to evolve into configurations where the core exceeds the Tolman-Oppenheimer-Volkoff limit for neutron star support, leading to black hole formation.[10] Stars below approximately 25 M⊙ tend to produce neutron stars or explode as ordinary supernovae, while those above 140 M⊙ may encounter pair-instability mechanisms that disrupt the star before collapse.[10]A low-metallicity environment is crucial for preserving the progenitor's mass and angular momentum, with metallicities Z ≲ 0.1 Z⊙ minimizing radiative mass loss through stellar winds. At higher metallicities, line-driven winds erode the envelope excessively, reducing the core's angular momentum and preventing the conditions for jet formation. This preference for low-Z progenitors explains the observed association of long gamma-ray bursts with metal-poor galaxies.[11]Rapid rotation is essential, with the specific angular momentum at the core needing to exceed j > 10^{16} cm² s⁻¹ to enable the formation of a centrifugally supported accretion disk around the black hole. This rotation can arise from single-star evolution with initial equatorial velocities near the critical limit or from binary interactions that spin up the star. Progenitors typically reach the Wolf-Rayet phase, where hydrogen and helium envelopes are stripped, leaving a carbon-oxygen (CO) core of 10 to 20 M⊙ that retains sufficient angular momentum for collapsar dynamics.[12][13]Magnetic fields play a key role in maintaining differential rotation throughout the progenitor's evolution by transporting angular momentum outward and stabilizing the core against excessive mixing. In models incorporating dynamo-generated fields, these magnetic torques help sustain the high core rotation needed for the collapsar engine, particularly in low-metallicity environments where other angular momentum loss mechanisms are reduced.[14]
Historical Development
Initial Proposal
The collapsar model was first proposed by S. E. Woosley in 1993 as a cosmological explanation for gamma-ray bursts (GRBs), positing that these events arise from the core collapse of rapidly rotating massive stars leading to failed Type Ib supernovae.[1] In this scenario, the progenitor is a Wolf-Rayet star with sufficient angular momentum to form a black hole surrounded by a massive accretion disk, rather than a successful explosion ejecting the envelope.[1] The term "collapsar" specifically refers to this black hole-accretion disk system resulting from the incomplete explosion.Prior to 1997, the distances to GRBs remained uncertain, with debates centering on whether they originated from Galactic neutron stars or more distant cosmological sources; Woosley's model was motivated by the need for ultra-relativistic outflows with Lorentz factors Γ > 100 to account for the observed isotropic energies and lack of Galactic counterparts.[1][15] The central engine relies on hyper-accretion rates of order 0.1–1 M⊙ s⁻¹ onto the black hole, powered by fallback of the stellar envelope, sustaining relativistic jets for durations of 10–100 seconds that match the observed lengths of long-duration GRBs.[1]The first empirical support for the collapsar model came in 1998 with the spatial and temporal coincidence of GRB 980425 and the Type Ic supernova SN 1998bw, a broad-lined event with high velocity ejecta indicative of core-collapse in a massive star. This association suggested an off-axis viewing of the GRB jet, as the event's isotropic-equivalent energy was orders of magnitude lower than typical cosmological GRBs, consistent with a beamed outflow from a collapsar not aligned with our line of sight.[16]
Key Advancements
Following the initial conceptualization of the collapsar model, significant progress was made in the late 1990s through computational simulations that explored the dynamics of jet formation. In 1999, MacFadyen and Woosley conducted the first two-dimensional (2D) hydrodynamics simulations of a rotating massive star undergoing core collapse, demonstrating that a black hole accretion disk could power relativistic jets along the rotation axis. These simulations revealed that energy deposition primarily occurs via neutrino-antineutrino annihilation in the polar regions above the disk, with supplementary contributions from magnetic processes in the accretion flow, leading to highly collimated outflows capable of penetrating the stellar envelope in seconds.[3] The total energy available for jet formation was estimated at 1-14 × 10^{51} ergs, sufficient to produce observed gamma-ray burst (GRB) luminosities if efficiently converted to relativistic motion.Into the 2000s, advancements incorporated three-dimensional (3D) effects and shifted emphasis toward magnetohydrodynamic (MHD) mechanisms for jet launching. Early 3D relativistic hydrodynamics simulations, such as those by Morsony et al. in 2007, illustrated how jets propagating through the stellar envelope develop instabilities, including kink modes that disrupt axisymmetry, while forming a hot, pressurized cocoon of shocked material that aids in collimating the jet. These simulations highlighted that the cocoon's expansion can influence jet directionality and energy distribution, with the jet head advancing at mildly relativistic speeds (~0.3c) before accelerating upon breakout.[17] Concurrently, MHD simulations by Proga et al. in 2003 demonstrated that magnetic fields in the accretion disk enable Blandford-Znajek-like extraction of rotational energy from the black hole, launching outflows with Lorentz factors up to 10 without relying heavily on neutrino heating, which proved inefficient for achieving ultra-relativistic speeds. This paradigm shift to MHD-dominated launching resolved limitations in neutrino-driven models, as magnetic torques more effectively accelerate and collimate polar material.[18]Integration of progenitor stellar structure with jet dynamics further refined the model. Heger, Langer, and Woosley in 2000 modeled the presupernova evolution of rotating massive stars (8-25 M_⊙), showing that differential rotation transports angular momentum outward, concentrating it in the core while allowing envelope loss via winds, which is crucial for forming the extended disks needed for jet collimation in collapsars. Their calculations indicated that progenitors with initial equatorial velocities of ~200-300 km/s retain sufficient core rotation (specific angular momentum j ~ 10^{16} cm²/s) to support bipolar outflows, linking stellar evolution parameters directly to the feasibility of successful jet production.[19]Early 2000s theoretical work also addressed observational challenges, such as the absence of radioactively powered supernova light curves in some GRBs. Rhoads in 2002 proposed that "choked" jets—those stalled within the stellar envelope—could explain this discrepancy by depositing energy into the cocoon without disrupting the star sufficiently for a bright Type Ib/c supernova, resulting in "failed" explosions with minimal ^{56}Ni production (~0.01 M_⊙ or less). In such scenarios, the choked jet's interaction with the envelope generates internal shocks and cocoon shocks that may produce detectable gamma rays but leave the outer layers intact, avoiding the characteristic radioactive decay signatures seen in successful GRB-supernova associations.[20]
Formation Mechanism
Stellar Evolution Leading to Collapse
Massive stars with initial masses exceeding approximately 25 solar masses (M⊙) spend their main sequence phase fusing hydrogen into helium in their cores via the CNO cycle, achieving surface rotation velocities up to 400 km s⁻¹ that induce efficient mixing and prevent the development of a distinct red giant branch.[21] These stars experience significant mass loss through radiatively driven winds, but in low-metallicity environments (Z ≲ 0.1 Z⊙), the reduced opacity leads to weaker winds (scaling roughly as Z^{0.5}), allowing retention of sufficient envelope mass and coreangular momentum essential for later collapsar conditions.[21]Corerotation remains rapid due to limited angular momentum transport during this phase, setting the stage for the necessary spin-up in subsequent evolution.As hydrogen exhaustion occurs, these stars progress through advanced nuclear burning stages: carbon burning lasts about 300 years, followed by neon (~1 year), oxygen (~0.5 years), and silicon (~1 day), collectively spanning roughly 100–1000 years while progressively building an iron core of ~1.5–2 M⊙.[21] During these phases, the envelope is stripped away, primarily through intensified winds in single-star scenarios or mass transfer in binaries, resulting in a compact Wolf-Rayet star with a carbon-oxygen (CO) or neon-oxygen (NeO) core. Angular momentum transport, mediated by convective instabilities and magnetic fields generated by the Tayler-Spruit dynamo, enforces near-rigid rotation and prevents excessive spin-down, preserving core angular velocities of order 10^6 rad s⁻¹ by countering viscous and magnetic braking. This mechanism ensures the core retains sufficient rotation for disk formation upon collapse, with major angular momentum redistribution occurring during carbon and oxygen shell burning.[21]A significant fraction of collapsar progenitors arise from binary evolution, where mass transfer from the primary strips the hydrogen envelope, yielding a rapidly rotating heliumstar of 8–12 M⊙ that evolves into a suitable Wolf-Rayet progenitor.[22] In such systems, accretion onto the companion can spin it up to near-breakup velocities, enhancing the likelihood of collapsar formation compared to single stars, which rely more heavily on intrinsic rotation and low-metallicity winds for envelope removal.[22]The pre-collapse endpoint is a Wolf-Rayet star with a final mass of ~10–15 M⊙, featuring a CO core (for helium core masses ~10–20 M⊙) or NeO core (for higher masses), a radius of approximately 10^{11} cm, and sufficient central angular momentum (j ~ 10^{16} cm² s⁻¹) to enable black hole formation with an accretion disk.[21] This configuration positions the star on the verge of dynamical instability, where the iron core's inability to support further fusion triggers collapse.[21]
Core Collapse and Black Hole Formation
In the collapsar model, the trigger for core collapse occurs when the iron core of a rapidly rotating massive star exceeds approximately 1.4 solar masses (M_⊙), approaching the Chandrasekhar limit of about 1.38 M_⊙, beyond which electron degeneracy pressure can no longer support the core against gravitational contraction.[23] At this stage, increasing densities promote electron capture reactions on protons and iron-group nuclei, which reduce the electron abundance and thus the degeneracy pressure, initiating a rapid implosion of the core at velocities reaching up to 0.25c.[23]During the collapse, the inner core, consisting of roughly 1.5 M_⊙ of iron-group material, compresses until nuclear densities are achieved, leading to an initial bounce that forms a hot proto-neutron star (PNS) with a radius of about 30 km.[23] However, in collapsar progenitors with cores exceeding 10 M_⊙ in helium mass, the high accretion rate—often surpassing 0.1 M_⊙ s⁻¹—causes the overlying material to rapidly fall back onto the PNS, overwhelming its Tolman-Oppenheimer-Volkoff limit and engulfing it within milliseconds to form a black hole of initial mass 2–3 M_⊙.[24][25]A significant portion of the gravitational binding energy released during this process, approximately 10⁵³ erg, is emitted primarily as neutrinos from the collapsing core and nascent PNS over the first few seconds, yet this outflow proves insufficient to revive the stalled shock and drive an explosion in these failed supernova scenarios.[23] Instead, most of the available energy, on the order of 10⁵¹ erg, remains trapped in the kinetic and thermal energy of the infalling envelope, which continues to accrete onto the black hole.[24]The resulting black hole features an event horizon with a Schwarzschild radius of roughly 10 km for masses in the 3–10 M_⊙ range, consistent with general relativistic expectations for stellar-mass black holes.[25] The spin parameter a, which quantifies the black hole's angular momentum (dimensionless, with a = 0 for non-spinning and a = 1 for extremal), typically starts at around 0.5 due to the progenitor's rotation and can reach 0.5–0.9 through subsequent accretion, enabling the formation of a centrifugally supported disk.[25]The entire process from core instability to black hole formation unfolds in less than 1 second, after which extended accretion from the stellar envelope sustains the system for seconds to tens of seconds, as detailed in progenitor evolution models leading to this phase.[24][25]
Jet Production and Dynamics
Accretion Disk and Jet Launching
In the collapsar model, following the formation of a black hole from the core collapse of a rapidly rotating massive star, fallback material with sufficient angular momentum accretes onto the black hole, forming a thick torus or accretion disk with a radius typically on the order of 100–1000 km. This disk arises from the rotational support of infalling gas that avoids immediate capture by the black hole, creating a centrifugally supported structure around the central engine.[26] The accretion process is primarily driven by viscosity, where magnetorotational instability (MRI) turbulence provides the dominant transport mechanism, enabling angular momentum redistribution and inward mass flow at rates of approximately 0.01–0.1 M_⊙ s⁻¹. These rates are sufficient to sustain high-energy outflows over the timescales relevant to long gamma-ray bursts (GRBs).[27]The gravitational binding energy released during accretion onto the black hole powers the central engine, with a total energy budget of around 10⁵² erg available from the fallback of several solar masses.[26] This energy is extracted through two primary mechanisms: magnetohydrodynamic (MHD) processes and neutrino-driven processes. In MHD scenarios, the Blandford-Znajek mechanism extracts rotational energy from the spinning black hole via magnetic field lines threading the event horizon, converting it into Poynting-flux-dominated outflows.[28] Neutrino processes involve the annihilation of neutrino-antineutrino pairs in the hot, dense disk and funnel regions above the poles, depositing thermal energy that can drive milder outflows or contribute to jet heating.[29]Jet launching in collapsars occurs primarily through these MHD and neutrino mechanisms, producing relativistic outflows with initial Lorentz factors up to Γ ≈ 100 and opening angles of about 5–10 degrees.[7] The Blandford-Znajek process is particularly efficient for rapidly spinning black holes (with dimensionless spin parameter a ≈ 0.5–1), generating magnetically accelerated jets that emerge from the disk-black hole interface.[30] Neutrino annihilation preferentially heats the low-density polar regions, potentially seeding or augmenting these jets by providing an initial thermal pressure gradient. Recent three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations as of 2025 confirm the MHD jet-launching process and reveal additional azimuthal and radial structure in the outflows.[31]Collimation of the jets is achieved through the pinching action of toroidal magnetic fields generated in the accretion disk, which confine the outflow into narrow, bipolar structures.[32] These fields, amplified by MRI-driven dynamo action, create a magnetic pressure gradient that accelerates the jet material along the axis, resulting in terminal Lorentz factors exceeding 300 in the innermost regions. This collimation ensures the jets remain relativistic and focused, facilitating their propagation through the stellar envelope.The accretion phase sustains the central engine for 10–100 seconds, aligning with the observed variability timescales of long GRBs.[33] During this period, the disk's viscous evolution and fallback supply maintain the necessary mass inflow, gradually tapering as outer material arrives or is ejected.[34]
Propagation Through Stellar Envelope
In the collapsar model, the relativistic jet launched from the central engine is initially confined by the surrounding stellar material, which shocks the envelope and forms a hot, pressurized cocoon of plasma.[35] This cocoon exerts sideways pressure on the jet, leading to recollimation shocks that can cause the jet to narrow and potentially precess due to instabilities or asymmetric energy deposition.[36] The interaction ensures that the jet remains collimated as it propagates, with the cocoon playing a key role in channeling the outflow along the rotation axis.[37]The head of the jet advances through the stellar envelope, typically spanning ~10^{11} cm, at velocities of approximately 0.1-0.3c, driven primarily by the pressure gradient from the expanding cocoon rather than the jet's internal energy alone.[20] This subrelativistic to mildly relativistic head speed allows the jet to drill through the dense material over timescales of several seconds, with the cocoon's thermal energy providing the necessary push against the ram pressure of the ambient medium.[37] During this phase, the jet-cocoon system evolves dynamically, with the cocoon's sound-crossing time influencing the overall propagation efficiency.[36]For the jet to successfully break out and produce an observable gamma-ray burst, its luminosity must exceed a critical value, typically L_j \gtrsim 10^{50} , \mathrm{erg , s^{-1}} for standard progenitors, enabling escape after ~10-20 seconds from launch.[38] Successful jets carve polar channels through the envelope, emerging with Lorentz factors of order 10-100, while in approximately 50% of cases, the jetsstall or choke within the star due to insufficient luminosity or dense progenitor structure, failing to produce GRBs but potentially generating X-ray transients from cocoon shocks.[39] These choked scenarios highlight the sensitivity of jet success to progenitor properties like rotation and envelopedensity.A significant fraction, around 10%, of the jet's total energy is deposited into heating the stellar envelope via cocoon expansion and shocks, contributing substantially to the luminosity of the associated core-collapse supernova by unbinding material and driving explosive outflows.[37] This energy transfer underscores the collapsar's dual role in producing both relativistic jets and broader explosive phenomena, with the heated envelope material potentially observable as early supernova light curves.[20]
Associated Astrophysical Phenomena
Long Gamma-Ray Bursts
In the collapsar model, long gamma-ray bursts (LGRBs) arise from the relativistic jets launched during the core collapse of massive, rapidly rotating stars, where the prompt emission phase is characterized by durations exceeding 2 seconds and isotropic-equivalent energies E_{\rm iso} typically ranging from $10^{52} to $10^{54} erg.[40] The prompt gamma-ray emission is primarily produced by internal shocks within the variable jet, where faster-moving ejecta catch up to slower material, dissipating energy through collisions that accelerate electrons to relativistic speeds. These electrons then radiate gamma rays via synchrotron processes in the amplified magnetic fields of the shocked regions, although photospheric emission from the jet's hot, optically thick base can also contribute significantly to the observed spectrum in some events.[41] In the collapsar framework, this emission occurs shortly after the jet breaks out of the stellar envelope, with the central black hole accretion disk powering the jet variability over timescales of seconds to minutes.Following the prompt phase, the LGRB afterglow emerges from external shocks as the decelerating jet interacts with the surrounding interstellar medium (ISM), generating synchrotron radiation across X-ray, optical, and radio wavelengths. The electron energy distribution in these shocks follows a power-law form N(\gamma_e) \propto \gamma_e^{-p} with a typical spectral index p \approx 2.2, leading to a broadbandspectrum that breaks at characteristic synchrotron frequencies and evolves predictably with time.[42] This model accounts for the observed flux decay and spectral softening, with the afterglow luminosity scaling as the jet sweeps up ISM mass, converting kinetic energy into radiation with efficiencies of about 10-20%.Jet collimation plays a crucial role in reconciling the immense apparent energies of LGRBs with more modest physical requirements, as the emission is beamed into a narrow cone with opening angle \theta_j \sim 5^\circ. The beaming factor, approximated as f_b \approx \theta_j^2 / 2, reduces the true radiated energy to around $10^{51} erg, aligning with the gravitational binding energy available from a massive star's core collapse and explaining the observed event rate without invoking extreme isotropic outputs. For observers positioned off-axis relative to the jetcone, such as in the case of GRB 980425, the emission appears weaker and more extended in duration due to the geometric dilution of the relativistic beaming, resulting in lower peak fluxes and delayed light curves compared to on-axis views.The rapid variability in LGRB prompt emission, featuring millisecond-scale pulses, stems from fluctuations in the central engine driven by instabilities in the hyperaccreting disk around the newborn black hole.[43] In collapsar disks, dead zones—regions of low ionization—can form due to insufficient MRI turbulence, leading to nonsteady accretion episodes that episodically inject material into the jet, producing the observed spiky light curves.[44] These disk instabilities, often triggered by gravitational fragmentation or magnetic reconnection, modulate the jet's Lorentz factor and luminosity on short timescales, directly imprinting the engine's dynamics onto the gamma-ray output.[45]
Core-Collapse Supernovae
In collapsar events, the supernova explosion arises from the energetic ejection of the stellar envelope following core collapse and black hole formation, primarily powered by the radioactive decay of ^{56}Ni synthesized during the explosion. Typical ^{56}Ni masses range from ~0.1 to 1 M_\sun, producing light curves that peak at 10–20 days after explosion with peak luminosities 10–100 times greater than those of standard core-collapse supernovae (CCSNe), reaching ~10^{43}–10^{44} erg s^{-1}. This enhanced brightness stems from both the higher nickel yield and the increased kinetic energy imparted to the ejecta, distinguishing these events from ordinary Type II or Ib/Ic supernovae.[46][47][48]These supernovae are identified as hypernovae, a subtype of broad-lined Type Ib/c events linked to stripped-envelope progenitors, where the explosion energy totals ~10^{52} erg—roughly an order of magnitude above the ~10^{51} erg of standard CCSNe. The broad spectral lines, with expansion velocities exceeding 30,000 km/s, reflect the acceleration of ejecta by relativistic jets originating from the central accretion disk around the newborn black hole. This jet interaction energizes the polar regions, leading to highly asymmetric outflow structures: bipolarejecta with polar velocities up to ~0.16c (~47,000 km/s) contrasted against slower equatorial expansion at ~10,000 km/s, imprinting distinct morphological and kinematic signatures.[3][49][50]Observationally, hypernovae exhibit reddened spectra due to line blanketing by heavy elements produced in the explosive nucleosynthesis, with light curves extending over ~100 days before fading. These traits are tied to progenitors that have undergone significant mass loss, shedding their hydrogen and much of their helium envelopes via winds or binary interactions, resulting in spectra dominated by oxygen, silicon, and iron-group lines without hydrogen features. Events such as SN 1998bw, associated with GRB 980425, serve as archetypal examples of these jet-driven hypernovae.[46]In certain collapsar scenarios, particularly those with high fallback rates where most ejecta accretes onto the black hole, no prominent supernova signature may emerge, yielding "failed" explosions invisible at optical wavelengths. However, faint emissions could still arise from low ^{56}Ni yields (~0.01 M_\sun or less) or dust obscuration in dense star-forming environments, potentially detectable in infrared or through associated gamma-ray burst afterglows if beaming aligns favorably.[51][52]
Observational Evidence
GRB-Supernova Associations
The association between gamma-ray bursts (GRBs) and supernovae (SNe) was first observationally established in 1998 with the underluminous GRB 980425, occurring at a redshift of z = 0.0085, and the broad-lined Type Ic supernova SN 1998bw, which exhibited unusual brightness and high expansion velocities consistent with a collapsar progenitor.[53] This event provided initial evidence linking long-duration GRBs to the deaths of massive stars, though its low luminosity raised questions about whether it represented a typical collapsar scenario.[54]Subsequent observations strengthened this connection through key events, including GRB 030329 at z \approx 0.169, whose afterglow revealed a spectroscopic match with the broad-lined Type Ic SN 2003dh, showing clear supernova features emerging from the fading GRB light approximately 10 days post-burst.[55] Another pivotal low-luminosity pair was GRB 060218 at z = 0.033, associated with the subluminous, rapidly evolving SN 2006aj, which displayed a wavelength-dependent early light curve component indicative of shock breakout and further supported the GRB-SN link in nearby events.[56]By 2025, approximately 61 confirmed GRB-SN associations have been identified, all involving long GRBs paired exclusively with broad-lined Type Ic (Ic-BL) supernovae, with no detections of Type II or Type Ia events in these systems.[57] This specificity underscores the collapsar model's prediction that only certain stripped-envelope core-collapse events produce relativistic jets capable of generating GRBs.[58]The scarcity of GRB detections among the broader population of Ic-BL supernovae—over 60 such events observed without associated GRBs—provides evidence for viewing angle dependence, where off-axis orientations relative to the collapsar jet may render the GRB emission undetectable while still producing an observable supernova.[59] For instance, low-luminosity cases like SN 2006aj highlight how marginal jet alignments can result in weaker GRB signatures.[56]Host galaxies of these GRB-SN associations are predominantly low-metallicity, star-forming dwarf galaxies, aligning with the environmental preferences for massive, Wolf-Rayet progenitors required by the collapsar mechanism.[60] This distribution reflects the need for reduced metal content to facilitate jet formation and escape.[61]
Recent Observations and Simulations
In the 2020s, advanced three-dimensional general-relativistic magnetohydrodynamic (3D GRMHD) simulations have provided deeper insights into collapsar dynamics, revealing multi-scale outflows that include powerful disk winds. These models, extending from black hole formation to jet breakout, demonstrate that accretion disks around the central black hole launch relativistic jets and milder sub-relativistic winds, with isotropic-equivalent energies in disk winds reaching approximately $10^{52} erg, contributing significantly to the overall energetics of gamma-ray burst (GRB) events.[31][62]Observations from the Einstein Probe mission have identified a class of relativistic supernova transients with multiband light curves that align with the predicted multicomponent outflows from collapsar models, featuring distinct prompt emission, a plateau phase, and extended afterglow. These transients, such as EP240414a, exhibit rapid X-ray rises followed by optical counterparts consistent with broad-lined Type Ic supernovae, supporting the structured outflow paradigm where jets and winds interact with the stellar envelope.[31][63]The supernova SN 2025kg, associated with the Einstein Probe fast X-ray transient EP 250108a, exemplifies a jet-driven explosion with early high-velocity spectral features indicative of a collapsar central engine. At a redshift of z = 0.176, this broad-lined Type Ic supernova displays rapid light-curve evolution and expansion velocities exceeding 30,000 km/s in its early spectra, suggesting a choked relativistic jet energizing the ejecta and interacting with circumstellar material.[64][63]Searches for heavy element production in GRB afterglows using the James Webb Space Telescope (JWST) have explored r-process signatures potentially originating from collapsar accretion disks, as predicted by models of neutrino-driven winds. For instance, JWST observations of the supernova linked to GRB 221009A revealed no detectable r-process enrichment in the near-infrared spectra, constraining the disk conditions required for rapid neutron capture in collapsar environments and highlighting the variability in outflow compositions.[65][66]Recent simulations have addressed longstanding challenges in collapsar jet production, particularly the baryon loading problem, by incorporating magnetic acceleration mechanisms that enable efficient relativistic outflows with reduced mass entrainment. High-magnetization disk models show that magneto-centrifugal winds and Blandford-Znajek processes launch jets with Lorentz factors \Gamma \gtrsim 100 while minimizing baryon pollution, resolving discrepancies between theoretical predictions and observed GRB prompt efficiencies.[31]
Theoretical Implications and Challenges
Role in Heavy Element Nucleosynthesis
In collapsar events, the rapid accretion of stellar material onto a central black hole forms a hyperaccretion disk that serves as a potential site for r-process nucleosynthesis, producing heavy elements beyond the iron peak through neutron capture on seed nuclei. Neutron-rich ejecta originate from disk winds and neutrino-irradiated regions, where the electron fraction Y_e ranges from approximately 0.1 to 0.4, providing conditions conducive to rapid neutron capture by maintaining a high neutron-to-proton ratio.[67][68] These outflows, driven by magnetic fields and neutrino heating, expel material with sufficient neutron excess to synthesize elements up to the third r-process peak, including lanthanides and actinides.[69]Recent simulations indicate that r-process yields vary significantly depending on disk viscosity, magnetic field strength, and neutrino interactions; while some models yield 0.01 to 0.1 solar masses (M_\odot) of lanthanide and actinide elements—comparable in mass to those from neutron star mergers—these occur under favorable conditions at a rate of roughly 1 per gigayear per galaxy, with many cases producing negligible heavy elements beyond the first r-process peak.[68][70][69] This production efficiency arises from the disk's high accretion rates (0.1–1 M_\odot s^{-1}) and the ejection of neutron-rich material via magnetically driven winds in select scenarios.[71] The cumulative output from collapsars can account for a potentially significant fraction of the r-process elements observed in the universe, particularly when integrated over cosmic time and assuming robust yields in low-metallicity environments, though recent models suggest contributions may be lower than previously estimated (e.g., <80%).[68][72]Neutrino physics plays a crucial role in enabling this nucleosynthesis, as the accretion disk emits high-luminosity neutrinos with total energies around $10^{52} erg s^{-1} .[73] These neutrinos drive collective oscillations above the disk, altering flavor content and reducing the electron fraction Y_e in the outflows, which enhances neutronization and supports robust r-process conditions. Without such oscillations, the ejecta would be proton-richer, limiting heavy element production.Observational evidence for collapsar-driven r-process includes blue spectral components in some GRB-associated supernovae, suggesting lanthanide-free inner ejecta from less neutron-rich disk regions, contrasting with redder, lanthanide-opaque outer layers.[74]Collapsars contribute significantly to galactic chemical enrichment, particularly in low-metallicity environments where they dominate r-process production and explain abundance patterns in very metal-poor stars, such as enhanced europium relative to iron ([Eu/Fe] trends).[70][75] This role aligns with observations of early universe metal enrichment, where collapsars provide prompt injection of heavy elements into the interstellar medium without requiring long delay times.[68]
Alternatives and Open Questions
One alternative to the collapsar model for long gamma-ray bursts (GRBs) is the supranova scenario, in which a neutron star forms following a core-collapse supernova and subsequently collapses into a black hole after accreting material over days to weeks, powering the GRB. This model predicts a temporal delay between the supernova and the GRB, but it has been disfavored by observations showing near-simultaneous occurrences in associated events, with no confirmed delayed GRBs detected.Another competing framework is the millisecond magnetar model, where a rapidly rotating, highly magnetized neutron star powers the GRB through magnetic dipole spin-down, injecting energy into a surrounding outflow.[76] While this mechanism can explain the energetics and durations of some GRBs and associated supernovae, it faces challenges in accounting for events requiring prolonged accretion or black hole signatures, such as those with isotropic energies exceeding 10^{52} erg, which favor collapsar-like black hole engines.[77]Several open questions persist within the collapsar framework. The success rate of jet breakout from the stellar envelope remains uncertain, estimated at 10-50% depending on progenitor rotation and magnetic field strength, with many simulations indicating that jets stall in a significant fraction of cases.[78] The role of three-dimensional hydrodynamic instabilities, such as the standing accretion shock instability and Kelvin-Helmholtz modes, in potentially reviving or choking jets during propagation is still under investigation through advanced simulations.[7] Additionally, the evolution of the central black hole's spin, influenced by accretion torques and jet launching, lacks precise observational constraints, complicating predictions of jet power and collimation.[79]Key challenges include the observed paucity of GRBs in high-metallicity galaxies at low redshift, attributed to metallicity-dependent angular momentum transport that hinders rapid rotation in progenitors, though the exact biasthreshold remains debated.[80] Furthermore, the existence of engine-driven supernovae lacking GRB counterparts, such as SN 2009bb with its relativistic outflow but no gamma-ray detection, suggests choked jets or off-axis viewing in collapsar events, implying a hidden population of failed GRBs.[81]Future multi-messenger observations offer promising tests. Gravitational waves from black hole formation and disk instabilities in collapsars could be detectable by LIGO/Virgo up to ~15 Mpc in optimistic models.[79] Complementary insights may come from joint detections using JWST for high-redshift afterglowspectroscopy and CTA for very-high-energy gamma rays, enabling constraints on jet composition and environment.[82]