Source author record

Masaru Shibata

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

ResearcherUnclaimed source record

gr-qc astro-ph.HE astro-ph.SR astro-ph astro-ph.CO hep-th hep-ph astro-ph.IM nucl-th

Catalog footprint

What is connected

80works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

A low-mass binary neutron star: long-term ejecta evolution and kilonovae with weak blue emission

We study the long-term evolution of ejecta formed in a binary neutron star (BNS) merger that results in a long-lived remnant NS by performing a hydrodynamics simulation with the outflow data of a numerical relativity simulation as the initial condition. At the homologously expanding phase, the total ejecta mass reaches $\approx0.1\,M_\odot$ with an average velocity of $\approx0.1\,c$ and lanthanide fraction of $\approx 0.005$. We further perform the radiative transfer simulation employing the obtained ejecta profile. We find that, contrary to a naive expectation from the large ejecta mass and low lanthanide fraction, the optical emission is not as bright as that in GW170817/AT2017gfo, while the infrared emission can be brighter. This light curve property is attributed to preferential diffusion of photons toward the equatorial direction due to the prolate ejecta morphology, large opacity contribution of Zr, Y, and lanthanides, and low specific heating rate of the ejecta. Our results suggest that these light curve features could be used as an indicator for the presence of a long-lived remnant NS. We also found that the bright optical emission broadly consistent with GW170817/AT2017gfo is realized for the case that the high-velocity ejecta components in the polar region are suppressed. These results suggest that the remnant in GW170817/AT2017gfo is unlikely to be a long-lived NS, but might have collapsed to a black hole within ${\cal O}(0.1)$ 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.

preprint2022arXiv

Electromagnetic counterparts of binary neutron star mergers leading to a strongly magnetized long-lived remnant neutron star

We explore the electromagnetic counterparts that will associate with binary neutron star mergers for the case that remnant massive neutron stars survive for $\gtrsim 0.5\,$s after the merger. For this study, we employ the outflow profiles obtained by long-term general-relativistic neutrino-radiation magneto-hydrodynamics simulations with a mean field dynamo effect. We show that a synchrotron afterglow with high luminosity can be associated with the merger event if the magnetic fields of the remnant neutron stars are significantly amplified by the dynamo effect. We also perform a radiative transfer calculation for kilonovae and find that for the highly amplified magnetic field cases, the kilonovae can be bright in the early epoch, while it shows the optical emission rapid declining in a few days and the long-lasting ($\sim 10\,{\rm d}$) emission very bright in the near-infrared wavelength. All these features have not been found in GW170817, indicating that the merger remnant neutron star formed in GW170817 might have collapsed to a black hole within several hundreds ms or magnetic-field amplification might be a minor effect.

preprint2022arXiv

Implementation of advanced Riemann solvers in a neutrino-radiation magnetohydrodynamics code in numerical relativity and its application to a binary neutron star merger

We implement advanced Riemann solvers HLLC and HLLD \cite{Mignone:2005ft,MUB:2009} together with an advanced constrained transport scheme \cite{Gardiner:2007nc} in a numerical-relativity neutrino-radiation magnetohydrodynamics code. We validate our implementation by performing a series of one- and multi-dimensional test problems for relativistic hydrodynamics and magnetohydrodynamics in both Minkowski spacetime and a static black hole spacetime. We find that the numerical solutions with the advanced Riemann solvers are more accurate than those with the HLLE solver \cite{DelZanna:2002rv}, which was originally implemented in our code. As an application to numerical relativity, we simulate an asymmetric binary neutron star merger leading to a short-lived massive neutron star both with and without magnetic fields. We find that the lifetime of the rotating massive neutron star formed after the merger and also the amount of the tidally-driven dynamical ejecta are overestimated when we employ the diffusive HLLE solver. We also find that the magnetorotational instability is less resolved when we employ the HLLE solver because of the solver's large numerical diffusivity. This causes a spurious enhancement both of magnetic winding resulting from large scale poloidal magnetic fields, and also of the energy of the outflow induced by magnetic pressure.

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.

preprint2020arXiv

Analytic properties of electromagnetic field of binary compact stars and electromagnetic precursors to gravitational waves

We analytically study the properties of the electromagnetic field in vacuum around close binary compact stars containing at least one neutron star. We show that the orbital motion of the neutron star induces high multipole modes of the electromagnetic field just before the merger. These modes are superimposed to form a spiral arm configuration and its edge is found to be a likely site for magnetic reconnection. These modes also enhance the total Poynting flux from neutron star binaries by a factor of 2--4. We also indicate that the electric field induced by the orbital motion lead to a magnetosphere around binaries and estimate its plasma density, which has a different parameter dependence than Goldreich-Julian density. With these properties, we discuss possible electromagnetic counterparts to gravitational wave events, and identify radio precursors, such as fast radio bursts, as the most promising observational targets.

preprint2020arXiv

Constraint on the ejecta mass for a black hole-neutron star merger event candidate S190814bv

We derive the upper limit to the ejecta mass of S190814bv, a black hole-neutron star merger candidate, through the radiative transfer simulations for kilonovae with the realistic ejecta density profile as well as the detailed opacity and heating rate models. The limits to the ejecta mass strongly depend on the viewing angle. For the face-on observations ($\le45^\circ$), the total ejecta mass should be smaller than $0.1\,M_\odot$ for the average distance of S190814bv ($D=267$ Mpc), while larger mass is allowed for the edge-on observations. We also derive the conservative upper limits of the dynamical ejecta mass to be $0.02\,M_\odot$, $0.03\,M_\odot$, and $0.05\,M_\odot$ for the viewing angle $\le 20^\circ$, $\le 50^\circ$, and for $\le 90^\circ$, respectively. We show that the {\it iz}-band observation deeper than $22$ mag within 2 d after the GW trigger is crucial to detect the kilonova with the total ejecta mass of $0.06\,M_\odot$ at the distance of $D=300$ Mpc. We also show that a strong constraint on the NS mass-radius relation can be obtained if the future observations put the upper limit of $0.03\,M_\odot$ to the dynamical ejecta mass for a BH-NS event with the chirp mass smaller than $\lesssim 3\,M_\odot$ and effective spin larger than $\gtrsim 0.5$.

