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.

The Simons Observatory (SO) is an upcoming cosmic microwave background (CMB) experiment located on Cerro Toco, Chile, that will map the microwave sky in temperature and polarization in six frequency bands spanning 27 to 285 GHz. SO will consist of one 6-meter Large Aperture Telescope (LAT) fielding $\sim$30,000 detectors and an array of three 0.42-meter Small Aperture Telescopes (SATs) fielding an additional 30,000 detectors. This synergy will allow for the extremely sensitive characterization of the CMB over angular scales ranging from an arcmin to tens of degrees, enabling a wide range of scientific output. Here we focus on the SATs targeting degree angular scales with successive dichroic instruments observing at Mid-Frequency (MF: 93/145 GHz), Ultra-High-Frequency (UHF: 225/285 GHz), and Low-Frequency (LF: 27/39 GHz). The three SATs will be able to map $\sim$10% of the sky to a noise level of 2 $\mu$K-arcmin when combining 93 and 145 GHz. The multiple frequency bands will allow the CMB to be separated from galactic foregrounds (primarily synchrotron and dust), with the primary science goal of characterizing the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma \left(r\right) \approx 0.003$.

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.

We established a production method of a good millimeter-wave absorber by using a 3D-printed mold. The mold has a periodic pyramid shape, and an absorptive material is filled into the mold. This shape reduces the surface reflection. The 3D-printed mold is made from a transparent material in the millimeter-wave range. Therefore, unmolding is not necessary. A significant benefit of this production method is easy prototyping with various shapes and various absorptive materials. We produced a test model and used a two-component epoxy encapsulant as the absorptive material. The test model achieved a low reflectance: $\sim 1\%$ at 100 GHz. The absorber is sometimes maintained at a low temperature condition for cases in which superconducting detectors are used. Therefore, cryogenic performance is required in terms of a mechanical strength for the thermal cycles, an adhesive strength, and a sufficient thermal conductivity. We confirmed the test-model strength by immersing the model into a liquid-nitrogen bath.

The gravitational wave event from the binary neutron star (BNS) merger GW170817 and the following multi-messenger observations present strong evidence for i) merger ejecta expanding with substantial velocities and ii) a relativistic jet which had to propagate through the merger ejecta. The ejecta's expansion velocity is not negligible for the jet head motion, which is a fundamental difference from the other systems like collapsars and active galactic nuclei. Here we present an analytic model of the jet propagation in an expanding medium. In particular, we notice a new term in the expression of the breakout time and velocity. In parallel, we perform a series of over a hundred 2D numerical simulations of jet propagation. The BNS merger ejecta is prepared based on numerical relativity simulations of a BNS merger with the highest-resolution to date. We show that our analytic results agree with numerical simulations over a wide parameter space. Then we apply our analytic model to GW170817, and obtain two solid constraints on: i) the central engine luminosity as $L_{iso,0} \sim 3\times10^{49}-2.5\times10^{52}$ erg s$^{-1}$, and on ii) the delay time between the merger and engine activation $t_0-t_m < 1.3$ s. The engine power implies that the apparently-faint \textit{short} gamma-ray burst (\textit{s}GRB) \textit{s}GRB 170817A is similar to typical \textit{s}GRBs if observed on-axis.

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.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

We revisit the constraint on the maximum mass of cold spherical neutron stars coming from the observational results of GW170817. We develop a new framework for the analysis by employing both energy and angular momentum conservation laws as well as solid results of latest numerical-relativity simulations and of neutron stars in equilibrium. The new analysis shows that the maximum mass of cold spherical neutron stars can be only weakly constrained as $M_{\rm max} \alt 2.3M_\odot$. Our present result illustrates that the merger remnant neutron star at the onset of collapse to a black hole is not necessarily rapidly rotating and shows that we have to take into account the angular momentum conservation law to impose the constraint on the maximum mass of neutron stars.

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.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

