Radiation GRMHD Models of Accretion onto Stellar-Mass Black Holes: II. Super-Eddington Accretion

This paper presents comprehensive GRMHD simulations with angle-discretized radiation transport showing that super-Eddington accretion onto stellar-mass black holes consistently forms radiation-pressure-supported thick disks with low efficiency due to outflows, where the presence of strong jets driven by high spin and magnetic flux can clear the funnel to enable beamed radiation escape, while weak jets result in obscured flows with distinct observational signatures applicable to ultraluminous X-ray sources and other transient systems.

The Gist

Imagine a black hole not as a cosmic vacuum cleaner that sucks everything in, but as a cosmic whirlpool that is so full of water (matter) it's actually overflowing. This paper is a detailed study of what happens when a black hole tries to eat more food than it possibly can handle—a state astronomers call "Super-Eddington Accretion."

Here is the story of the black hole's messy, high-energy dinner, explained simply.

1. The Overstuffed Table (The Disk)

Normally, matter falling into a black hole forms a flat, thin disk, like a vinyl record spinning around a turntable. But in this study, the black hole is eating so fast that the "food" piles up.

• The Analogy: Imagine trying to pour a firehose into a teacup. The water doesn't just sit flat; it swells up into a giant, puffy, three-dimensional blob.
• The Result: The disk becomes geometrically thick. It's no longer a flat pancake; it's a puffy, radiation-filled cloud supported by the pressure of light itself. The light generated by the friction of the swirling gas is so intense that it pushes back against gravity, keeping the disk puffed up.

2. The Traffic Jam and the Escape Tunnel (The Funnel)

Because the black hole is eating so fast, the light (photons) generated in the center gets trapped. It's like a traffic jam where the cars (light particles) can't get out because the road is too crowded.

• The Analogy: Think of a crowded concert hall where everyone is trying to leave through the main doors. The crowd is so thick that no one can move. However, if you clear a path straight up to the ceiling, people can escape.
• The Result: The simulations show that the black hole naturally clears a conical tunnel (a funnel) straight up from its poles. This is the only place where light can escape. Because most light is trapped in the "traffic jam" of the disk, the black hole appears surprisingly dim (low efficiency) to an outside observer, even though it's generating massive amounts of energy.

3. The Two Types of Jets (The Sprinklers)

Black holes often shoot out powerful beams of energy called jets. This paper found that the strength of these jets depends on the black hole's "magnetic hairdo" (magnetic field) and how fast it spins.

• The Strong Jet (The Firehose):
  • Scenario: If the black hole spins fast and has a strong, organized magnetic field (like a single loop of wire), it acts like a powerful pump.
  • The Result: It blasts a super-fast jet that clears the funnel completely. This opens the "ceiling" of the traffic jam, allowing a huge burst of light to shoot out. This is like a firehose clearing a path through a crowd.
• The Weak Jet (The Garden Hose):
  • Scenario: If the magnetic field is messy (double loops) or the spin is slow, the jet is weak.
  • The Result: It can't clear the funnel. Instead, the "crowd" of gas and light stays trapped. The black hole is still eating, but it's choking on its own light. The energy escapes slowly as a slow, puffy wind rather than a focused beam.

4. The Invisible Drivers (Magnetic Stress)

How does the matter actually fall in? It needs to lose its "spin" (angular momentum) to drop toward the center.

• The Analogy: Imagine a spinning ice skater. To spin faster, they pull their arms in. To slow down and fall inward, they need to push their arms out against something.
• The Result: The paper found that magnetic fields act like invisible tethers. They grab the spinning gas and drag the "spin" outward, allowing the gas to fall in. It's not the gas rubbing against itself (friction) that does the work; it's the magnetic "ropes" pulling the energy away.

5. The Spiral Waves (The Plunging Region)

Right before the matter crosses the point of no return (the event horizon), the paper found something cool: Spiral waves.

• The Analogy: Think of a whirlpool in a bathtub. As the water spirals down the drain, you see spiral arms forming.
• The Result: The matter doesn't just fall straight down; it forms spiral density waves, like a spiral staircase of gas, as it plunges into the black hole. These waves help transport the remaining spin energy outward one last time.

