Collapsars

Gamma-ray bursts and explosions in "failed supernovae"

A. I. MacFadyen, S. E. Woosley

    Research output: Contribution to journalArticle

    Abstract

    Using a two-dimensional hydrodynamics code (PROMETHEUS), we explore the continued evolution of rotating helium stars, Mα ≳ 10 M⊙, in which iron-core collapse does not produce a successful outgoing shock but instead forms a black hole of 2-3 M⊙. The model explored in greatest detail is the 14 M helium core of a 35 M main-sequence star. The outcome is sensitive to the angular momentum. For J16 ≡ j/(1016 cm2 s-1) ≲ 3, material falls into the black hole almost uninhibited. No outflows are expected. For j16 ≳ 20, the infalling matter is halted by centrifugal force outside 1000 km where neutrino losses are negligible. The equatorial accretion rate is very low, and explosive oxygen burning may power a weak equatorial explosion. For 3 ≲ j16 ≲ 20, however, a reasonable value for such stars, a compact disk forms at a radius at which the gravitational binding energy can be efficiently radiated as neutrinos or converted to beamed outflow by magnetohydrodynamical (MHD) processes. These are the best candidates for producing gamma-ray bursts (GRBs). Here we study the formation of such a disk, the associated flow patterns, and the accretion rate for disk viscosity parameter α ≈ 0.001 and 0.1. Infall along the rotational axis is initially uninhibited, and an evacuated channel opens during the first few seconds. Meanwhile the black hole is spun up by the accretion (to a ≈ 0.9), and energy is dissipated in the disk by MHD processes and radiated by neutrinos. For the α = 0.1 model, appreciable energetic outflows develop between polar angles of 30° and 45°. These outflows, powered by viscous dissipation in the disk, have an energy of up to a few times 1051 ergs and a mass ∼ 1 M⊙ and are rich in 56Ni. They constitute a supernova-like explosion by themselves. Meanwhile accretion through the disk is maintained for approximately 10-20 s but is time variable (±30%) because of hydrodynamical instabilities at the outer edge in a region where nuclei are experiencing photodisintegration. Because the efficiency of neutrino energy deposition is sensitive to the accretion rate, this instability leads to highly variable energy deposition in the polar regions. Some of this variability, which has significant power at 50 ms and overtones, may persist in the time structure of the burst. During the time followed, the average accretion rate for the standard α = 0.1 and j16 = 10 model is 0.07 M⊙ s-1. The total energy deposited along the rotational axes by neutrino annihilation is (1-14) × 1051 ergs, depending upon the evolution of the Kerr parameter and uncertain neutrino efficiencies. Simulated deposition of energy in the polar regions, at a constant rate of 5 × 1050 ergs s-1 per pole, results in strong relativistic outflow jets beamed to about 1% of the sky. These jets may be additionally modulated by instabilities in the sides of the "nozzle" through which they flow. The jets blow aside the accreting material, remain highly focused, and are capable of penetrating the star in ≃ 10 s. After the jet breaks through the surface of the star, highly relativistic flow can emerge. Because of the sensitivity of the mass ejection and jets to accretion rate, angular momentum, and disk viscosity, and the variation of observational consequences with viewing angle, a large range of outcomes is possible, ranging from bright GRBs like GRB 971214 to faint GRB-supernovae like SN 1998bw. X-ray precursors are also possible as the jet first breaks out of the star. While only a small fraction of supernovae make GRBs, we predict that collapsars will always make supernovae similar to SN 1998bw. However, hard, energetic GRBs shorter than a few seconds will be difficult to produce in this model and may require merging neutron stars and black holes for their explanation.

    Original languageEnglish (US)
    Pages (from-to)262-289
    Number of pages28
    JournalAstrophysical Journal Letters
    Volume524
    Issue number1 PART 1
    StatePublished - Oct 10 1999

    Fingerprint

    gamma ray bursts
    supernovae
    explosions
    explosion
    accretion
    outflow
    neutrinos
    energy
    erg
    polar region
    angular momentum
    helium
    stars
    energetics
    viscosity
    polar regions
    gravitational binding energy
    flow pattern
    explosive
    rate

    Keywords

    • Accretion, accretion disks
    • Black hole physics
    • Gamma rays: bursts
    • Supernovae: general

    ASJC Scopus subject areas

    • Space and Planetary Science

    Cite this

    MacFadyen, A. I., & Woosley, S. E. (1999). Collapsars: Gamma-ray bursts and explosions in "failed supernovae". Astrophysical Journal Letters, 524(1 PART 1), 262-289.

    Collapsars : Gamma-ray bursts and explosions in "failed supernovae". / MacFadyen, A. I.; Woosley, S. E.

    In: Astrophysical Journal Letters, Vol. 524, No. 1 PART 1, 10.10.1999, p. 262-289.

    Research output: Contribution to journalArticle

    MacFadyen, AI & Woosley, SE 1999, 'Collapsars: Gamma-ray bursts and explosions in "failed supernovae"', Astrophysical Journal Letters, vol. 524, no. 1 PART 1, pp. 262-289.
    MacFadyen AI, Woosley SE. Collapsars: Gamma-ray bursts and explosions in "failed supernovae". Astrophysical Journal Letters. 1999 Oct 10;524(1 PART 1):262-289.
    MacFadyen, A. I. ; Woosley, S. E. / Collapsars : Gamma-ray bursts and explosions in "failed supernovae". In: Astrophysical Journal Letters. 1999 ; Vol. 524, No. 1 PART 1. pp. 262-289.
    @article{2f1fbf22d5d242c9bf4c16a7f9ddf013,
    title = "Collapsars: Gamma-ray bursts and explosions in {"}failed supernovae{"}",
    abstract = "Using a two-dimensional hydrodynamics code (PROMETHEUS), we explore the continued evolution of rotating helium stars, Mα ≳ 10 M⊙, in which iron-core collapse does not produce a successful outgoing shock but instead forms a black hole of 2-3 M⊙. The model explored in greatest detail is the 14 M⊙ helium core of a 35 M⊙ main-sequence star. The outcome is sensitive to the angular momentum. For J16 ≡ j/(1016 cm2 s-1) ≲ 3, material falls into the black hole almost uninhibited. No outflows are expected. For j16 ≳ 20, the infalling matter is halted by centrifugal force outside 1000 km where neutrino losses are negligible. The equatorial accretion rate is very low, and explosive oxygen burning may power a weak equatorial explosion. For 3 ≲ j16 ≲ 20, however, a reasonable value for such stars, a compact disk forms at a radius at which the gravitational binding energy can be efficiently radiated as neutrinos or converted to beamed outflow by magnetohydrodynamical (MHD) processes. These are the best candidates for producing gamma-ray bursts (GRBs). Here we study the formation of such a disk, the associated flow patterns, and the accretion rate for disk viscosity parameter α ≈ 0.001 and 0.1. Infall along the rotational axis is initially uninhibited, and an evacuated channel opens during the first few seconds. Meanwhile the black hole is spun up by the accretion (to a ≈ 0.9), and energy is dissipated in the disk by MHD processes and radiated by neutrinos. For the α = 0.1 model, appreciable energetic outflows develop between polar angles of 30° and 45°. These outflows, powered by viscous dissipation in the disk, have an energy of up to a few times 1051 ergs and a mass ∼ 1 M⊙ and are rich in 56Ni. They constitute a supernova-like explosion by themselves. Meanwhile accretion through the disk is maintained for approximately 10-20 s but is time variable (±30{\%}) because of hydrodynamical instabilities at the outer edge in a region where nuclei are experiencing photodisintegration. Because the efficiency of neutrino energy deposition is sensitive to the accretion rate, this instability leads to highly variable energy deposition in the polar regions. Some of this variability, which has significant power at 50 ms and overtones, may persist in the time structure of the burst. During the time followed, the average accretion rate for the standard α = 0.1 and j16 = 10 model is 0.07 M⊙ s-1. The total energy deposited along the rotational axes by neutrino annihilation is (1-14) × 1051 ergs, depending upon the evolution of the Kerr parameter and uncertain neutrino efficiencies. Simulated deposition of energy in the polar regions, at a constant rate of 5 × 1050 ergs s-1 per pole, results in strong relativistic outflow jets beamed to about 1{\%} of the sky. These jets may be additionally modulated by instabilities in the sides of the {"}nozzle{"} through which they flow. The jets blow aside the accreting material, remain highly focused, and are capable of penetrating the star in ≃ 10 s. After the jet breaks through the surface of the star, highly relativistic flow can emerge. Because of the sensitivity of the mass ejection and jets to accretion rate, angular momentum, and disk viscosity, and the variation of observational consequences with viewing angle, a large range of outcomes is possible, ranging from bright GRBs like GRB 971214 to faint GRB-supernovae like SN 1998bw. X-ray precursors are also possible as the jet first breaks out of the star. While only a small fraction of supernovae make GRBs, we predict that collapsars will always make supernovae similar to SN 1998bw. However, hard, energetic GRBs shorter than a few seconds will be difficult to produce in this model and may require merging neutron stars and black holes for their explanation.",
    keywords = "Accretion, accretion disks, Black hole physics, Gamma rays: bursts, Supernovae: general",
    author = "MacFadyen, {A. I.} and Woosley, {S. E.}",
    year = "1999",
    month = "10",
    day = "10",
    language = "English (US)",
    volume = "524",
    pages = "262--289",
    journal = "Astrophysical Journal Letters",
    issn = "2041-8205",
    publisher = "IOP Publishing Ltd.",
    number = "1 PART 1",

    }

    TY - JOUR

    T1 - Collapsars

    T2 - Gamma-ray bursts and explosions in "failed supernovae"

    AU - MacFadyen, A. I.

    AU - Woosley, S. E.

    PY - 1999/10/10

    Y1 - 1999/10/10

    N2 - Using a two-dimensional hydrodynamics code (PROMETHEUS), we explore the continued evolution of rotating helium stars, Mα ≳ 10 M⊙, in which iron-core collapse does not produce a successful outgoing shock but instead forms a black hole of 2-3 M⊙. The model explored in greatest detail is the 14 M⊙ helium core of a 35 M⊙ main-sequence star. The outcome is sensitive to the angular momentum. For J16 ≡ j/(1016 cm2 s-1) ≲ 3, material falls into the black hole almost uninhibited. No outflows are expected. For j16 ≳ 20, the infalling matter is halted by centrifugal force outside 1000 km where neutrino losses are negligible. The equatorial accretion rate is very low, and explosive oxygen burning may power a weak equatorial explosion. For 3 ≲ j16 ≲ 20, however, a reasonable value for such stars, a compact disk forms at a radius at which the gravitational binding energy can be efficiently radiated as neutrinos or converted to beamed outflow by magnetohydrodynamical (MHD) processes. These are the best candidates for producing gamma-ray bursts (GRBs). Here we study the formation of such a disk, the associated flow patterns, and the accretion rate for disk viscosity parameter α ≈ 0.001 and 0.1. Infall along the rotational axis is initially uninhibited, and an evacuated channel opens during the first few seconds. Meanwhile the black hole is spun up by the accretion (to a ≈ 0.9), and energy is dissipated in the disk by MHD processes and radiated by neutrinos. For the α = 0.1 model, appreciable energetic outflows develop between polar angles of 30° and 45°. These outflows, powered by viscous dissipation in the disk, have an energy of up to a few times 1051 ergs and a mass ∼ 1 M⊙ and are rich in 56Ni. They constitute a supernova-like explosion by themselves. Meanwhile accretion through the disk is maintained for approximately 10-20 s but is time variable (±30%) because of hydrodynamical instabilities at the outer edge in a region where nuclei are experiencing photodisintegration. Because the efficiency of neutrino energy deposition is sensitive to the accretion rate, this instability leads to highly variable energy deposition in the polar regions. Some of this variability, which has significant power at 50 ms and overtones, may persist in the time structure of the burst. During the time followed, the average accretion rate for the standard α = 0.1 and j16 = 10 model is 0.07 M⊙ s-1. The total energy deposited along the rotational axes by neutrino annihilation is (1-14) × 1051 ergs, depending upon the evolution of the Kerr parameter and uncertain neutrino efficiencies. Simulated deposition of energy in the polar regions, at a constant rate of 5 × 1050 ergs s-1 per pole, results in strong relativistic outflow jets beamed to about 1% of the sky. These jets may be additionally modulated by instabilities in the sides of the "nozzle" through which they flow. The jets blow aside the accreting material, remain highly focused, and are capable of penetrating the star in ≃ 10 s. After the jet breaks through the surface of the star, highly relativistic flow can emerge. Because of the sensitivity of the mass ejection and jets to accretion rate, angular momentum, and disk viscosity, and the variation of observational consequences with viewing angle, a large range of outcomes is possible, ranging from bright GRBs like GRB 971214 to faint GRB-supernovae like SN 1998bw. X-ray precursors are also possible as the jet first breaks out of the star. While only a small fraction of supernovae make GRBs, we predict that collapsars will always make supernovae similar to SN 1998bw. However, hard, energetic GRBs shorter than a few seconds will be difficult to produce in this model and may require merging neutron stars and black holes for their explanation.

    AB - Using a two-dimensional hydrodynamics code (PROMETHEUS), we explore the continued evolution of rotating helium stars, Mα ≳ 10 M⊙, in which iron-core collapse does not produce a successful outgoing shock but instead forms a black hole of 2-3 M⊙. The model explored in greatest detail is the 14 M⊙ helium core of a 35 M⊙ main-sequence star. The outcome is sensitive to the angular momentum. For J16 ≡ j/(1016 cm2 s-1) ≲ 3, material falls into the black hole almost uninhibited. No outflows are expected. For j16 ≳ 20, the infalling matter is halted by centrifugal force outside 1000 km where neutrino losses are negligible. The equatorial accretion rate is very low, and explosive oxygen burning may power a weak equatorial explosion. For 3 ≲ j16 ≲ 20, however, a reasonable value for such stars, a compact disk forms at a radius at which the gravitational binding energy can be efficiently radiated as neutrinos or converted to beamed outflow by magnetohydrodynamical (MHD) processes. These are the best candidates for producing gamma-ray bursts (GRBs). Here we study the formation of such a disk, the associated flow patterns, and the accretion rate for disk viscosity parameter α ≈ 0.001 and 0.1. Infall along the rotational axis is initially uninhibited, and an evacuated channel opens during the first few seconds. Meanwhile the black hole is spun up by the accretion (to a ≈ 0.9), and energy is dissipated in the disk by MHD processes and radiated by neutrinos. For the α = 0.1 model, appreciable energetic outflows develop between polar angles of 30° and 45°. These outflows, powered by viscous dissipation in the disk, have an energy of up to a few times 1051 ergs and a mass ∼ 1 M⊙ and are rich in 56Ni. They constitute a supernova-like explosion by themselves. Meanwhile accretion through the disk is maintained for approximately 10-20 s but is time variable (±30%) because of hydrodynamical instabilities at the outer edge in a region where nuclei are experiencing photodisintegration. Because the efficiency of neutrino energy deposition is sensitive to the accretion rate, this instability leads to highly variable energy deposition in the polar regions. Some of this variability, which has significant power at 50 ms and overtones, may persist in the time structure of the burst. During the time followed, the average accretion rate for the standard α = 0.1 and j16 = 10 model is 0.07 M⊙ s-1. The total energy deposited along the rotational axes by neutrino annihilation is (1-14) × 1051 ergs, depending upon the evolution of the Kerr parameter and uncertain neutrino efficiencies. Simulated deposition of energy in the polar regions, at a constant rate of 5 × 1050 ergs s-1 per pole, results in strong relativistic outflow jets beamed to about 1% of the sky. These jets may be additionally modulated by instabilities in the sides of the "nozzle" through which they flow. The jets blow aside the accreting material, remain highly focused, and are capable of penetrating the star in ≃ 10 s. After the jet breaks through the surface of the star, highly relativistic flow can emerge. Because of the sensitivity of the mass ejection and jets to accretion rate, angular momentum, and disk viscosity, and the variation of observational consequences with viewing angle, a large range of outcomes is possible, ranging from bright GRBs like GRB 971214 to faint GRB-supernovae like SN 1998bw. X-ray precursors are also possible as the jet first breaks out of the star. While only a small fraction of supernovae make GRBs, we predict that collapsars will always make supernovae similar to SN 1998bw. However, hard, energetic GRBs shorter than a few seconds will be difficult to produce in this model and may require merging neutron stars and black holes for their explanation.

    KW - Accretion, accretion disks

    KW - Black hole physics

    KW - Gamma rays: bursts

    KW - Supernovae: general

    UR - http://www.scopus.com/inward/record.url?scp=0033543855&partnerID=8YFLogxK

    UR - http://www.scopus.com/inward/citedby.url?scp=0033543855&partnerID=8YFLogxK

    M3 - Article

    VL - 524

    SP - 262

    EP - 289

    JO - Astrophysical Journal Letters

    JF - Astrophysical Journal Letters

    SN - 2041-8205

    IS - 1 PART 1

    ER -