The observed nonthermal afterglow spectrum of the binary neutron star (BNS) merger GW170817 from radio to X-ray are consistent with synchrotron radiation by shock-accelerated electrons. However, previous afterglow modeling studies were based on a simplified assumption that the acceleration efficiency is extremely high, i.e. all electrons in the shock are accelerated as a nonthermal population. This affects the estimate of the minimum electron energy and hence $\nu_m$, the peak frequency of the afterglow spectrum. Here we present Bayesian fitting to the observed data with a more natural electron energy distribution, in which the acceleration efficiency is a free parameter. Interestingly, the maximum likelihood solutions are found with radio flux below $\nu_m$ in the early phase, in contrast to previous studies that found the radio frequency always above $\nu_m$. Therefore the $\nu_m$ passage through the radio band could have been clearly detected for GW170817, if sufficient low-frequency radio data had been taken in early time. In the new solutions, the lowest energy of electrons is found close to equipartition with the post shock protons, but only a small fraction ($<$10\%) of electrons are accelerated as nonthermal particles. The jet energy and interstellar medium density are increased by 1--2 orders of magnitude from the conventional modeling, though these are still consistent with other constraints. We encourage to take densely sampled low-frequency radio data in the early phase for future BNS merger events, which would potentially detect $\nu_m$ passage and give a strong constraint on electron energy distribution and particle acceleration efficiency.

We perform long-term general relativistic neutrino radiation hydrodynamics simulations (in axisymmetry) for a massive neutron star (MNS) surrounded by a torus, which is a canonical remnant formed after the binary neutron star merger. We take into account the effects of viscosity, which is likely to arise in the merger remnant due to magnetohydrodynamical turbulence. As the initial condition, we employ the azimuthally averaged data of the MNS-torus system derived in a three-dimensional, numerical-relativity simulation for the binary neutron star merger. The viscous effect plays key roles for the remnant evolution and mass ejection from it in two phases of the evolution. In the first $t\lesssim10$ ms, a differential rotation state of the MNS is changed to a rigidly rotating state, and as a result, a sound wave, which subsequently becomes a shock wave, is formed in the vicinity of the MNS due to the variation of the quasi-equilibrium state of the MNS. The shock wave induces significant mass ejection of mass $\sim(0.5-2.0)\times 10^{-2}M_\odot$ for the alpha viscosity parameter of $0.01-0.04$. For the longer-term evolution with $\sim 0.1-10$ s, a significant fraction of the torus material is ejected. The ejecta mass is likely to be of order $10^{-2}M_\odot$, so that the total mass of the viscosity-driven ejecta could dominate that of the dynamical ejecta of mass $\lesssim 10^{-2}M_\odot$. The electron fraction, $Y_e$, of the ejecta is always high enough ($Y_e\gtrsim0.25$) that this post-merger ejecta is lanthanide-poor; hence, the opacity of the ejecta is likely to be $\sim 10-100$ times lower than that of the dynamical ejecta. This indicates that the electromagnetic signal from the ejecta would be rapidly evolving, bright, and blue if it is observed from a small viewing angle ($\lesssim 45^\circ$) for which the effect of the dynamical ejecta is minor.

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 find, using high resolution numerical relativistic simulations, that the tail of the dynamical ejecta of neutron star mergers extends to mildly relativistic velocities faster than $0.7c$. The kinetic energy of this fast tail is $\sim 10^{47}$--$10^{49}$ erg, depending on the neutron star equation of state and on the binary masses. The synchrotron flare arising from the interaction of this fast tail with the surrounding ISM can power the observed non-thermal emission that followed GW170817, provided that the ISM density is $\sim 10^{-2}\,{\rm cm^{-3}}$, the two neutron stars had roughly equal masses and the neutron star equation of state is soft (small neutron star radii). One of the generic predictions of this scenario is that the cooling frequency crosses the X-ray band on a time scale of a few months to a year, leading to a cooling break in the X-ray light curve. If this dynamical ejecta scenario is correct, we expect that the synchrotron radio flare from the ejecta that have produced the macronova/kilonova emission will be observable on time scales of $10^3$ to $10^5$ days. Further multi-frequency observations will confirm or rule out this dynamical ejecta scenario.

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

Most of fast radio bursts (FRB) do not show evidence for repetition, and such non-repeating FRBs may be produced at the time of a merger of binary neutron stars (BNS), provided that the BNS merger rate is close to the high end of the currently possible range. However, the merger environment is polluted by dynamical ejecta, which may prohibit the radio signal to propagate. We examine this by using a general-relativistic simulation of a BNS merger, and show that the ejecta appears about 1 ms after the rotation speed of the merged star becomes the maximum. Therefore there is a time window in which an FRB signal can reach outside, and the short duration of non-repeating FRBs can be explained by screening after ejecta formation. A fraction of BNS mergers may leave a rapidly rotating and stable neutron star, and such objects may be the origin of repeating FRBs like FRB 121102. We show that a merger remnant would appear as a repeating FRB in a time scale of about 1-10 yrs, and expected properties are consistent with the observations of FRB 121102. We construct an FRB rate evolution model including these two populations of repeating and non-repeating FRBs from BNS mergers, and show that the detection rate of repeating FRBs relative to non-repeating ones rapidly increases with improving search sensitivity. This may explain that the only repeating FRB 121102 was discovered by the most sensitive FRB search with Arecibo. Several predictions are made, including appearance of a repeating FRB 1-10 years after a BNS merger that is localized by gravitational wave and subsequent electromagnetic radiation.

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.

