We study the post-merger mass ejection of low-mass binary neutron stars (NSs) with the system mass of $2.5\, M_\odot$, and subsequent nucleosynthesis by performing general-relativistic, neutrino-radiation viscous-hydrodynamics simulations in axial symmetry. We find that the merger remnants are long-lived massive NSs surviving more than several seconds, irrespective of the nuclear equations of state (EOSs) adopted. The ejecta masses of our fiducial models are $\sim 0.06$-$0.1\, M_\odot$ (depending on the EOS), being $\sim 30\%$ of the initial disk masses ($\sim 0.15$-$0.3\, M_\odot$). Post-processing nucleosynthesis calculations indicate that the ejecta is composed mainly of light $r$-process nuclei with small amounts of lanthanides (mass fraction $\sim 0.002$-$0.004$) and heavier species due to the modest average electron fraction ($\sim 0.32$-$0.34$) for a reasonable value of the viscous coefficient. Such abundance distributions are incompatible with the solar $r$-process-like abundance patterns found in all measured $r$-process-enhanced metal-poor stars. Therefore, low-mass binary NS mergers should be rare. If such low-mass NS mergers occur, their electromagnetic counterparts, kilonovae, will be characterized by an early bright blue emission because of the large ejecta mass as well as the small lanthanide fraction. We also show, however, that if the effective turbulent viscosity is very high, or there is an efficient mass ejection working in the early post-merger phase, the electron fraction of the ejecta could be low enough that the solar $r$-process-like abundance pattern is reproduced and the lanthanide fraction becomes so high that the kilonova would be characterized by early bright blue and late bright red emissions.

We reanalyze gravitational waves from binary-neutron-star mergers GW170817 and GW190425 using a numerical-relativity (NR) calibrated waveform model, the TF2+_Kyoto model, which includes nonlinear tidal terms. For GW170817, by imposing a uniform prior on the binary tidal deformability $\tilde{\Lambda}$, the symmetric $90\%$ credible interval of $\tilde{\Lambda}$ is estimated to be $481^{+436}_{-359}$ ($402^{+465}_{-279}$) for the case of $f_\mathrm{max}=1000~\mathrm{Hz}$ ($2048~\mathrm{Hz}$), where $f_\mathrm{max}$ is the maximum frequency in the analysis. We also reanalyze the event with other waveform models: two post-Newtonian waveform models (TF2_PNTidal and TF2+_PNTidal), the TF2+_NRTidal model that is another NR calibrated waveform model, and its upgrade, the TF2+_NRTidalv2 model. While estimates of parameters other than $\tilde{\Lambda}$ are broadly consistent among various waveform models, our results indicate that estimates of $\tilde{\Lambda}$ depend on waveform models. However, the difference is smaller than the statistical error. For GW190425, we can only obtain little information on the binary tidal deformability. The systematic error among the NR calibrated waveform models will become significant to measure $\tilde{\Lambda}$ as the number of detectors and events increase and sensitivities of detectors are improved.

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.

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.

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.

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

We find that the Hanford and Livingston detectors of Advanced LIGO derive a distinct posterior probability distribution of binary tidal deformability tilde{Lambda} of the first binary-neutron-star merger GW170817. By analyzing public data of GW170817 with a nested-sampling engine and the default TaylorF2 waveform provided by the LALInference package, the probability distribution of the binary tidal deformability derived by the LIGO-Virgo detector network turns out to be determined dominantly by the Hanford detector. Specifically, by imposing the flat prior on tidal deformability of individual stars, symmetric 90% credible intervals of tilde{Lambda} are estimated to be 527^{+619}_{-345} with the Hanford detector, 927^{+522}_{-619} with the Livingston detector, and 455^{+668}_{-281} with the LIGO-Virgo detector network. Furthermore, the distribution derived by the Livingston detector changes irregularly when we vary the maximum frequency of the data used in the analysis. This feature is not observed for the Hanford detector. While they are all consistent, the discrepancy and irregular behavior suggest that an in-depth study of noise properties might improve our understanding of GW170817 and future events.

