BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Koutarou Kyutoku

Koutarou Kyutoku appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-ph.SRastro-ph.COhep-phnucl-thastro-ph.GA

Catalog footprint

What is connected

42works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2022arXiv

Comprehensive study of mass ejection and nucleosynthesis in binary neutron star mergers leaving short-lived massive neutron stars

By performing general relativistic hydrodynamics simulations with an approximate neutrino-radiation transfer, the properties of ejecta in dynamical and post-merger phases are investigated for the cases in which the remnant massive neutron star collapses into a black hole in $\lesssim 20$ ms after the onset of the merger. The dynamical mass ejection is investigated in three-dimensional simulations. The post-merger mass ejection is investigated in two-dimensional axisymmetric simulations with viscosity using the three-dimensional post-merger systems as the initial conditions. We show that the typical neutron-richness of the dynamical ejecta is higher for the merger of more asymmetric binaries; hence, heavier $r$-process nuclei are dominantly synthesized. The post-merger ejecta are shown to have only a mild neutron-richness, which results in the production of lighter $r$-process nuclei, irrespective of binary mass ratios. Because of the larger disk mass, the post-merger ejecta mass is larger for more asymmetric binary mergers. Thus, the post-merger ejecta can compensate for the underproduced lighter $r$-process nuclei for asymmetric merger cases. As a result, by summing up both ejecta components, the solar residual $r$-process pattern is reproduced within the average deviation of a factor of three, irrespective of the binary mass ratio. Our result also indicates that the (about a factor of a few) light-to-heavy abundance scatter observed in $r$-process-enhanced stars can be attributed to variation in the binary mass ratio and total mass. Implications of our results associated with the mass distribution of compact neutron star binaries and the magnetar scenario of short gamma-ray bursts are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

How many extra-Galactic stellar-mass binary black holes will be detected by space gravitational-wave interferometers?

On the basis of GWTC-3, we discuss the detection prospect of extra-Galactic binary black holes (BBHs) by space gravitational-wave interferometers. In particular, targeting BBHs with component masses around 5-100$M_\odot$, we directly incorporate the chirp mass distribution of the 62 BBHs detected at high significance. We find that, due to the reduction of both the comoving merger rate and a weighted average of chirp masses, the expected detection numbers are generally much smaller than the results obtained by the same authors immediately after the report of GW150914. For LISA, the total BBH detections are estimated to be $N_{\rm tot}\sim 2 (T/4{\rm yr})^{3/2}(ρ_{\rm thr}/10)^{-3}$, dominated by nearly monochromatic BBHs ($ρ_{\rm thr}$: the detection threshold, $T$: the observational period). TianQin will have a total detection number $N_{\rm tot}$ similar to LISA. Meanwhile, TianQin has potential to find $N_{\rm mer}\sim0.6 (T/4{\rm yr})^{7/4}(ρ_{\rm thr}/10)^{-3}$ BBHs that merge in the observational period. This number for merging BBHs is 4-5 times larger than that of LISA, because of the difference between the optimal bands. We also investigate prospects for joint operations of multiple detectors, finding that concurrent observations will be more advantageous than sequential ones.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Reducing orbital eccentricity in initial data of black hole--neutron star binaries in the puncture framework

We develop a method to compute low-eccentricity initial data of black hole--neutron star binaries in the puncture framework extending previous work on other types of compact binaries. In addition to adjusting the orbital angular velocity of the binary, the approaching velocity of a neutron star is incorporated by modifying the helical Killing vector used to derive equations of the hydrostationary equilibrium. The approaching velocity of the black hole is then induced by requiring the vanishing of the total linear momentum of the system, differently from the case of binary black holes in the puncture framework where the linear momentum of each black hole is specified explicitly. We successfully reduce the orbital eccentricity to <~0.001 by modifying the parameters iteratively using simulations of ~3 orbits both for nonprecessing and precessing configurations. We find that empirical formulas for binary black holes derived in the excision framework do not reduce the orbital eccentricity to ~0.001 for black hole--neutron star binaries in the puncture framework, although they work for binary neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the possibility of GW190425 being a black hole--neutron star binary merger

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Properties of the remnant disk and the dynamical ejecta produced in low-mass black hole-neutron star mergers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sub-radian-accuracy gravitational waves from coalescing binary neutron stars II: Systematic study on the equation of state, binary mass, and mass ratio

We report results of numerical relativity simulations for {\it new} 26 non-spinning binary neutron star systems with 6 grid resolutions using an adaptive mesh refinement numerical re\ lativity code {\tt SACRA-MPI}. The finest grid spacing is $\approx 64$--$85$ m, depending on the systems. First, we derive long-term high-precision inspiral gravitational waveforms and show that the accumulated gravitational-wave phase error due to the finite grid resolution is less than $0.5$ rad during more than $200$ rad phase evolution irrespective of the systems. We also find that the gravitational-wave phase error for a binary system with a tabulated equation of state (EOS) is comparable to that for a piecewise polytropic EOS. Then we validate the SACRA inspiral gravitational waveform template, which will be used to extract tidal deformability from gravitational wave observation, and find that accuracy of \ our waveform modeling is $\lesssim 0.1$ rad in the gravitational-wave phase and $\lesssim 20 \%$ in the gravitational-wave amplitude up to the gravitational-wave frequency $1000$ Hz.\ Finally, we calibrate the proposed universal relations between a post-merger gravitational wave signal and tidal deformability/neutron star radius in the literature and show that th\ ey suffer from systematics and many relations proposed as universal are not very universal. Improved fitting formulae are also proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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 $χ\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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