We investigate current and future prospects for coincident detection of high-energy neutrinos and gravitational waves (GWs). Short gamma-ray bursts (SGRBs) are believed to originate from mergers of compact star binaries involving neutron stars. We estimate high-energy neutrino fluences from prompt emission, extended emission, X-ray flares, and plateau emission, and show that neutrino signals associated with the extended emission are the most promising. Assuming that the cosmic-ray loading factor is $\sim10$ and the Lorentz factor distribution is lognormal, we calculate the probability of neutrino detection from extended emission by current and future neutrino detectors, and find that the quasi-simultaneous detection of high-energy neutrinos, gamma rays, and GWs is possible with future instruments or even with current instruments for nearby SGRBs having extended emission. We also discuss stacking analyses that will also be useful with future experiments such as IceCube-Gen2.

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 performed general relativistic, long-term, axisymmetric neutrino radiation hydrodynamics simulations for the remnant formed after the binary neutron star merger, which consist of a massive neutron star and a torus surrounding it. As an initial condition, we employ the result derived in a three-dimensional, numerical relativity simulation for the binary neutron star merger. We investigate the properties of neutrino-driven ejecta. Due to the pair-annihilation heating, the dynamics of the neutrino-driven ejecta is significantly modified. The kinetic energy of the ejecta is about two times larger than that in the absence of the pair-annihilation heating. This suggests that the pair-annihilation heating plays an important role in the evolution of the merger remnants. The relativistic outflow, which is required for driving gamma-ray bursts, is not observed because the specific heating rate around the rotational axis is not sufficiently high due to the baryon loading caused by the neutrino-driven ejecta from the massive neutron star. We discuss the condition for launching the relativistic outflow and the nucleosynthesis in the ejecta.

Employing a simplified version of the Israel-Stewart formalism of general-relativistic shear-viscous hydrodynamics, we explore the evolution of a remnant massive neutron star of binary neutron star merger and pay special attention to the resulting gravitational waveforms. We find that for the plausible values of the so-called viscous alpha parameter of the order $10^{-2}$, the degree of the differential rotation in the remnant massive neutron star is significantly reduced in the viscous timescale, $\alt 5$\,ms. Associated with this, the degree of non-axisymmetric deformation is also reduced quickly, and as a consequence, the amplitude of quasi-periodic gravitational waves emitted also decays in the viscous timescale. Our results indicate that for modeling the evolution of the merger remnants of binary neutron stars, we would have to take into account magnetohydrodynamics effects, which in nature could provide the viscous effects.

Employing a simplified version of the Israel-Stewart formalism for general-relativistic shear-viscous hydrodynamics, we perform axisymmetric general-relativistic simulations for a rotating neutron star surrounded by a massive torus, which can be formed from differentially rotating stars. We show that with our choice of a shear-viscous hydrodynamics formalism, the simulations can be stably performed for a long time scale. We also demonstrate that with a possibly high shear-viscous coefficient, not only viscous angular momentum transport works but also an outflow could be driven from a hot envelope around the neutron star for a time scale $\gtrsim 100\,$ms with the ejecta mass $\gtrsim 10^{-2}M_\odot$ which is comparable to the typical mass for dynamical ejecta of binary neutron star mergers. This suggests that massive neutron stars surrounded by a massive torus, which are typical outcomes formed after the merger of binary neutron stars, could be the dominant source for providing neutron-rich ejecta, if the effective shear viscosity is sufficiently high, i.e., if the viscous $\alpha$ parameter is $\gtrsim 10^{-2}$. The present numerical result indicates the importance of a future high-resolution magnetohydrodynamics simulation that is the unique approach to clarify the viscous effect in the merger remnants of binary neutron stars by the first-principle manner.

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.

