Jets from Scratch: A 3D Dynamo Origin of Long Gamma-Ray Burst Jets

This paper presents the first 3D general-relativistic magnetohydrodynamic simulations demonstrating that an in situ accretion-disk dynamo can convert toroidal magnetic fields into coherent poloidal structures, thereby launching highly variable, wobbling relativistic jets with sufficient power to explain long gamma-ray bursts without requiring a pre-existing large-scale poloidal field.

The Gist

Imagine a massive star, much larger than our Sun, running out of fuel. Instead of gently fading away, its core collapses under its own weight, creating a tiny, incredibly dense object called a black hole. Usually, this event is a quiet implosion. But sometimes, it explodes with the energy of a billion suns, shooting out two beams of light so powerful they can be seen across the entire universe. These are called Long Gamma-Ray Bursts (LGRBs).

For decades, scientists have been puzzled by a specific question: How does the black hole know to shoot these beams?

To fire a beam like a laser, the black hole needs a specific kind of magnetic "wire" (a large-scale magnetic field) threaded through it. The problem is that the dying star mostly creates "twisted" magnetic fields (like a tangled rubber band), not the straight, organized wires needed to launch the beam.

This paper, "Jets from Scratch," solves that mystery by showing that the black hole's surrounding disk of gas acts like a giant, self-organizing machine that untangles the mess and builds the necessary wires from scratch.

Here is how the authors explain this process using simple analogies:

1. The Tangled Rubber Band (The Problem)

Before the star dies, its rotation creates magnetic fields that are mostly toroidal. Imagine a rubber band stretched around a ball; it goes around the equator but doesn't go from the North Pole to the South Pole.

• The Issue: To launch a jet, you need a field that goes from pole to pole (poloidal).
• The Old Theory: Scientists thought the star might have had a strong pole-to-pole field hidden inside that survived the collapse. But calculations show that by the time the black hole forms, that field is usually too weak or too messy to work.

2. The Kitchen Blender (The Solution: The Dynamo)

The authors ran the most detailed 3D computer simulations ever done for this scenario. They started with the "tangled rubber band" (the weak, twisted magnetic field) and watched what happened as gas swirled around the new black hole.

They found that the swirling gas acts like a kitchen blender:

• The Spin: As the gas spins around the black hole, it stretches and twists the magnetic "rubber bands."
• The Dynamo Effect: This stretching and twisting creates a feedback loop (a dynamo). Just like a bicycle dynamo generates electricity by spinning a magnet, the spinning gas generates a new, organized magnetic field that points from the North Pole to the South Pole.
• The Result: Within a few seconds, this "blender" creates strong, straight magnetic loops out of the chaos.

3. The Garden Hose and the Kink (The Jet Launch)

Once these new magnetic loops form, they get pulled toward the black hole.

• The Connection: These loops connect to the spinning black hole.
• The Launch: The black hole acts like a spinning top. Because the magnetic "hose" is attached to it, the spinning motion twists the hose, shooting out a powerful stream of energy (the jet) along the poles.
• The Wobble: The paper notes that these jets aren't perfectly straight. Because the gas falling in comes from random directions, it pushes the disk around. This makes the jet wobble like a garden hose that hasn't been clamped down tight. This wobble explains why the light from these bursts flickers and changes so rapidly.

4. The "Striped" Pattern

The simulations showed something fascinating: the magnetic field doesn't just stay one color (one direction). It flips back and forth.

• The Analogy: Imagine a zebra. The jet isn't just one solid beam; it's a striped jet with alternating magnetic directions.
• The Implication: These stripes might be the reason why the light curves of these bursts look the way they do, with rapid spikes and dips.

Why This Matters

The paper proves that you don't need the star to have a perfect magnetic field to begin with. Even if the star starts with a weak, messy magnetic field, the accretion disk (the swirling gas around the black hole) can generate the necessary power on its own.

• Robustness: This means that almost any rapidly spinning massive star that forms a disk has the potential to create a gamma-ray burst. It doesn't rely on a "lucky" initial magnetic setup.
• Timing: The process happens quickly, launching the jet within a few seconds of the black hole's birth.

In short, the universe doesn't need a pre-made laser pointer to create a gamma-ray burst. It just needs a spinning black hole and a messy disk of gas, which naturally organizes itself into a powerful engine to shoot light across the cosmos.

Technical Summary

Technical Summary: Jets from Scratch: A 3D Dynamo Origin of Long Gamma-Ray Burst Jets

Problem Statement
The origin of the large-scale, coherent poloidal magnetic field required to power relativistic jets in collapsars (collapsing massive stars) remains uncertain. While the Blandford–Znajek (BZ) mechanism is the standard model for jet launching, it requires a dynamically important poloidal field threading the black hole (BH). Previous hypotheses suggested this field is inherited from the progenitor star or amplified in the proto-neutron star (PNS) phase. However, stellar evolution models indicate progenitors primarily generate toroidal fields via the Tayler-Spruit dynamo, with poloidal components being small-scale and prone to cancellation upon collapse. Furthermore, the "PNS inheritance" scenario faces uncertainties regarding flux retention during the transition from PNS to BH and the formation of a centrifugally supported disk. Consequently, it is unclear whether an in situ magnetohydrodynamic (MHD) dynamo within the collapsar accretion disk can efficiently convert the abundant toroidal fields into the necessary large-scale poloidal fields, particularly given the weak seed fields ($\sim 10^{12}$ G) and high mass accretion rates characteristic of collapsars.

