Three-dimensional General Relativistic Magnetohydrodynamic Simulations Research Articles

Multizone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales

Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. Elaborating on H. Cho et al., we describe and test a “multizone” technique, which is designed to tackle this difficult problem in three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. While short-timescale variability should be interpreted with caution, the method is demonstrated to be well-suited for finding dynamical steady states over a wide dynamic range. We simulate accretion on a nonspinning SMBH (a * = 0) using initial conditions and the external galactic potential from a large-scale galaxy simulation and achieve a steady state over eight decades in radius. As found in H. Cho et al., the density scales with radius as ρ ∝ r −1 inside the Bondi radius R B , which is located at R B = 2 × 105 r g (≈60 pc for M87), where r g is the gravitational radius of the SMBH; the plasma-β is ∼ unity, indicating an extended magnetically arrested state; the mass accretion rate Ṁ is ≈1% of the analytical Bondi accretion rate ṀB; and there is continuous energy feedback out to ≈100R B (or beyond > kpc) at a rate ≈0.02Ṁc2 . Surprisingly, no ordered rotation in the external medium survives as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multizone method, we simulate GRMHD accretion over a wide range of Bondi radii, R B ∼ 102−107 r g, and find that Ṁ/ṀB≈(RB/6rg)−0.5 .

Shock-induced Partial Alignment in Geometrically Thick Tilted Accretion Disks Around Black Holes

We carry out idealized three-dimensional general-relativistic magnetohydrodynamic simulations of prograde, weakly magnetized, and geometrically thick accretion flows where the gas distribution is misaligned from the black hole (BH) spin axis. We evolve the disk for three BH spins: a = 0.5, 0.75, and 0.9375, and we contrast them with a standard aligned disk simulation with a = 0.9375. The tilted disks achieve a warped and twisted steady-state structure, with the outer disk misaligning further away from the BH and surpassing the initial 24° misalignment. However, closer to the BH, there is evidence of partial alignment, as the inclination angle decreases with radius in this regime. Standing shocks also emerged in proximity to the BH, roughly at ∼6 gravitational radii. We show that these shocks act to partially align the inner disk with the BH spin. The rate of alignment increases with increasing BH spin magnitude, but in all cases is insufficient to fully align the gas before it accretes. Additionally, we present a toy model of orbit crowding that can predict the location of the shocks in moderate-to-fast rotating BHs, illustrating a potential physical origin for the behavior seen in simulations—with possible applications in determining the positions of shocks in real misaligned astrophysical systems.

Dynamics and emission properties of flux ropes from two-temperature GRMHD simulations with multiple magnetic loops

Context. Flux ropes erupting in the vicinity of a black hole are thought to be a potential model for the flares observed in Sagittarius A*. Aims. In this study, we examine the radiative properties of flux ropes that emerged in the vicinity of a black hole. Methods. We performed three-dimensional two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of magnetized accretion flows with alternating multiple magnetic loops and general relativistic radiation transfer (GRRT) calculations. In the GRMHD simulations, we implemented two different sizes of initial magnetic loops. Results. In the small loop case, magnetic dissipation leads to a weaker excitement of magneto-rotational instability inside the torus, which generates a lower accretion rate compared to the large loop case. However, the small loop case generates more flux ropes due to frequent reconnection by magnetic loops with different polarities. By calculating the thermal synchrotron emission, we found that the variability of light curves and the emitting region are tightly related. At 230 GHz and higher, the emission from the flux ropes is relatively stronger compared to the background, which is responsible for the filamentary structure in the images. At lower frequencies (e.g., 43 GHz), emission comes from more extended regions, which have a less filamentary structure in the image. Conclusions. Our study shows that self-consistent electron temperature models are essential for the calculation of thermal synchrotron radiation and the morphology of the GRRT images. Flux ropes contribute considerable emission at frequencies ≳230 GHz.

Magnetorotational dynamo can generate large-scale vertical magnetic fields in 3D GRMHD simulations of accreting black holes

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $t_{\rm gen}$, and one for the large-scale scale advection, $t_{\rm adv}$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $\alpha$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $\alpha$-dynamos.