We revisit the lower bound on binary tidal deformability tilde{Lambda} imposed by a luminous kilonova/macronova, AT 2017gfo, by numerical-relativity simulations of models that are consistent with gravitational waves from the binary neutron star merger GW170817. Contrary to the claim made in the literature, we find that binaries with tilde{Lambda}<~400 can explain the luminosity of AT 2017gfo, as long as moderate mass ejection from the remnant is assumed as had been done in previous work. The reason is that the maximum mass of a neutron star is not strongly correlated with the tidal deformability of neutron stars with a typical mass of ~1.4M_sun. If the maximum mass is so large that the binary does not collapse into a black hole immediately after merger, the mass of the ejecta can be sufficiently large irrespective of the binary tidal deformability. We present models of binary mergers with tilde{Lambda} down to 242 that satisfy the requirement on the mass of the ejecta from the luminosity of AT 2017gfo. We further find that the luminosity of AT 2017gfo could be explained by models that do not experience bounce after merger. We conclude that the luminosity of AT 2017gfo is not very useful for constraining the binary tidal deformability. Accurate estimation of the mass ratio will be necessary to establish a lower bound using electromagnetic counterparts in the future. We also caution that merger simulations that employ a limited class of tabulated equations of state could be severely biased due to the lack of generality.

We discuss a multimessenger strategy to detect radio pulses from Galactic binary neutron stars in a very tight orbit with the period shorter than 10 min. On one hand, all-sky surveys by radio instruments are inefficient for detecting faint pulsars in very tight binaries due partly to the rarity of targets and primarily to the need of correction for severe Doppler smearing. On the other hand, the Laser Interferometer Space Antenna (LISA) will detect these binaries with a very large signal-to-noise ratio and determine the orbital frequency, binary parameters, and sky location to high accuracy. The information provided by LISA will reduce the number of required pointings by two to six orders of magnitude and that of required trials for the corrections by about nine orders of magnitude, increasing the chance of discovering radio pulsars. For making full use of this strategy, it is desirable to operate high-sensitivity radio instruments such as Square Kilometer Array Phase 2 simultaneously with LISA.

The Gravitational Wave (GW) event GW 170817 was generated by the coalescence of two neutron stars (NS) and produced an electromagnetic transient, labelled AT 2017gfo, that was target of a massive observational campaign. Polarimetry, a powerful diagnostic tool for probing the geometry and emission processes of unresolved sources, was obtained for this event. The observed linear polarization was consistent with being mostly induced by intervening dust, suggesting that the intrinsic emission was weakly polarized ($P < 0.4-0.5$ %). In this paper, we present and discuss a detailed analysis of the linear polarization expected from a merging NS binary system by means of 3D Monte Carlo radiative transfer simulations assuming a range of possible configurations, wavelengths, epochs and viewing angles. We find that polarization originates from the non-homogeneous opacity distribution within the ejecta and can reach levels of $P\sim1$ % at early times (1$-$2 days after the merger) and in the optical R band. Smaller polarization signals are expected at later epochs and/or different wavelengths. From the viewing-angle dependence of the polarimetric signal, we constrain the observer orientation of AT 2017gfo within $\sim$65$^\circ$ from the polar direction. The detection of non-zero polarization in future events will unambiguously reveal the presence of a lanthanide-free ejecta component and unveil its spatial and angular distribution.

We investigate the effect of inhomogeneity on primordial black hole formation in the matter dominated era. In the gravitational collapse of an inhomogeneous density distribution, a black hole forms if apparent horizon prevents information of the central region of the configuration from leaking. Since information cannot propagate faster than the speed of light, we identify the threshold of the black hole formation by considering the finite speed for propagation of information. We show that the production probability $\beta_{inhom}(\sigma)$ of primordial black holes, where $\sigma$ is density fluctuation at horizon entry, is significantly enhanced from that derived in previous work in which the speed of propagation was effectively regarded as infinite. For $\sigma \ll 1$, we obtain $\beta_{inhom}\simeq 3.70 \sigma^{3/2}$, which is larger by about an order of magnitude than the probability derived in earlier work by assuming instantaneous propagation of information.