preprint2020arXiv

Current status of space gravitational wave antenna DECIGO and B-DECIGO

Deci-hertz Interferometer Gravitational Wave Observatory (DECIGO) is the future Japanese space mission with a frequency band of 0.1 Hz to 10 Hz. DECIGO aims at the detection of primordial gravitational waves, which could be produced during the inflationary period right after the birth of the universe. There are many other scientific objectives of DECIGO, including the direct measurement of the acceleration of the expansion of the universe, and reliable and accurate predictions of the timing and locations of neutron star/black hole binary coalescences. DECIGO consists of four clusters of observatories placed in the heliocentric orbit. Each cluster consists of three spacecraft, which form three Fabry-Perot Michelson interferometers with an arm length of 1,000 km. Three clusters of DECIGO will be placed far from each other, and the fourth cluster will be placed in the same position as one of the three clusters to obtain the correlation signals for the detection of the primordial gravitational waves. We plan to launch B-DECIGO, which is a scientific pathfinder of DECIGO, before DECIGO in the 2030s to demonstrate the technologies required for DECIGO, as well as to obtain fruitful scientific results to further expand the multi-messenger astronomy.

preprint2020arXiv

Extreme mass ratio inspirals on the equatorial plane in the adiabatic order

We compute gravitational waves from inspiraling stellar-mass compact objects on the equatorial plane of a massive spinning black hole (BH). Our inspiral orbits are computed by taking into account the adiabatic change of orbital parameters due to gravitational radiation in the lowest order in mass ratio. We employ an interpolation method to compute the adiabatic change at arbitrary points inside the region of orbital parameter space computed in advance. Using the obtained inspiral orbits and associated gravitational waves, we compute power spectra of gravitational waves and the signal-to-noise ratio (SNR) for several values of the BH spin, the masses of the binary, and the initial orbital eccentricity during a hypothetical three-yrs LISA observation before final plunge. We find that (i) the SNR increases as the BH spin and the mass of the compact object increase for the BH mass $M \agt 10^6M_\odot$, (ii) the SNR has a maximum for $M \approx 10^6M_\odot$, and (iii) the SNR increases as the initial eccentricity increases for $M=10^6M_\odot$. We also show that incorporating the contribution from the higher multipole modes of gravitational waves is crucial for enhancing the detection rate.

preprint2020arXiv

Magnetosphere of an orbiting neutron star

We conduct force-free simulations of a single neutron star undergoing orbital motion in flat spacetime, mimicking the trajectory of the star about the center of mass on a compact binary system. Our attention is focused on the kinetic energy being extracted from the orbit by the acceleration of the magnetic dipole moment of the neutron star, and particularly, on how this energy gets distributed within its surrounding magnetosphere. A detailed study of the resulting magnetospheric configurations in our setting is presented, incorporating as well the effects due to neutron star spin and the misalignment of the magnetic and orbital axes. We find many features resembling those of pulsar magnetospheres for the orbiting neutron star --even in the absence of spin--, being of particular interest the development of a spiral current sheet that extends beyond the light cylinder. Then, we use recent advances in pulsar theory to estimate electromagnetic emissions produced at the reconnection regions of such current sheets.

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.

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.

preprint2020arXiv

Properties of neutrino transfer in a deformed remnant of neutron star merger

We study properties of neutrino transfer in a remnant of neutron star merger, consisting of a massive neutron star and a surrounding torus. We perform numerical simulations of the neutrino transfer by solving the Boltzmann equation with momentum-space angles and energies of neutrinos for snapshots of the merger remnant having elongated shapes. The evaluation of the neutrino distributions in the multi-dimensions enable us to provide the detailed information of angle and energy spectra and neutrino reaction rates. We demonstrate features of asymmetric neutrino fluxes from the deformed remnant and investigate the neutrino emission region by determining the neutrinosphere for each energy. We examine the emission and absorption of neutrinos to identify important ingredients of heating rates through neutrino irradiation. We show that the contributions of $μ$- and $τ$-types neutrinos are important for the heating in the region above the massive neutron star. We also examine the angle moments and the Eddington tensor calculated directly by the neutrino distribution functions and compare them with those obtained by a moment closure approach, which is often used in the study of neutrino-radiation hydrodynamics. We show that the components of the Eddington tensor have non-monotonic behaviors and the approximation of the closure relation may become inaccurate for high energy neutrinos, whose fluxes are highly aspherical due to the extended merger remnant.

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.

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.

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.

preprint2019arXiv

Diversity of kilonova light curves

We perform radiative transfer simulations for kilonova in various situations, including the cases of prompt collapse to a black hole from neutron-star mergers, high-velocity ejecta possibly accelerated by magnetars, and a black hole-neutron star merger. Our calculations are done employing ejecta profiles predicted by numerical-relativity simulations and a new line list for all the r-process elements. We found that (i) the optical emission for binary neutron stars promptly collapsing to a black hole would be fainter by $\gtrsim1$-$2\,{\rm mag}$ than that found in GW170817, while the infrared emission could be as bright as that in GW170817 if the post-merger ejecta is as massive as $\approx0.01\,M_\odot$; (ii) the kilonova would be brighter than that observed in GW170817 for the case that the ejecta is highly accelerated by the electromagnetic energy injection from the remnant, but it would decline rapidly and the magnitude would become fainter than in GW170817 within a few days; (iii) the optical emission from a black hole-neutron star merger ejecta could be as bright as that observed in GW170817 for the case that sufficiently large amount of matter is ejected ($\gtrsim0.02\,M_\odot$), while the infrared brightness would be brighter by $1$-$2\,{\rm mag}$ at the same time. We show that the difference in the ejecta properties would be imprinted in the differences in the peak brightness and time of peak. This indicates that we may be able to infer the type of the central engine for kilonovae by observation of the peak in the multiple band.

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.

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.

preprint2016arXiv

