Simulations In Full General Relativity Research Articles

Outflow energy and black-hole spin evolution in collapsar scenarios

We explore the collapsar scenario for long gamma-ray bursts by performing axisymmetric neutrino-radiation magnetohydrodynamics simulations in full general relativity for the first time. In this paper, we pay particular attention to the outflow energy and the evolution of the black-hole spin. We show that for a strong magnetic field with an aligned field configuration initially given, a jet is launched by magnetohydrodynamical effects before the formation of a disk and a torus, and after the jet launch, the matter accretion onto the black hole is halted by the strong magnetic pressure, leading to the spin-down of the black hole due to the Blandford-Znajek mechanism. The spin-down timescale depends strongly on the magnetic-field strength initially given because the magnetic-field strength on the black-hole horizon, which is determined by the mass infall rate at the jet launch, depends strongly on the initial condition, although the total jet-outflow energy appears to be huge >1053 erg depending only weakly on the initial field strength and configuration. For the models in which the magnetic-field configuration is not suitable for quick jet launch, a torus is formed, and after a long-term magnetic-field amplification, a jet can be launched. For this case, the matter accretion onto the black hole continues even after the jet launch, and black-hole spin-down is not found. We also find that the jet launch is often accompanied by the powerful explosion of the entire star with the explosion energy of order 1052 erg by magnetohydrodynamical effects. We discuss an issue of the overproduced energy for the early-jet-launch models. Published by the American Physical Society 2024

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.

Minidisc influence on flow variability in accreting spinning black hole binaries: simulations in full general relativity

ABSTRACT We perform magnetohydrodynamic simulations of accreting, equal-mass binary black holes in full general relativity focusing on the effect of spin and minidiscs on the accretion rate and Poynting luminosity variability. We report on the structure of the minidiscs and periodicities in the mass of the minidiscs, mass accretion rates, and Poynting luminosity. The accretion rate exhibits a quasi-periodic behaviour related to the orbital frequency of the binary in all systems that we study, but the amplitude of this modulation is dependent on the existence of persistent minidiscs. In particular, systems that are found to produce persistent minidiscs have a much weaker modulation of the mass accretion rate, indicating that minidiscs can increase the inflow time of matter on to the black holes, and dampen out the quasi-periodic behaviour. This finding has potential consequences for binaries at greater separations where minidiscs can be much larger and may dampen out the periodicities significantly.

Fate of twin stars on the unstable branch: Implications for the formation of twin stars

Hybrid hadron-quark equations of state that give rise a third family of stable compact stars have been shown to be compatible with the LIGO-Virgo event GW170817. Stable configurations in the third family are called hybrid hadron-quark stars. The equilibrium stable hybrid hadron-quark star branch is separated by the stable neutron star branch with a branch of unstable hybrid hadron-quark stars. The end-state of these unstable configurations has not been studied, yet, and it could have implications for the formation and existence of twin stars -- hybrid stars with the same mass as neutron stars but different radii. We modify existing hybrid hadron-quark equations of state with a first-order phase transition in order to guarantee a well-posed initial value problem of the equations of general relativistic hydrodynamics, and study the dynamics of non-rotating or rotating unstable twin stars via 3-dimensional simulations in full general relativity. We find that unstable twin stars naturally migrate toward the hadronic branch. Before settling into the hadronic regime, these stars undergo (quasi)radial oscillations on a dynamical timescale while the core bounces between the two phases. Our study suggests that it may be difficult to form stable twin stars if the phase transition is sustained over a large jump in energy density, and hence it may be more likely that astrophysical hybrid hadron-quark stars have masses above the twin star regime. We also study the minimum-mass instability for hybrid stars, and find that these configurations do not explode, unlike the minimum-mass instability for neutron stars. Additionally, our results suggest that oscillations between the two Quantum Chromodynamic phases could provide gravitational wave signals associated with such phase transitions in core-collapse supernovae and white dwarf-neutron star mergers.

Magnetohydrodynamic Simulations of Self-Consistent Rotating Neutron Stars with Mixed Poloidal and Toroidal Magnetic Fields.