Revisiting flares in Sagittarius A* based on general relativistic magnetohydrodynamic numerical simulations of black hole accretion

ABSTRACT High-resolution observations with GRAVITY-Very Large Telescope Interferometer (VLTI) instrument have provided abundant information about the flares in Sgr A*, the supermassive black hole in our Galactic centre, including the time-dependent location of the centroid (a ‘hotspot’), the light curve, and polarization. Yuan et al. (2009) proposed a ‘coronal mass ejection’ model to explain the flares and their association with the plasma ejection. The key idea is that magnetic reconnection in the accretion flow produces the flares and results in the formation and ejection of flux ropes. The dynamical process proposed in the model has been confirmed by three-dimensional general-relativistic magnetohydrodynamic simulations in a later work. Based on this scenario, in our previous works the radiation of the flux rope has been calculated analytically and compared to the observations. In the present paper, we develop the model by directly using numerical simulation data to interpret observations. We first identify flux ropes formed due to reconnection from the data. By assuming that electrons are accelerated in the reconnection current sheet and flow into the flux rope and emit their radiation there, we have calculated the time-dependent energy distribution of electrons after phenomenologically considering their injection due to reconnection acceleration, radiative and adiabatic cooling. The radiation of these electrons is calculated using the ray-tracing approach. The trajectory of the hotspot, the radiation light curve during the flare, and the polarization are calculated. These results are compared with the GRAVITY observations and good consistencies are found.

Enhanced Blandford Znajek jet in loop quantum black hole

The Blandford-Znajek (BZ) process powers energetic jets by extracting the rotating energy of a Kerr black hole. It is important to understand this process in non-Kerr black hole spacetimes. In this study, we conduct two-dimensional and three-dimensional two-temperature General Relativistic Magnetohydrodynamic (GRMHD) simulations of magnetized accretion flows onto a rotating Loop-Quantum black hole (LQBH). Our investigation focuses on the accretion flow structure and jet launching dynamics from our simulations.We observe that the loop quantum effects increase the black hole angular frequency for spinning black holes.This phenomenon intensifies the frame-dragging effect, leading to an amplification of the toroidal magnetic field within the funnel region and enhancement of the launching jet power.It is possible to fit the jet power following a similar fitting formula of the black hole angular frequency as seen in the Kerr black hole.Based on the General Relativistic Radiation Transfer (GRRT) calculation, we find that the jet image from LQBH has a wider opening angle and an extended structure than the Kerr BH.

The relationship between simulated sub-millimeter and near-infrared images of sagittarius A* from a magnetically arrested black hole accretion flow

Abstract Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, undergoes large-amplitude near-infrared (NIR) flares that can coincide with the continuous rotation of the NIR emission region. One promising explanation for this observed NIR behavior is a magnetic flux eruption, which occurs in three-dimensional General Relativistic Magneto-Hydrodynamic (3D GRMHD) simulations of magnetically arrested accretion flows. After running two-temperature 3D GRMHD simulations, where the electron temperature is evolved self-consistently along with the gas temperature, it is possible to calculate ray-traced images of the synchotron emission from thermal electrons in the accretion flow. Changes in the gas dominated (σ = b2/2ρ < 1) regions of the accretion flow during a magnetic flux eruption reproduce the NIR flaring and NIR emission region rotation of Sgr A* with durations consistent with observation. In this paper, we demonstrate that these models also predict that large (1.5x - 2x) size increases of the sub-millimeter (sub-mm) and millimeter (mm) emission region follow most NIR flares by 20 - 50 minutes. These size increases occur across a wide parameter space of black hole spin (a = 0.3, 0.5, −0.5, 0.9375) and initial tilt angle between the accretion flow and black hole spin axes θ0 (θ0 = 0○, 16○, 30○). We also calculate the sub-mm polarization angle rotation and the shift of the sub-mm spectral index from zero to -0.8 during a prominent NIR flare in our high spin (a = 0.9375) simulation. We show that, during a magnetic flux eruption, a large (∼10rg), magnetically dominated (σ > 1), low density, and high temperature “bubble” forms in the accretion flow. The drop in density inside the bubble and additional electron heating in accretion flow between 15 rg- 25 rg leads to a sub-mm size increase in corresponding images.