Gravitational waves from supermassive stars collapsing to a supermassive black hole

We derive the gravitational waveform from the collapse of a rapidly rotating supermassive star (SMS) core leading directly to a seed of a supermassive black hole (SMBH) in axisymmetric numerical-relativity simulations. We find that the peak strain amplitude of gravitational waves emitted during the black-hole formation is $\approx 5 \times 10^{-21}$ at the frequency $f \approx 5$\,mHz for an event at the cosmological redshift $z=3$, if the collapsing SMS core is in the hydrogen-burning phase. Such gravitational waves will be detectable by space laser interferometric detectors like eLISA with signal-to-noise ratio $\approx 10$, if the sensitivity is as high as LISA for $f=1$--10\,mHz. The detection of the gravitational-wave signal will provide a potential opportunity for testing the direct-collapse scenario for the formation of a seed of SMBHs.

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.

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.

preprint2016arXiv

Stability of rigidly rotating supermassive stars against gravitational collapse

We revisit secular stability against quasi-radial collapse for rigidly rotating supermassive stars (SMSs) in general relativity. We suppose that the SMSs are in a nuclear-burning phase and can be modeled by polytropic equations of state with the polytropic index $n_p$ slightly smaller than $3$. The stability is determined in terms of the turning-point method. We find a fitting formula of the stability condition for the plausible range of $n_p$ ($2.95 \alt n_p \alt 3$) for SMSs. This condition reconfirms that, while non-rotating SMSs with mass $\sim 10^5M_\odot$--$10^6M_\odot$ may undergo a general-relativistically induced quasi-radial collapse, rigidly rotating SMSs with a ratio of rotational to gravitational potential energy ($β$) of $\sim 10^{-2}$ are likely to be stable against collapse unless they are able to accrete $\sim 5$ times more mass during the (relatively brief) hydrogen-burning phase of their evolution. We discuss implications of our results.

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.

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.].

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.

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.

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.

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.

preprint2015arXiv

Neutrino-driven explosions of ultra-stripped type Ic supernovae generating binary neutron stars

We study explosion characteristics of ultra-stripped supernovae (SNe), which are candidates of SNe generating binary neutron stars (NSs). As a first step, we perform stellar evolutionary simulations of bare carbon-oxygen cores of mass from 1.45 to 2.0 $M_\odot$ until the iron cores become unstable and start collapsing. We then perform axisymmetric hydrodynamics simulations with spectral neutrino transport using these stellar evolution outcomes as initial conditions. All models exhibit successful explosions driven by neutrino heating. The diagnostic explosion energy, ejecta mass, Ni mass, and NS mass are typically $\sim 10^{50}$ erg, $\sim 0.1 M_\odot$, $\sim 0.01M_\odot$, and $\approx 1.3 M_\odot$, which are compatible with observations of rapidly-evolving and luminous transient such as SN 2005ek. We also find that the ultra-stripped SN is a candidate for producing the secondary low-mass NS in the observed compact binary NSs like PSR J0737-3039.

preprint2015arXiv

Quasiequilibrium sequences of binary neutron stars undergoing dynamical scalarization

We calculate quasiequilibrium sequences of equal-mass, irrotational binary neutron stars (BNSs) in a scalar-tensor (ST) theory of gravity that admits dynamical scalarization. We model neutron stars with realistic equations of state (notably through piecewise polytropic equations of state). Using these quasiequilibrium sequences we compute the binary's scalar charge and binding energy versus orbital angular frequency. We find that the absolute value of the binding energy is smaller than in general relativity (GR), differing at most by ~14% at high frequencies for the cases considered. We use the newly computed binding energy and the balance equation to estimate the number of gravitational-wave (GW) cycles during the adiabatic, quasicircular inspiral stage up to the end of the sequence, which is the last stable orbit or the mass-shedding point, depending on which comes first. We find that, depending on the ST parameters, the number of GW cycles can be substantially smaller than in GR. In particular, we obtain that when dynamical scalarization sets in around a GW frequency of ~130 Hz, the sole inclusion of the ST binding energy causes a reduction of GW cycles from ~120 Hz up to the end of the sequence (~1200 Hz) of ~11% with respect to the GR case. We estimate that when the ST energy flux is also included the reduction in GW cycles becomes of ~24%. Quite interestingly, dynamical scalarization can produce a difference in the number of GW cycles with respect to the GR point-particle case that is much larger than the effect due to tidal interactions, which is on the order of only a few GW cycles. These results further clarify and confirm recent studies that have evolved BNSs either in full numerical relativity or in post-Newtonian theory, and point out the importance of developing accurate ST-theory waveforms for systems composed of strongly self-gravitating objects, such as BNSs.

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.

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.

preprint2014arXiv

Higher dimensional Numerical Relativity: code comparison

The nonlinear behavior of higher dimensional black hole spacetimes is of interest in several contexts, ranging from an understanding of cosmic censorship to black hole production in high-energy collisions. However, nonlinear numerical evolutions of higher dimensional black hole spacetimes are tremendously complex, involving different diagnostic tools and "dimensional reduction methods". In this work we compare two different successful codes to evolve Einstein's equations in higher dimensions, and show that the results of such different procedures agree to numerical precision, when applied to the collision from rest of two equal-mass black holes. We calculate the total radiated energy to be E/M=9x10^{-4} in five dimensions and E/M=8.1x10^{-4} in six dimensions.

preprint2014arXiv

Jet collimation in the ejecta of double neutron star merger: new canonical picture of short gamma-ray bursts

The observations of jet breaks in the afterglows of short gamma-ray bursts (SGRBs) indicate that the jet has a small opening angle of < 10°. The collimation mechanism of the jet is a longstanding theoretical problem. We numerically analyze the jet propagation in the material ejected by double neutron star merger, and demonstrate that if the ejecta mass is > 10^{-2} M_{sun}, the jet is well confined by the cocoon and emerges from the ejecta with the required collimation angle. Our results also suggest that there are some populations of choked (failed) SGRBs or low-luminous new types of event. By constructing a model for SGRB 130603B, which is associated with the first kilonova/macronova can- didate, we infer that the equation-of-state of neutron stars would be soft enough to provide sufficient ejecta to collimate the jet, if this event was associated with a double neutron star merger.

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.

