Search

We perform fully general-relativistic hydrodynamics simulations of binary neutron star mergers over $100\,\rm ms$ post-merger to investigate the dynamics of remnant massive neutron stars (NSs). Our focus is mainly on the analysis of convective stability and mode characteristics of the massive NSs. We derive stability criteria for hot, differentially rotating relativistic stars that account for both buoyant and rotational restoring forces, and apply them for the first time to the post-merger massive NSs. Our results show no evidence of large-scale convective instability, as both angle-averaged specific entropy and specific angular momentum increase outward within the massive NSs. Rotational effects significantly enhance stability for local regions that would be otherwise unstable by the Schwarzschild criterion. Additionally, our mode analysis of matter fields and gravitational waves reveals no excitation of inertial modes after the damping of quadrupolar $f$-modes in the massive NSs, contrasting with previous studies. As in many previous works, we observe the excitation of an $m=1$ one-armed mode. However, we also find that the growth of the $m=1$ mode amplitude after the merger may correlate strongly with the violation of linear momentum conservation, indicating that we cannot reject the possibility that the excitation of the one-armed mode has a numerical origin., Comment: 22 pages, 18 figures; submitted to Phys. Rev. D

We performed the longest numerical-relativity neutrino-radiation magnetohydrodynamics simulation for a binary neutron star merger that extends to $\approx1.5\mathrm{\,s}$ after the merger. We consider the binary model that undergoes the prompt collapse to a black hole after the merger with asymmetric mass 1.25$\,M_{\odot}$ and 1.65$\,M_{\odot}$ and SFHo equation of state. We find the Poynting flux-driven collimated outflow as well as the gravitational wave emission, neutrino emission, dynamical mass ejection, and post-merger mass ejection facilitated by magnetorotational instability-driven turbulent viscosity in a single self-consistent binary neutron star merger simulation. A magnetosphere dominated by the aligned global magnetic field penetrating the black hole develops along the black-hole spin axis after the turbulence in the remnant disk is enhanced. A jet with the Poynting flux with isotropic-equivalent luminosity of $\sim10^{49}\mathrm{\,erg/s}$ is launched, and the duration of the high luminosity is expected to be $O(1)\mathrm{\,s}$., Comment: 14 pages, 9 figures

We study the merger of black hole-neutron star (BH-NS) binaries in numerical relativity, focusing on the properties of the remnant disk and the ejecta, varying the mass of compactness of the NS and the mass and spin of the BH. We find that within the precision of our numerical simulations, the remnant disk mass and ejecta mass normalized by the NS baryon mass ($\hat{M}_{\rm{rem}}$ and $\hat{M}_{\rm{eje}}$, respectively), and the cutoff frequency $f_{\rm{cut}}$ normalized by the initial total gravitational mass of the system at infinite separation approximately agree among the models with the same NS compactness $C_{\rm{NS}}=M_{\rm{NS}}/R_{\rm{NS}}$, mass ratio $Q=M_{\rm{BH}}/M_{\rm{NS}}$, and dimensionless BH spin $\chi_{\rm{BH}}$ irrespective of the NS mass $M_{\rm{NS}}$ in the range of $1.092$--$1.691\,M_\odot$. This result shows that the merger outcome depends sensitively on $Q$, $\chi_{\rm BH}$, and $C_{\rm{NS}}$ but only weekly on $M_{\rm{NS}}$. This justifies the approach of studying the dependence of NS tidal disruptions on the NS compactness by fixing the NS mass but changing the EOS. We further perform simulations with massive NSs of $M_{\rm{NS}}=1.8M_{\odot}$, and compare our results of $\hat{M}_{\rm{rem}}$ and $\hat{M}_{\rm{eje}}$ with those given by existing fitting formulas to test their robustness for more compact NSs. We find that the fitting formulas obtained in the previous studies are accurate within the numerical errors assumed, while our results also suggest that further improvement is possible by systematically performing more precise numerical simulations., Comment: 19 pages, 14 figures