Fast-spinning strongly magnetized newborn neutron stars, including nascent magnetars, are popularly implemented as the engine of luminous stellar explosions. Here, we consider the scenario that they power various stripped-envelope supernovae, not only super-luminous supernovae Ic but also broad-line supernovae Ibc and possibly some ordinary supernovae Ibc. This scenario is also motivated by the hypothesis that Galactic magnetars largely originate from fast-spinning neutron stars as remnants of stripped-envelope supernovae. By consistently modeling the energy injection from magnetized wind and Ni decay, we show that proto-neutron stars with >~ 10 ms rotation and B_dip >~ 5 x 10^14 G can be harbored in ordinary supernovae Ibc. On the other hand, millisecond proto-neuton stars can solely power broad-line supernovae Ibc if they are born with poloidal magnetic field of B_dip >~ 5 x 10^14 G, and superluminous supernovae Ic with B_dip >~ 10^13 G. Then, we study how multi-messenger emission can be used to discriminate such pulsar-driven supernova models from other competitive scenarios. First, high-energy x-ray and gamma-ray emission from embryonic pulsar wind nebulae is a promising smoking gun of the underlying newborn pulsar wind. Follow-up observations of stripped-envelope supernovae using NuSTAR ~ 50-100 days after the explosion is strongly encouraged for nearby objects. We also discuss possible effects of gravitational-waves on the spin-down of proto-neutron stars. If millisecond proto-neutron stars with B_dip <~ a few x 10^13 G emit gravitational waves through e.g., non-axisymmetric rotation deformed by the inner toroidal fields of B_t >~ 10^16 G, the gravitational wave signal can be detectable from ordinary supernova Ibc in the Virgo cluster by Advanced LIGO, Advanced Virgo, and KAGRA.

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.

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.

It has been suggested that some classes of luminous supernovae (SNe) and gamma-ray bursts (GRBs) are driven by newborn magnetars. Fast-rotating proto-neutron stars have also been of interest as potential sources of gravitational waves (GWs). We show that for a range of rotation periods and magnetic fields, hard X rays and GeV gamma rays provide us with a promising probe of pulsar-aided SNe. It is observationally known that young pulsar wind nebulae (PWNe) in the Milky Way are very efficient lepton accelerators. We argue that, if embryonic PWNe satisfy similar conditions at early stages of SNe (in ~1-10 months after the explosion), external inverse-Compton emission via upscatterings of SN photons is naturally expected in the GeV range as well as broadband synchrotron emission. To fully take into account the Klein-Nishina effect and two-photon annihilation process that are important at early times, we perform detailed calculations including electromagnetic cascades. Our results suggest that hard X-ray telescopes such as NuSTAR can observe such early PWN emission by followup observations in months-to-years. GeV gamma rays may also be detected by Fermi for nearby SNe, which serve as counterparts of these GW sources. Detecting the signals will give us an interesting probe of particle acceleration at early times of PWNe, as well as clues to driving mechanisms of luminous SNe and GRBs. Since the Bethe-Heitler cross section is lower than the Thomson cross section, gamma rays would allow us to study subphotospheric dissipation. We encourage searches for high-energy emission from nearby SNe, especially Type Ibc SNe including super-luminous objects.

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