preprint2013arXiv

Coalescence of binary neutron stars in a scalar-tensor theory of gravity

We carry out numerical-relativity simulations of coalescing binary neutron stars in a scalar-tensor theory that admits spontaneous scalarization. We model neutron stars with realistic equations of state. We choose the free parameters of the theory taking into account the constraints imposed by the latest observations of neutron-star-- white-dwarf binaries with pulsar timing. We show that even within those severe constraints, scalarization can still affect the evolution of the binary neutron stars not only during the late inspiral, but also during the merger stage. We also confirm that even when both neutron stars have quite small scalar charge at large separations, they can be strongly scalarized dynamically during the final stages of the inspiral. In particular, we identify the binary parameters for which scalarization occurs either during the late inspiral or only after the onset of the merger when a remnant, supramassive or hypermassive neutron star is formed. We also discuss how those results can impact the extraction of physical information on gravitational waves once they are detected.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

preprint2012arXiv

NR/HEP: roadmap for the future

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

preprint2012arXiv

Stably stratified magnetized stars in general relativity

We construct magnetized stars composed of a fluid stably stratified by entropy gradients in the framework of general relativity, assuming ideal magnetohydrodynamics and employing a barotropic equation of state. We first revisit basic equations for describing stably-stratified stationary axisymmetric stars containing both poloidal and toroidal magnetic fields. As sample models, the magnetized stars considered by Ioka and Sasaki (2004), inside which the magnetic fields are confined, are modified to the ones stably stratified. The magnetized stars newly constructed in this study are believed to be more stable than the existing relativistic models because they have both poloidal and toroidal magnetic fields with comparable strength, and magnetic buoyancy instabilities near the surface of the star, which can be stabilized by the stratification, are suppressed.

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.

preprint2011arXiv

Afterglow of binary neutron star merger

The merger of two neutron stars results often in a rapidly and differentially rotating hypermassive neutron star (HMNS). We show by numerical-relativity simulation that the magnetic-field profile around such HMNS is dynamically varied during its subsequent evolution, and as a result, electromagnetic radiation with a large luminosity ~ 0.1 B^2 R^3 Omega is emitted with baryon (B, R, and Omega are poloidal magnetic-field strength at stellar surface, stellar radius, and angular velocity of a HMNS). The predicted luminosity of electromagnetic radiation, which is primarily emitted along the magnetic-dipole direction, is ~ 10^{47} (B/10^{13} G)^2(R/10 km)^3(Omega/10^4 rad/s) ergs/s, that is comparable to the luminosity of quasars.

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.

preprint2011arXiv

Effects of hyperons in binary neutron star mergers

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.

preprint2011arXiv

Formation of black hole and accretion disk in a massive high-entropy stellar core collapse

We present the first numerical result of fully general relativistic axisymmetric simulations for the collapse of a rotating high-entropy stellar core to a black hole and an accretion disk. The simulations are performed taking into account the relevant microphysics. We adopt as initial condition a spherical core with constant electron fraction ($Y_e = 0.5$) and entropy per baryon $s$ = 8 $k_B$, and angular velocity is superimposed. In the early phase, the core collapses in a homologous manner. Then, it experiences a weak bounce due to the gas pressure of free nucleons. Because the bounce is weak, the core collapses eventually to a black hole. Subsequent evolution depends on initial angular velocity. When the rotation is not fast, a geometrically thin (but optically thick) accretion disk is formed, and shock waves are formed in the inner part of the disk. For the moderately rotating case, the thin accretion disk expands eventually to be a geometrically thick torus after sufficient accumulation of the thermal energy generated at the shocks. Furthermore, convection occurs inside the torus. Neutrino luminosities vary violently with time because of the convective motion. For the rapidly rotating case, by contrast, a geometrically thick torus is formed soon after the black hole formation, and convective activity is weak due to the presence of epicyclic mode.

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.

preprint2011arXiv

Gravitational waves from the Papaloizou-Pringle instability in black hole-torus systems

Black hole (BH)--torus systems are promising candidates for the central engine of gamma-ray bursts (GRBs), and also possible outcomes of the collapse of supermassive stars to supermassive black holes (SMBHs). By three-dimensional general relativistic numerical simulations, we show that an $m=1$ nonaxisymmetric instability grows for a wide range of self-gravitating tori orbiting BHs. The resulting nonaxisymmetric structure persists for a timescale much longer than the dynamical one, becoming a strong emitter of large amplitude, quasiperiodic gravitational waves. Our results indicate that both, the central engine of GRBs and newly formed SMBHs, can be strong gravitational wave sources observable by forthcoming ground-based and spacecraft detectors.

preprint2011arXiv

Is super-Planckian physics visible? -- Scattering of black holes in 5 dimensions

It may be widely believed that probing short-distance physics is limited by the presence of the Planck energy scale above which scale any information is cloaked behind a horizon. If this hypothesis is correct, we could observe quantum behavior of gravity only through a black hole of Planck mass. We numerically show that in a scattering of two black holes in the 5-dimensional spacetime, a visible domain, whose curvature radius is much shorter than the Planck length, can be formed. Our result indicates that super-Planckian phenomena may be observed without an obstruction by horizon formation in particle accelerators.

preprint2011arXiv

Neutron star equation of state via gravitational wave observations

Gravitational wave observations can potentially measure properties of neutron star equations of state by measuring departures from the point-particle limit of the gravitational waveform produced in the late inspiral of a neutron star binary. Numerical simulations of inspiraling neutron star binaries computed for equations of state with varying stiffness are compared. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly if the neutron stars are more compact. This suggests that gravitational wave observations at frequencies around 1 kHz will be able to measure a compactness parameter and place stringent bounds on possible neutron star equations of state. Advanced laser interferometric gravitational wave observatories will be able to tune their frequency band to optimize sensitivity in the required frequency range to make sensitive measures of the late-inspiral phase of the coalescence.