r-Process enrichment in the Galactic halo characterized by nucleosynthesis variation in the ejecta of coalescing neutron star binaries

A large star-to-star variation in the abundances of r-process elements, as seen in the [Eu/Fe] ratio for Galactic halo stars, is a prominent feature that is distinguishable from other heavy elements. It is, in part, caused by the presence of highly r-process enriched stars, classified as r-II stars ([Eu/Fe]>+1). In parallel, halo stars show that the ratio of a light r-process element (Y) to Eu is tightly correlated with [Eu/Fe], giving the lowest [Y/Eu] ratio that levels off at r-II stars. On the other hand, recent hydrodynamical simulations of coalescing double neutron stars (cNSNSs) have suggested that r-process sites may be separated into two classes providing different electron-fraction distributions: tidally-driven dynamical ejecta and (dynamical or postmerger) non-tidal ejecta. Here, we show that a widely spanning feature of [Eu/Fe] can be reproduced by models that consider the different masses of tidally-driven dynamical ejecta from both cNSNSs and coalescing black hole/neutron star binaries (cBHNSs). In addition, the observed [Y/Eu] trend is explained by the combined nucleosynthesis in two kinds of ejecta with varying mass asymmetry in double NS systems. Our scenario suggests that massive tidally-driven dynamical ejecta accompanied by massive non-tidal part from cNSNSs or cBHNSs could alone accommodate r-II abundances, including an actinide boost in some cases. The event rate for cNSNSs estimated from our study agrees with the latest result of ~1000 (90% confidence interval of 110-3840) Gpc$^{-3}$yr$^{-1}$ by gravitational-wave detection, and a few events per Gpc$^3$ per year of cBHNSs associated with r-process production are predicted to emerge.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Concise estimate of the expected number of detections for stellar-mass binary black holes by eLISA

We study prospects for detecting extragalactic binary black holes similar to GW150914 by evolved Laser Interferometer Space Antenna (eLISA). We find that the majority of detected binary black holes will not merge within reasonable observation periods of eLISA in any configuration. While long-arm detectors are highly desired for promoting multiband gravitational-wave astronomy by increasing the detections of merging binaries, the number of total detections can be increased also by improving the acceleration noise. A monochromatic approximation works well to derive semiquantitative features of observational prospects for non-merging binaries with clearly indicating the parameter dependence. Our estimate also suggests that the number of galaxies in the error volume is so small that the host galaxy may be determined uniquely with high confidence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical mass ejection from the merger of asymmetric binary neutron stars: Radiation-hydrodynamics study in general relativity

We perform neutrino radiation-hydrodynamics simulations for the merger of asymmetric binary neutron stars in numerical relativity. Neutron stars are modeled by soft and moderately stiff finite-temperature equations of state (EOS). We find that the properties of the dynamical ejecta such as the total mass, neutron richness profile, and specific entropy profile depend on the mass ratio of the binary systems for a given EOS in a unique manner. For the soft EOS (SFHo), the total ejecta mass depends weakly on the mass ratio, but the average of electron number per baryon ($Y_e$) and specific entropy ($s$) of the ejecta decreases significantly with the increase of the degree of mass asymmetry. For the stiff EOS (DD2), with the increase of the mass asymmetry degree, the total ejecta mass significantly increases while the average of $Y_e$ and $s$ moderately decreases. We find again that only for the soft EOS (SFHo), the total ejecta mass exceeds $0.01M_\odot$ irrespective of the mass ratio chosen in this paper. The ejecta have a variety of electron number per baryon with its average approximately between $Y_e \sim 0.2$ and $\sim 0.3$ irrespective of the EOS employed, which is well-suited for the production of the r-process heavy elements (second and third peaks), although its averaged value decreases with the increase of the degree of mass asymmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effects of neutron-star dynamic tides on gravitational waveforms within the effective-one-body approach

Extracting the unique information on ultradense nuclear matter from the gravitational waves emitted by merging, neutron-star binaries requires robust theoretical models of the signal. We develop a novel effective-one-body waveform model that includes, for the first time, dynamic (instead of only adiabatic) tides of the neutron star, as well as the merger signal for neutron-star--black-hole binaries. We demonstrate the importance of the dynamic tides by comparing our model against new numerical-relativity simulations of nonspinning neutron-star--black-hole binaries spanning more than 24 gravitational-wave cycles, and to other existing numerical simulations for double neutron-star systems. Furthermore, we derive an effective description that makes explicit the dependence of matter effects on two key parameters: tidal deformability and fundamental oscillation frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Measurability of the tidal deformability by gravitational waves from coalescing binary neutron stars

Combining new gravitational waveforms derived by long-term (14--16 orbits) numerical-relativity simulations with waveforms by an effective-one-body (EOB) formalism for coalescing binary neutron stars, we construct hybrid waveforms and estimate the measurability for the dimensionless tidal deformability of the neutron stars, $Λ$, by advanced gravitational-wave detectors. We focus on the equal-mass case with the total mass $2.7M_\odot$. We find that for an event at a hypothetical effective distance of $D_{\rm eff}=200$ Mpc, the distinguishable difference in the dimensionless tidal deformability will be $\approx 100$, 400, and 800 at 1-$σ$, 2-$σ$, and 3-$σ$ levels, respectively, for advanced LIGO. If the true equation of state is stiff and the typical neutron-star radius is $R \gtrsim 13 $ km, our analysis suggests that the radius will be constrained within $\approx 1$ km at 2-$σ$ level for an event at $D_{\rm eff}=200$ Mpc. On the other hand, if the true equation of state is soft and the typical neutron-star radius is $R\lesssim 12$ km , it will be difficult to narrow down the equation of state among many soft ones, although it is still possible to discriminate the true one from stiff equations of state with $R\gtrsim 13$ km. We also find that gravitational waves from binary neutron stars will be distinguished from those from spinless binary black holes at more than 2-$σ$ level for an event at $D_{\rm eff}=200$ Mpc. The validity of the EOB formalism, Taylor-T4, and Taylor-F2 approximants as the inspiral waveform model is also examined.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Models of Kilonova/macronova emission from black hole-neutron star mergers