The merger of two neutron stars leaves behind a rapidly spinning hypermassive object whose survival is believed to depend on the maximum mass supported by the nuclear equation of state, angular momentum redistribution by (magneto-)rotational instabilities, and spindown by gravitational waves. The high temperatures (~5-40 MeV) prevailing in the merger remnant may provide thermal pressure support that could increase its maximum mass and, thus, its life on a neutrino-cooling timescale. We investigate the role of thermal pressure support in hypermassive merger remnants by computing sequences of spherically-symmetric and axisymmetric uniformly and differentially rotating equilibrium solutions to the general-relativistic stellar structure equations. Using a set of finite-temperature nuclear equations of state, we find that hot maximum-mass critically spinning configurations generally do not support larger baryonic masses than their cold counterparts. However, subcritically spinning configurations with mean density of less than a few times nuclear saturation density yield a significantly thermally enhanced mass. Even without decreasing the maximum mass, cooling and other forms of energy loss can drive the remnant to an unstable state. We infer secular instability by identifying approximate energy turning points in equilibrium sequences of constant baryonic mass parametrized by maximum density. Energy loss carries the remnant along the direction of decreasing gravitational mass and higher density until instability triggers collapse. Since configurations with more thermal pressure support are less compact and thus begin their evolution at a lower maximum density, they remain stable for longer periods after merger.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We investigate by two-dimensional axisymmetric relativistic hydrodynamical simulations (1) jet propagations through an envelope of a rapidly rotating and collapsing massive star, which is supposed to be a progenitor of long duration gamma ray bursts (GRBs), (2) breakouts and subsequent expansions into stellar winds and (3) accompanying photospheric emissions. We find that if the envelope rotates uniformly almost at the mass shedding limit, its outer part stops contracting eventually when the centrifugal force becomes large enough. Then another shock wave is formed, propagates outwards and breaks out of the envelope into the stellar wind. Which breaks out earlier, the jet or the centrifugal bounce-induced shock, depends on the timing of jet injection. If the shock breakout occurs earlier owing to a later injection, the jet propagation and subsequent photospheric emissions are affected substantially. We pay particular attention to observational consequences of the difference in the timing of jet injection. We calculate optical depths to find the location of photospheres, extracting densities and temperatures at appropriate retarded times from the hydrodynamical data. We show that the luminosity and observed temperature of the photospheric emissions are both much lower than those reported in previous studies. Although luminosities are still high enough for GRBs, the observed temperature are lower than the energy at the spectral peak expected by the Yonetoku-relation. This may imply that energy exchanges between photons and matter are terminated deeper inside or some non-thermal processes are operating to boost photon energies.

We present a new formula to numerically construct configurations in rotational equilibrium, which consist of multiple layers. Each layer rotates uniformly or differentially according to cylindrical rotation-laws that are different from layer to layer. Assuming a different barotropic equation of state (EOS) for each layer, we solve the Bernoulli equation in each layer separately and combine the solutions by imposing continuity of the pressure at each boundary of the layers. It is confirmed that a single continuous barotropic EOS is incompatible with the junction condition. Identifying appropriate variables to be solved, we construct a convergent iteration scheme. For demonstration, we obtain two-layered configurations, each layer of which rotates rapidly with either an "$\Omega$-constant law" or a "$j$-constant law" or a "$v$-constant law". Other rotation laws and/or a larger number of layers can be treated similarly. We hope that this formula will be useful in studying the stellar evolution in multi-dimension with the non-spherical configuration induced by rotation being fully taken into account.

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 $\gamma$-ray bursts, which are likely to be composed of a black hole and surrounding disk, therefore could be constrained by GW observation.

We report the effects of quark-hadron phase transition on the structures of general relativistic stars with purely toroidal magnetic field. For the mixed phase, we take into account of the finite-size effects, which lead to non-uniform "Pasta" structures. Our study is based on axisymetric and stationary formalism including purely toroidal magnetic field. For hybrid stars, we find the characteristic distribution of magnetic field, which has a discontinuity originated in the quark-hadron mixed phase. These distributions of magnetic field will change astrophysical phenomena, such as cooling processes.

We investigate structures of hybrid stars, which feature quark core surrounded by a hadronic matter mantle, with super-strong toroidal magnetic fields in full general relativity. Modeling the equation of state (EOS) with a first order transition by bridging the MIT bag model for the description of quark matter and the nuclear EOS by Shen et al., we numerically construct thousands of the equilibrium configurations to study the effects of the phase transition. It is found that the appearance of the quark phase can affect distributions of the magnetic fields inside the hybrid stars, making the maximum field strength about up to 30 % larger than for the normal neutron stars. Using the equilibrium configurations, we explore the possible evolutionary paths to the formation of hybrid stars due to the spin-down of magnetized rotating neutron stars. We find that the energy release by the phase transition to the hybrid stars is quite large ($\la 10^{52} \rm erg$) even for super strongly magnetized compact stars. Our results suggest that the strong gravitational-wave emission and the sudden spin-up signature could be observable signals of the QCD phase transition, possibly for a source out to Megaparsec distances.