preprint2011arXiv

Non-axisymmetric instabilities of neutron star with toroidal magnetic fields

The aim of this paper is to clarify the stabilities of neutron stars with strong toroidal magnetic fields against non-axisymmetric perturbation. The motivation comes from the fact that super magnetized neutron stars of $\sim 10^{15}$G, magnetars, and magnetized proto-neutron stars born after the magnetically-driven supernovae are likely to have such strong toroidal magnetic fields. Long-term, three-dimensional general relativistic magneto-hydrodynamic simulations are performed, preparing isentropic neutron stars with toroidal magnetic fields in equilibrium as initial conditions. To explore the effects of rotations on the stability, simulations are done for both non-rotating and rigidly rotating models. We find the emergence of the Parker and/or Tayler instabilities in both the non-rotating and rotating models. For both non-rotating and rotating models, the Parker instability is the primary instability as predicted by the local linear perturbation analysis. The interchange instability also appears in the rotating models. It is found that rapid rotation is not enough to suppress the Parker instability, and this finding does not agree with the perturbation analysis. The reason for this is that rigidly and rapidly rotating stars are marginally stable, and hence, in the presence of stellar pulsations by which the rotational profile is deformed, unstable regions with negative gradient of angular momentum profile is developed. After the onset of the instabilities, a turbulence is excited. Contrary to the axisymmetric case, the magnetic fields never reach an equilibrium state after the development of the turbulence. This conclusion suggests that three-dimensional simulation is indispensable for exploring the formation of magnetars or prominence activities of magnetars such as giant flares.

preprint2011arXiv

Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity

A truncated moment formalism for general relativistic radiation hydrodynamics, based on the Thorne's moment formalism, is derived. The fluid rest frame is chosen to be the fiducial frame for defining the radiation moments. Then, zeroth-, first-, and second-rank radiation moments are defined from the distribution function with a physically reasonable assumption for it in the optically thin and thick limits. The source terms are written, focusing specifically on the neutrino transfer and neglecting higher harmonic angular dependence of the reaction angle. Finally, basic equations for a truncated moment formalism for general relativistic radiation hydrodynamics in a closed covariant form are derived assuming a closure relation among the radiation stress tensor, energy density, and energy flux, and a variable Eddington factor, which works well.

preprint2010arXiv

Bar-mode instability of rapidly spinning black hole in higher dimensions: Numerical simulation in general relativity

Numerical-relativity simulation is performed for rapidly spinning black holes (BHs) in a higher-dimensional spacetime of special symmetries for the dimensionality $6 \leq d \leq 8$. We find that higher-dimensional BHs, spinning rapidly enough, are dynamically unstable against nonaxisymmetric bar-mode deformation and spontaneously emit gravitational waves, irrespective of $d$ as in the case $d=5$ \cite{SY09}. The critical values of a nondimensional spin parameter for the onset of the instability are $q:=a/μ^{1/(d-3)} \approx 0.74$ for $d=6$, $\approx 0.73$ for $d=7$, and $\approx 0.77$ for $d=8$ where $μ$ and $a$ are mass and spin parameters. Black holes with a spin smaller than these critical values ($q_{\rm crit}$) appear to be dynamically stable for any perturbation. Longterm simulations for the unstable BHs are also performed for $d=6$ and 7. We find that they spin down as a result of gravitational-wave emission and subsequently settle to a stable stationary BH of a spin smaller than $q_{\rm crit}$. For more rapidly spinning unstable BHs, the timescale, for which the new state is reached, is shorter and fraction of the spin-down is larger. Our findings imply that a highly rapidly spinning BH with $q > q_{\rm crit}$ cannot be a stationary product in the particle accelerators, even if it would be formed as a consequence of a TeV-gravity hypothesis. Its implications for the phenomenology of a mini BH are discussed.

preprint2010arXiv

Binary Neutron Stars in Quasi-equilibrium

Quasi-equilibrium sequences of binary neutron stars are constructed for a variety of equations of state in general relativity. Einstein's constraint equations in the Isenberg-Wilson-Mathews approximation are solved together with the relativistic equations of hydrostationary equilibrium under the assumption of irrotational flow. We focus on unequal-mass sequences as well as equal-mass sequences, and compare those results. We investigate the behavior of the binding energy and total angular momentum along a quasi-equilibrium sequence, the endpoint of sequences, and the orbital angular velocity as a function of time, changing the mass ratio, the total mass of the binary system, and the equation of state of a neutron star. It is found that the orbital angular velocity at the mass-shedding limit can be determined by an empirical formula derived from an analytic estimation. We also provide tables for 160 sequences which will be useful as a guideline of numerical simulations for the inspiral and merger performed in the near future.

preprint2010arXiv

Binary neutron-star mergers with Whisky and SACRA: First quantitative comparison of results from independent general-relativistic hydrodynamics codes

We present the first quantitative comparison of two independent general-relativistic hydrodynamics codes, the Whisky code and the SACRA code. We compare the output of simulations starting from the same initial data and carried out with the configuration (numerical methods, grid setup, resolution, gauges) which for each code has been found to give consistent and sufficiently accurate results, in particular in terms of cleanness of gravitational waveforms. We focus on the quantities that should be conserved during the evolution (rest mass, total mass energy, and total angular momentum) and on the gravitational-wave amplitude and frequency. We find that the results produced by the two codes agree at a reasonable level, with variations in the different quantities but always at better than about 10%.

preprint2010arXiv

Exploring binary-neutron-star-merger scenario of short-gamma-ray bursts by gravitational-wave observation

We elucidate the feature of gravitational waves (GWs) from binary neutron star merger collapsing to a black hole by general relativistic simulation. We show that GW spectrum imprints the coalescence dynamics, formation process of disk, equation of state for neutron stars, total masses, and mass ratio. A formation mechanism of the central engine of short $γ$-ray bursts, which are likely to be composed of a black hole and surrounding disk, therefore could be constrained by GW observation.

preprint2010arXiv

Influence of self-gravity on the runaway instability of black hole-torus systems