Black hole-neutron star mergers are among the promising gravitational-wave sources for ground-based detectors, and gravitational waves from black hole-neutron star mergers are expected to be detected in the next few years. Simultaneous detection of electromagnetic counterparts with gravitational-wave detection provides rich information about the merger events. Among the possible electromagnetic counterparts from the black hole-neutron star merger, the emission powered by the decay of radioactive r-process nuclei, so called kilonova/macronova, is one of the best targets for follow-up observation. We derive fitting formulas for the mass and the velocity of ejecta from a generic black hole-neutron star merger based on recently performed numerical relativity simulations. We combined these fitting formulas with a new semi-analytic model for a black hole-neutron star kilonova/macronova lightcurve which reproduces the results of radiation-transfer simulations. Specifically, the semi-analytic model reproduces the result of each band magnitude obtained by the radiation transfer simulations within ~1 mag. By using this semi-analytic model, we found that, at 400 Mpc, the kilonova/macronova is as bright as 22-24 mag for the cases with a small chirp mass and a high black hole spin, and >28 mag for a large chirp mass and a low black hole spin. We also apply our model to GRB130603B as an illustration, and show that a black hole-neutron star merger with a rapidly spinning black hole and a large neutron star radius is favored.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The unreasonable weakness of r-process cosmic rays in the neutron-star-merger nucleosynthesis scenario

We reach the robust conclusion that, by combining the observed cosmic rays of r-process elements with the fact that the velocity of the neutron-star-merger ejecta is much higher than that of the supernova ejecta, either (1) the reverse shock in the neutron-star-merger ejecta is a very inefficient accelerator that converts less than 0.003% of the ejecta kinetic energy to the cosmic-ray energy or (2) the neutron star merger is not the origin of the Galactic r-process elements. We also find that the acceleration efficiency should be less than 0.1% for the reverse shock of the supernova ejecta with the observed cosmic rays lighter than the iron.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Aligned spin neutron star-black hole mergers: a gravitational waveform amplitude model

The gravitational radiation emitted during the merger of a black hole with a neutron star is rather similar to the radiation from the merger of two black holes when the neutron star is not tidally disrupted. When tidal disruption occurs, gravitational waveforms can be broadly classified in two groups, depending on the spatial extent of the disrupted material. Extending previous work by some of us, here we present a phenomenological model for the gravitational waveform amplitude in the frequency domain encompassing the three possible outcomes of the merger: no tidal disruption, "mild" and "strong" tidal disruption. The model is calibrated to 134 general-relativistic numerical simulations of binaries where the black hole spin is either aligned or antialigned with the orbital angular momentum. All simulations were produced using the SACRA code and piecewise polytropic neutron star equations of state. The present model can be used to determine when black-hole binary waveforms are sufficient for gravitational-wave detection, to extract information on the equation of state from future gravitational-wave observations, to obtain more accurate estimates of black hole-neutron star merger event rates, and to determine the conditions under which these systems are plausible candidates as central engines of gamma-ray bursts, macronovae and kilonovae.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Black hole-neutron star binary merger: Dependence on black hole spin orientation and equations of state

We systematically performed numerical-relativity simulations for black hole (BH) - neutron star (NS) binary mergers with a variety of the BH spin orientation and equations of state (EOS) of the NS. The initial misalignment angles of the BH spin are chosen in the range of i_tilt,0 = 30--90[deg.]. We employed four models of NS EOS with which the compactness of the NS is in the range of C = M_NS/R_NS = 0.138--0.180, where M_NS and R_NS are the mass and the radius of the NS, respectively. The mass ratio of the BH to the NS, Q = M_BH/M_NS, and the dimensionless spin parameter of the BH, chi, are chosen to be Q = 5 and chi = 0.75, together with M_NS = 1.35 M_sun. We obtain the following results: (i) The inclination angle of i_tilt,0 < 70[deg.] and i_tilt,0 < 50[deg.] are required for the formation of a remnant disk with its mass larger than 0.1 M_sun for the case C = 0.140 and C = 0.160, respectively, while the disk mass is always smaller than 0.1M_sun for C = 0.175. The ejecta with its mass larger than 0.01 M_sun is obtained for i_tilt,0 < 85[deg.] with C = 0.140, for i_tilt,0 < 65[deg.] with C = 0.160, and for i_tilt,0 < 30[deg.] with C = 0.175. (ii) The rotational axis of the dense part of the remnant disk is approximately aligned with the remnant BH spin for i_tilt,0 = 30[deg.]. On the other hand, the disk axis is misaligned initially with ~ 30[deg.] for i_tilt,0 = 60[deg.], and the alignment with the remnant BH spin is achieved at ~ 50--60 ms after the onset of merger. The accretion time scale of the remnant disk is typically ~ 100 ms and depends only weakly on the misalignment angle and the EOS. (iii) The ejecta velocity is typically ~ 0.2--0.3c and depends only weakly on i_tilt,0 and the EOS of the NS, while the morphology of the ejecta depends on its mass. (iv) The gravitational-wave spectra contains the information of the NS compactness in the cutoff frequency for i_tilt,0 < 60[deg.].

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamical mass ejection from black hole-neutron star binaries

