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

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

Major construction and initial-phase operation of a second-generation gravitational-wave detector KAGRA has been completed. The entire 3-km detector is installed underground in a mine in order to be isolated from background seismic vibrations on the surface. This allows us to achieve a good sensitivity at low frequencies and high stability of the detector. Bare-bones equipment for the interferometer operation has been installed and the first test run was accomplished in March and April of 2016 with a rather simple configuration. The initial configuration of KAGRA is named {\it iKAGRA}. In this paper, we summarize the construction of KAGRA, including the study of the advantages and challenges of building an underground detector and the operation of the iKAGRA interferometer together with the geophysics interferometer that has been constructed in the same tunnel.

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

We investigate the nucleosynthesis process in high-entropy ($s/k_{\rm B}\gtrsim100$) and very fast-expanding ($\tau_{\rm exp}\sim10^{-3}\ {\rm s}$) materials. In such a material with the electron fraction near 0.5, an interesting nucleosynthesis process occurs. In this process, the abundance distribution of heavy-nuclei of $A>100$ achieve quasi-statistical equilibrium (QSE) at high temperature and the abundances are frozen at the end of the nucleosynthesis. We explain this abundance distribution using the "alpha-constrained QSE" abundances formulated in this paper. We demonstrate that this nucleosynthesis would occur in neutrino-driven winds from massive proto-neutron stars in hypernovae, where $A\sim140$ $p$-nuclei are synthesized.

We investigate the nucleosynthesis in the neutrino-driven winds blown off from a $3M_\odot$ massive proto-neutron star (mPNS) temporarily formed during the collapse of a $100M_\odot$ star. Such mPNSs would be formed in hypernovae. We construct steady and spherically symmetric wind solutions. We set large neutrino luminosities of $\sim 10^{53}\ {\rm erg\ s^{-1}}$ and average energies of electron neutrinos and antineutrinos in the ranges of $\epsilon_{\nu_e}=9-16\ {\rm MeV}$ and $\epsilon_{\bar{\nu}_e}=11-18\ {\rm MeV}$ based on a recent numerical relativity simulation. The wind solutions indicate much shorter temperature-decrease timescale than that of the winds from ordinary PNSs and, depending on $\epsilon_\nu$, the winds can be both neutron-rich and proton-rich. In the neutron-rich wind, the $r$-process occurs and the abundance distribution of a fiducial wind model of the mPNS gives an approximate agreement with the abundance pattern of metal-poor weak $r$ star HD~122563, although the third-peak elements are produced only when the $\bar{\nu}_e$ energy is much larger than the $\nu_e$ energy. In the proton-rich wind, the strong $\nu p$-process occurs and $A>100$ nuclides are synthesized. The synthesized nuclei can be neutron-rich in some cases because the large neutrino luminosity of the mPNS supplies a sufficient amount of neutrons.

The cosmological scenario of the dark matter generated in the hidden gauge theory based on the grand unification is discussed. It is found that the stability of the dark matter halo of our Galaxy and the cosmic ray observation constrain, respectively, the dark matter mass and the unification scale between the standard model and the hidden gauge theory sectors. To obtain a phenomenologically consistent thermal evolution, the entropy of the standard model sector needs to be increased. We therefore propose a scenario where the mini-inflation is induced from the potential coupled to the Standard model sector, in particular the Higgs sector. This scenario makes consistent the current dark matter density as well as the baryon-to-photon ratio for the case of pion dark matter. For the glueball or heavy pion of hidden gauge theory, an additional mini-inflation in the standard model sector before the leptogenesis is required. We also propose the possibility to confirm this scenario by known prospective experimental approaches.