Neutrino-cooled Accretion Disks Research Articles

Fast neutrino conversion in hydrodynamic simulations of neutrino-cooled accretion disks

The outflows from neutrino-cooled black-hole (BH) accretion disks formed in neutron-star mergers or cores of collapsing stars are expected to be neutron-rich enough to explain a large fraction of elements created by the rapid neutron-capture (r-) process, but their precise chemical composition remains elusive. Here, we investigate the role of fast neutrino flavor conversion, motivated by the findings of our post-processing analysis that shows evidence of electron-neutrino lepton-number (ELN) crossings deep inside the disk, hence suggesting possibly non-trivial effects due to neutrino flavor mixing. We implement a parametric, dynamically self-consistent treatment of fast conversion in time-dependent simulations and examine the impact on the disk and its outflows. By activating the, otherwise inefficient, emission of heavy-lepton neutrinos, fast conversions enhance the disk cooling rates and reduce the absorption rates of electron-type neutrinos, causing a reduction of the electron fraction in the disk by 0.03-0.06 and in the ejected material by 0.01-0.03. The r-process yields are enhanced by typically no more than a factor of two, rendering the overall impact of fast conversions modest. The kilonova is prolonged as a net result of increased lanthanide opacities and enhanced radioactive heating rates. We observe only mild sensitivity to the disk mass, the condition for the onset of flavor conversion, and to the considered cases of flavor mixing. Remarkably, parametric models of flavor mixing that conserve the lepton numbers per family result in an overall smaller impact than models invoking three-flavor equipartition, often assumed in previous works.

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.

Evolution of the magnetized, neutrino-cooled accretion disk in the aftermath of a black hole-neutron star binary merger

Black hole-torus systems from compact binary mergers are possible engines for gamma-ray bursts (GRBs). During the early evolution of the post-merger remnant, the state of the torus is determined by a combination of neutrino cooling and magnetically-driven heating processes, so realistic models must include both effects. In this paper, we study the post-merger evolution of a magnetized black hole-neutron star binary system using the Spectral Einstein Code (SpEC) from an initial post-merger state provided by previous numerical relativity simulations. We use a finite-temperature nuclear equation of state and incorporate neutrino effects in a leakage approximation. To achieve the needed accuracy, we introduce improvements to SpEC's implementation of general-relativistic magnetohydrodynamics (MHD), including the use of cubed-sphere multipatch grids and an improved method for dealing with supersonic accretion flows where primitive variable recovery is difficult. We find that a seed magnetic field triggers a sustained source of heating, but its thermal effects are largely cancelled by the accretion and spreading of the torus from MHD-related angular momentum transport. The neutrino luminosity peaks at the start of the simulation, and then drops significantly over the first 20\,ms but in roughly the same way for magnetized and nonmagnetized disks. The heating rate and disk's luminosity decrease much more slowly thereafter. These features of the evolution are insensitive to grid structure and resolution, formulation of the MHD equations, and seed field strength, although turbulent effects are not fully converged

Three-Dimensional General-Relativistic Magnetohydrodynamic Simulations of Remnant Accretion Disks from Neutron Star Mergers: Outflows and r-Process Nucleosynthesis.

The merger of binary neutron stars, or of a neutron star and a stellar-mass black hole, can result in the formation of a massive rotating torus around a spinning black hole. In addition to providing collimating media for γ-ray burst jets, unbound outflows from these disks are an important source of mass ejection and rapid neutron capture (r-process) nucleosynthesis. We present the first three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations of neutrino-cooled accretion disks in neutron star mergers, including a realistic equation of state valid at low densities and temperatures, self-consistent evolution of the electron fraction, and neutrino cooling through an approximate leakage scheme. After initial magnetic field amplification by magnetic winding, we witness the vigorous onset of turbulence driven by the magnetorotational instability (MRI). The disk quickly reaches a balance between heating from MRI-driven turbulence and neutrino cooling, which regulates the midplane electron fraction to a low equilibrium value Y_{e}≈0.1. Over the 380-ms duration of the simulation, we find that a fraction ≈20% of the initial torus mass is unbound in powerful outflows with asymptotic velocities v≈0.1c and electron fractions Y_{e}≈0.1-0.25. Postprocessing the outflows through a nuclear reaction network shows the production of a robust second- and third-peak r process. Though broadly consistent with the results of previous axisymmetric hydrodynamical simulations, extrapolation of our results to late times suggests that the total ejecta mass from GRMHD disks is significantly higher. Our results provide strong evidence that postmerger disk outflows are an important site for the r process.