Future Prospects for Constraining Black Hole Spacetime: Horizon-scale Variability of Astrophysical Jets

The Event Horizon Telescope (EHT) Collaboration has recently published the first horizon-scale images of the supermassive black holes M87* and Sgr A* and provided some first information on the physical conditions in their vicinity. The comparison between the observations and the three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations has enabled the EHT to set initial constraints on the properties of these black hole spacetimes. However, accurately distinguishing the properties of the accretion flow from those of the spacetime, most notably, the black hole mass and spin, remains challenging because of the degeneracies the emitted radiation suffers when varying the properties of the plasma and those of the spacetime. The next-generation EHT (ngEHT) observations are expected to remove some of these degeneracies by exploring the complex interplay between the disk–jet dynamics, which represents one of the most promising tools for extracting information on the black hole spin. By using GRMHD simulations of magnetically arrested disks and general relativistic radiative transfer (GRRT) calculations of the emitted radiation, we have studied the properties of the jet and the accretion disk dynamics on spatial scales that are comparable with the horizon. In this way, we are able to highlight that the radial and azimuthal dynamics of the jet are well correlated with the black hole spin. Based on the resolution and image reconstruction capabilities of the ngEHT observations of M87*, we can assess the detectability and associated uncertainty of this correlation. Overall, our results serve to assess the prospects for constraining the black hole spin with future EHT observations.

Three-dimensional GRMHD Simulations of Neutron Star Jets

Neutron stars and black holes in X-ray binaries are observed to host strong collimated jets in the hard spectral state. Numerical simulations can act as a valuable tool in understanding the mechanisms behind jet formation and its properties. Although there have been significant efforts in understanding black hole jets from general relativistic magnetohydrodynamic (GRMHD) simulations in recent years, neutron star jets still remain poorly explored. We present the results from three-dimensional GRMHD simulations of accreting neutron stars with oblique magnetospheres for the very first time. The jets in our simulations are produced due to the anchored magnetic field of the rotating star in analogy with the Blandford–Znajek process. We find that for accreting stars, the star–disk magnetic field interaction plays a significant role, and as a result, the jet power becomes directly proportional to , where Φjet is the open flux in the jet. The jet power decreases with increasing stellar magnetic inclination, and finally, for an orthogonal magnetosphere, it reduces by a factor of ≃2.95 compared to the aligned case. We also find that in the strong propeller regime, with a highly oblique magnetosphere, the disk-induced collimation of the open stellar flux preserves parts of the striped wind, resulting in a striped jet.

Jets from Neutron-Star Merger Remnants and Massive Blue Kilonovae.

We perform three-dimensional general-relativistic magnetohydrodynamic simulations with weak interactions of binary neutron-star (BNS) mergers resulting in a long-lived remnant neutron star, with properties typical of galactic BNS and consistent with those inferred for the first observed BNS merger GW170817. We demonstrate self-consistently that within ≲30 ms postmerger magnetized (σ∼5-10) incipient jets emerge with asymptotic Lorentz factor Γ∼5-10, which successfully break out from the merger debris within ≲20 ms. A fast (v≲0.6c), magnetized (σ∼0.1) wind surrounds the jet core and generates a UV/blue kilonova precursor on timescales of hours, similar to the precursor signal due to free neutron decay in fast dynamical ejecta. Postmerger ejecta are quickly dominated by magnetohydrodynamically driven outflows from an accretion disk. We demonstrate that, within only 50ms postmerger, ≳2×10^{-2}M_{⊙} of lanthanide-free, quasispherical ejecta with velocities ∼0.1-0.2c is launched, yielding a kilonova signal consistent with GW170817 on timescales of ≲5 d.

Three-dimensional General-relativistic Simulations of Neutrino-driven Winds from Magnetized Proto–Neutron Stars