Methodology
The authors present the first three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulations of collapsars initialized with physically motivated magnetic field profiles.

• Initial Conditions: The simulations utilize radial profiles of density, angular velocity, and magnetic fields derived from an ensemble average of 120 MESA stellar evolution models. These models incorporate the Tayler-Spruit dynamo for angular momentum transport. The initial magnetic configuration is a weak toroidal field ($B_\phi$) closely following pre-collapse stellar models, with negligible initial poloidal fields.
• Simulation Setup: The authors employ the GPU-accelerated code H-AMR to solve the ideal GRMHD equations. The computational domain spans from $1.16 GM/c^2$ (inside the event horizon) to $10^5 GM/c^2$.
• Parameter Space: The study varies the specific angular momentum ($j_{shell}$) to test different circularization radii ($r_{circ} = 18, 36, 72 GM/c^2$) and the initial toroidal field strength ($B_0$), defining "Strong," "Weak," and "Very Weak" magnetic field models. The simulations run for several seconds in physical time.
• Resolution: The grid employs static mesh refinement (SMR) with a base resolution of $224 \times 144 \times 128$ and two levels of refinement in the disk region ($4 \le r \le 500 GM/c^2$) to ensure the magnetorotational instability (MRI) is resolved.

Key Results

1. Dynamo Operation and Jet Launching: The simulations demonstrate that an MHD dynamo operates self-consistently in collapsar disks initialized with weak toroidal fields. As the toroidal field becomes dynamically important, it seeds the dynamo, generating coherent poloidal magnetic loops at radii of $\sim O(100) GM/c^2$. These loops are advected inward, threading the BH and launching relativistic jets.
2. Jet Power and Variability: The resulting jets sustain powers of $\gtrsim 10^{50}$ erg s$^{-1}$, comparable to observed long gamma-ray bursts (LGRBs). The jets are highly variable and exhibit "wobbling" (misalignment with the BH spin axis) due to the stochastic nature of the dynamo and perturbations from infalling gas.
3. MAD State and Flux Inversion: In successful models, the system transitions into a Magnetically Arrested Disk (MAD) state, characterized by a saturated dimensionless magnetic flux ($\Phi_{MAD} \gtrsim 10$). The authors identify stochastic magnetic-flux inversions driven by the dynamo, where loops of opposite polarity are advected inward and reconnect with existing fields. This leads to the formation of "striped jets" and transient jet shutdowns or sign reversals in the magnetic coherence parameter.
4. Dependence on Initial Conditions: The timing of jet onset depends on the initial magnetic field strength and the circularization radius.
   • Optimal Regime: The model with $r_{circ} = 36 GM/c^2$ and a strong seed field launches the most robust and earliest jet ($\sim 1.25$ s). This radius balances the need for a large enough disk to generate large-scale loops against the need for sufficient initial magnetization to amplify the field quickly.
   • Failure Cases: Models with very weak seed fields (e.g., $rc\ 36\ B\ vw$) fail to generate sufficient poloidal flux to launch a jet, remaining in a low-flux state ($\Phi_{MAD} < 5$).
5. Morphology: Unlike steady-state torus simulations, the collapsar jets are not strictly aligned with the spin axis. The continuous perturbation by infalling gas misaligns the disk and the threading magnetic field lines, resulting in jets that propagate out of the midplane.

Significance and Claims
The paper claims to provide the first demonstration of self-consistent, 3D collapsar jet formation driven entirely by an in situ accretion-disk dynamo, without relying on inherited magnetic flux from a PNS.

• Robustness: The results suggest that the accretion-disk dynamo provides a robust pathway for jet production across a broad range of progenitors, potentially making rotation alone a sufficient condition for LGRB production, provided an accretion disk forms.
• Resolution of Uncertainties: This work addresses the uncertainty regarding the origin of the poloidal field by showing that even weak, toroidal-dominated seed fields can be amplified to the levels required for jet launching ($\sim 10^{15}$ G effective fields) within seconds.
• Observational Implications: The stochastic nature of the dynamo and the resulting flux inversions offer a potential explanation for the variability and "striped" structure observed in LGRB light curves.
• Comparison to Prior Work: The authors note that their results differ from previous axisymmetric studies (e.g., Shibata et al. 2025), which required parameterized dynamo prescriptions and showed delayed or less coherent jet launching. The 3D nature of this simulation allows for the self-consistent growth of large-scale fields and captures the complex, wobbling morphology of the jets.

The authors acknowledge limitations, including the neglect of neutrino cooling and the uncertainty regarding the resolution dependence of the dynamo in very weak-field regimes, suggesting these as areas for future investigation.