Magnetorotational Instability Research Articles

Parametric Survey of Nonaxisymmetric Accretion Disk Instabilities: Magnetorotational Instability to Super-Alfvénic Rotational Instability

Accretion disks are highly unstable to magnetic instabilities driven by shear flow, where classically, the axisymmetric, weak-field magnetorotational instability (MRI) has received much attention through local WKB approximations. In contrast, discrete nonaxisymmetric counterparts require a more involved analysis through a full global approach to deal with the influence of the nearby magnetohydrodynamic (MHD) continua. Recently, rigorous MHD spectroscopy identified a new type of ultralocalized, nonaxisymmetric instability in global disks with super-Alfvénic flow. These super-Alfvénic rotational instabilities (SARIs) fill vast unstable regions in the complex eigenfrequency plane with (near eigen)modes that corotate at the local Doppler velocity and are radially localized between Alfvénic resonances. Unlike discrete modes, they are utterly insensitive to the radial disk boundaries. In this work, we independently confirm the existence of these unprecedented modes using our novel spectral MHD code Legolas, reproducing and extending our earlier study with detailed eigenspectra and eigenfunctions. We calculate the growth rates of SARIs and MRI in a variety of disk equilibria, highlighting the impact of field strength and orientation, and find correspondence with analytical predictions for thin, weakly magnetized disks. We show that nonaxisymmetric modes can significantly extend instability regimes at high mode numbers, with maximal growth rates comparable to those of the MRI. Furthermore, we explicitly show a region filled with quasi-modes whose eigenfunctions are extremely localized in all directions. These modes must be ubiquitous in accretion disks, and play a role in local shearing box simulations. Finally, we revisit recent dispersion relations in the appendix, highlighting their relation to our global framework.

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

Flux ropes erupting in the vicinity of a black hole are thought to be a potential model for the flares observed in Sagittarius\,A$^*$. In this study, we examine the radiative properties of flux ropes that emerged in the vicinity of a black hole. 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. 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. 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 $ GHz$.

The solar dynamo begins near the surface

The magnetic dynamo cycle of the Sun features a distinct pattern: a propagating region of sunspot emergence appears around 30° latitude and vanishes near the equator every 11 years (ref. 1). Moreover, longitudinal flows called torsional oscillations closely shadow sunspot migration, undoubtedly sharing a common cause2. Contrary to theories suggesting deep origins of these phenomena, helioseismology pinpoints low-latitude torsional oscillations to the outer 5–10% of the Sun, the near-surface shear layer3,4. Within this zone, inwardly increasing differential rotation coupled with a poloidal magnetic field strongly implicates the magneto-rotational instability5,6, prominent in accretion-disk theory and observed in laboratory experiments7. Together, these two facts prompt the general question: whether the solar dynamo is possibly a near-surface instability. Here we report strong affirmative evidence in stark contrast to traditional models8 focusing on the deeper tachocline. Simple analytic estimates show that the near-surface magneto-rotational instability better explains the spatiotemporal scales of the torsional oscillations and inferred subsurface magnetic field amplitudes9. State-of-the-art numerical simulations corroborate these estimates and reproduce hemispherical magnetic current helicity laws10. The dynamo resulting from a well-understood near-surface phenomenon improves prospects for accurate predictions of full magnetic cycles and space weather, affecting the electromagnetic infrastructure of Earth.

Magnetic amplification in premerger neutron stars through resonance-induced magnetorotational instabilities

Magnetic amplification in premerger neutron stars through resonance-induced magnetorotational instabilities

Large Fluctuations within 1 au in Protoplanetary Disks

Protoplanetary disks are often assumed to change slowly and smoothly during planet formation. Here, we investigate the time evolution of isolated disks subject to viscosity and a disk wind. The viscosity is assumed to increase rapidly at around 900 K due to thermal ionization of alkali metals, or thermionic and ion emission from dust, and the onset of magnetorotational instability (MRI). The disks generally undergo large, rapid fluctuations for a wide range of time-averaged mass accretion rates. Fluctuations involve coupled waves in temperature and surface density that move radially in either direction through the inner 1.5 au of the disk. Two types of waves are seen with radial speeds of roughly 50 and 1000 cm s−1, respectively. The pattern of waves repeats with a period of roughly 10,000 yr that depends weakly on the average mass accretion rate. Viscous transport due to MRI is confined to the inner disk. This region is resupplied by mass flux from the outer disk driven by the disk wind. Interior to 1 au, the temperature and surface density can vary by a factor of 2–10 on timescales of years to kiloyears. The stellar mass accretion rate varies by 3 orders of magnitude on a similar timescale. This behavior lasts for at least 1 Myr for initial disks comparable to the minimum-mass solar nebula.