We investigate properties of material ejected dynamically in the merger of black hole-neutron star binaries by numerical-relativity simulations. We systematically study the dependence of ejecta properties on the mass ratio of the binary, spin of the black hole, and equation of state of the neutron-star matter. Dynamical mass ejection is driven primarily by tidal torque, and the ejecta is much more anisotropic than that from binary neutron star mergers. In particular, the dynamical ejecta is concentrated around the orbital plane with a half opening angle of 10--20deg and often sweeps out only a half of the plane. The ejecta mass can be as large as ~0.1M_sun, and the velocity is subrelativistic with ~0.2--0.3c for typical cases. The ratio of the ejecta mass to the bound mass (disk and fallback components) is larger, and the ejecta velocity is larger, for larger values of the binary mass ratio, i.e., for larger values of the black-hole mass. The remnant black hole-disk system receives a kick velocity of O(100)km/s due to the ejecta linear momentum, and this easily dominates the kick velocity due to gravitational radiation. Structures of postmerger material, velocity distribution of the dynamical ejecta, fallback rates, and gravitational waves are also investigated. We also discuss the effect of ejecta anisotropy on electromagnetic counterparts, specifically a macronova/kilonova and synchrotron radio emission, developing analytic models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exploring tidal effects of coalescing binary neutron stars in numerical relativity II: Longterm simulations

We perform new longterm (15-16 orbits) simulations of coalescing binary neutron stars in numerical relativity using an updated Einstein's equation solver, employing low-eccentricity initial data, and modeling the neutron stars by a piecewise polytropic equation of state. A convergence study shows that our new results converge more rapidly than the third order and using the determined convergence order, we construct an extrapolated waveform for which the estimated total phase error should be less than 1 radian. We then compare the extrapolated waveforms with those calculated by the latest effective-one-body (EOB) formalism in which the so-called tidal deformability, higher post-Newtonian corrections, and gravitational self-force effects are taken into account. We show that for a binary of compact neutron stars with their radius 11.1 km, the waveform by the EOB formalism agrees quite well with the numerical waveform so that the total phase error is smaller than 1 radian for the total phase of $\sim 200$ radian up to the merger. By contrast, for a binary of less compact neutron stars with their radius 13.6 km, the EOB and numerical waveforms disagree with each other in the last few wave cycles, resulting in the total phase error of $\sim 3$ radian.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gravitational-wave cutoff frequencies of tidally disruptive neutron star-black hole binary mergers

Tidal disruption has a dramatic impact on the outcome of neutron star-black hole mergers. The phenomenology of these systems can be divided in three classes: nondisruptive, mildly disruptive or disruptive. The cutoff frequency of the gravitational radiation produced during the merger (which is potentially measurable by interferometric detectors) is very different in each regime, and when the merger is disuptive it carries information on the neutron star equation of state. Here we use semianalytical tools to derive a formula for the critical binary mass ratio $Q=M_{\rm BH}/M_{\rm NS}$ below which mergers are disruptive as a function of the stellar compactness $\mathcal{C}=M_{\rm NS}/R_{\rm NS}$ and the dimensionless black hole spin $χ$. We then employ a new gravitational waveform amplitude model, calibrated to $134$ general relativistic numerical simulations of binaries with black hole spin (anti-)aligned with the orbital angular momentum, to obtain a fit to the gravitational-wave cutoff frequency in the disruptive regime as a function of $\mathcal{C}$, $Q$ and $χ$. Our findings are important to build gravitational wave template banks, to determine whether neutron star-black hole mergers can emit electromagnetic radiation (thus helping multimessenger searches), and to improve event rate calculations for these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-resolution magnetohydrodynamics simulation of black hole-neutron star merger: Mass ejection and short gamma-ray burst

We report results of a high-resolution numerical-relativity simulation for the merger of black hole-magnetized neutron star binaries on Japanese supercomputer "K". We focus on a binary that is subject to tidal disruption and subsequent formation of a massive accretion torus. We find the launch of thermally driven torus wind, subsequent formation of a funnel wall above the torus and a magnetosphere with collimated poloidal magnetic field, and high Blandford-Znajek luminosity. We show for the first time this picture in a self-consistent simulation. The turbulence-like motion induced by the non-axisymmetric magnetorotational instability as well as the Kelvin-Helmholtz instability inside the accretion torus works as an agent to drive the mass accretion and converts the accretion energy to thermal energy, which results in the generation of a strong wind. By an in-depth resolution study, we reveal that high resolution is essential to draw such a picture. We also discuss the implication for the r-process nucleosynthesis, the radioactively-powered transient emission, and short gamma-ray bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The dynamical mass ejection from binary neutron star mergers: Radiation-hydrodynamics study in general relativity

We perform radiation-hydrodynamics simulations of binary neutron star mergers in numerical relativity on the Japanese "K" supercomputer, taking into account neutrino cooling and heating by an updated leakage-plus-transfer scheme for the first time. Neutron stars are modeled by three modern finite-temperature equations of state (EOS) developed by Hempel and his collaborators. We find that the properties of the dynamical ejecta of the merger such as total mass, average electron fraction, and thermal energy depend strongly on the EOS. Only for a soft EOS (the so-called SFHo), the ejecta mass exceeds $0.01M_{\odot}$. In this case, the distribution of the electron fraction of the ejecta becomes broad due to the shock heating during the merger. These properties are well-suited for the production of the solar-like $r$-process abundance. For the other stiff EOS (DD2 and TM1), for which a long-lived massive neutron star is formed after the merger, the ejecta mass is smaller than $0.01M_{\odot}$, although broad electron-fraction distributions are achieved by the positron capture and the neutrino heating.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effective No-Hair Relations for Neutron Stars and Quark Stars: Relativistic Results