We identify various contributors of systematic effects in the measurement of the neutron star (NS) tidal deformability and quantify their magnitude for several types of neutron star - black hole (NSBH) binaries. Gravitational waves from NSBH mergers contain information about the components' masses and spins as well as the NS equation of state. Extracting this information requires comparison of the signal in noisy detector data with theoretical templates derived from some combination of post-Newtonian (PN) approximants, effective one-body (EOB) models and %analytic fits to numerical relativity (NR) simulations. The accuracy of these templates is limited by errors in the NR simulations, by the approximate nature of the PN/EOB waveforms, and by the hybridization procedure used to combine them. In this paper, we estimate the impact of these errors by constructing and comparing a set of PN-NR hybrid waveforms, for the first time with NR waveforms from two different codes, namely, SpEC and SACRA, for such systems. We then attempt to recover the parameters of the binary using two non-precessing template approximants. We find that systematic errors are too large for tidal effects to be accurately characterized for any realistic NS equation of state model. We conclude that NSBH waveform models must be significantly improved if they are to be useful for the extraction of NS equation of state information or even for distinguishing NSBH systems from binary black holes.

We perform a general-relativistic magnetohydrodynamics simulation for $\approx 30$ ms after merger of a binary neutron star to a remnant massive neutron star (RMNS) with a high spatial resolution of the finest grid resolution $12.5$ m. First, we estimate that the Kelvin-Helmholtz instability at merger could amplify the magnetic-field energy up to $\sim 1\%$ of the thermal energy. Second, we find that the magnetorotational instability in the RMNS envelope and torus with $\rho < 10^{13}~{\rm g~cm^{-3}}$ sustains magneto-turbulent state and the effective viscous parameter in these regions is likely to converge to $\approx 0.01$--$0.02$ with respect to the grid resolution. Third, the current grid resolution is not still fine enough to sustain magneto-turbulent state in the RMNS with $\rho \ge 10^{13}~{\rm g~cm^{-3}}$.

We propose a long-term strategy for detecting thermal neutrinos from the remnant of binary-neutron-star mergers with a future M-ton water-Cherenkov detector such as Hyper-Kamiokande. Monitoring >~2500 mergers within <~200 Mpc, we may be able to detect a single neutrino with a human-time-scale operation of ~80 Mt years for the merger rate of 1 Mpc^{-3} Myr^{-1}, which is slightly lower than the median value derived by the LIGO-Virgo collaboration with GW170817. Although the number of neutrino events is minimal, contamination from other sources of neutrinos can be reduced efficiently to ~0.03 by analyzing only ~1 s after each merger identified with gravitational-wave detectors if Gadolinium is dissolved in the water. The contamination may be reduced further to ~0.01 if we allow the increase of waiting time by a factor of ~1.7. The detection of even a single neutrino can pin down the energy scale of thermal neutrino emission from binary-neutron-star mergers and could strongly support or disfavor formation of remnant massive neutron stars. Because the dispersion relation of gravitational waves are now securely constrained to that of massless particles with a corresponding limit on the graviton mass of <~10^{-22} eV/c^2 by binary-black-hole mergers, the time delay of a neutrino from gravitational waves can be used to put an upper limit of <~O(10) meV/c^2 on the absolute neutrino mass in the lightest eigenstate. Large neutrino detectors will enhance the detectability, and in particular, 5 Mt Deep-TITAND and 10 Mt MICA planned in the future will allow us to detect thermal neutrinos every ~16 and 8 years, respectively, increasing the significance.

We develop a model for frequency-domain gravitational waveforms from inspiraling binary neutron stars. Our waveform model is calibrated by comparison with hybrid waveforms constructed from our latest high-precision numerical-relativity waveforms and the SEOBNRv2T waveforms in the frequency range of $10$--$1000\,{\rm Hz}$. We show that the phase difference between our waveform model and the hybrid waveforms is always smaller than $0.1\, {\rm rad}$ for the binary tidal deformability, ${\tilde \Lambda}$, in the range $300\lesssim{\tilde \Lambda}\lesssim1900$ and for the mass ratio between 0.73 and 1. We show that, for $10$--$1000\,{\rm Hz}$, the distinguishability for the signal-to-noise ratio $\lesssim50$ and the mismatch between our waveform model and the hybrid waveforms are always smaller than 0.25 and $1.1\times10^{-5}$, respectively. The systematic error of our waveform model in the measurement of ${\tilde \Lambda}$ is always smaller than $20$ with respect to the hybrid waveforms for $300\lesssim{\tilde \Lambda}\lesssim1900$. The statistical error in the measurement of binary parameters is computed employing our waveform model, and we obtain results consistent with the previous studies. We show that the systematic error of our waveform model is always smaller than $20\%$ (typically smaller than $10\%$) of the statistical error for events with the signal-to-noise ratio of $50$.