The High Mass Accretion in the Innermost Regions of a Viscously Evolved Protoplanetary Disk

In this paper, we investigate the mass accretion properties in the innermost regions of a viscously evolved protoplanetary disk and try to find some clues to the outburst events. In our newly developed one-dimensional time-dependent disk model based on the diffusion equation for surface density, we take into account the following physical effects: the gravitational collapse of the parent molecular cloud core, the irradiation from the central star to the disk, the effect of the photoevaporation mechanism, the viscosity due to the magnetorotational instability (MRI) and the gravitational instability (GI), and the thermal ionization mechanism in the inner regions. We find that the mass accretion rate M·disk in the innermost regions is statistically high enough to generate outbursts, although there are regions where the accretion rate is low. Additionally, we find that there is a weak correlation between the high mass accretion rate M·disk and the molecular cloud core’s properties (angular velocity ω and mass Mcd), as well as a strong correlation with the minimum viscosity parameter αmin. In general, there are two regions of outburst, the inner Region I and outer Region II. The outburst of Region I is caused by the MRI mechanism and thermal instability, while neither the MRI, the GI, nor the thermal instability causes the outburst of Region II. Our analysis suggests that the outer Region II is dominated by, or largely related to, the Rosseland mean opacity κR and the αmin parameter.

Viscous torque in turbulent magnetized active galactic nucleus accretion disks and its effects on the gravitational waves of extreme mass ratio inspirals

The merger of supermassive black holes produces millihertz gravitational waves (GWs), which are potentially detectable by the future Laser Interferometer Space Antenna (LISA). Such binary systems are usually embedded in an accretion disk environment at the center of an active galactic nucleus (AGN). Recent studies suggest the plasma environment imposes measurable imprints on the GW signal if the mass ratio of the binary is around q $ $. The effect of the gaseous environment on the GW signal is strongly dependent on the disk’s parameters; therefore, it is believed that future low-frequency GW detections will provide us with precious information about the physics of AGN accretion disks. We investigated this effect by measuring the viscous torque via modeling of the evolution of magnetized tori around the primary massive black hole. Using the general relativistic magnetohydrodynamic HARM-COOL code, we performed 2D and 3D simulations of weakly magnetized, thin accretion disks, with a possible truncation and transition to advection-dominated accretion flow. We studied the angular momentum transport and turbulence generated by the magnetorotational instability. We quantified the disk’s effective alpha viscosity and its evolution over time. We applied our numerical results to quantify the relativistic viscous torque on a hypothetical low-mass secondary black hole via a 1D analytical approach, and we estimated the GW phase shift due to the gas environment.

Dynamics of cold circumstellar gas in debris disks

Abstract Mounting observational evidence indicates that cold circumstellar gas is present in debris disk systems. This work focuses on various dynamical processes that debris-disk gas may undergo. We review five mechanisms that can transport angular momentum and their applications to debris disks. These include molecular viscosity, hydrodynamic turbulence, magnetohydrodynamic turbulence, magnetized disk winds, and laminar magnetic stress. We find that molecular viscosity can result in α as high as ≲ 0.1 for sufficiently low densities, while the Rossby wave instability is a possible source of hydrodynamic turbulence and structure formation. We argue that the vertical shear instability is unlikely due to the long cooling times. The onset of the magnetorotational instability (MRI) is dichotomous: for low density disks the MRI can be excited at the midplane, while for high mass disks it may only be operating at z > 2 − 3H, if at all. The MHD wind and laminar magnetic stress mechanisms rely on the configuration and strength of any background large-scale magnetic field, the existence of which is uncertain and possibly unlikely. We conclude that the dominant mechanism and its efficiency in transporting angular momentum varies from one system to the other, depending especially closely on the gas density. More detailed analyses shall be performed in the future focusing on representative, nearby debris disks.

Shedding light on the MRI-driven dynamo in a stratified shearing box