Astrophysical charge-free black holes are known to satisfy no-hair relations through which all multipole moments can be specified in terms of just their mass and spin angular momentum. We here investigate the possible existence of no-hair-like relations among multipole moments for neutron stars and quark stars that are independent of their equation of state. We calculate the multipole moments of these stars up to hexadecapole order by constructing uniformly-rotating and unmagnetized stellar solutions to the Einstein equations. For slowly-rotating stars, we construct stellar solutions to quartic order in spin in a slow-rotation expansion, while for rapidly-rotating stars, we solve the Einstein equations numerically with the LORENE and RNS codes. We find that the multipole moments extracted from these numerical solutions are consistent with each other. We confirm that the current-dipole is related to the mass-quadrupole in an approximately equation of state independent fashion, which does not break for rapidly rotating neutron stars or quark stars. We further find that the current-octupole and the mass-hexadecapole moments are related to the mass-quadrupole in an approximately equation of state independent way to $\sim 10%$, worsening in the hexadecapole case. All of our findings are in good agreement with previous work that considered stellar solutions to leading-order in a weak-field expansion. The quartic in spin, slowly-rotating solutions found here allow us to estimate the systematic errors in the measurement of the neutron star's mass and radius with future X-ray observations, such as NICER and LOFT. We find that the effect of these quartic-in-spin terms on the quadrupole and hexadecapole moments and stellar eccentricity may dominate the error budget for very rapidly-rotating neutron stars. The new universal relations found here should help to reduce such systematic errors.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High resolution numerical-relativity simulations for the merger of binary magnetized neutron stars

We perform high-resolution magnetohydrodynamics simulations of binary neutron star mergers in numerical relativity on the Japanese supercomputer K. The neutron stars and merger remnants are covered by a grid spacing of 70\,m, which yields the highest-resolution results among those derived so far. By an in-depth resolution study, we clarify several amplification mechanisms of magnetic fields during the binary neutron star merger for the first time. First, the Kelvin-Helmholtz instability developed in the shear layer at the onset of the merger significantly amplifies the magnetic fields. A hypermassive neutron star (HMNS) formed after the merger is then subject to the nonaxisymmetric magnetorotational instability, which amplifies the magnetic field in the HMNS. These two amplification mechanisms cannot be found with insufficient-resolution runs. We also show that the HMNS eventually collapses to a black hole surrounded by an accretion torus which is strongly magnetized at birth.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Production of all the r-process nuclides in the dynamical ejecta of neutron star mergers

Recent studies suggest that binary neutron star (NS-NS) mergers robustly produce the heavy r-process nuclei above the atomic mass number A ~ 130 because of their ejecta consisting of almost pure neutrons (electron fraction of Y_e < 0.1). However, little production of the lighter r-process nuclei (A = 90-120) conflicts with the spectroscopic results of r-process-enhanced Galactic halo stars. We present, for the first time, the result of nucleosynthesis calculations based on the fully general-relativistic simulation of a NS-NS merger with approximate neutrino transport. It is found that the bulk of the dynamical ejecta are appreciably shock-heated and neutrino-processed, resulting in a wide range of Y_e (= 0.09-0.45). The mass-averaged abundance distribution of calculated nucleosynthesis yields is in reasonable agreement with the full-mass range (A = 90-240) of the solar r-process curve. This implies, if our model is representative of such events, that the dynamical ejecta of NS-NS mergers can be the origin of the Galactic r-process nuclei. Our result also shows that the radioactive heating after ~ 1 day from the merging, giving rise to r-process-powered transient emission, is dominated by the beta-decays of several species close to stability with precisely measured half-lives. This implies that the total radioactive heating rate for such an event can be well constrained within about a factor of two if the ejected material has a solar-like r-process pattern.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exploring tidal effects of coalescing binary neutron stars in numerical relativity

We study gravitational waves emitted in the late inspiral stage of binary neutron stars by analyzing the waveform obtained in numerical-relativity simulations. For deriving the physical gravitational waveforms from the numerical results, the resolution extrapolation plays an essential role for our simulations. The extrapolated gravitational-wave phases are compared with those calculated in the post-Newtonian (PN) and effective-one-body (EOB) formalisms including corrections of tidal effects. We show that the extrapolated gravitational-wave phases in numerical relativity agree well with those by the PN and EOB calculations for most of the inspiral stage except for a tidally-dominated, final inspiral stage, in which the PN and EOB results underestimate the tidal effects. Nevertheless, the accumulated phase difference between our extrapolated results and the results by the PN/EOB calculations is at most 1--3 radian in the last 15 cycles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Extracting equation of state parameters from black hole-neutron star mergers: aligned-spin black holes and a preliminary waveform model

Information about the neutron-star equation of state is encoded in the waveform of a black hole-neutron star system through tidal interactions and the possible tidal disruption of the neutron star. During the inspiral this information depends on the tidal deformability Lambda of the neutron star, and we find that Lambda is the best measured parameter during the merger and ringdown as well. We performed 134 simulations where we systematically varied the equation of state as well as the mass ratio, neutron star mass, and aligned spin of the black hole. Using these simulations we have developed an analytic representation of the full inspiral-merger-ringdown waveform calibrated to these numerical waveforms, and we use this analytic waveform to estimate the accuracy to which Lambda can be measured with gravitational-wave detectors. We find that although the inspiral tidal signal is small, coherently combining this signal with the merger-ringdown matter effect improves the measurability of Lambda by a factor of ~3 over using just the merger-ringdown matter effect alone. However, incorporating correlations between all the waveform parameters then decreases the measurability of Lambda by a factor of ~3. The uncertainty in Lambda increases with the mass ratio, but decreases as the black hole spin increases. Overall, a single Advanced LIGO detector can measure Lambda for mass ratios Q = 2--5, black hole spins J_BH/M_BH^2 = -0.5--0.75, neutron star masses M_NS = 1.2M_sun--1.45M_sun, and an optimally oriented distance of 100Mpc to an uncertainty of ~10%--100%. For the proposed Einstein Telescope, the uncertainty in Lambda is an order of magnitude smaller.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High-Energy Radiation from Remnants of Neutron Star Binary Mergers