Why Does This Matter?

This isn't just about math; it explains real things we see in the universe:

• ULXs (Ultraluminous X-ray Sources): These are incredibly bright objects. Our model suggests they are black holes with strong jets that have cleared the funnel, letting the light beam directly at us.
• Little Red Dots: These are faint, red objects. Our model suggests they are black holes with weak jets, where the light is trapped in a puffy, obscuring cloud, making them look dim and red.
• Tidal Disruption Events: When a star gets eaten by a black hole, this model helps us understand the flash of light and the wind of gas that follows.

The Bottom Line

The universe is messy. When a black hole eats too fast, it doesn't just suck everything in neatly. It puffs up, gets trapped in its own light, and shoots out jets that act like traffic cops, either clearing the road for a blinding beam of light or failing to clear it, leaving the black hole hidden in a fog of its own making. This paper gives us the blueprint for how that chaotic dinner party plays out.

Technical Summary

Here is a detailed technical summary of the paper "Radiation GRMHD Models of Accretion onto Stellar-Mass Black Holes: II. Super-Eddington Accretion" by Zhang et al.

1. Problem Statement

Super-Eddington accretion (accretion rates exceeding the Eddington limit, $\dot{M} > \dot{M}_{\text{Edd}}$) is a critical process powering various high-energy astrophysical phenomena, including Ultraluminous X-ray sources (ULXs), Little Red Dots (LRDs), Tidal Disruption Events (TDEs), and stellar-mass black hole binaries (BHBs). While theoretical frameworks like the "slim disk" model exist, they rely on simplifying assumptions (e.g., geometrically thin disks, constant viscosity $\alpha$) that may not hold in the presence of strong radiation pressure, massive outflows, and complex magnetic field topologies.

There is a need for a comprehensive, self-consistent numerical framework that solves the General Relativistic Magnetohydrodynamics (GRMHD) equations coupled with full radiation transport to understand:

• The global structure and stability of super-Eddington flows.
• The mechanisms driving angular momentum transport and energy redistribution.
• The formation and properties of relativistic jets and radiation-driven outflows.
• The observational signatures of these systems under varying magnetic and spin conditions.

2. Methodology

The authors utilize the ATHENAK code to perform a suite of 12 radiation GRMHD simulations.

• Physics Framework:

  • Equations: Solves the GRMHD equations coupled with angle-discretized radiation transport (frequency-integrated).
  • Coordinates: Cartesian Kerr-Schild coordinates.
  • Radiation Treatment: Uses a full transport method (not the M1 closure approximation) to accurately capture anisotropy and photon trapping.
  • Numerical Fixes: Implements a specialized density reduction technique in low-density funnel regions to prevent artificial Compton cooling caused by the numerical density floor interacting with high-density units.

• Parameter Space:

  • Black Hole Spin ($a$): Two values tested: $0.3$(low spin) and$0.9375$ (high spin).
  • Magnetic Topology: Single-loop (standard) and Double-loop configurations.
  • Accretion Rates: A range of super-Eddington rates (normalized to Eddington, $\dot{M}/\dot{M}_{\text{Edd}} \approx 9$ to $150$).
  • Resolutions: Intermediate-resolution (2048x2048x1152) and Low-resolution (1152x1152x640) meshes, plus a resolution study on radiation angular grid.
  • Comparisons: Includes non-radiative GRMHD runs and comparisons with the slim disk analytical model.

• Analysis Routine:

  • Time-averaging over steady-state periods ($\sim 10,000 r_g/c$).
  • Decomposition of angular momentum flux (Reynolds, Maxwell, Radiation diffusion).
  • Four-force analysis to dissect the balance between gravity, pressure, and radiation forces.
  • Energetic partitioning of heating/cooling mechanisms in disk, wind, and jet regions.

3. Key Contributions