General relativistic simulations for the merger of binary neutron stars are performed as an extension of a previous work\cite{Shibata:2006nm}. We prepare binary neutron stars with a large initial orbital separation and employ the moving-puncture formulation, which enables to follow merger and ringdown phases for a long time, even after black hole formation. For modeling inspiraling neutron stars, which should be composed of cold neutron stars, the Akmal-Pandhalipande-Ravenhall (APR) equation of state (EOS) is adopted. After the onset of the merger, the hybrid-type EOS is used; i.e., the cold and thermal parts are given by the APR and $\Gamma$-law EOSs, respectively. Three equal-mass binaries each with mass $1.4M_\odot,1.45M_\odot,1.5M_\odot$ and two unequal-mass binaries with mass 1.3--$1.6M_\odot$, 1.35--$1.65M_\odot$ are prepared. We focus primarily on the black hole formation case, and explore mass and spin of the black hole, mass of disks which surround the black hole, and gravitational waves emitted during the black hole formation. We find that (i) for the systems of $m_0=2.9$--$3.0M_\odot$ and of mass ratio $\approx 0.8$, the mass of disks which surround the formed black hole is 0.006--$0.02M_{\odot}$; (ii) the spin of the formed black hole is $0.78 \pm 0.02$ when a black hole is formed after the merger in the dynamical time scale. This value depends weakly on the total mass and mass ratio, and is about 0.1 larger than that of a black hole formed from nonspinning binary black holes; (iii) for the black-hole formation case, Fourier spectrum shape of gravitational waves emitted in the merger and ringdown phases has a universal qualitative feature irrespective of the total mass and mass ratio, but quantitatively, the spectrum reflects the parameters of the binary neutron stars.

We study structures of general relativistic compact stars with exotic matter. Our study is based on axisymetric and stationary formalism including purely toroidal magnetic field. We also study the finite size effects of quark-hadron mixed phase on structures of magnetars. For hybrid stars, we find a characteristic distribution of magnetic field, which has a discontinuity originated in the quark-hadron mixed phase. These distributions of magnetic field will change astrophysical phenomena, such as cooling processes.

We investigate equilibrium sequences of relativistic stars containing purely toroidal magnetic fields with four kinds of realistic equations of state (EOSs) of SLy (Douchin et al.), FPS (Pandharipande et al.), Shen (Shen et al.), and LS (Lattimer & Swesty). We numerically construct thousands of equilibrium configurations. Particularly we pay attention to the equilibrium sequences of constant baryon mass and/or constant magnetic flux, which model evolutions of an isolated neutron star. Important properties obtained in this study are summarized as follows ; (1) The dependence of the mass-shedding angular velocity on the EOSs is determined from that of the non-magnetized case. The stars with Shen(FPS) EOS reach the mass-shedding limit at the smallest(largest) angular velocity, while the stars with SLy or Lattimer-Swesty EOSs take the moderate values. (2) For the supramassive sequences, the equilibrium configurations are found to be generally oblate for the realistic EOSs in sharp contrast to the polytropic stars. For FPS(LS) EOS, the parameter region which permits the prolately deformed stars is widest(narrowest). For SLy and Shen EOS, it is in medium. Furthermore, the angular velocities $\Omega_{\rm up}$, above which the stars start to spin up as they lose angular momentum, are found to depend sharply on the realistic EOSs. Our analysis indicates that the hierarchy of this spin up angular velocity is $\Omega_{\rm up,SLy} > \Omega_{\rm up,FPS} > \Omega_{\rm up,LS}>\Omega_{\rm up,Shen}$ and this relation holds even if the sequences have strong magnetic fields. Our results suggest the EOSs within the relativistic stars containing purely toroidal magnetic fields can be constrained by observing the angular velocity, the gravitational wave, and the signature of the spin up.

We study the stability of neutron stars with toroidal magnetic fields by magnetohydrodynamic simulation in full general relativity under assumption of axial symmetry. Nonrotating and rigidly rotating neutron stars are prepared for a variety of magnetic field configuration. For modeling the neutron stars, the polytropic equation of state with the adiabatic index $\Gamma=2$ is used for simplicity. It is found that nonrotating neutron stars are dynamically unstable for the case that toroidal magnetic field strength varies $\propto \varpi^{2k-1}$ with $k\geq 2$ (here $\varpi$ is the cylindrical radius), whereas for $k=1$ the neutron stars are stable. After the onset of the instability, unstable modes grow approximately in the Alfv\'en time scale and, as a result, a convective motion is excited to change the magnetic field profile until a new state, which is stable against axisymmetric perturbation, is reached. We also find that rotation plays a role in stabilization, although the instability still sets in in the Alfv\'en time scale when the ratio of magnetic energy to rotational kinetic energy is larger than a critical value $\sim 0.2$. Implication for the evolution of magnetized protoneutron stars is discussed.