We study high-energy emission from the mergers of neutron star binaries as electromagnetic counterparts to gravitational waves aside from short gamma-ray bursts. The mergers entail significant mass ejection, which interacts with the surrounding medium to produce similar but brighter remnants than supernova remnants in a few years. We show that electrons accelerated in the remnants can produce synchrotron radiation in X-rays detectable at $\sim 100$ Mpc by current generation telescopes and inverse Compton emission in gamma rays detectable by the \emph{Fermi} Large Area Telescopes and the Cherenkov Telescope Array under favorable conditions. The remnants may have already appeared in high-energy surveys such as the Monitor of All-sky X-ray Image and the \emph{Fermi} Large Area Telescope as unidentified sources. We also suggest that the merger remnants could be the origin of ultra-high-energy cosmic rays beyond the knee energy, $\sim 10^{15}$ eV, in the cosmic-ray spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Matter effects on binary neutron star waveforms

Using an extended set of equations of state and a multiple-group multiple-code collaborative effort to generate waveforms, we improve numerical-relativity-based data-analysis estimates of the measurability of matter effects in neutron-star binaries. We vary two parameters of a parameterized piecewise-polytropic equation of state (EOS) to analyze the measurability of EOS properties, via a parameter Λ that characterizes the quadrupole deformability of an isolated neutron star. We find that, to within the accuracy of the simulations, the departure of the waveform from point-particle (or spinless double black-hole binary) inspiral increases monotonically with Λ, and changes in the EOS that did not change Λ are not measurable. We estimate with two methods the minimal and expected measurability of Λ in second- and third- generation gravitational-wave detectors. The first estimate, using numerical waveforms alone, shows two EOS which vary in radius by 1.3km are distinguishable in mergers at 100Mpc. The second estimate relies on the construction of hybrid waveforms by matching to post-Newtonian inspiral, and estimates that the same EOS are distinguishable in mergers at 300Mpc. We calculate systematic errors arising from numerical uncertainties and hybrid construction, and we estimate the frequency at which such effects would interfere with template-based searches.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonspinning black hole-neutron star mergers: a model for the amplitude of gravitational waveforms

Black hole-neutron star binary mergers display a much richer phenomenology than black hole-black hole mergers, even in the relatively simple case - considered in this paper - in which both the black hole and the neutron star are nonspinning. When the neutron star is tidally disrupted, the gravitational wave emission is radically different from the black hole-black hole case and it can be broadly classified in two groups, depending on the spatial extent of the disrupted material. We present a phenomenological model for the gravitational waveform amplitude in the frequency domain that encompasses the three possible outcomes of the merger: no tidal disruption, "mild", and "strong" tidal disruption. The model is calibrated to general relativistic numerical simulations using piecewise polytropic neutron star equations of state. It should prove useful to extract information on the nuclear equation of state from future gravitational-wave observations, and also to obtain more accurate estimates of black hole-neutron star merger event rates in second- and third-generation interferometric gravitational-wave detectors. We plan to extend and improve the model as longer and more accurate gravitational waveforms become available, and we will make it publicly available online as a Mathematica package. We also present in appendix analytical fits of the projected KAGRA noise spectral density, that should be useful in data analysis applications.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Progenitor models of the electromagnetic transient associated with the short GRB 130603B

An electromagnetic transient powered by the radioactive decay of r-process elements, a so-called kilonova/macronova, is one of the possible observable consequences of compact binary mergers including at east one neutron star. Recent observations strongly suggest the first discovery of the electromagnetic transient, which is associated with the short GRB 130603B. We explore a possible progenitor of this event combining the numerical-relativity simulations and radiative transfer simulations of the dynamical ejecta of binary neutron star and black hole - neutron star mergers. We show that the ejecta models within a realistic parameter range consistently reproduce the observed near-infrared excess. We also show that the soft equation of state models for binary neutron star mergers and the stiff equation of state models for black hole - neutron star mergers are favored to reproduce the observed luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radioactively Powered Emission from Black Hole-Neutron Star Mergers

Detection of electromagnetic counterparts of gravitational wave (GW) sources is important to unveil the nature of compact binary coalescences. We perform three-dimensional, time-dependent, multi-frequency radiative transfer simulations for radioactively powered emission from the ejecta of black hole (BH) - neutron star (NS) mergers. Depending on the BH to NS mass ratio, spin of the BH, and equations of state of dense matter, BH-NS mergers can eject more material than NS-NS mergers. In such cases, radioactively powered emission from the BH-NS merger ejecta can be more luminous than that from NS-NS mergers. We show that, in spite of the expected larger distances to BH-NS merger events, observed brightness of BH-NS mergers can be comparable to or even higher than that of NS-NS mergers. We find that, when the tidally disrupted BH-NS merger ejecta are confined to a small solid angle, the emission from BH-NS merger ejecta tends to be bluer than that from NS-NS merger ejecta for a given total luminosity. Thanks to this property, we might be able to distinguish BH-NS merger events from NS-NS merger events by multi-band observations of the radioactively powered emission. In addition to the GW observations, such electromagnetic observations can potentially provide independent information on the nature of compact binary coalescences.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Remnant massive neutron stars of binary neutron star mergers: Evolution process and gravitational waveform

