Magnetically-Driven Neutron-Rich Ejecta Unleashed: Global 3D Neutrino-GRMHD Simulations of Collapsars Reveal the Conditions for r-process Nucleosynthesis

Collapsars - rapidly rotating stellar cores that form black holes (BHs) - can power gamma-ray bursts (GRBs) and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process ($r$-process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included magnetic fields or the BH, both of which are essential for launching mass outflows. We present $\nu$H-AMR, a novel neutrino-transport general relativistic magnetohydrodynamic ($\nu$GRMHD) code, which we use to perform the first 3D $\nu$GRMHD collapsar simulations. We find a self-consistent formation of a disk with initially weak magnetic flux, resulting in a low accretion speed and leaving sufficient time for the disk to neutronize. However, once substantial magnetic flux accumulates near the BH, it becomes dynamically important, leading to a magnetically arrested disk that unbinds some of the neutron-rich material. The strong flux also accelerates the accretion speed, preventing further disk neutronization. The neutron-rich disk ejecta collides with the infalling stellar gas, generating a shocked cocoon with an electron fraction, $Y_\text{e}\gtrsim0.2$. Continuous mixing between the cocoon and neutron-poor stellar gas incrementally raises the outflow $Y_\text{e}$, but the final $r$-process yield is determined earlier at the point of neutron capture freeze-out. Our models require extreme magnetic fluxes and mass accretion rates to eject neutron-rich material ($Y_\text{e}\lesssim0.3$), implying very high $r$-process ejecta masses $M_\text{ej}\lesssim{}M_\odot$. Future work will explore under what conditions more typical collapsar engines become $r$-process factories.