ABSTRACT We study the magneto-rotational instability (MRI) dynamo in a geometrically thin disc (H/R ≪ 1) using stratified zero net (vertical) flux shearing box simulations. We find that mean fields and electromotive forces (EMFs) oscillate with a primary frequency fdyn = 0.017 Ω (approximately nine orbital period), but also have higher harmonics at 3fdyn. Correspondingly, the current helicity has two frequencies 2fdyn and 4fdyn, which appear to be the beat frequencies of mean fields and EMFs, respectively, as expected from the magnetic helicity density evolution equation. Further, we adopt a novel inversion algorithm called the ‘Iterative Removal Of Sources’, to extract the turbulent dynamo coefficients in the mean-field closure using the mean magnetic fields and EMFs obtained from the shearing box simulation. We show that an α-effect (αyy) is predominantly responsible for the creation of the poloidal field from the toroidal field, while shear generates back a toroidal field from the poloidal field, indicating that an α–Ω-type dynamo is operative in MRI-driven accretion discs. We also find that both strong outflow ($\bar{v}_z$) and turbulent pumping (γz) transport mean fields away from the mid-plane. Instead of turbulent diffusivity, they are the principal sink terms in the mean magnetic energy evolution equation. We find encouraging evidence that a generative helicity flux is responsible for the effective α-effect. Finally, we point out potential limitations of horizontal (x − y) averaging in defining the ‘mean’ on the extraction of dynamo coefficients and their physical interpretations.

Particle-in-cell Simulations of the Magnetorotational Instability in Stratified Shearing Boxes

Abstract The magnetorotational instability (MRI) plays a crucial role in regulating the accretion efficiency in astrophysical accretion disks. In low-luminosity disks around black holes, such as Sgr A* and M87, Coulomb collisions are infrequent, making the MRI physics effectively collisionless. The collisionless MRI gives rise to kinetic plasma effects that can potentially affect its dynamic and thermodynamic properties. We present 2D and 3D particle-in-cell (PIC) plasma simulations of the collisionless MRI in stratified disks using shearing boxes with net vertical field. We use pair plasmas, with initial β = 100 and concentrate on sub-relativistic plasma temperatures (kBT ≲ mc2). Our 2D and 3D runs show disk expansion, particle and magnetic field outflows, and a dynamo-like process. They also produce magnetic pressure dominated disks with (Maxwell stress dominated) viscosity parameter α ∼ 0.5 − 1. By the end of the simulations, the dynamo-like magnetic field tends to dominate the magnetic energy and the viscosity in the disks. Our 2D and 3D runs produce fairly similar results, and are also consistent with previous 3D MHD simulations. Our simulations also show nonthermal particle acceleration, approximately characterized by power-law tails with temperature dependent spectral indices −p. For temperatures kBT ∼ 0.05 − 0.3 mc2, we find p ≈ 2.2 − 1.9. The maximum accelerated particle energy depends on the scale separation between MHD and Larmor-scale plasma phenomena in a way consistent with previous PIC results of magnetic reconnection-driven acceleration. Our study constitutes a first step towards modeling from first principles potentially observable stratified MRI effects in low-luminosity accretion disks around black holes.

MRI turbulence in vertically stratified accretion discs at large magnetic Prandtl numbers

ABSTRACT The discovery of the first binary neutron star merger, GW170817, has spawned a plethora of global numerical relativity simulations. These simulations are often ideal (with dissipation determined by the grid) and/or axisymmetric (invoking ad hoc mean-field dynamos). However, binary neutron star mergers (similar to X-ray binaries and active galactic nuclei inner discs) are characterized by large magnetic Prandtl numbers, $\rm Pm$, (the ratio of viscosity to resistivity). $\rm Pm$ is a key parameter determining dynamo action and dissipation but it is ill-defined (and likely of order unity) in ideal simulations. To bridge this gap, we investigate the magnetorotational instability (MRI) and associated dynamo at large magnetic Prandtl numbers using fully compressible, three-dimensional, vertically stratified, isothermal simulations of a local patch of a disc. We find that, within the bulk of the disc (z ≲ 2H, where H is the scale-height), the turbulent intensity (parametrized by the stress-to-thermal-pressure ratio α), and the saturated magnetic field energy density, Emag, produced by the MRI dynamo, both scale as a power with Pm at moderate Pm (4 ≲ Pm ≲ 32): Emag ∼ Pm0.74 and α ∼ Pm0.71, respectively. At larger Pm (≳ 32), we find deviations from power-law scaling and the onset of a plateau. Compared to our recent unstratified study, this scaling with Pm becomes weaker further away from the disc mid-plane, where the Parker instability dominates. We perform a thorough spectral analysis to understand the underlying dynamics of small-scale MRI-driven turbulence in the mid-plane and of large-scale Parker-unstable structures in the atmosphere.

