GRMHD Simulations of Accreting Proto-Magnetars I. Implications for Gamma-Ray Burst Jets and Energetic Explosions

Abstract

Newly formed, rapidly rotating, strongly magnetized neutron stars ("millisecond proto-magnetars") are promising central engines for gamma-ray bursts (GRBs) and luminous supernovae. Although often modeled in isolation, they can be born surrounded by accretion disks in stellar collapse, neutron-star mergers, or accretion-induced collapse. We present axisymmetric GRMHD simulations of hyperaccretion onto such objects, including a physical equation of state and charged-current weak interactions. Holding the weakly magnetized accretion torus fixed, we vary the stellar dipole field strength to span crushed-magnetosphere, magnetically channeled accretion, and centrifugal-propeller regimes, and compare with an otherwise similar accreting black hole. Accretion compresses the stellar magnetosphere and opens additional magnetic flux, producing relativistic jet powers that exceed isolated-dipole spin-down estimates by factors of a few to ~10. Even while the magnetosphere remains compressed against the stellar surface, stronger fields increasingly impede accretion and enhance outflows. Channeled-accretion models show strong jet variability driven by plasmoid eruptions and intermittent magnetospheric accretion, whereas the propeller model produces a steadier, more powerful jet and rapid spin-down. The disk-magnetosphere interaction also regulates how efficiently the neutron star grows and whether it spins up or down; near spin equilibrium, inefficient accretion can delay collapse to a black hole relative to estimates based on the external mass-supply rate. Accreting proto-magnetars can therefore power relativistic jets and baryon-rich outflows with energetics comparable to those inferred for long GRBs and GRB-supernovae. A companion paper explores implications for neutron-rich ejecta and r-process nucleosynthesis.