We perform the first magnetohydrodynamic simulations in full general relativity of self-consistent rotating neutron stars (NSs) with ultrastrong mixed poloidal and toroidal magnetic fields. The initial uniformly rotating NS models are computed assuming perfect conductivity, stationarity, and axisymmetry. Although the specific geometry of the mixed field configuration can delay or accelerate the development of various instabilities known from analytic perturbative studies, all our models finally succumb to them. Differential rotation is developed spontaneously in the cores of our magnetars which, after sufficient time, is converted back to uniform rotation. The rapidly rotating magnetars show a significant amount of ejecta, which can be responsible for transient kilonova signatures. However, no highly collimated, helical magnetic fields or incipient jets, which are necessary for γ-ray bursts, arise at the poles of these magnetars by the time our simulations are terminated.

Impact of the nuclear symmetry energy on the post-merger phase of a binary neutron star coalescence

The nuclear symmetry energy plays a key role in determining the equation of state of dense, neutron-rich matter, which governs the properties of both terrestrial nuclear matter as well as astrophysical neutron stars. A recent measurement of the neutron skin thickness from the PREX Collaboration has lead to new constraints on the slope of the nuclear symmetry energy, $L$, which can be directly compared to inferences from gravitational wave observations of the first binary neutron star merger inspiral, GW170817. In this paper, we explore a new regime for potentially constraining the slope, $L$, of the nuclear symmetry energy with future gravitational wave events: the post-merger phase of a binary neutron star coalescence. In particular, we go beyond the inspiral phase, where imprints of the slope parameter $L$ may be inferred from measurements of the tidal deformability, to consider imprints on the post-merger dynamics, gravitational wave emission, and dynamical mass ejection. To this end, we perform a set of targeted neutron star merger simulations in full general relativity using new finite-temperature equations of state, which systematically vary $L$ while keeping the magnitude of the symmetry energy at the saturation density, $S$, fixed. We find that the post-merger dynamics and gravitational wave emission are mostly insensitive to the slope of the nuclear symmetry energy. In contrast, we find that dynamical mass ejection contains a weak imprint of $L$, with large values of $L$ leading to systematically enhanced ejecta.

Ultra-delayed Neutrino-driven Explosion of Rotating Massive-star Collapse

Abstract Long-term neutrino-radiation hydrodynamics simulations in full general relativity are performed for the collapse of rotating massive stars that are evolved from He-stars with initial masses of 20 and 32 M ⊙. It is shown that if the collapsing stellar core has sufficient angular momentum, the rotationally supported proto-neutron star (PNS) survives for seconds accompanying the formation of a massive torus of mass larger than 1 M ⊙. Subsequent mass accretion onto the central region produces a massive and compact central object, and eventually enhances the neutrino luminosity beyond 1053 erg s−1, resulting in a very delayed neutrino-driven explosion, in particular toward the polar direction. The kinetic energy of the explosion can be appreciably higher than 1052 erg for a massive progenitor star and compatible with that of energetic supernovae like broad-line type-Ic supernovae. By the subsequent accretion, the massive PNS collapses eventually into a rapidly spinning black hole, which could be a central engine for gamma-ray bursts if a massive torus surrounds it.

Long-term evolution of neutron-star merger remnants in general relativistic resistive magnetohydrodynamics with a mean-field dynamo term

Long-term neutrino-radiation resistive-magnetohydrodynamics simulations in full general relativity are performed for a system composed of a massive neutron star and a torus formed as a remnant of binary neutron-star mergers. The simulation is performed in axial symmetry incorporating a mean-field dynamo term for a hypothetical amplification of the magnetic-field strength. We first calibrate the mean-field dynamo parameters by comparing the results for the evolution of black hole--disk systems with viscous hydrodynamics results. We then perform simulations for the system of a remnant massive neutron star and a torus. As in the viscous hydrodynamics case, the mass ejection occurs primarily from the torus surrounding the massive neutron star. The total ejecta mass and electron fraction in the new simulation are similar to those in the viscous hydrodynamics case. However, the velocity of the ejecta can be significantly enhanced by magnetohydrodynamics effects caused by global magnetic fields.

Accretion onto a small black hole at the center of a neutron star.

We revisit the system consisting of a neutron star that harbors a small, possibly primordial, black hole at its center, focusing on a nonspinning black hole embedded in a nonrotating neutron star. Extending earlier treatments, we provide an analytical treatment describing the rate of secular accretion of the neutron star matter onto the black hole, adopting the relativistic Bondi accretion formalism for stiff equations of state that we presented elsewhere. We use these accretion rates to sketch the evolution of the system analytically until the neutron star is completely consumed. We also perform numerical simulations in full general relativity for black holes with masses up to nine orders of magnitude smaller than the neutron star mass, including a simulation of the entire evolution through collapse for the largest black hole mass. We construct relativistic initial data for these simulations by generalizing the black hole puncture method to allow for the presence of matter, and evolve these data with a code that is optimally designed to resolve the vastly different length scales present in this problem. We compare our analytic and numerical results, and provide expressions for the lifetime of neutron stars harboring such endoparasitic black holes.