The Impact of Ambipolar Diffusion on the Rossby Wave Instability in a Protoplanetary Disk

Recent observational and simulation studies have revealed that ambipolar diffusion is an important phenomenon in the outer regions of a protoplanetary disk (PPD). However, numerous simulation studies have found that ambipolar diffusion suppresses the turbulence caused by the magnetorotational instability (MRI) in these regions of a PPD. The aim of this study is to investigate the impact of ambipolar diffusion on the Rossby wave instability (RWI) at large radii of a PPD. To accomplish this, we examine the occurrence of the RWI in a PPD threaded by the magnetic field in the presence of ambipolar diffusion. Additionally, we scale the ambipolar diffusivity with respect to both the toroidal magnetic field and an important parameter known as the Elsässer number. We obtain the growth rate of unstable RWI modes in the outer regions of a PPD using linear perturbation analysis. In our nonaxisymmetric perturbation analysis, we find that the amplitude of the toroidal field oscillates in various modes for small values of the Elsässer numbers. For small Elsässer numbers, the growth rate of unstable modes associated with the RWI decreases. In other words, ambipolar diffusion suppresses the perturbation caused by the RWI. This effect is similar to the effect of ambipolar diffusion on the MRI. In contrast to the MRI, where there is a range around 1 for the Elsässer numbers, the existence of ambipolar diffusion supports the occurrence of RWI. Finally, we compare our findings with those of simulation studies to emphasize the importance of RWI in the outer regions of a PPD.

High turbulence in the IM Lup protoplanetary disk

Context. Constraining turbulence in disks is key to understanding their evolution via the transport of angular momentum. Measurements of high turbulence remain elusive, and methods for estimating turbulence mostly rely on complex radiative transfer models of the data. Using the disk emission from IM Lup, a source proposed to be undergoing magneto-rotational instabilities (MRIs) and to possibly have high turbulence values in the upper disk layers, we present a new way of directly measuring turbulence without the need of radiative transfer or thermochemical models. Aims. Through the characterization of the CN and C2H emission in IM Lup, we aim to connect the information on the vertical and thermal structure of a particular disk region to derive the turbulence at that location. By using an optically thin tracer, it is possible to directly measure turbulence from the nonthermal broadening of the line. Methods. The vertical layers of the CN and C2H emission were traced directly from the channel maps using ALFAHOR. By comparing their position to that of optically thick CO observations, we were able to characterize the kinetic temperature of the emitting region. Using a simple parametric model of the line intensity with DISCMINER, we accurately measured the emission linewidth and separated the thermal and nonthermal components. Assuming that the nonthermal component is fully turbulent, we were able to directly estimate the turbulent motions at the studied radial and vertical location of CN emission. Results. IM Lup shows a high turbulence of Mach 0.4–0.6 at z/r ~ 0.25. Considering previous estimates of low turbulence near the midplane, this may indicate a vertical gradient in the disk turbulence, which is a key prediction in MRI studies. CN and C2H are both emitting from a localized upper disk region at z/r = 0.2–0.3, in agreement with thermochemical models.

Simulating X-Ray Reverberation in the Ultraviolet-emitting Regions of Active Galactic Nuclei Accretion Disks with Three-dimensional Multifrequency Radiation Magnetohydrodynamic Simulations

Active galactic nuclei (AGN) light curves observed with different wave bands show that the variability in longer wavelength bands lags the variability in shorter wavelength bands. Measuring these lags, or reverberation mapping, is used to measure the radial temperature profile and extent of AGN disks, typically with a reprocessing model that assumes X-rays are the main driver of the variability in other wavelength bands. To demonstrate how this reprocessing works with realistic accretion disk structures, we use 3D local shearing box multifrequency radiation magnetohydrodynamic simulations to model the UV-emitting region of an AGN disk, which is unstable to the magnetorotational instability and convection. At the same time, we inject hard X-rays (>1 keV) into the simulation box to study the effects of X-ray irradiation on the local properties of the turbulence and the resulting variability of the emitted UV light curve. We find that disk turbulence is sufficient to drive intrinsic variability in emitted UV light curves and that a damped random walk model is a good fit to this UV light curve for timescales >5 days. Meanwhile, X-ray irradiation has negligible impact on the power spectrum of the emitted UV light curve. Furthermore, the injected X-ray and emitted UV light curves are only correlated if there is X-ray variability on timescales >1 day, in which case we find a correlation coefficient r = 0.34. These results suggest that if the opacity for hard X-rays is scattering dominated as in the standard disk model, hard X-rays are not the main driver of reverberation signals.