We study long-term evolution of the matter ejected in a black-hole neutron-star (BH-NS) merger employing the results of a long-term numerical-relativity simulation and nucleosynthesis calculation, in which both dynamical and post-merger ejecta formation is consistently followed. In particular, we employ the results for the merger of a $1.35\,M_\odot$ NS and a $5.4\,M_\odot$ BH with the dimensionless spin of 0.75. We confirm the finding in the previous studies that thermal pressure induced by radioactive heating in the ejecta significantly modifies the morphology of the ejecta. We then compute the kilonova (KN) light curves employing the ejecta profile obtained by the long-term evolution. We find that our present BH-NS model results in a KN light curve that is fainter yet more enduring than that observed in AT2017gfo. This is due to the fact that the emission is primarily powered by the lanthanide-rich dynamical ejecta, in which a long photon diffusion time scale is realized by the large mass and high opacity. While the peak brightness of the KN emission in both the optical and near-infrared bands is fainter than or comparable to those of binary NS models, the time-scale maintaining the peak brightness is much longer in the near-infrared band for the BH-NS KN model. Our result indicates that a BH-NS merger with massive ejecta can observationally be identified by the long lasting ($>$two weeks) near-infrared emission., Comment: 21 pages, 13 figures, Accepted for publication in MNRAS

We examine nucleosynthesis in the ejecta of black hole-neutron star mergers based on the results of long-term neutrino-radiation-magnetohydrodynamics simulations for the first time. We find that the combination of dynamical and post-merger ejecta reproduces a solar-like r-process pattern. Moreover, the enhancement level of actinides is highly sensitive to the distribution of both electron fraction and the velocity of the dynamical ejecta. Our result implies that the mean electron fraction of dynamical ejecta should be >~ 0.05 in order to reconcile the nucleosynthetic abundances with those in r-process-enhanced, actinide-boost stars. Since the tidal ejecta preserve the neutron-richness in the inner crust of pre-merging neutron stars, this result provides an important constraint for nuclear equations of state, if black hole-neutron star mergers are responsible for actinide-boost stars., Comment: 8 pages, 6 figures (including Supplemental Material), accepted for publication in PRL

As a follow-up study of our previous work, numerical-relativity simulations for seconds-long black hole-neutron star mergers are performed for a variety of setups. Irrespective of the initial and symmetry conditions, we find qualitatively universal evolution processes: The dynamical mass ejection takes place together with a massive accretion disk formation after the neutron star is tidally disrupted; Subsequently, the magnetic field in the accretion disk is amplified by the magnetic winding, Kelvin-Helmholtz instability, and magnetorotational instability, which establish a turbulent state inducing the dynamo and angular momentum transport; The post-merger mass ejection by the effective viscous processes stemming from the magnetohydrodynamics turbulence sets in at $\sim300$-$500$ ms after the merger and continues for several hundred ms; A magnetosphere near the black-hole spin axis is developed and the collimated strong Poynting flux is generated with its lifetime of $\sim0.5$-$2$ s. We have newly found that the model of no equatorial-plane symmetry shows the reverse of the magnetic-field polarity in the magnetosphere, which is caused by the dynamo associated with the magnetorotational instability in the accretion disk. The model with initially toroidal fields shows the tilt of the disk and magnetosphere in the late post-merger stage because of the anisotropic post-merger mass ejection. These effects could terminate the strong Poynting-luminosity stage within the timescale of $\sim0.5$-$2$ s., Comment: 28 pages, 20 figures, and 2 tables. arXiv admin note: text overlap with arXiv:2111.04621