Formed in the aftermath of a core-collapse supernova or neutron star merger, a hot proto–neutron star (PNS) launches an outflow driven by neutrino heating lasting for up to tens of seconds. Though such winds are considered potential sites for the nucleosynthesis of heavy elements via the rapid neutron capture process (r-process), previous work has shown that unmagnetized PNS winds fail to achieve the necessary combination of high entropy and/or short dynamical timescale in the seed nucleus formation region. We present three-dimensional general-relativistic magnetohydrodynamical simulations of PNS winds which include the effects of a dynamically strong (B ≳ 1015 G) dipole magnetic field. After initializing the magnetic field, the wind quickly develops a helmet-streamer configuration, characterized by outflows along open polar magnetic field lines and a “closed” zone of trapped plasma at lower latitudes. Neutrino heating within the closed zone causes the thermal pressure of the trapped material to rise in time compared to the polar outflow regions, ultimately leading to the expulsion of this matter from the closed zone on a timescale of ∼60 ms, consistent with the predictions of Thompson. The high entropies of these transient ejecta are still growing at the end of our simulations and are sufficient to enable a successful second-peak r-process in at least a modest ≳1% of the equatorial wind ejecta.

R-process nucleosynthesis and kilonovae from hypermassive neutron star post-merger remnants

ABSTRACT We investigate r-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network SkyNet. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction Ye, peaking between ∼0.25–0.4 depending on the neutrino luminosity employed. We find that the resulting r-process abundance patterns differ from solar, with no significant production of material beyond the second r-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code SNEC and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the r-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.

Dynamical Unification of Tidal Disruption Events

The ∼100 tidal disruption events (TDEs) observed so far exhibit a wide range of emission properties both at peak and over their lifetimes. Some TDEs radiate predominantly at X-ray energies, while others radiate chiefly at UV and optical wavelengths. While the peak luminosities across TDEs show distinct properties, the evolutionary behavior can also vary between TDEs with similar peak emission properties. In particular, for optical TDEs, while their UV and optical emissions decline somewhat following the fallback pattern, some events can greatly rebrighten in X-rays at late time. In this Letter, we conduct three-dimensional general relativistic radiation magnetohydrodynamics simulations of TDE accretion disks at varying accretion rates in the regime of super-Eddington accretion. We make use of Monte Carlo radiative transfer simulations to calculate the reprocessed spectra at various inclinations and at different evolutionary stages. We confirm the unified model proposed by Dai et al., which predicts that the observed emission largely depends on the viewing angle of the observer with respect to the disk orientation. Furthermore, we find that disks with higher accretion rates have elevated wind and disk densities, which increases the reprocessing of the high-energy radiation and thus generally augments the optical-to-X-ray flux ratio along a particular viewing angle. This implies that at later times, as the accretion level declines, we expect that more X-rays will leak out along intermediate viewing angles. Such dynamical model for TDEs can provide a natural explanation for the diversity in the emission properties observed in TDEs at peak and along their temporal evolution.

Modeling the Gamma-Ray Burst Jet Properties with 3D General Relativistic Simulations of Magnetically Arrested Accretion Flows

We investigate the dependence of the gamma-ray burst (GRB) jet structure and its evolution on the properties of the accreting torus in the central engine. Our models numerically evolve the accretion disk around a Kerr black hole using three-dimensional general relativistic magnetohydrodynamic simulations. We use two different analytical hydrodynamical models of the accretion disk, based on the Fishbone–Moncrief and Chakrabarti solutions, as our initial states for the structure of the collapsar disk and the remnant after a binary neutron star (BNS) merger, respectively. We impose poloidal magnetic fields of two different geometries upon the initial stable solutions. We study the formation and evolution of the magnetically arrested disk state and its effect on the properties of the emitted jet. The jets produced in our models are structured and have a relatively hollow core and reach higher Lorentz factors at an angle ≳9° from the axis. The jet in our short GRB model has an opening angle of up to ∼25° while our long GRB engine produces a narrower jet, of up to ∼11°. We also study the time variability of the jets and provide an estimate of the minimum variability timescale in our models. The application of our models to the GRB jets in the BNS postmerger system and to the ultrarelativistic jets launched from collapsing stars are briefly discussed.