Angular Momentum Transport in Binary Star Formation: The Enhancement of Magnetorotational Instability and Role of Outflows

The formation of binary stars is highly influenced by magnetic fields, which play a crucial role in transporting angular momentum. We conducted 3D numerical simulations of binary star accretion via a circumbinary disk, taking into account a magnetic field perpendicular to the disk and an infalling envelope. Our simulations reproduce the following phenomena: (1) the spiral arms associated with circumstellar disks; (2) the turbulence in the circumbinary disk, induced by magnetorotational instability (MRI); (3) a fast outflow launched from each circumstellar disk; and (4) a slow outflow from the circumbinary disk. The binary models exhibit a higher α-parameter than the corresponding single-star models, indicating that the binary stars enhance the MRI turbulence. Moreover, an infalling envelope also enhances the turbulence, leading to a high α-parameter. While the spiral arms promote radial flow, causing the transfer of mass and angular momentum within the circumbinary disk, the MRI turbulence and outflows are the main drivers of angular momentum transfer to reduce the specific angular momentum of the system.

Nonaxisymmetric modes of magnetorotational and possible hydrodynamical instabilities in the upcoming DRESDYN-MRI experiments: Linear and nonlinear dynamics

Nonaxisymmetric modes of magnetorotational and possible hydrodynamical instabilities in the upcoming DRESDYN-MRI experiments: Linear and nonlinear dynamics

Angular momentum and chemical transport by azimuthal magnetorotational instability in radiative stellar interiors

Context. The transport of angular momentum and chemical elements within evolving stars remains poorly understood. Asteroseismic and spectroscopic observations of low-mass main sequence stars and red giants reveal that their radiative cores rotate orders of magnitude slower than classical predictions from stellar evolution models and that the abundances of their surface light elements are too small. Magnetohydrodynamic (MHD) turbulence is considered a primary mechanism to enhance the transport in radiative stellar interiors but its efficiency is still largely uncertain. Aims. We explore the transport of angular momentum and chemical elements due to azimuthal magnetorotational instability, one of the dominant instabilities expected in differentially rotating radiative stellar interiors. Methods. We employed 3D MHD direct numerical simulations in a spherical shell of unstratified and stably stratified flows under the Boussinesq approximation. The background differential rotation was maintained by a volumetric body force. We examined the transport of chemical elements using a passive scalar. Results. We provide evidence of magnetorotational instability for purely azimuthal magnetic fields in the parameter regime expected from local and global linear stability analyses. Without stratification and when the Reynolds number Re and the background azimuthal field strength are large enough, we observed dynamo action driven by the instability at values of the magnetic Prandtl number Pm in the range 0.6 − 1, which is the smallest ever reported in a global setup. When considering stable stratification at Pm = 1, the turbulence is transitional and becomes less homogeneous and isotropic upon increasing buoyancy effects. The transport of angular momentum occurs radially outward and is dominated by the Maxwell stresses when stratification is large enough. We find that the turbulent viscosity decreases when buoyancy effects strengthen and scales with the square root of the ratio of the reference rotation rate Ωa to the Brunt–Väisälä frequency N. The chemical turbulent diffusion coefficient scales with stratification similarly to the turbulent viscosity, but is lower in amplitude so that the transport of chemicals is slower than the one of angular momentum, in agreement with recent stellar evolution models of low-mass stars. Conclusions. We show that the transport induced by azimuthal magnetorotational instability scales somewhat slowly with stratification and may enforce rigid rotations of red giant cores on a timescale of a few thousand years. In agreement with recent stellar evolution models of low-mass stars, the instability transports chemical elements less efficiently than angular momentum.

A large-scale magnetic field produced by a solar-like dynamo in binary neutron star mergers