Massive (hypermassive and supramassive) neutron stars are likely to be often formed after the merger of binary neutron stars. We explore the evolution process of the remnant massive neutron stars and gravitational waves emitted by them, based on numerical-relativity simulations for binary neutron star mergers employing a variety of equations of state and choosing a plausible range of the neutron-star mass of binaries. We show that the lifetime of remnant hypermassive neutron stars depends strongly on the total binary mass and also on the equations of state. Gravitational waves emitted by the remnant massive neutron stars universally have a quasiperiodic nature of an approximately constant frequency although the frequency varies with time. We also show that the frequency and time-variation feature of gravitational waves depend strongly on the equations of state. We derive a fitting formula for the quasiperiodic gravitational waveforms, which may be used for the data analysis of a gravitational-wave signal.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Current Status of Numerical-Relativity Simulations in Kyoto

We describe the current status of our numerical simulations for the collapse of a massive stellar core to a BH and the BNS mergers, performed in the framework of full general relativity incorporating finite-temperature EOS and neutrino cooling. For the stellar core collapse simulation, we present the latest numerical results. We employed a purely nucleonic EOS (Shen-EOS). As an initial condition, we adopted a 100 $M_{\odot}$ presupernova model calculated by Umeda and Nomoto. Changing the degree of rotation for the initial condition, we clarify the strong dependence of the outcome of the collapse on this. When the rotation is rapid enough, the shock wave formed at the core bounce is deformed to be a torus-like shape. Then, the infalling matter is accumulated in the central region due to the oblique shock at the torus surface, hitting the PNS and dissipating the kinetic energy there. As a result, outflows can be launched. The PNS eventually collapses to a BH and an accretion torus is formed around it. We also found that the evolution of the BH and torus depends strongly on the rotation initially given. In the BNS merger simulations, we in addition employ an EOS incorporating a degree of freedom for hyperons. The numerical simulations show that for the purely nucleonic EOS, a HMNS with a long lifetime ($\gg 10$ ms) is the outcome for the total mass $M \lesssim 3.0M_{\odot}$. By contrast, the formed HMNS collapses to a BH in a shorter time scale with the hyperonic EOS for $M \gtrsim 2.7M_{\odot}$. It is shown that the typical total neutrino luminosity of the HMNS is $\sim 3$--$10\times 10^{53}$ ergs/s and the effective amplitude of gravitational waves from the HMNS is 2--$6 \times 10^{-22}$ at $f\approx 2$--2.5 kHz for a source distance of 100 Mpc.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Geometrical Aspects on Parameter estimation of stochastic gravitational wave background: beyond the Fisher analysis

The maximum likelihood method is often used for parameter estimation in gravitational wave astronomy. Recently, an interesting approach was proposed by Vallisneri to evaluate the distributions of parameter estimation errors expected for the method. This approach is to statistically analyze the local peaks of the likelihood surface, and works efficiently even for signals with low signal-to-noise ratios. Focusing special attention to geometric structure of the likelihood surface, we follow the proposed approach and derive formulae for a simplified model of data analysis where the target signal has only one intrinsic parameter, along with its overall amplitude. Then we apply our formulae to correlation analysis of stochastic gravitational wave background with a power-law spectrum. We report qualitative trends of the formulae using numerical results specifically obtained for correlation analysis with two Advanced-LIGO detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational waves from black hole-neutron star binaries I: Classification of waveforms

Using our new numerical-relativity code SACRA, long-term simulations for inspiral and merger of black hole (BH)-neutron star (NS) binaries are performed, focusing particularly on gravitational waveforms. As the initial conditions, BH-NS binaries in a quasiequilibrium state are prepared in a modified version of the moving-puncture approach. The BH is modeled by a nonspinning moving puncture and for the NS, a polytropic equation of state with $Γ=2$ and the irrotational velocity field are employed. The mass ratio of the BH to the NS, $Q=M_{\rm BH}/M_{\rm NS}$, is chosen in the range between 1.5 and 5. The compactness of the NS, defined by ${\cal C}=GM_{\rm NS}/c^2R_{\rm NS}$, is chosen to be between 0.145 and 0.178. For a large value of $Q$ for which the NS is not tidally disrupted and is simply swallowed by the BH, gravitational waves are characterized by inspiral, merger, and ringdown waveforms. In this case, the waveforms are qualitatively the same as that from BH-BH binaries. For a sufficiently small value of $Q \alt 2$, the NS may be tidally disrupted before it is swallowed by the BH. In this case, the amplitude of the merger and ringdown waveforms is very low, and thus, gravitational waves are characterized by the inspiral waveform and subsequent quick damping. The difference in the merger and ringdown waveforms is clearly reflected in the spectrum shape and in the "cut-off" frequency above which the spectrum amplitude steeply decreases. When an NS is not tidally disrupted (e.g., for Q=5), kick velocity, induced by asymmetric gravitational wave emission, agrees approximately with that derived for the merger of BH-BH binaries, whereas for the case that the tidal disruption occurs, the kick velocity is significantly suppressed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational waves, neutrino emissions, and effects of hyperons in binary neutron star mergers

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating both nucleonic and hyperonic finite-temperature equations of state (EOS) and neutrino cooling. It is found that for the nucleonic and hyperonic EOS, a hyper massive neutron star (HMNS) with a long lifetime $(t_{\rm life}\gtrsim 10 {\rm ms})$ is the outcome for the total mass $\approx 2.7 M_\odot$. For the total mass $\approx 3 M_\odot$, a long-lived (short-lived with $t_{\rm life}\approx 3 {\rm ms}$) HMNS is the outcome for the nucleonic (hyperonic) EOS. It is shown that the typical total neutrino luminosity of the HMNS is $\sim 3$ -- $6 \times 10^{53} {\rm erg /s}$ and the effective amplitude of gravitational waves from the HMNS is 1 -- $4\times 10^{-22}$ at $f\approx 2$ -- $3.2 {\rm kHz}$ for a source of distance of 100 Mpc. During the HMNS phase, characteristic frequencies of gravitational waves shift to a higher frequency for the hyperonic EOS in contrast to the nucleonic EOS in which they remain constant approximately. Our finding suggests that the effects of hyperons are well imprinted in gravitational wave and its detection will give us a potential opportunity to explore the composition of the neutron star matter. We present the neutrino luminosity curve when a black hole is formed as well.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Three dimensional evolution of differentially rotating magnetized neutron stars