Minidisk Dynamics in Accreting, Spinning Black Hole Binaries: Simulations in Full General Relativity

Abstract We perform magnetohydrodynamic simulations of accreting, equal-mass binary black holes in full general relativity focusing on the impact of black hole spin on the dynamical formation and evolution of minidisks. We find that during the late inspiral the sizes of minidisks are primarily determined by the interplay between the tidal field and the effective innermost stable orbit around each black hole. Our calculations support that a minidisk forms when the Hill sphere around each black hole is significantly larger than the black hole’s effective innermost stable orbit. As the binary inspirals, the radius of the Hill sphere decreases, and minidisks consequently shrink in size. As a result, electromagnetic signatures associated with minidisks may be expected to gradually disappear prior to merger when there are no more stable orbits within the Hill sphere. In particular, a gradual disappearance of a hard electromagnetic component in the spectrum of such systems could provide a characteristic signature of merging black hole binaries. For a binary of given total mass, the timescale to minidisk “evaporation” should therefore depend on the black hole spins and the mass ratio. We also demonstrate that accreting binary black holes with spin have a higher efficiency for converting accretion power to jet luminosity. These results could provide new ways to estimate black hole spins in the future.

Long-term evolution of a merger-remnant neutron star in general relativistic magnetohydrodynamics: Effect of magnetic winding

Long-term ideal and resistive magnetohydrodynamics (MHD) simulations in full general relativity are performed for a massive neutron star formed as a remnant of binary neutron star mergers. Neutrino radiation transport effects are taken into account as in our previous papers. The simulation is performed in axial symmetry and without considering dynamo effects as a first step. In the ideal MHD, the differential rotation of the remnant neutron star amplifies the magnetic-field strength by the winding in the presence of a seed poloidal field until the electromagnetic energy reaches $\sim 10\%$ of the rotational kinetic energy, $E_{\rm kin}$, of the neutron star. The timescale until the maximum electromagnetic energy is reached depends on the initial magnetic-field strength and it is $\sim 1$ s for the case that the initial maximum magnetic-field strength is $\sim 10^{15}$ G. After a significant amplification of the magnetic-field strength by the winding, the magnetic braking enforces the initially differentially rotating state approximately to a rigidly rotating state. In the presence of the resistivity, the amplification is continued only for the resistive timescale, and if the maximum electromagnetic energy reached is smaller than $\sim 3\%$ of $E_{\rm kin}$, the initial differential rotation state is approximately preserved. In the present context, the post-merger mass ejection is induced primarily by the neutrino irradiation/heating and the magnetic winding effect plays only a minor role for the mass ejection.

Gravitational waves from disks around spinning black holes: Simulations in full general relativity.

We present fully general-relativistic numerical evolutions of self-gravitating tori around spinning black holes with dimensionless spin a/M = 0.7 parallel or antiparallel to the disk angular momentum. The initial disks are unstable to the hydrodynamic Papaloizou-Pringle instability which causes them to grow persistent orbiting matter clumps. The effect of black hole spin on the growth and saturation of the instability is assessed. We find that the instability behaves similarly to prior simulations with nonspinning black holes, with a shift in frequency due to spin-induced changes in disk orbital period. Copious gravitational waves are generated by these systems, and we analyze their detectability by current and future gravitational wave observatories for a large range of masses. We find that systems of 10 M ⊙-relevant for black hole-neutron star mergers-are detectable by Cosmic Explorer out to ~300 Mpc, while DECIGO (LISA) will be able to detect systems of 1000 M ⊙ (105 M ⊙)-relevant for disks forming in collapsing supermassive stars-out to cosmological redshift of z ~ 5 (z ~ 1). Computing the accretion rate of these systems we find that these systems may also be promising sources of coincident electromagnetic signals.

Viscous evolution of a massive disk surrounding stellar-mass black holes in full general relativity