The merger of two neutron stars launches a relativistic jet, which must be driven by a strong large-scale magnetic field. However, the magnetohydrodynamical mechanism required to build up this magnetic field remains uncertain. By performing an ab initio super-high-resolution neutrino-radiation magnetohydrodynamics merger simulation in full general relativity, we show that the αΩ dynamo mechanism, driven by the magnetorotational instability, builds up the large-scale magnetic field inside the long-lived remnant of the binary neutron star merger. As a result, the magnetic field induces a Poynting-flux-dominated relativistic outflow with an isotropic equivalent luminosity of ~1052 erg s−1 and a magnetically driven post-merger mass ejection of ~0.1 M⊙. Therefore, the magnetar hypothesis, in which an ultra-strongly magnetized neutron star drives a relativistic jet in binary neutron star mergers, is possible. Magnetars can be the engines of short, hard gamma-ray bursts, and they should be associated with very bright kilonovae, which current telescopes could observe. Therefore, this scenario is testable in future observations.

The gap-size influence on the excitation of magnetorotational instability in cylindricTaylor–Couette flows

The excitation conditions of the magnetorotational instability (MRI) are studied for axially unbounded Taylor–Couette (TC) flows of various gap widths between the cylinders. The cylinders are considered as made from both perfect-conducting or insulating material and the conducting fluid with a finite but small magnetic Prandtl number rotates with a quasi-Keplerian velocity profile. The solutions are optimized with respect to the wavenumber and the Reynolds number of the rotation of the inner cylinder. For the axisymmetric modes, we find the critical Lundquist number of the applied axial magnetic field: the lower, the wider the gap between the cylinders. A similar result is obtained for the induced cell structure: the wider the gap, the more spherical the cells are. The marginal rotation rate of the inner cylinder – for a fixed size of the outer cylinder – always possesses a minimum for not too wide and not too narrow gap widths. For perfect-conducting walls the minimum lies at $r_{{\rm in}}\simeq 0.4$ , where $r_{{\rm in}}$ is the ratio of the radii of the two rotating cylinders. The lowest magnetic field amplitudes to excite the instability are required for TC flows between perfect-conducting cylinders with gaps corresponding to $r_{{\rm in}}\simeq ~0.2$ . For even wider and also for very thin gaps the needed magnetic fields and rotation frequencies are shown to become rather huge. Also the non-axisymmetric modes with $|m|=1$ have been considered. Their excitation generally requires stronger magnetic fields and higher magnetic Reynolds numbers in comparison with those for the axisymmetric modes. If TC experiments with too slow rotation for the applied magnetic fields yield unstable modes of any azimuthal symmetry, such as the currently reported Princeton experiment (Wang et al., Phys. Rev. Lett., vol. 129, 115001), then also other players, including axial boundary effects, than the MRI-typical linear combination of current-free fields and differential rotation should be in the game.

Ab-initio General-relativistic Neutrino-radiation Hydrodynamics Simulations of Long-lived Neutron Star Merger Remnants to Neutrino Cooling Timescales

We perform the first 3D ab-initio general-relativistic neutrino-radiation hydrodynamics of a long-lived neutron star merger remnant spanning a fraction of its cooling timescale. We find that neutrino cooling becomes the dominant energy loss mechanism after the gravitational-wave dominated phase (∼20 ms postmerger). Electron flavor antineutrino luminosity dominates over electron flavor neutrino luminosity at early times, resulting in a secular increase of the electron fraction in the outer layers of the remnant. However, the two luminosities become comparable ∼20–40 ms postmerger. A dense gas of electron antineutrinos is formed in the outer core of the remnant at densities ∼1014.5 g cm−3, corresponding to temperature hot spots. The neutrinos account for ∼10% of the lepton number in this region. Despite the negative radial temperature gradient, the radial entropy gradient remains positive, and the remnant is stably stratified according to the Ledoux criterion for convection. A massive accretion disk is formed from the material squeezed out of the collisional interface between the stars. The disk carries a large fraction of the angular momentum of the system, allowing the remnant massive neutron star to settle to a quasi-steady equilibrium within the region of possible, stable, rigidly rotating configurations. The remnant is differentially rotating, but it is stable against the magnetorotational instability. Other MHD mechanisms operating on longer timescales are likely responsible for the removal of the differential rotation. Our results indicate the remnant massive neutron star is thus qualitatively different from a protoneutron stars formed in core-collapse supernovae.