NEUTRINO-COOLED ACCRETION MODEL WITH MAGNETIC COUPLING FOR X-RAY FLARES IN GAMMA-RAY BURSTS

The neutrino-cooled accretion disk, which was proposed to work as the central engine of gamma-ray bursts, encounters difficulty in interpreting the X-ray flares after the prompt gamma-ray emission. In this paper, the magnetic coupling (MC) between the inner disk and the central black hole (BH) is taken into consideration. For mass accretion rates around 0.001 ∼ 0.1 M☉ s−1, our results show that the luminosity of neutrino annihilation can be significantly enhanced due to the coupling effects. As a consequence, after the gamma-ray emission, a remnant disk with mass Mdisk ≲ 0.5 M☉ may power most of the observed X-ray flares with the rest frame duration less than 100 s. In addition, a comparison between the MC process and the Blandford–Znajek mechanism is shown on the extraction of BH rotational energy.

Electron-positron energy deposition rate from neutrino pair annihilation on the rotation axis of neutron and quark stars

We investigate the deposition of energy due to the annihilations of neutrinos and antineutrinos on the rotation axis of rotating neutron and quark stars, respectively. The source of the neutrinos is assumed to be a neutrino-cooled accretion disk around the compact object. Under the assumption of the separability of the neutrino null geodesic equation of motion we obtain the general relativistic expression of the energy deposition rate for arbitrary stationary and axisymmetric space-times. The neutrino trajectories are obtained by using a ray tracing algorithm, based on numerically solving the Hamilton-Jacobi equation for neutrinos by reversing the proper time evolution. We obtain the energy deposition rates for several classes of rotating neutron stars, described by different equations of state of the neutron matter, and for quark stars, described by the MIT bag model equation of state and in the CFL (Color-Flavor-Locked) phase, respectively. The electron-positron energy deposition rate on the rotation axis of rotating neutron and quark stars is studied for two accretion disk models (isothermal disk and accretion disk in thermodynamical equilibrium). Rotation and general relativistic effects modify the total annihilation rate of the neutrino-antineutrino pairs on the rotation axis of compact stellar, as measured by an observer at infinity. The differences in the equations of state for neutron and quark matter also have important effects on the spatial distribution of the energy deposition rate by neutrino-antineutrino annihilation.

Annihilation luminosity of a neutrino-cooled accretion disk in a gamma-ray burst

We discuss how the annihilation luminosity of a neutrino-cooled accretion disk in a gamma-ray burst, $$ L_{V\bar V} $$ , is determined by the disk’s fundamental parameters, namely, the mass of the central black hole M, the mass accretion rate Ṁ, and the viscosity parameter α. It is shown that $$ L_{V\bar V} $$ depends mainly on Ṁ in evidence, and decreases with increasing M, but is almost independent of α. This result argues additionally that the central black hole in a gamma-ray burst must be with a stellar mass.

Precession of neutrino-cooled accretion disks in gamma-ray burst engines

We study the precession of accretion disks in the context of gamma-ray burst inner engines. With an accretion disk model that allows for neutrino cooling, we evaluate the possible periods of disk precession and nutation due to the Lense-Thirring effect. Assuming jet ejection perpendicular to the disk midplane and a typical intrinsic time dependence for the burst, we find possible gamma-ray light curves with temporal microstructure similar to what is observed in some subsamples. We conclude that the precession and nutation of a neutrino-cooled accretion disk in the burst engine might be responsible for some events, especially those with a slow rise and a fast decay.