Long-term viscous neutrino-radiation hydrodynamics simulations in full general relativity are performed for a massive disk surrounding spinning stellar-mass black holes with mass $M_{\rm BH}=4$, $6$, and $10M_\odot$ and initial dimensionless spin $\chi \approx 0.8$. The initial disk is chosen to have mass $M_{\rm disk}\approx 0.1$ or $3M_\odot$ as plausible models of the remnants for the merger of black hole-neutron star binaries or the stellar core collapse from a rapidly rotating progenitor, respectively. For $M_{\rm disk} \approx 0.1M_\odot$ with the outer disk edge initially located at $r_{\rm out} \sim 200$ km, we find that $15$%-$20$% of $M_{\rm disk}$ is ejected and the average electron fraction of the ejecta is $\langle Y_e \rangle = 0.30$-$0.35$ as found in the previous study. For $M_{\rm disk} \approx 3M_\odot$, we find that $\approx 10$%-$20$% of $M_{\rm disk}$ is ejected for $r_{\rm out}\approx 200$-$1000$ km. In addition, $\langle Y_e \rangle$ of the ejecta can be enhanced to be $\gtrsim 0.4$ because the electron fraction is increased significantly during the long-term viscous expansion of the disk with high neutrino luminosity until the mass ejection sets in. Our results suggest that not heavy $r$-process elements but light trans-iron elements would be synthesized in the matter ejected from a massive torus surrounding stellar-mass black holes. We also find that the outcomes of the viscous evolution for the high-mass disk case is composed of a rapidly spinning black hole surrounded by a torus with a narrow funnel, which appears to be suitable for generating gamma-ray bursts.

Hybrid equation of state approach in binary neutron-star merger simulations

We investigate the use of hybrid equations of state in binary neutron-star simulations in full general relativity, where thermal effects are included in an approximate way through the adiabatic index $\Gamma_{\rm{th}}$. We employ a newly developed finite-temperature equation of state derived in the Brueckner-Hartree-Fock approach and carry out comparisons with the corresponding hybrid versions of the same equation of state, investigating how different choices of $\Gamma_{\rm{th}}$ affect the gravitational-wave signal and the hydrodynamical properties of the remnant. We also perform comparisons with the widely used SFHo equation of state, detailing the differences between the two cases. Overall, we determine that when using a hybrid equation of state in binary neutron-star simulations, the value of thermal adiabatic index $\Gamma_{\rm{th}} \approx 1.7$ best approximates the dynamical and thermodynamical behaviour of matter computed using complete, finite-temperature equations of state.

Mass ejection from neutron-star mergers

AbstractMerger of binary neutron stars and black hole-neutron star binaries is the promising source of short-hard gamma-ray bursts, the most promising site for the r-process nucleosynthesis, and the source of kilonovae. To theoretically predict the merger and mass ejection processes and resulting electromagnetic emission, numerical simulation in full general relativity (numerical relativity) is the unique approach. We summarize our current understanding for the processes of neutron-star mergers and subsequent mass ejection based on the results of long-term numerical-relativity simulations. We pay particular attention to the electron fraction of the ejecta.

Mass ejection from disks surrounding a low-mass black hole: Viscous neutrino-radiation hydrodynamics simulation in full general relativity

New viscous neutrino-radiation hydrodynamics simulations are performed for accretion disks surrounding a spinning black hole with low mass $3M_\odot$ and dimensionless spin 0.8 or 0.6 in full general relativity, aiming at modeling the evolution of a merger remnant of massive binary neutron stars or low-mass black hole-neutron star binaries. We reconfirm the following results found by previous studies of other groups: 15-30% of the disk mass is ejected from the system with the average velocity of $\sim $5-10% of the speed of light for the plausible profile of the disk as merger remnants. In addition, we find that for the not extremely high viscous coefficient case, the neutron richness of the ejecta does not become very high, because weak interaction processes enhance the electron fraction during the viscous expansion of the disk before the onset of the mass ejection, resulting in the suppression of the lanthanide synthesis. For high-mass disks, the viscous expansion timescale is increased by a longer-term neutrino emission, and hence, the electron fraction of the ejecta becomes even higher. We also confirm that the mass distribution of the electron fraction depends strongly on the magnitude of the given viscous coefficient. This demonstrates that a first-principle magnetohydrodynamics simulation is necessary for black hole-disk systems with sufficient grid resolution and with sufficiently long timescale (longer than seconds) to clarify the nucleosynthesis and electromagnetic signals from them.

Merger and Mass Ejection of Neutron Star Binaries