We construct a new three-dimensional general relativistic magnetohydrodynamics code, in which a fixed mesh refinement technique is implemented. To ensure the divergence-free condition as well as the magnetic flux conservation, we employ the method by Balsara (2001). Using this new code, we evolve differentially rotating magnetized neutron stars, and find that a magnetically driven outflow is launched from the star exhibiting a kink instability. The matter ejection rate and Poynting flux are still consistent with our previous finding (Shibata et al., 2011) obtained in axisymmetric simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Binary Neutron Star Mergers: Dependence on the Nuclear Equation of State

We perform a numerical-relativity simulation for the merger of binary neutron stars with 6 nuclear-theory-based equations of state (EOSs) described by piecewise polytropes. Our purpose is to explore the dependence of the dynamical behavior of the binary neutron star merger and resulting gravitational waveforms on the EOS of the supernuclear-density matter. The numerical results show that the merger process and the first outcome are classified into three types; (i) a black hole is promptly formed, (ii) a short-lived hypermassive neutron star (HMNS) is formed, (iii) a long-lived HMNS is formed. The type of the merger depends strongly on the EOS and on the total mass of the binaries. For the EOS with which the maximum mass is larger than 2Msun, the lifetime of the HMNS is longer than 10 ms for a total mass m_0=2.7Msun. A recent radio observation suggests that the maximum mass of spherical neutron stars is M_max \geq 1.97\pm 0.04Msun in one σlevel. This fact and our results support the possible existence of a HMNS soon after the onset of the merger for a typical binary neutron star with m_0=2.7Msun. We also show that the torus mass surrounding the remnant black hole is correlated with the type of the merger process; the torus mass could be large, \geq 0.1Msun, in the case that a long-lived HMNS is formed. We also show that gravitational waves carry information of the merger process, the remnant, and the torus mass surrounding a black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Effects of hyperons in binary neutron star mergers

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating both nucleonic and hyperonic finite-temperature equations of state (EOS) and neutrino cooling for the first time. It is found that even for the hyperonic EOS, a hypermassive neutron star is first formed after the merger for the typical total mass $\approx$ 2.7M\bigodot, and subsequently collapses to a black hole (BH). It is shown that hyperons play a substantial role in the post-merger dynamics, torus formation around the BH, and emission of gravitational waves (GWs). In particular, the existence of hyperons is imprinted in GWs. Therefore, GW observations will provide a potential opportunity to explore the composition of the neutron star matter.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Extracting equation of state parameters from black hole-neutron star mergers. I. Nonspinning black holes

The late inspiral, merger, and ringdown of a black hole-neutron star (BHNS) system can provide information about the neutron-star equation of state (EOS). Candidate EOSs can be approximated by a parametrized piecewise-polytropic EOS above nuclear density, matched to a fixed low-density EOS; and we report results from a large set of BHNS inspiral simulations that systematically vary two parameters. To within the accuracy of the simulations, we find that, apart from the neutron-star mass, a single physical parameter Lambda, describing its deformability, can be extracted from the late inspiral, merger, and ringdown waveform. This parameter is related to the radius, mass, and l=2 Love number, k_2, of the neutron star by Lambda = 2k_2 R^5/3M_{NS}^5, and it is the same parameter that determines the departure from point-particle dynamics during the early inspiral. Observations of gravitational waves from BHNS inspiral thus restrict the EOS to a surface of constant Lambda in the parameter space, thickened by the measurement error. Using various configurations of a single Advanced LIGO detector, we find that neutron stars are distinguishable from black holes of the same mass and that Lambda^{1/5} or equivalently R can be extracted to 10-40% accuracy from single events for mass ratios of Q=2 and 3 at a distance of 100 Mpc, while with the proposed Einstein Telescope, EOS parameters can be extracted to accuracy an order of magnitude better.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational waves and neutrino emission from the merger of binary neutron stars

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating a finite-temperature (Shen's) equation of state (EOS) and neutrino cooling for the first time. It is found that for this stiff EOS, a hypermassive neutron star (HMNS) with a long lifetime ($\gg 10$ ms) is the outcome for the total mass $\alt 3.0M_{\odot}$. It is shown that the typical total neutrino luminosity of the HMNS is $\sim 3$--$8\times 10^{53}$ ergs/s and the effective amplitude of gravitational waves from the HMNS is 4--$6 \times 10^{-22}$ at $f=2.1$--2.5 kHz for a source distance of 100 Mpc. We also present the neutrino luminosity curve when a black hole is formed for the first time.

0 citations0 reviewsNo code linkedOpen access