Simulating the Magnetorotational Collapse of Supermassive Stars: Incorporating Gas Pressure Perturbations and Different Rotation Profiles

TL;DR

This study uses GRMHD simulations to explore how gas pressure perturbations and different rotation profiles influence the collapse of supermassive stars, leading to black hole formation and jet launching, relevant for understanding SMBH seeds and gamma-ray bursts.

Contribution

It extends previous models by incorporating gas pressure effects and varied rotation profiles in supermassive star collapse simulations, revealing their impact on black hole mass and jet properties.

Findings

Black hole mass depends sharply on gas pressure and rotation profile.

Jets are launched within 1-2×10^3 seconds after black hole formation.

Jet luminosities are consistent with the Blandford-Znajek mechanism and a universal accretion model.

Abstract

Collapsing supermassive stars (SMSs) with masses $M \gtrsim 10^{4-6}M_\odot$ have long been speculated to be the seeds that can grow and become supermassive black holes (SMBHs). We previously performed GRMHD simulations of marginally stable magnetized $\Gamma = 4/3$ polytropes uniformly rotating at the mass-shedding limit to model the direct collapse of SMSs. These configurations are supported entirely by thermal radiation pressure and model SMSs with $M \gtrsim 10^{6}M_\odot$. We found that around $90\%$ of the initial stellar mass forms a spinning black hole (BH) surrounded by a massive, hot, magnetized torus, which eventually launches an incipient jet. Here we perform GRMHD simulations of $\Gamma \gtrsim 4/3$, polytropes to account for the perturbative role of gas pressure in SMSs with $M \lesssim 10^{6}M_\odot$. We also consider different initial stellar rotation profiles. The stars are initially seeded with a dynamically weak dipole magnetic field that is either confined to the stellar interior or extended from its interior into the stellar exterior. We find that the mass of the BH remnant is $90\%-99\%$ of the initial stellar mass, depending sharply on $\Gamma -4/3$ as well as on the initial stellar rotation profile. After $t\sim 250-550M\approx 1-2\times 10^3(M/10^6M_\odot)$s following the BH formation, a jet is launched and it lasts for $\sim 10^4-10^5(M/10^6M_\odot)$s, consistent with the duration of long gamma-ray bursts. Our results suggest that the Blandford-Znajek mechanism powers the jet. They are also in agreement with our proposed universal model that estimates accretion rates and luminosities that characterize magnetized BH-disk remnant systems that launch a jet. This model helps explain why the outgoing luminosities for vastly different BH-disk formation scenarios all reside within a narrow range ($\sim 10^{52 \pm 1} \rm erg/s$), roughly independent of $M$.