Results from the first fully general relativistic numerical simulations in axisymmetry of a system formed by a black hole surrounded by a self-gravitating torus in equilibrium are presented, aiming to assess the influence of the torus self-gravity on the onset of the runaway instability. We consider several models with varying torus-to-black hole mass ratio and angular momentum distribution orbiting in equilibrium around a non-rotating black hole. The tori are perturbed to induce the mass transfer towards the black hole. Our numerical simulations show that all models exhibit a persistent phase of axisymmetric oscillations around their equilibria for several dynamical timescales without the appearance of the runaway instability, indicating that the self-gravity of the torus does not play a critical role favoring the onset of the instability, at least during the first few dynamical timescales.

preprint2009arXiv

Maximal slicing of D-dimensional spherically-symmetric vacuum spacetime

We study the foliation of a $D$-dimensional spherically symmetric black-hole spacetime with $D\ge 5$ by two kinds of one-parameter family of maximal hypersurfaces: a reflection-symmetric foliation with respect to the wormhole slot and a stationary foliation that has an infinitely long trumpet-like shape. As in the four-dimensional case, the foliations by the maximal hypersurfaces have the singularity avoidance nature irrespective of dimensionality. This indicates that the maximal slicing condition will be useful for simulating higher-dimensional black-hole spacetimes in numerical relativity. For the case of D=5, we present analytic solutions of the intrinsic metric, the extrinsic curvature, the lapse function, and the shift vector for the foliation by the stationary maximal hypersurfaces. This data will be useful for checking five-dimensional numerical relativity codes based on the moving puncture approach.

preprint2009arXiv

Measuring the neutron star equation of state with gravitational wave observations

We report the results of a first study that uses numerical simulations to estimate the accuracy with which one can use gravitational wave observations of double neutron star inspiral to measure parameters of the neutron-star equation of state. The simulations use the evolution and initial-data codes of Shibata and Uryu to compute the last several orbits and the merger of neutron stars, with matter described by a parametrized equation of state. Previous work suggested the use of an effective cutoff frequency to place constraints on the equation of state. We find, however, that greater accuracy is obtained by measuring departures from the point-particle limit of the gravitational waveform produced during the late inspiral. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly for smaller values of the neutron star compactness (the ratio of the mass of the neutron star to its radius). We estimate that realistic equations of state will lead to gravitational waveforms that are distinguishable from point particle inspirals at an effective distance (the distance to an optimally oriented and located system that would produce an equivalent waveform amplitude) of 100 Mpc or less. As Lattimer and Prakash observed, neutron-star radius is closely tied to the pressure at density not far above nuclear. Our results suggest that broadband gravitational wave observations at frequencies between 500 and 1000 Hz will constrain this pressure, and we estimate the accuracy with which it can be measured. Related first estimates of radius measurability show that the radius can be determined to an accuracy of ~1 km at 100 Mpc.

preprint2009arXiv

Non-conformally flat initial data for binary compact objects

A new method is described for constructing initial data for a binary neutron-star (BNS) system in quasi-equilibrium circular orbit. Two formulations for non-conformally flat data, waveless (WL) and near-zone helically symmetric (NHS), are introduced; in each formulation, the Einstein-Euler system, written in 3+1 form on an asymptotically flat spacelike hypersurface, is exactly solved for all metric components, including the spatially non-conformally flat potentials, and for irrotational flow. A numerical method applicable to both formulations is explained with an emphasis on the imposition of a spatial gauge condition. Results are shown for solution sequences of irrotational BNS with matter approximated by parametrized equations of state that use a few segments of polytropic equations of state. The binding energy and total angular momentum of solution sequences computed within the conformally flat -- Isenberg-Wilson-Mathews (IWM) -- formulation are closer to those of the third post-Newtonian (3PN) two point particles up to the closest orbits, for the more compact stars, whereas sequences resulting from the WL/NHS formulations deviate from the 3PN curve even more for the sequences with larger compactness. We think it likely that this correction reflects an overestimation in the IWM formulation as well as in the 3PN formula, by $\sim 1$ cycle in the gravitational wave phase during the last several orbits. The work suggests that imposing spatial conformal flatness results in an underestimate of the quadrupole deformation of the components of binary neutron-star systems in the last few orbits prior to merger.

preprint2009arXiv

Nonaxisymmetric instability of rapidly rotating black hole in five dimensions

We present results from numerical solution of Einstein's equation in five dimensions describing evolution of rapidly rotating black holes. We show, for the first time, that the rapidly rotating black holes in higher dimensions are unstable against nonaxisymmetric deformation; for the five-dimensional case, the critical value of spin parameter for onset of the instability is $\approx 0.87$.

preprint2007arXiv

The Final Fate of Binary Neutron Stars: What Happens After the Merger?

The merger of two neutron stars usually produces a remnant with a mass significantly above the single (nonrotating) neutron star maximum mass. In some cases, the remnant will be stabilized against collapse by rapid, differential rotation. MHD-driven angular momentum transport eventually leads to the collapse of the remnant's core, resulting in a black hole surrounded by a massive accretion torus. Here we present simulations of this process. The plausibility of generating short duration gamma ray bursts through this scenario is discussed.

preprint2004arXiv

Deriving formulations for numerical computation of binary neutron stars in quasicircular orbits

Two relations, the virial relation $M_{\rm ADM}=M_{\rm K}$ and the first law in the form $δM_{\rm ADM}=ΩδJ$, should be satisfied by a solution and a sequence of solutions describing binary compact objects in quasiequilibrium circular orbits. Here, $M_{\rm ADM}$, $M_{\rm K}$, $J$, and $Ω$ are the ADM mass, Komar mass, angular momentum, and orbital angular velocity, respectively. $δ$ denotes an Eulerian variation. These two conditions restrict the allowed formulations that we may adopt. First, we derive relations between $M_{\rm ADM}$ and $M_{\rm K}$ and between $δM_{\rm ADM}$ and $ΩδJ$ for general asymptotically flat spacetimes. Then, to obtain solutions that satisfy the virial relation and sequences of solutions that satisfy the first law at least approximately, we propose a formulation for computation of quasiequilibrium binary neutron stars in general relativity. In contrast to previous approaches in which a part of the Einstein equation is solved, in the new formulation, the full Einstein equation is solved with maximal slicing and in a transverse gauge for the conformal three-metric. Helical symmetry is imposed in the near zone, while in the distant zone, a waveless condition is assumed. We expect the solutions obtained in this formulation to be excellent quasiequilibria as well as initial data for numerical simulations of binary neutron star mergers.