We investigate the effects of the purely toroidal magnetic field on the equilibrium structures of the relativistic stars. The master equations for obtaining equilibrium solutions of relativistic rotating stars containing purely toroidal magnetic fields are derived for the first time. To solve these master equations numerically, we extend the Cook-Shapiro-Teukolsky scheme for calculating relativistic rotating stars containing no magnetic field to incorporate the effects of the purely toroidal magnetic fields. By using the numerical scheme, we then calculate a large number of the equilibrium configurations for a particular distribution of the magnetic field in order to explore the equilibrium properties. We also construct the equilibrium sequences of the constant baryon mass and/or the constant magnetic flux, which model the evolution of an isolated neutron star as it loses angular momentum via the gravitational waves. Important properties of the equilibrium configurations of the magnetized stars obtained in this study are summarized as follows ; (1) For the non-rotating stars, the matter distribution of the stars is prolately distorted due to the toroidal magnetic fields. (2) For the rapidly rotating stars, the shape of the stellar surface becomes oblate because of the centrifugal force. But, the matter distribution deep inside the star is sufficiently prolate for the mean matter distribution of the star to be prolate. (3) The stronger toroidal magnetic fields lead to the mass-shedding of the stars at the lower angular velocity. (4) For some equilibrium sequences of the constant baryon mass and magnetic flux, the stars can spin up as they lose angular momentum.

We present our numerical comparisons between the BSSN formulation widely used in numerical relativity today and its adjusted versions using constraints. We performed three testbeds: gauge-wave, linear wave, and Gowdy-wave tests, proposed by the Mexico workshop on the formulation problem of the Einstein equations. We tried three kinds of adjustments, which were previously proposed from the analysis of the constraint propagation equations, and investigated how they improve the accuracy and stability of evolutions. We observed that the signature of the proposed Lagrange multipliers are always right and the adjustments improve the convergence and stability of the simulations. When the original BSSN system already shows satisfactory good evolutions (e.g., linear wave test), the adjusted versions also coincide with those evolutions; while in some cases (e.g., gauge-wave or Gowdy-wave tests) the simulations using the adjusted systems last 10 times as long as those using the original BSSN equations. Our demonstrations imply a potential to construct a robust evolution system against constraint violations even in highly dynamical situations.

We investigate equilibrium sequences of magnetized rotating stars with four kinds of realistic equations of state (EOSs) of SLy (Douchin et al.), FPS (Pandharipande et al.), Shen (Shen et al.), and LS (Lattimer & Swesty). Employing the Tomimura-Eriguchi scheme to construct the equilibrium configurations. we study the basic physical properties of the sequences in the framework of Newton gravity. In addition we newly take into account a general relativistic effect to the magnetized rotating configurations. With these computations, we find that the properties of the Newtonian magnetized stars, e.g., structure of magnetic field, highly depends on the EOSs. The toroidal magnetic fields concentrate rather near the surface for Shen and LS EOSs than those for SLy and FPS EOSs. The poloidal fields are also affected by the toroidal configurations. Paying attention to the stiffness of the EOSs, we analyze this tendency in detail. In the general relativistic stars, we find that the difference due to the EOSs becomes small because all the employed EOSs become sufficiently stiff for the large maximum density, typically greater than $10^{15}\rm{g} \rm{cm}^{-3}$. The maximum baryon mass of the magnetized stars with axis ratio $q\sim 0.7$ increases about up to twenty percents for that of spherical stars. We furthermore compute equilibrium sequences at finite temperature, which should serve as an initial condition for the hydrodynamic study of newly-born magnetars. Our results suggest that we may obtain information about the EOSs from the observation of the masses of magnetars.

We study gravitational waves from a particle moving around a system of a point mass with a disk in Newtonian gravitational theory. A particle motion in this system can be chaotic when the gravitational contribution from a surface density of a disk is comparable with that from a point mass. In such an orbit, we sometimes find that there appears a phase of the orbit in which particle motion becomes to be nearly regular (the so-called ``stagnant motion'') for a finite time interval between more strongly chaotic phases. To study how these different chaotic behaviours affect on observation of gravitational waves, we investigate a correlation of the particle motion and the waves. We find that such a difference in chaotic motions reflects on the wave forms and energy spectra. The character of the waves in the stagnant motion is quite different from that either in a regular motion or in a more strongly chaotic motion. This suggests that we may make a distinction between different chaotic behaviours of the orbit via the gravitational waves.