Linear and circular polarization of a 1D relativistic jet model

Context. Polarimetric observations of black holes allow us to probe structures of magnetic fields and plasmas in strong gravity. Aims. We present a study of the polarimetric properties of a synchrotron spectrum emitted from a relativistic jet using a low-dimensional model. Methods. A novel numerical scheme is used to integrate relativistic polarized radiative transfer equations in a slab geometry where the plasma conditions change along the integration path. Results. We find that the simple model of a non-uniform jet can recover basic observational characteristics of some astrophysical sources with a relativistic jet, such as extremely high rotation measures. Our models incorporate a time-dependent component. A small fluctuation in density or temperature of the plasma along the jet produces significant amounts of fluctuations not only in the fractional linear and circular polarizations, but also in the jet internal rotation measures. Conclusions. The low-dimensional models presented here are developed within the same computational framework as the complex three-dimensional general relativistic magnetohydrodynamics simulations of black hole disks and jets, and they offer guidance when interpreting the results from more complex polarization models. The models presented here are scalable to stationary and transient polarized radio emissions produced by relativistic plasma ejected from around compact objects, in both stellar-mass and supermassive black hole systems.

Predicting the X-Ray Spectra of Stellar-mass Black Holes from Simulations

Abstract We describe results from a new technique for the prediction of complete, self-consistent X-ray spectra from three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations of black hole accretion flows. Density and cooling rate data from a harm3d GRMHD simulation are post-processed by an improved version of the Monte Carlo radiation transport code pandurata (in the corona) and the Feautrier solver ptransx (in the disk), with xstar subroutines. The codes are run in a sequential, iterative fashion to achieve globally energy-conserving and self-consistent radiation fields, temperature maps, and photoionization equilibria. The output is the X-ray spectrum as seen by a distant observer, including features, such as the Fe Kα emission line and corresponding K-edge absorption trough, due to disk reprocessing of coronal power. For the example cases we consider—a non-rotating 10 M ⊙ black hole with solar abundances, accreting at 0.01, 0.03, 0.1, or 0.3 Eddington—we find spectra resembling actual observations of stellar-mass black holes in the soft or steep power-law state: broad thermal peaks (at 1–3 keV), steep power laws extending to high energy (Γ = 2.7–4.5), and prominent, asymmetric Fe Kα emission lines with equivalent widths in the range 40–400 eV (larger EW at lower accretion rates). By starting with simulation data, we obviate the need for parameterized descriptions of the accretion flow geometry—no a priori specification of the corona’s shape or flux, or the disk temperature or density, etc., is needed. Instead, we apply the relevant physical principles to simulation output using appropriate numerical techniques; this procedure allows us to calculate inclination-dependent spectra after choosing only a small number of physically meaningful parameters: black hole mass and spin, accretion rate, and elemental abundances.

Observing supermassive black holes in virtual reality

We present a 360∘ (i.e., 4π steradian) general-relativistic ray-tracing and radiative transfer calculations of accreting supermassive black holes. We perform state-of-the-art three-dimensional general-relativistic magnetohydrodynamical simulations using the BHAC code, subsequently post-processing this data with the radiative transfer code RAPTOR. All relativistic and general-relativistic effects, such as Doppler boosting and gravitational redshift, as well as geometrical effects due to the local gravitational field and the observer’s changing position and state of motion, are therefore calculated self-consistently. Synthetic images at four astronomically-relevant observing frequencies are generated from the perspective of an observer with a full 360∘ view inside the accretion flow, who is advected with the flow as it evolves. As an example we calculated images based on recent best-fit models of observations of Sagittarius A*. These images are combined to generate a complete 360∘ Virtual Reality movie of the surrounding environment of the black hole and its event horizon. Our approach also enables the calculation of the local luminosity received at a given fluid element in the accretion flow, providing important applications in, e.g., radiation feedback calculations onto black hole accretion flows. In addition to scientific applications, the 360∘ Virtual Reality movies we present also represent a new medium through which to interactively communicate black hole physics to a wider audience, serving as a powerful educational tool.