preprint2004arXiv

Gravitational waves from axisymmetric rotating stellar core collapse to a neutron star in full general relativity

Axisymmetric numerical simulations of rotating stellar core collapse to a neutron star are performed in the framework of full general relativity. The so-called Cartoon method, in which the Einstein field equations are solved in the Cartesian coordinates and the axisymmetric condition is imposed around the $y=0$ plane, is adopted. The hydrodynamic equations are solved in the cylindrical coordinates (on the $y=0$ plane in the Cartesian coordinates) using a high-resolution shock-capturing scheme with the maximum grid size $(2500,2500)$. A parametric equation of state is adopted to model collapsing stellar cores and neutron stars following Dimmelmeier et al. It is found that the evolution of central density during the collapse, bounce, and formation of protoneutron stars agree well with those in the work of Dimmelmeier et al. in which an approximate general relativistic formulation is adopted. This indicates that such approximation is appropriate for following axisymmetric stellar core collapses and subsequent formation of protoneutron stars. Gravitational waves are computed using a quadrupole formula. It is found that the waveforms are qualitatively in good agreement with those by Dimmelmeier et al. However, quantitatively, two waveforms do not agree well. Possible reasons for the disagreement are discussed.

preprint2004arXiv

Numerical evolution of secular bar-mode instability induced by the gravitational radiation reaction in rapidly rotating neutron stars

The evolution of a nonaxisymmetric bar-mode perturbation of rapidly rotating stars due to a secular instability induced by gravitational wave emission is studied in post-Newtonian simulations taking into account gravitational radiation reaction. A polytropic equation of state with the polytropic index $n=1$ is adopted. The ratio of the rotational kinetic energy to the gravitational potential energy $T/|W|$ is chosen in the range between 0.2 and 0.26. Numerical simulations were performed until the perturbation grows to the nonlinear regime, and illustrate that the outcome after the secular instability sets in is an ellipsoidal star of a moderately large ellipticity $\agt 0.7$. A rapidly rotating protoneutron star may form such an ellipsoid, which is a candidate for strong emitter of gravitational waves for ground-based laser interferometric detectors. A possibility that effects of magnetic fields neglected in this work may modify the growth of the secular instability is also mentioned.

preprint2003arXiv

Computing gravitational waves from slightly nonspherical stellar collapse to black hole: Odd-parity perturbation

Nonspherical stellar collapse to a black hole is one of the most promising gravitational wave sources for gravitational wave detectors. We numerically study gravitational waves from a slightly nonspherical stellar collapse to a black hole in linearized Einstein theory. We adopt a spherically collapsing star as the zeroth-order solution and gravitational waves are computed using perturbation theory on the spherical background. In this paper we focus on the perturbation of odd-parity modes. Using the polytropic equations of state with polytropic indices $n_p=1$ and 3, we qualitatively study gravitational waves emitted during the collapse of neutron stars and supermassive stars to black holes from a marginally stable equilibrium configuration. Since the matter perturbation profiles can be chosen arbitrarily, we provide a few types for them. For $n_p=1$, the gravitational waveforms are mainly characterized by a black hole quasinormal mode ringing, irrespective of perturbation profiles given initially. However, for $n_p=3$, the waveforms depend strongly on the initial perturbation profiles. In other words, the gravitational waveforms strongly depend on the stellar configuration and, in turn, on the ad hoc choice of the functional form of the perturbation in the case of supermassive stars.

preprint2003arXiv

Merger of binary neutron stars of unequal mass in full general relativity

We present results of three dimensional numerical simulations of the merger of unequal-mass binary neutron stars in full general relativity. A $Γ$-law equation of state $P=(Γ-1)ρε$ is adopted, where $P$, $ρ$, $\varep$, and $Γ$ are the pressure, rest mass density, specific internal energy, and the adiabatic constant, respectively. We take $Γ=2$ and the baryon rest-mass ratio $Q_M$ to be in the range 0.85--1. The typical grid size is $(633,633,317)$ for $(x,y,z)$ . We improve several implementations since the latest work. In the present code, the radiation reaction of gravitational waves is taken into account with a good accuracy. This fact enables us to follow the coalescence all the way from the late inspiral phase through the merger phase for which the transition is triggered by the radiation reaction. It is found that if the total rest-mass of the system is more than $\sim 1.7$ times of the maximum allowed rest-mass of spherical neutron stars, a black hole is formed after the merger irrespective of the mass ratios. The gravitational waveforms and outcomes in the merger of unequal-mass binaries are compared with those in equal-mass binaries. It is found that the disk mass around the so formed black holes increases with decreasing rest-mass ratios and decreases with increasing compactness of neutron stars. The merger process and the gravitational waveforms also depend strongly on the rest-mass ratios even for the range $Q_M= 0.85$--1.

preprint2002arXiv

Dynamical instability of differentially rotating stars

We study the dynamical instability against bar-mode deformation of differentially rotating stars. We performed numerical simulation and linear perturbation analysis adopting polytropic equations of state with the polytropic index $n=1$. It is found that rotating stars of a high degree of differential rotation are dynamically unstable even for the ratio of the kinetic energy to the gravitational potential energy of $O(0.01)$. Gravitational waves from the final nonaxisymmetric quasistationary states are calculated in the quadrupole formula. For rotating stars of mass $1.4M_{\odot}$ and radius several 10 km, gravitational waves have frequency several 100 Hz and effective amplitude $\sim 5 \times 10^{-22}$ at a distance of $\sim 100$ Mpc.

preprint2001arXiv

Thermodynamics of binary black holes and neutron stars