Mergers of binary neutron stars and black hole–neutron star binaries are among the most promising sources for ground-based gravitational-wave (GW) detectors and are also high-energy astrophysical phenomena, as illustrated by the observations of GWs and electromagnetic (EM) waves in the event of GW170817. Mergers of these neutron star binaries are also the most promising sites for r-process nucleosynthesis. Numerical simulation in full general relativity (numerical relativity) is a unique approach to the theoretical prediction of the merger process, GWs emitted, mass ejection process, and resulting EM emission. We summarize the current understanding of the processes of neutron star mergers and subsequent mass ejection based on the results of the latest numerical-relativity simulations. We emphasize that the predictions of the numerical-relativity simulations agree broadly with the optical and IR observations of GW170817.

Dynamical stability of quasitoroidal differentially rotating neutron stars

We investigate the dynamical stability of relativistic, differentially rotating, quasi-toroidal models of neutron stars through hydrodynamical simulations in full general relativity. We find that all quasi-toroidal configurations studied in this work are dynamically unstable against the growth of non-axisymmetric modes. Both one-arm and bar mode instabilities grow during their evolution. We find that very high rest mass configurations collapse to form black holes. Our calculations suggest that configurations whose rest mass is less than the binary neutron star threshold mass for prompt collapse to black hole transition dynamically to spheroidal, differentially rotating stars that are dynamically stable, but secularly unstable. Our study shows that the existence of extreme quasi-toroidal neutron star equilibrium solutions does not imply that long-lived binary neutron star merger remnants can be much more massive than previously found. Finally, we find models that are initially supra-Kerr ($J/M^2>1$) and undergo catastrophic collapse on a dynamical timescale, in contrast to what was found in earlier works. However, cosmic censorship is respected in all of our cases. Our work explicitly demonstrates that exceeding the Kerr bound in rotating neutron star models does not imply dynamical stability.

Are fast radio bursts the most likely electromagnetic counterpart of neutron star mergers resulting in prompt collapse?

Inspiraling and merging binary neutron stars (BNSs) are important sources of both gravitational waves and coincident electromagnetic counterparts. If the BNS total mass is larger than a threshold value, a black hole ensues promptly after merger. Through a statistical study in conjunction with recent LIGO/Virgo constraints on the nuclear equation of state, we estimate that up to $\sim 25\%$ of BNS mergers may result in prompt collapse. Moreover, we find that most models of the BNS mass function we study here predict that the majority of prompt-collapse BNS mergers have $q\gtrsim 0.8$. Prompt-collapse BNS mergers with mass ratio $q \gtrsim 0.8$ may not be accompanied by detectable kilonovae or short gamma-ray bursts, because they unbind a negligible amount of mass and form negligibly small accretion disks onto the remnant black hole. We call such BNS mergers "orphan". However, recent studies have found that ${10^{41-43}(B_p/10^{12}\rm G)^2 erg\, s^{-1}}$ electromagnetic signals can be powered by magnetospheric interactions several milliseconds prior to merger. Moreover, the energy stored in the magnetosphere of an orphan BNS merger remnant will be radiated away in ${\mathcal O}(1\ \rm ms)$. Through simulations in full general relativity of BNSs endowed with an initial dipole magnetosphere, we find that the energy in the magnetosphere following black hole formation is $E_B \sim 10^{39-41} (B_p/10^{12}\rm G)^2$ erg. Radiating $\sim 1\%$ of $E_B$ in 1 ms, as has been found in previous studies, matches the premerger magnetospheric luminosity. These magnetospheric signals are not beamed, and their duration and power agrees with those of non-repeating fast radio bursts (FRBs). These results combined with our statistical study suggest that a non-repeating FRB may be the most likely electromagnetic counterpart of prompt-collapse BNSs.

The fate of dense scalar stars

Long-lived pseudo-solitonic objects, known as oscillons/oscillatons, which we collectively call real scalar stars, are ubiquitous in early Universe cosmology of scalar field theories. Typical examples are axions stars and moduli stars. Using numerical simulations in full general relativity to include the effects of gravity, we study the fate of real scalar stars and find that depending on the scalar potential they are either meta-stable or collapse to black holes. In particular we find that for KKLT potentials the configurations are meta-stable despite the asymmetry of the potential, consistently with the results from lattice simulations that do not include gravitational effects. For α-attractor potentials collapse to black holes is possible in a region of the parameter space where scalar stars would instead seem to be meta-stable or even disperse without including gravity. Each case gives rise to different cosmological implications which may affect the stochastic spectrum of gravitational waves.