We study the merger of black hole-neutron star binaries by fully general-relativistic neutrino-radiation-hydrodynamics simulations throughout the coalescence, particularly focusing on the role of neutrino irradiation in dynamical mass ejection. Neutrino transport is incorporated by an approximate transfer scheme based on the truncated moment formalism. While we fix the mass ratio of the black hole to the neutron star to be 4 and the dimensionless spin parameter of the black hole to be 0.75, the equations of state for finite-temperature neutron-star matter are varied. The hot accretion disk formed after tidal disruption of the neutron star emits a copious amount of neutrinos with the peak total luminosity ~1--3x10^53 erg s^(-1) via thermal pair production and subsequent electron/positron captures on free nucleons. Nevertheless, the neutrino irradiation does not modify significantly the electron fraction of the dynamical ejecta from the neutrinoless beta-equilibrium value at zero temperature of initial neutron stars. The mass of the wind component driven by neutrinos from the remnant disk is negligible compared to the very neutron-rich dynamical component, throughout our simulations performed until a few tens milliseconds after the onset of merger, for the models considered in this study. These facts suggest that the ejecta from black hole-neutron star binaries are very neutron rich and are expected to accommodate strong r-process nucleosynthesis, unless magnetic or viscous processes contribute substantially to the mass ejection from the disk. We also find that the peak neutrino luminosity does not necessarily increase as the disk mass increases, because tidal disruption of a compact neutron star can result in a remnant disk with a small mass but high temperature.