R-process Nucleosynthesis from Three-dimensional Magnetorotational Core-collapse Supernovae

Abstract We investigate r-process nucleosynthesis in 3D general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. The simulations include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission and absorption. We track the composition of the ejected material using the nuclear reaction network SkyNet. Our results show that the 3D dynamics of magnetorotational core-collapse supernovae (CCSN) are important for their nucleosynthetic signature. We find that production of r-process material beyond the second peak is reduced by a factor of 100 when the magnetorotational jets produced by the rapidly rotating core undergo a kink instability. Our results indicate that 3D magnetorotationally powered CCSNe are robust r-process sources only if they are obtained by the collapse of cores with unrealistically large precollapse magnetic fields of the order of 1013 G. Additionally, a comparison simulation that we restrict to axisymmetry results in overly optimistic r-process production for lower magnetic field strengths.

A Unified Model for Tidal Disruption Events

Abstract In the past few years wide-field optical and UV transient surveys and X-ray telescopes have allowed us to identify a few dozen candidate tidal disruption events (TDEs). While in theory the physical processes in TDEs are ubiquitous, a few distinct classes of TDEs have been observed. Some TDEs radiate mainly in NUV/optical, while others produce prominent X-rays. Moreover, relativistic jets have been observed in only a handful of TDEs. This diversity might be related to the details of the super-Eddington accretion and emission physics relevant to TDE disks. In this Letter, we utilize novel three-dimensional general relativistic radiation magnetohydrodynamics simulations to study the super-Eddington compact disk phase expected in TDEs. Consistent with previous studies, geometrically thick disks, wide-angle optically thick fast outflows, and relativistic jets are produced. The outflow density and velocity depend sensitively on the inclination angle, and hence so does the reprocessing of emission produced from the inner disk. We then use Monte Carlo radiative transfer to calculate the reprocessed spectra and find that that the observed ratio of optical to X-ray fluxes increases with increasing inclination angle. This naturally leads to a unified model for different classes of TDEs in which the spectral properties of the TDE depend mainly on the viewing angle of the observer with respect to the orientation of the disk.

Three-dimensional GRMHD Simulations of Neutrino-cooled Accretion Disks from Neutron Star Mergers

Abstract Merging binaries consisting of two neutron stars (NSs) or an NS and a stellar-mass black hole typically form a massive accretion torus around the remnant black hole or long-lived NS. Outflows from these neutrino-cooled accretion disks represent an important site for r-process nucleosynthesis and the generation of kilonovae. We present the first three-dimensional, general-relativistic magnetohydrodynamic (GRMHD) simulations including weak interactions and a realistic equation of state of such accretion disks over viscous timescales (380 ms). We witness the emergence of steady-state MHD turbulence, a magnetic dynamo with an ∼20 ms cycle, and the generation of a “hot” disk corona that launches powerful thermal outflows aided by the energy released as free nucleons recombine into α-particles. We identify a self-regulation mechanism that keeps the midplane electron fraction low (Y e ∼ 0.1) over viscous timescales. This neutron-rich reservoir, in turn, feeds outflows that retain a sufficiently low value of Y e ≈ 0.2 to robustly synthesize third-peak r-process elements. The quasi-spherical outflows are projected to unbind 40% of the initial disk mass with typical asymptotic escape velocities of 0.1c and may thus represent the dominant mass ejection mechanism in NS–NS mergers. Including neutrino absorption, our findings agree with previous hydrodynamical α-disk simulations that the entire range of r-process nuclei from the first to the third r-process peak can be synthesized in the outflows, in good agreement with observed solar system abundances. The asymptotic escape velocities and quantity of ejecta, when extrapolated to moderately higher disk masses, are consistent with those needed to explain the red kilonova emission following the NS merger GW170817.