We consider compact binary systems, modeled in general relativity as vacuum or perfect-fluid spacetimes with a helical Killing vector k^α, heuristically, the generator of time-translations in a corotating frame. Systems that are stationary in this sense are not asymptotically flat, but have asymptotic behavior corresponding to equal amounts of ingoing and outgoing radiation. For black-hole binaries, a rigidity theorem implies that the Killing vector lies along the horizon's generators, and from this one can deduce the zeroth law (constant surface gravity of the horizon). Remarkably, although the mass and angular momentum of such a system are not defined, there is an exact first law, relating the change in the asymptotic Noether charge to the changes in the vorticity, baryon mass, and entropy of the fluid, and in the area of black holes. Binary systems with MΩsmall have an approximate asymptopia in which one can write the first law in terms of the asymptotic mass and angular momentum. Asymptotic flatness is precise in two classes of solutions used to model binary systems: spacetimes satisfying the post-Newtonian equations, and solutions to a modified set of field equations that have a spatially conformally flat metric. (The spatial conformal flatness formalism with helical symmetry, however, is consistent with maximal slicing only if replaces the extrinsic curvature in the field equations by an artificially tracefree expression in terms of the shift vector.) For these spacetimes, nearby equilibria whose stars have the same vorticity obey the relation δM = ΩδJ, from which one can obtain a turning point criterion that governs the stability of orbits.

preprint1999arXiv

On the Maximum Mass of Differentially Rotating Neutron Stars

We construct relativistic equilibrium models of differentially rotating neutron stars and show that they can support significantly more mass than their nonrotating or uniformly rotating counterparts. We dynamically evolve such ``hypermassive'' models in full general relativity and show that there do exist configurations which are dynamically stable against radial collapse and bar formation. Our results suggest that the remnant of binary neutron star coalescence may be temporarily stabilized by differential rotation, leading to delayed collapse and a delayed gravitational wave burst.

preprint1995arXiv

Hubble Parameter in Void Universe

We investigate the distance-redshift relation in the simple void model. As discussed by Moffat and Tatarski, if the observer stays at the center of the void, the observed Hubble parameter is not so different from the background Hubble parameter. However, if the position of observer is off center of the void, we must consider the peculiar velocity correction which is measured by the observed dipole anisotropy of cosmic microwave background. This peculiar velocity correction for the redshift is crucial to determine the Hubble parameter and we shall discuss this effect. Further the results of Turner et al by the N-body simulation will be also considered.

Masaru Shibata

What is connected

Connect this record

See the researcher in context

Building this map preview

80 published item(s)

A low-mass binary neutron star: long-term ejecta evolution and kilonovae with weak blue emission

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

Electromagnetic counterparts of binary neutron star mergers leading to a strongly magnetized long-lived remnant neutron star

Implementation of advanced Riemann solvers in a neutrino-radiation magnetohydrodynamics code in numerical relativity and its application to a binary neutron star merger

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

Analytic properties of electromagnetic field of binary compact stars and electromagnetic precursors to gravitational waves

Constraint on the ejecta mass for a black hole-neutron star merger event candidate S190814bv

Current status of space gravitational wave antenna DECIGO and B-DECIGO

Extreme mass ratio inspirals on the equatorial plane in the adiabatic order

Magnetosphere of an orbiting neutron star

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

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

Properties of neutrino transfer in a deformed remnant of neutron star merger

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

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

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

Diversity of kilonova light curves

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

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

Gravitational waves from supermassive stars collapsing to a supermassive black hole

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

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

Stability of rigidly rotating supermassive stars against gravitational collapse

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

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

Dynamical mass ejection from black hole-neutron star binaries

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

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

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

Neutrino-driven explosions of ultra-stripped type Ic supernovae generating binary neutron stars

Quasiequilibrium sequences of binary neutron stars undergoing dynamical scalarization

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

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

Higher dimensional Numerical Relativity: code comparison

Jet collimation in the ejecta of double neutron star merger: new canonical picture of short gamma-ray bursts

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

Coalescence of binary neutron stars in a scalar-tensor theory of gravity

Exploring tidal effects of coalescing binary neutron stars in numerical relativity

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

Matter effects on binary neutron star waveforms

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

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

Radioactively Powered Emission from Black Hole-Neutron Star Mergers

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

Current Status of Numerical-Relativity Simulations in Kyoto

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

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

NR/HEP: roadmap for the future

Stably stratified magnetized stars in general relativity

Three dimensional evolution of differentially rotating magnetized neutron stars

Afterglow of binary neutron star merger

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

Effects of hyperons in binary neutron star mergers

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

Formation of black hole and accretion disk in a massive high-entropy stellar core collapse

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

Gravitational waves from the Papaloizou-Pringle instability in black hole-torus systems

Is super-Planckian physics visible? -- Scattering of black holes in 5 dimensions

Neutron star equation of state via gravitational wave observations

Non-axisymmetric instabilities of neutron star with toroidal magnetic fields

Truncated Moment Formalism for Radiation Hydrodynamics in Numerical Relativity

Bar-mode instability of rapidly spinning black hole in higher dimensions: Numerical simulation in general relativity

Binary Neutron Stars in Quasi-equilibrium

Binary neutron-star mergers with Whisky and SACRA: First quantitative comparison of results from independent general-relativistic hydrodynamics codes

Exploring binary-neutron-star-merger scenario of short-gamma-ray bursts by gravitational-wave observation

Influence of self-gravity on the runaway instability of black hole-torus systems

Maximal slicing of D-dimensional spherically-symmetric vacuum spacetime

Measuring the neutron star equation of state with gravitational wave observations

Non-conformally flat initial data for binary compact objects

Nonaxisymmetric instability of rapidly rotating black hole in five dimensions

The Final Fate of Binary Neutron Stars: What Happens After the Merger?

Deriving formulations for numerical computation of binary neutron stars in quasicircular orbits

Gravitational waves from axisymmetric rotating stellar core collapse to a neutron star in full general relativity

Numerical evolution of secular bar-mode instability induced by the gravitational radiation reaction in rapidly rotating neutron stars