We perform a general-relativistic neutrino-radiation magnetohydrodynamic simulation of a one second-long binary neutron star merger on Japanese supercomputer Fugaku using about $72$ million CPU hours with $20,736$ CPUs. We consider an asymmetric binary neutron star merger with masses of $1.2$ and $1.5M_\odot$ and a `soft' equation of state SFHo. It results in a short-lived remnant with the lifetime of $\approx 0.017$\,s, and subsequent massive torus formation with the mass of $\approx 0.05M_\odot$ after the remnant collapses to a black hole. For the first time, we confirm that after the dynamical mass ejection, which drives the fast tail and mildly relativistic components, the post-merger mass ejection from the massive torus takes place due to the magnetorotational instability-driven turbulent viscosity and the two ejecta components are seen in the distributions of the electron fraction and velocity with distinct features., Comment: Accepted in PRL. 8 pages, 5 figure, Supplement Material is https://www2.yukawa.kyoto-u.ac.jp/~kenta.kiuchi/anime/FUGAKU/FUGAKU2022_Supplement_Material.pdf

Seconds-long numerical-relativity simulations for black hole-neutron star mergers are performed for the first time to obtain a self-consistent picture of the merger and post-merger evolution processes. To investigate the case that tidal disruption takes place, we choose the initial mass of the black hole to be $5.4M_\odot$ or $8.1M_\odot$ with the dimensionless spin of 0.75. The neutron-star mass is fixed to be $1.35M_\odot$. We find that after the tidal disruption, dynamical mass ejection takes place spending $\lesssim 10$ ms together with the formation of a massive accretion disk. Subsequently, the magnetic field in the disk is amplified by the magnetic winding and magnetorotational instability, establishing a turbulent state and inducing the angular momentum transport. The post-merger mass ejection by the magnetically-induced viscous effect sets in at $\sim 300$-500\,ms after the tidal disruption, at which the neutrino luminosity drops below $\sim 10^{51.5}\,{\rm erg/s}$, and continues for several hundreds ms. A magnetosphere near the rotational axis of the black hole is developed after the matter and magnetic flux fall into the black hole from the accretion disk, and high-intensity Poynting flux generation sets in at a few hundreds ms after the tidal disruption. The intensity of the Poynting flux becomes low after the significant post-merger mass ejection, because the opening angle of the magnetosphere increases. The lifetime for the stage with the strong Poynting flux is $1$-2 s, which agrees with the typical duration of short-hard gamma-ray bursts., Comment: 25 pages, 17 figures, and 2 tables

The third observing run of the LIGO-Virgo collaboration has resulted in about hundred gravitational-wave triggers including the binary neutron star merger GW190425. However, none of these events have been accompanied with an electromagnetic transient found during extensive follow-up searches. In this article, we perform new numerical-relativity simulations of binary neutron star and black hole - neutron star systems that have a chirp mass consistent with GW190425. Assuming that the GW190425's sky location was covered with sufficient accuracy during the electromagnetic follow-up searches, we investigate whether the non-detection of the kilonova is compatible with the source parameters estimated through the gravitational-wave analysis and how one can use this information to place constraints on the properties of the system. Our simulations suggest that GW190425 is incompatible with an unequal mass binary neutron star merger with a mass ratio $q<0.8$ when considering stiff or moderately stiff equations of state if the binary was face-on and covered by the observation. Our analysis shows that a detailed observational result for kilonovae will be useful to constrain the mass ratio of binary neutron stars in future events., Comment: 10 pages, 7 figures

We systematically perform numerical-relativity simulations for low-mass black hole-neutron star mergers for the models with seven mass ratios $Q=M_{\rm BH}/M_{\rm NS}$ ranging from 1.5 to 4.4, and three neutron-star equations of state, focusing on properties of matter remaining outside the black hole and ejected dynamically during the merger. We pay particular attention to the dependence on the mass ratio of the binaries. It is found that the rest mass remaining outside the apparent horizon after the merger depends only weakly on the mass ratio for the models with low mass ratios. It is also clarified that the rest mass of the ejecta has a peak at $Q \sim 3$, and decreases steeply as the mass ratio decreases for the low mass-ratio case. We present a novel analysis method for the behavior of matter during the merger, focusing on the matter distribution in the phase space of specific energy and specific angular momentum. Then we model the matter distribution during and after the merger. Using the result of the analysis, we discuss the properties of the ejecta., Comment: 16 pages, 9 figures

We argue that the kilonova/macronova associated with the gravitational-wave event GW190425 could have been bright enough to be detected if it was caused by the merger of a low-mass black hole and a neutron star. Although tidal disruption occurs for such a low-mass black hole as is generally expected, the masses of the dynamical ejecta are limited to <~10^{-3}M_sun, which is consistent with previous work in the literature. The remnant disk could be as massive as 0.05--0.1M_sun, and the disk outflow of ~0.01--0.03M_sun is likely to be driven by viscous or magnetohydrodynamic effects. The disk outflow may not be neutron-rich enough to synthesize an abundance of lanthanide elements, even in the absence of strong neutrino emitter, if the ejection is driven on the viscous time scale of >~0.3s. If this is the case, the opacity of the disk outflow is kept moderate, and a kilonova/macronova at the distance of GW190425 reaches a detectable brightness of 20--21mag at 1day after merger for most viewing angles. If some disk activity ejects the mass within ~0.1s, instead, lanthanide-rich outflows would be launched and the detection of emission becomes challenging. Future possible detections of kilonovae/macronovae from GW190425-like systems will disfavor the prompt collapse of binary neutron stars and a non-disruptive low-mass black hole--neutron star binary associated with a small neutron-star radius, whose mass ejection is negligible. The host-galaxy distance will constrain the viewing angle and deliver further information about the mass ejection., Comment: 8 pages, 1 figure, matched to the published version with an additional reference