• Unified Analysis Framework: Establishes a systematic routine for analyzing super-Eddington GRMHD simulations, applied consistently across different parameter regimes.
• Resolution Studies: Provides rigorous convergence tests for spatial and angular resolution, demonstrating that while small-scale turbulence is under-resolved, the primary drivers of accretion (Maxwell stress) are well-converged.
• Density Floor Mitigation: Introduces a novel method to suppress artificial Compton cooling in the low-density jet funnel, ensuring physical fidelity in regions where density approaches numerical limits.
• Jet-Flow Coupling: Explicitly links the strength of the magnetic field topology to the ability of jets to clear the polar funnel, thereby determining the system's radiative efficiency and observational appearance.

4. Key Results

A. Disk Structure and Dynamics

• Geometry: All super-Eddington models develop geometrically thick disks supported primarily by radiation pressure, regardless of magnetic topology.
• Photon Trapping: Radiation generated in the inner disk is largely trapped by the optically thick inflow. The thermal energy transport is dominated by radiation advection rather than diffusion.
• Angular Momentum Transport:
  • Maxwell Stress: The dominant mechanism for outward angular momentum transport (contributing ~60–90%).
  • Reynolds Stress: Subdominant (~10–40%).
  • Radiation Stress: Negligible contribution to angular momentum transport.
  • Viscosity: The effective viscosity ($\alpha$) is not constant; it follows a power-law dependence on radius, varying by a factor of two across the inflow equilibrium region.

B. Outflows and Jets

• Wind: Radiation pressure drives significant, mildly relativistic outflows from the disk surface. In high accretion rate models, these winds are optically thick.
• Strong Jets (Single-Loop):
  • Formed when sufficient net vertical magnetic flux accumulates (SANE regime).
  • Powered by the Blandford-Znajek process.
  • Effect: Successfully evacuate the polar funnel, creating a low-density channel. This allows radiation to escape via strong geometric beaming.
  • Power: Jet power can equal or exceed radiation luminosity.
• Weak Jets (Double-Loop):
  • Intermittent and weak; fail to clear the funnel.
  • Effect: The funnel remains filled with optically thick, radiation-driven outflows. High-energy photons are trapped, leading to very low radiation efficiency and obscured inner disks.
  • Power: Negligible compared to radiation.

C. Plunging Region

• Spiral structures (density waves) are observed in the plunging region near the ISCO.
• These structures exhibit a $90^\circ$ phase offset between density and velocity divergence, behaving like classic density waves that extract angular momentum via Maxwell stress.

D. Energetics and Efficiency

• Radiation Efficiency: Generally very low ($\eta < 1\%$) due to photon trapping and funnel obscuration. Efficiency decreases as accretion rate increases.
• Cooling Mechanisms:
  • Disk: Dominated by advection.
  • Wind: Kinetic energy and Poynting flux dominate cooling at large radii.
  • Jet: Poynting flux dominates near the horizon; kinetic energy dominates at large radii.

5. Significance and Observational Implications

• Breaking Degeneracies: The models suggest that systems with strong jets (cleared funnels) will appear as steady, highly variable sources with hard spectra, while weak-jet systems (obscured funnels) will show high variability, soft spectra, and potential X-ray under-detection (relevant for LRDs).
• ULXs and TDEs: The conical funnel structure cleared by strong jets provides a physical explanation for high X-ray polarization in ULXs (e.g., Cyg X-3) and the spectral pivoting observed in these systems.
• Slim Disk Model Validation: The numerical results agree well with the slim disk model only when outflows are negligible. When outflows are significant, the slim disk model's assumption of constant $\alpha$ and thin geometry fails to capture the disk's true profile.
• MAD vs. SANE: Unlike previous M1-based simulations that often reach the Magnetically Arrested Disk (MAD) regime with high efficiency, these SANE (Stable, Non-Arrested) models show that without sufficient initial vertical flux, super-Eddington flows remain inefficient and geometrically thick.

In conclusion, this paper provides a robust physical picture of super-Eddington accretion, highlighting the critical role of magnetic topology in determining whether a system produces a powerful, observable jet or remains an obscured, low-efficiency accretor.