Gravitational-wave observation together with a large number of electromagnetic observations shows that the source of the latest gravitational-wave event, GW170817, detected primarily by advanced LIGO, is the merger of a binary neutron star. We attempt to interpret this observational event based on our results of numerical-relativity simulations performed so far paying particular attention to the optical and infra-red observations. We finally reach a conclusion that this event is described consistently by the presence of a long-lived hypermassive or supramassive neutron star as the merger remnant, because (i) significant contamination by lanthanide elements along our line of sight to this source can be avoided by the strong neutrino irradiation from it and (ii) it could play a crucial role to produce an ejecta component of appreciable mass with fast motion in the post-merger phase. We also point out that (I) the neutron-star equation of state has to be sufficiently stiff (i.e., the maximum mass of cold spherical neutron stars, M_max, has to be appreciably higher than 2M_sun in order that a long-lived massive neutron star can be formed as the merger remnant for the binary systems of GW170817, for which the initial total mass is >~ 2.73M_sun and (II) no detection of relativistic optical counterpart suggests a not-extremely high value of M_max approximately as 2.15-2.25M_sun.

Following the detection of the GW170817 signal and its associated electromagnetic emissions, we discuss the prospects of the local Hubble parameter measurement using double neutron stars (DNSs). The kilonova emissions of GW170817 are genuinely unique in terms of the rapid evolution of color and magnitude and we expect that, for a good fraction $\gtrsim 50\%$ of the DNS events within $\sim 200$Mpc, we could identify their host galaxies, using their kilonovae. At present, the estimated DNS merger rate $(1.5^{+3.2}_{-1.2})\times 10^{-6} {\rm Mpc^{-3} yr^{-1}}$ has a large uncertainty. But, if it is at the high end, we could measure the local Hubble parameter $H_L$ with the level of $\Delta H_L/H_L\sim 0.042$ ($1\sigma$ level), after the third observational run (O3). This accuracy is four times better than that obtained from GW170817 alone, and we will be able to examine the Hubble tension at $2.1\sigma$ level.

The location of the cutoff in the gravitational-wave spectrum of black hole-neutron star mergers is correlated strongly with the neutron-star radius for the case that the neutron star is disrupted by the black hole during the merger. However, the modulation which appears in the spectrum due to the mode mixing makes it difficult to measure the cutoff frequency if gravitational waves are observed from inclined direction or the binary is precessing. In this letter, we show that the cutoff frequency can be measured even in such situations with a method we have recently proposed to reconstruct the face-on waveforms only from the strain observed from a particular direction. We show that the systematic error in the measurement of the neutron-star radius can be reduced to $\lesssim 5\%$ for the case that tidal disruption of the neutron star occurs significantly.

Extending our previous studies, we perform high-resolution simulations of inspiraling binary neutron stars in numerical relativity. We thoroughly carry through a convergence study in our currently available computational resources with the smallest grid spacing of $\approx 63$--86~meter for the neutron-star radius 10.9--13.7\,km. The estimated total error in the gravitational-wave phase is of order 0.1~rad for the total phase of $\gtrsim 210$\,rad in the last $\sim 15$--16 inspiral orbits. We then compare the waveforms (without resolution extrapolation) with those calculated by the latest effective-one-body formalism (tidal SEOBv2 model referred to as TEOB model). We find that for any of our models of binary neutron stars, the waveforms calculated by the TEOB formalism agree with the numerical-relativity waveforms up to $\approx 3$\,ms before the peak of the gravitational-wave amplitude is reached: For this late inspiral stage, the total phase error is $\lesssim 0.1$\,rad. Although the gravitational waveforms have an inspiral-type feature for the last $\sim 3$\,ms, this stage cannot be well reproduced by the current TEOB formalism, in particular, for neutron stars with large tidal deformability (i.e., lager radius). The reason for this is described.

We propose that stellar-mass binary black holes like GW150914 will become a tool to explore the local Universe within ~100Mpc in the era of the Laser Interferometer Space Antenna (LISA). High calibration accuracy and annual motion of LISA could enable us to localize up to ~60 binaries more accurately than the error volume of ~100Mpc^3 without electromagnetic counterparts under moderately optimistic assumptions. This accuracy will give us a fair chance to determine the host object solely by gravitational waves. By combining the luminosity distance extracted from gravitational waves with the cosmological redshift determined from the host, the local value of the Hubble parameter will be determined up to a few % without relying on the empirically constructed distance ladder. Gravitational-wave cosmography would pave the way for resolution of the disputed Hubble tension, where the local and global measurements disagree in the value of the Hubble parameter at 3.4sigma level, which amounts to ~9%.

We discuss the gravitational wave emission and the orbital evolution of a hierarchical triple system composed of an inner binary black hole (BBH) and an outer tertiary. Depending on the kick velocity at the merger, the merged BBH could tidally disrupt the tertiary. Even though the fraction of BBH mergers accompanied by such disruptions is expected to be much smaller than unity, the existence of a tertiary and its basic parameters (e.g. semimajor axis, projected mass) can be examined for more than 1000 BBHs with the space GW detector LISA and its follow-on missions. This allows us to efficiently prescreen the targets for the follow-up searches for the tidal disruption events (TDEs). The TDE probability would be significantly higher for triple systems with aligned orbital- and spin-angular momenta, compared with random configurations.

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.

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

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

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.

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.

We explore magnetic-field amplification due to the Kelvin-Helmholtz instability during binary neutron star mergers. By performing high-resolution general relativistic magnetohydrodynamics simulations with a resolution of $17.5$ m for $4$--$5$ ms after the onset of the merger on the Japanese supercomputer "K", we find that an initial magnetic field of moderate maximum strength $10^{13}$ G is amplified at least by a factor of $\approx 10^3$. We also explore the saturation of the magnetic-field energy and our result shows that it is likely to be $\gtrsim 4 \times 10^{50}$ erg, which is $\gtrsim 0.1\%$ of the bulk kinetic energy of the merging binary neutron stars.

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

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.

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.

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.

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

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.

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.

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

We develop a method to compute low-eccentricity initial data of binary neutron stars required to perform realistic simulations in numerical relativity. The orbital eccentricity is controlled by adjusting the orbital angular velocity of a binary and incorporating an approaching relative velocity of the neutron stars. These modifications improve the solution primarily through the hydrostatic equilibrium equation for the binary initial data. The orbital angular velocity and approaching velocity of initial data are updated iteratively by performing time evolutions over ~3 orbits. We find that the eccentricity can be reduced by an order of magnitude compared to standard quasicircular initial data, specifically from ~0.01 to <~0.001, by three successive iterations for equal-mass binaries leaving ~10 orbits before the merger.

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.

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.

We study the possibility of pre-merger localization of eccentric compact binary coalescences by second-generation gravitational-wave detector networks. Gravitational waves from eccentric binaries can be regarded as a sequence of pulses, which are composed of various higher harmonic modes than ones with twice the orbital frequency. The higher harmonic modes from a very early inspiral phase will not only contribute to the signal-to-noise ratio, but also allow us to localize the gravitational-wave source before the merger sets in. This is due to the fact that high-frequency gravitational waves are essential for the source localization via triangulation by ground-based detector networks. We found that the single-detector signal-to-noise ratio exceeds 5 at 10 min before the merger for a 1.4--1.4 $M_\odot$ eccentric binary neutron stars at 100 Mpc in optimal cases, and it can be localized up to 10 deg$^2$ at half a minute before the merger by a four-detector network. We will even be able to achieve $\sim$ 10 deg$^2$ at 10 min before the merger for a face-on eccentric compact binary by a five-detector network.

We propose a possibility of ultrarelativistic electromagnetic counterparts to gravitational waves from binary neutron star mergers at nearly all the viewing angles. Our proposed mechanism relies on the merger-shock propagation accelerating a smaller mass in the outer parts of the neutron star crust to a larger Lorentz factor $\Gamma$ with smaller energy $\sim 10^{47} \Gamma^{-1}$ erg. This mechanism is difficult to resolve by current 3D numerical simulations. The outflows emit synchrotron flares for seconds to days by shocking the ambient medium. Ultrarelativistic flares shine at an early time and in high-energy bands, potentially detectable by current X-ray to radio instruments, such as Swift XRT and Pan-STARRS, and even in low ambient density $\sim 10^{-2}$ cm$^{-3}$ by EVLA. The flares probe the merger position and time, and the merger types as black hole--neutron star outflows would be non-/mildly relativistic.

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.

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.

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.

The merger of black hole-neutron star binaries can eject substantial material with the mass ~0.01-0.1M_sun when the neutron star is disrupted prior to the merger. The ejecta shows significant anisotropy, and travels in a particular direction with the bulk velocity ~0.2c. This is drastically different from the binary neutron star merger, for which ejecta is nearly isotropic. Anisotropic ejecta brings electromagnetic-counterpart diversity which is unique to black hole-neutron star binaries, such as viewing-angle dependence, polarization, and proper motion. The kick velocity of the black hole, gravitational-wave memory emission, and cosmic-ray acceleration are also discussed.

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

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.

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 {\Lambda} 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 {\Lambda}, and changes in the EOS that did not change {\Lambda} are not measurable. We estimate with two methods the minimal and expected measurability of {\Lambda} 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.

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.

Numerical-relativity simulations for the merger of binary neutron stars are performed for a variety of equations of state (EOSs) and for a plausible range of the neutron-star mass, focusing primarily on the properties of the material ejected from the system. We find that a fraction of the material is ejected as a mildly relativistic and mildly anisotropic outflow with the typical and maximum velocities $\sim 0.15$ -- $0.25c$ and $\sim 0.5$ -- $0.8c$ (where $c$ is the speed of light), respectively, and that the total ejected rest mass is in a wide range $10^{-4}$ -- $10^{-2}M_{\odot}$, which depends strongly on the EOS, the total mass, and the mass ratio. The total kinetic energy ejected is also in a wide range between $10^{49}$ and $10^{51} {\rm ergs}$. The numerical results suggest that for a binary of canonical total mass $2.7M_{\odot}$, the outflow could generate an electromagnetic signal observable by the planned telescopes through the production of heavy-element unstable nuclei via the $r$-process or through the formation of blast waves during the interaction with the interstellar matter, if the EOS and mass of the binary are favorable ones.

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.

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 $\Gamma=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.

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.

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

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.

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

We study the merger of black hole (BH)-neutron star (NS) binaries with a variety of BH spins aligned or anti-aligned with the orbital angular momentum, and with the mass ratio in the range MBH/MNS = 2--5, where MBH and MNS are the mass of the BH and NS, respectively. We model NS matter by systematically parametrized piecewise polytropic equations of state. The initial condition is computed in the puncture framework adopting an isolated horizon framework to estimate the BH spin and assuming an irrotational velocity field for the fluid inside the NS. Dynamical simulations are performed in full general relativity by an adaptive mesh refinement code, SACRA. The treatment of hydrodynamic equations and estimation of the disk mass are improved. We find that the NS is tidally disrupted irrespective of the mass ratio when the BH has a moderately large prograde spin, whereas only binaries with low mass ratios, MBH/MNS <~ 3 or small compactnesses of the NSs, bring the tidal disruption when the BH spin is zero or retrograde. The mass of the remnant disk is accordingly large as >~ 0.1 Msun, which is required by central engines of short gamma-ray bursts, if the BH spin is prograde. Information of the tidal disruption is reflected in a clear relation between the compactness of the NS and an appropriately defined "cutoff frequency" in the gravitational-wave spectrum, above which the spectrum damps exponentially. We find that the tidal disruption of the NS and excitation of the quasinormal mode of the remnant BH occur in a compatible manner in high mass-ratio binaries with the prograde BH spin. The correlation between the compactness and the cutoff frequency still holds for such cases. It is also suggested by extrapolation that the merger of an extremely spinning BH and an irrotational NS binary does not lead to formation of an overspinning BH.

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.

We report results of a numerical-relativity simulation for the merger of a black hole-neutron star binary with a variety of equations of state (EOSs) modeled by piecewise polytropes. We focus in particular on the dependence of the gravitational waveform at the merger stage on the EOSs. The initial conditions are computed in the moving-puncture framework, assuming that the black hole is nonspinning and the neutron star has an irrotational velocity field. For a small mass ratio of the binaries (e.g., MBH/MNS = 2 where MBH and MNS are the masses of the black hole and neutron star, respectively), the neutron star is tidally disrupted before it is swallowed by the black hole irrespective of the EOS. Especially for less-compact neutron stars, the tidal disruption occurs at a more distant orbit. The tidal disruption is reflected in a cutoff frequency of the gravitational-wave spectrum, above which the spectrum amplitude exponentially decreases. A clear relation is found between the cutoff frequency of the gravitational-wave spectrum and the compactness of the neutron star. This relation also depends weakly on the stiffness of the EOS in the core region of the neutron star, suggesting that not only the compactness but also the EOS at high density is reflected in gravitational waveforms. The mass of the disk formed after the merger shows a similar correlation with the EOS, whereas the spin of the remnant black hole depends primarily on the mass ratio of the binary, and only weakly on the EOS. Properties of the remnant disks are also analyzed.

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.

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 \sigma 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.

General relativistic quasiequilibrium states of black hole-neutron star binaries are computed in the moving-puncture framework. We propose three conditions for determining the quasiequilibrium states and compare the numerical results with those obtained in the excision framework. We find that the results obtained in the moving-puncture framework agree with those in the excision framework and with those in the third post-Newtonian approximation for the cases that (i) the mass ratio of the binary is close to unity irrespective of the orbital separation, and (ii) the orbital separation is large enough ($m_0\Omega \alt 0.02$ where $m_0$ and $\Omega$ are the total mass and the orbital angular velocity, respectively) irrespective of the mass ratio. For $m_0 \Omega \agt 0.03$, both of the results in the moving-puncture and excision frameworks deviate, more or less, from those in the third post-Newtonian approximation. Thus the numerical results do not provide a quasicircular state, rather they seem to have a nonnegligible eccentricity of order 0.01--0.1. We show by numerical simulation that a method in the moving-puncture framework can provide approximately quasicircular states in which the eccentricity is by a factor of $\sim 2$ smaller than those in quasiequilibrium given by other approaches.