LIGO has detected gravitational waves from massive binary black hole mergers. In order to explain the origin of such massive stellar-mass black holes, extreme metal poor stars including first stars have been invoked. However, black holes do not carry information of the metallicity. In order to check the metallicity dependence of the black hole formation, we focus on galactic black hole-main sequence binaries (BH-MSs). Using a binary population synthesis method, we find that $\gaia$ can detect $\sim200-400$ BH-MSs whose metallicity is $\zsun$ and $\sim70-400$ BH-MSs whose metallicity is $0.1\zsun$. With the spectroscopic observation on 4-m class telescopes, we can check the metallicity of BH-MSs. The metallicity dependence of the black hole formation might be checked by the astrometry and spectroscopic observations.

We apply the independent component analysis (ICA) to the real data from a gravitational wave detector for the first time. Specifically we use the iKAGRA data taken in April 2016, and calculate the correlations between the gravitational wave strain channel and 35 physical environmental channels. Using a couple of seismic channels which are found to be strongly correlated with the strain, we perform ICA. Injecting a sinusoidal continuous signal in the strain channel, we find that ICA recovers correct parameters with enhanced signal-to-noise ratio, which demonstrates usefulness of this method. Among the two implementations of ICA used here, we find the correlation method yields the optimal result for the case environmental noises act on the strain channel linearly.

We have devised fitting formulae for evolution tracks of massive stars with $8 \lesssim M/M_\odot \lesssim 160$ under extreme metal poor (EMP) environments for $\log (Z/Z_\odot) = -2, -4, -5, -6$, and $-8$, where $M_\odot$ and $Z_\odot$ are the solar mass and metallicity, respectively. Our fitting formulae are based on reference stellar models which we have newly obtained by simulating the time evolutions of EMP stars. Our fitting formulae take into account stars ending with blue supergiant (BSG) stars, and stars skipping Hertzsprung gap (HG) phases and blue loops, which are characteristics of massive EMP stars. In our fitting formulae, stars may remain BSG stars when they finish their core Helium burning (CHeB) phases. Our fitting formulae are in good agreement with our stellar evolution models. We can use these fitting formulae on the SSE, BSE, NBODY4, and NBODY6 codes, which are widely used for population synthesis calculations and star cluster simulations. These fitting formulae should be useful to make theoretical templates of binary black holes formed under EMP environments.

We performed Population III (Pop III) binary evolution by using population synthesis simulations for seven different models. We found that Pop III binaries tend to be binary black holes (BBHs) with chirp mass $M_{\rm chirp} \sim 30~M_{\odot}$ and they can merge at present day due to long merger time. The merger rate densities of Pop III BBHs at $z=0$ ranges 3.34--21.2 $\rm /yr/Gpc^3$ which is consistent with the aLIGO/aVIRGO result of 9.7--101 $\rm /yr/Gpc^3$. These Pop III binaries might contribute to some part of the massive BBH gravitational wave (GW) sources detected by aLIGO/aVIRGO. We also calculated the redshift dependence of Pop III BBH mergers. We found that Pop III low spin BBHs tend to merge at low redshift, while Pop III high spin BBHs do at high redshift, which can be confirmed by future GW detectors such as ET, CE, and DECIGO. These detectors can also check the redshift dependence of BBH merger rate and spin distribution. Our results show that except for one model, the mean effective spin $\left\langle \chi_{\rm eff} \right\rangle$ at $z=0$ ranges $0.02$--$0.3$ while at $z=10$ it does $0.16$--$0.64$. Therefore, massive stellar-mass BBH detection by GWs will be a key for the stellar evolution study in the early universe.

We estimate formation rates of LB-1-like systems through dynamical interactions in the framework of the theory of stellar evolution before the discovery of the LB-1 system. The LB-1 system contains $\sim 70M_\odot$ black hole (BH), so-called pair instability (PI)-gap BH, and B-type star with solar metallicity, and has nearly zero eccentricity. The most efficient formation mechanism is as follows. In an open cluster, a naked helium (He) star (with $\sim 20M_\odot$) collides with a heavy main-sequence (MS) star (with $\sim 50M_\odot$) which has a B-type companion. The collision results in a binary consisting of the collision product and B-type star with a high eccentricity. The binary can be circularized through the dynamical tide with radiative damping of the collision-product envelope. Finally, the collision product collapses to a PI-gap BH, avoiding pulsational pair instability and pair instability supernovae because its He core is as massive as the pre-colliding naked He star. We find that the number of LB-1-like systems in the Milky Way galaxy is $\sim 0.01 (\rho_{\rm oc} / 10^4 M_\odot \mbox{pc}^{-3})$, where $\rho_{\rm oc}$ is the initial mass densities of open clusters. If we take into account LB-1-like systems with O-type companion stars, the number increases to $\sim 0.03 (\rho_{\rm oc} / 10^4 M_\odot \mbox{pc}^{-3})$. This mechanism can form LB-1-like systems at least 10 times more efficiently than the other mechanisms: captures of B-type stars by PI-gap BHs, stellar collisions between other type stars, and stellar mergers in hierarchical triple systems. We conclude that no dynamical mechanism can explain the presence of the LB-1 system.

It is generally believed that Type Ia supernovae are thermonuclear explosions of carbon-oxygen white dwarfs (WDs). However, there is currently no consensus regarding the events leading to the explosion. A binary WD (WD-WD) merger is a possible progenitor of Type Ia supernovae. Space-based gravitational wave (GW) detectors with great sensitivity in the decihertz range like DECIGO can observe WD-WD mergers directly. Therefore, access to the deci-Hz band of GWs would enable multi-messenger observations of Type Ia supernovae to constrain their progenitor and explosion mechanism. In this paper, we consider the event rate of WD-WD mergers and minimum detection range to observe one WD-WD merger per year, using nearby galaxy catalog and the relation between the Ia supernova and the host galaxy. Furthermore, we calculate the DECIGO's ability to localize WD-WD mergers and to determine the masses of binary mergers. We estimate that if the deci-Hz GW observatory can detect the GW whose amplitude is $h\sim10^{-20}[\rm Hz^{-1/2}]$ at 0.1 Hz, 1000 times higher than the detection limit of DECIGO. In fact, DECIGO is expected to detect WD-WD ($1M_{\odot}-1M_{\odot}$) mergers within $z=0.115$,corresponding to the detection rate of $\sim20000\,\rm yr^{-1}$, and identify the host galaxy of WD-WD mergers for $\sim8000$ WD-WDs only by the GW detection.

Double neutron star (DNS) systems are produced from massive binaries. A supernova (SN) explosion of an extremely stripped star is expected to occur at the final stage of DNS formation. This type of SN is called ultra-stripped SN (USSN). A recent research reveals that a type Ic SN, iPTF 14gqr (SN 2014ft) has low ejecta mass ($\approx0.2~\mathrm{M}_\odot$) and its progenitor has a helium envelope with its mass $\sim0.01~\mathrm{M}_\odot$. Then this SN is interpreted as an USSN, and this is the first discovery of the USSN. Furthermore, the observation of iPTF 14gqr provides us with some information about its formation history. Here, we perform population synthesis calculations so as to estimate the detection rate of iPTF 14gqr like USSNe by optical transient surveys, intermediate Palomar Transient Factory (iPTF), Zwicky Transient Facility (ZTF), and Large Synoptic Survey Telescope (LSST). We find that iPTF, ZTF, and LSST can observe iPTF 14gqr like USSNe at the rates of 1, 40, and 5 $\mathrm{yr}^{-1}$, respectively. The iPTF may detect 4 iPTF 14gqr like USSNe during four years observation. Additionally, we investigate effects of some binary physics on the orbital parameters of DNS systems and formation channels.

KAGRA is a new gravitational wave detector which aims to begin joint observation with Advanced LIGO and Advanced Virgo from late 2019. Here, we present KAGRA's possible upgrade plans to improve the sensitivity in the decade ahead. Unlike other state-of-the-art detectors, KAGRA requires different investigations for the upgrade since it is the only detector which employs cryogenic cooling of the test mass mirrors. In this paper, investigations on the upgrade plans which can be realized by changing the input laser power, increasing the mirror mass, and injecting frequency dependent squeezed vacuum are presented. We show how each upgrade affects to the detector frequency bands and also discuss impacts on gravitational-wave science. We then propose an effective progression of upgrades based on technical feasibility and scientific scenarios.

Future missions for long gammma-ray burst (GRB) observations at high redshift such as HiZ-GUNDAM and THESEUS will provide clue to the star formation history in our universe. In this paper focusing on high redshift (z>8) GRBs, we calculate the detection rate of long GRBs by future observations, considering both Population (Pop) I&II stars and Pop III stars as GRB progenitors. For the Pop I&II star formation rate (SFR), we adopt an up-to-date model of high-redshift SFR based on the halo mass function and dark matter accretion rate obtained from cosmological simulations. We show that the Pop I&II GRB rate steeply decreases with redshift. This would rather enable us to detect the different type of GRBs, Pop III GRBs, at very high redshift. If 10% or more Pop III stars die as an ultra-long GRB, the future missions would detect such GRBs in one year in spite of their low fluence. More luminous GRBs are expected from massive compact Pop III stars produced via the binary merger. In our conventional case, the detection rate of such luminous GRBs is 3-20 /yr (z>8). Those future observations contribute to revealing of the Pop III star formation history.

KAGRA is a second-generation interferometric gravitational-wave detector with 3-km arms constructed at Kamioka, Gifu in Japan. It is now in its final installation phase, which we call bKAGRA (baseline KAGRA), with scientific observations expected to begin in late 2019. One of the advantages of KAGRA is its underground location of at least 200 m below the ground surface, which brings small seismic motion at low frequencies and high stability of the detector. Another advantage is that it cools down the sapphire test mass mirrors to cryogenic temperatures to reduce thermal noise. In April-May 2018, we have operated a 3-km Michelson interferometer with a cryogenic test mass for 10 days, which was the first time that km-scale interferometer was operated at cryogenic temperatures. In this article, we report the results of this "bKAGRA Phase 1" operation. We have demonstrated the feasibility of 3-km interferometer alignment and control with cryogenic mirrors.

The recent detections of gravitational waves (GWs) reported by LIGO/Virgo collaborations have made significant impact on physics and astronomy. A global network of GW detectors will play a key role to solve the unknown nature of the sources in coordinated observations with astronomical telescopes and detectors. Here we introduce KAGRA (former name LCGT; Large-scale Cryogenic Gravitational wave Telescope), a new GW detector with two 3-km baseline arms arranged in the shape of an "L", located inside the Mt. Ikenoyama, Kamioka, Gifu, Japan. KAGRA's design is similar to those of the second generations such as Advanced LIGO/Virgo, but it will be operating at the cryogenic temperature with sapphire mirrors. This low temperature feature is advantageous for improving the sensitivity around 100 Hz and is considered as an important feature for the third generation GW detector concept (e.g. Einstein Telescope of Europe or Cosmic Explorer of USA). Hence, KAGRA is often called as a 2.5 generation GW detector based on laser interferometry. The installation and commissioning of KAGRA is underway and its cryogenic systems have been successfully tested in May, 2018. KAGRA's first observation run is scheduled in late 2019, aiming to join the third observation run (O3) of the advanced LIGO/Virgo network. In this work, we describe a brief history of KAGRA and highlights of main feature. We also discuss the prospects of GW observation with KAGRA in the era of O3. When operating along with the existing GW detectors, KAGRA will be helpful to locate a GW source more accurately and to determine the source parameters with higher precision, providing information for follow-up observations of a GW trigger candidate.

The long gamma ray bursts (GRBs) may arise from the core collapse of massive stars. However, the long GRB rate does not follow the star formation rate (SFR) at high redshifts. In this Letter, we focus on the binary merger model and consider the high spin helium stars after the merger as the progenitor of long GRBs. With this scenario, we estimate the GRB rate by the population synthesis method with the metallicity evolution. Low metallicity binaries are easier to become long GRB progenitors than those for solar metallicity due to the weak wind mass loss and the difference in the stellar evolution. In our results, the long GRB rate roughly agrees with the observed rate, and shows a similar behavior to the observed redshift evolution.

KAGRA is a 3-km interferometric gravitational wave telescope located in the Kamioka mine in Japan. It is the first km-class gravitational wave telescope constructed underground to reduce seismic noise, and the first km-class telescope to use cryogenic cooling of test masses to reduce thermal noise. The construction of the infrastructure to house the interferometer in the tunnel, and the initial phase operation of the interferometer with a simple 3-km Michelson configuration have been completed. The first cryogenic operation is expected in 2018, and the observing runs with a full interferometer are expected in 2020s. The basic interferometer configuration and the current status of KAGRA are described.

We propose a method for confirmation of the existence of Population III (Pop III) stars with massive black hole binaries as GW150914 in gravitational wave (GW) observation. When we get enough number of events, we want to determine which model is closer to reality, with and without Pop III stars. We need to prepare various "Pop I/II models" and various "Pop I/II/III models" and investigate which model is consistent with the events. To demonstrate our analysis, we simulate detections of GW events for some examples of population synthesis models with and without Pop III stars. We calculate the likelihood ratio with the realistic number of events and evaluate the probability of identifying the existence of Pop III stars. In typical cases, our analysis can distinguish between Pop I/II model and Pop I/II/III model with 90% probability by 22 GW signals from black hole-black hole binary mergers.

We study formation of stellar mass binary black holes (BBHs) originating from Population III (PopIII) stars, performing stellar evolution simulations for PopIII binaries with MESA. We find that a significant fraction of PopIII binaries form massive BBHs through stable mass transfer between two stars in a binary, without experiencing common envelope phases. We investigate necessary conditions required for PopIII binaries to form BBHs coalescing within the Hubble time with a semi-analytical model calibrated by the stellar evolution simulations. The formation efficiency of coalescing PopIII BBHs is estimated for two different initial conditions for PopIII binaries with large and small separations, respectively. Consequently, in both models, $\sim 10\%$ of the total PopIII binaries form BBHs only through stable mass transfer and $\sim 10\%$ of these BBHs merge due to gravitational wave emission within the Hubble time. Furthermore, the chirp mass of merging BBHs has a flat distribution over $15\lesssim M_{\rm chirp}/M_\odot \lesssim 35$. This formation pathway of PopIII BBHs is presumably robust because stable mass transfer is less uncertain than common envelope evolution, which is the main formation channel for Population II BBHs. We also test the hypothesis that the BBH mergers detected by LIGO originate from PopIII stars using our result and the total number of PopIII stars formed in the early universe as inferred from the optical depth measured by Planck. We conclude that the PopIII BBH formation scenario can explain the mass-weighted merger rate of the LIGO's O1 events with the maximal PopIII formation efficiency inferred from the Planck measurement, even without BBHs formed by unstable mass transfer or common envelope phases.

In the population synthesis simulations of Pop III stars, many BH (Black Hole)-BH binaries with merger time less than the age of the Universe $(\tau_{\rm H})$ are formed, while NS (Neutron Star)-BH binaries are not. The reason is that Pop III stars have no metal so that no mass loss is expected. Then, in the final supernova explosion to NS, much mass is lost so that the semi major axis becomes too large for Pop III NS-BH binaries to merge within $\tau_{\rm H}$. However it is almost established that the kick velocity of the order of $200-500{\rm~ km~s^{-1}}$ exists for NS from the observation of the proper motion of the pulsar. Therefore, the semi major axis of the half of NS-BH binaries can be smaller than that of the previous argument for Pop III NS-BH binaries to decrease the merging time. We perform population synthesis Monte Carlo simulations of Pop III NS-BH binaries including the kick of NS and find that the event rate of Pop III NS-BH merger rate is $\sim 1 {\rm Gpc^{-3} yr^{-1}}$. This suggests that there is a good chance of the detection of Pop III NS-BH mergers in O2 of Advanced LIGO and Advanced Virgo from this autumn.

Pre-DECIGO consists of three spacecraft arranged in an equilateral triangle with 100km arm lengths orbiting 2000km above the surface of the earth. It is hoped that the launch date will be in the late 2020s. Pre-DECIGO has one clear target: binary black holes (BBHs) like GW150914 and GW151226. Pre-DECIGO can detect $\sim 30M_\odot-30M_\odot$ BBH mergers up to redshift $z\sim 30$. The cumulative event rate is $\sim 1.8\times 10^{5}\,{\rm events~yr^{-1}}$ in the Pop III origin model of BBHs like GW150914, and it saturates at $z\sim 10$, while in the primordial BBH (PBBH) model, the cumulative event rate is $ \sim 3\times 10^{4}\,{\rm events~ yr^{-1}}$ at $z=30$ even if only $0.1\%$ of the dark matter consists of PBHs, and it is still increasing at $z=30$. In the Pop I/II model of BBHs, the cumulative event rate is $(3-10)\times10^{5}\,{\rm events~ yr^{-1}}$ and it saturates at $z \sim 6$. We present the requirements on orbit accuracy, drag free techniques, laser power, frequency stability, and interferometer test mass. For BBHs like GW150914 at 1Gpc, SNR$\sim 90$ is achieved with the definition of Pre-DECIGO in the $0.01-100$Hz band. Pre-DECIGO can measure the mass spectrum and the $z$-dependence of the merger rate to distinguish various models of BBHs like GW150914. Pre-DECIGO can also predict the direction of BBHs at $z=0.1$ with an accuracy of $\sim 0.3\,{\rm deg}^2$ and a merging time accuracy of $\sim 1$s at about a day before the merger so that ground-based GW detectors further developed at that time as well as electromagnetic follow-up observations can prepare for the detection of merger in advance. For intermediate mass BBHs at a large redshift $z > 10$, the QNM frequency after the merger can be within the Pre-DECIGO band so that the ringing tail can also be detectable to confirm the Einstein theory of general relativity with SNR$\sim 35$. [abridged]

Focusing on the remnant black holes after merging binary black holes, we show that ringdown gravitational waves of Population III binary black holes mergers can be detected with the rate of $5.9-500~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$ for various parameters and functions. This rate is estimated for the events with SNR$>8$ for the second generation gravitational wave detectors such as KAGRA. Here, ${\rm SFR_p}$ and ${\rm f_b}$ are the peak value of the Population III star formation rate and the fraction of binaries, respectively. When we consider only the events with SNR$>35$, the event rate becomes $0.046-4.21~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$. This suggest that for remnant black hole's spin $q_f>0.95$ we have the event rate with SNR$>35$ less than $0.037~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$, while it is $3-30~{\rm events~yr^{-1}}~({\rm SFR_p}/ (10^{-2.5}~M_\odot~{\rm yr^{-1}~Mpc^{-3}})) \cdot ({\rm [f_b/(1+f_b)]/0.33})$ for the third generation detectors such as Einstein Telescope. If we detect many Population III binary black holes merger, it may be possible to constrain the Population III binary evolution paths not only by the mass distribution but also by the spin distribution.

The existence of the ergoregion of the Kerr space-time has not been confirmed observationally yet. We show that the confirmation would be possible by observing the quasinormal mode in gravitational waves. As an example, using the recent population synthesis results of Pop III binary black holes, we find that the peak of the final merger mass ($M_f$) is about $50~\rm M_{\odot}$, while the fraction of the final spin $q_f = a_f/M_f > 0.7$ needed for the confirmation of a part of ergoregion is $\sim 77\%$. To confirm the frequency of the quasinormal mode, ${\rm SNR} > 35$ is needed. The standard model of Pop III population synthesis tells us that the event rate for the confirmation of more than $50\%$ of the ergoregion by the second generation gravitational wave detectors is $\sim 2.3$ ${\rm events\ yr^{-1}\ (SFR_p/(10^{-2.5}\ M_\odot yr^{-1}\ Mpc^{-3}))} \cdot (\rm [f_b/(1+f_b)]/0.33)$ where ${\rm SFR_p}$ and ${\rm f_b}$ are the peak value of the Pop III star formation rate and the fraction of binaries, respectively.

Using our population synthesis code, we found that the typical chirp mass of Pop III BH-BHs is $\sim30~\msun$ with the total mass of $\sim60~\msun$ so that the inspiral chirp signal as well as quasi normal mode (QNM) of the merging BH are interesting targets of KAGRA. The detection rate of the coalescing Pop III BH-BHs is $\sim$180 $\rm events~yr^{-1}$$(\rm SFR_p/(10^{-2.5}~\msun \rm~yr^{-1}~Mpc^{-3}))\cdot([f_b/(1+f_b)]/0.33)\cdot Err_{sys}$ in our standard model where $\rm SFR_{p},~f_b$ and $\rm Err_{sys}$ are the peak value of the Pop III star formation rate, the binary fraction and the systematic error with $\rm Err_{sys}=1$ for our standard model, respectively. To evaluate the robustness of chirp mass distribution and the range of $\rm Err_{sys}$, we examine the dependence of the results on the unknown parameters and the distribution functions. We found that the chirp mass has a peak at $\sim 30 ~\msun$ in most of parameters and distribution functions as well as $\rm Err_{sys}$ ranges from 0.046 to 4. The minimum rate corresponds to the worst model which we think unlikely so that unless $ {\rm ~(SFR_p/(10^{-2.5}~\msun~yr^{-1}~Mpc^{-3}))\cdot([f_b/(1+f_b)]/0.33) \ll 0.1}$, we expect the Pop III BH-BHs merger rate of at least one event per year by KAGRA. Nakano, Tanaka \& Nakamura (2015) show that if S/N of QNM is larger than 35, we can confirm or refute the General Relativity (GR) more than 5 sigma level. In our standard model, the detection rate of Pop III BH-BHs whose S/N is larger than 35 is $3.2~\rm events~yr^{-1}$$(\rm SFR_p/(10^{-2.5}~\msun \rm~yr^{-1}~Mpc^{-3}))\cdot([f_b/(1+f_b)]/0.33)\cdot Err_{sys}$. Thus, there is a good chance to check whether GR is correct or not in the strong gravity region.

We perform population synthesis simulations for Population III (Pop III) coalescing binary neutron stars (NS-NSs), neutron star - black hole binaries (NS-BHs), and binary black holes (BH-BHs) which merge within the age of the universe. We found that the typical mass of Pop III BH-BHs is $\sim 30\ \rm{M}_{\odot}$ so that the inspiral chirp signal of gravitational waves can be detected up to z=0.28 by KAGRA, Adv. LIGO, Adv. Virgo and GEO network. Our simulations suggest that the detection rate of the coalescing Pop III BH-BHs is $140 (68) \cdot ({\rm SFR}_{\rm p}/10^{-2.5} \rm{M}_{\odot} {\rm yr}^{-1} {\rm Mpc}^{-3}) \cdot {\rm Err}_{\rm sys} ~{\rm events} ~{\rm yr}^{-1}$ for the flat (Salpeter) initial mass function (IMF), respectively, where $\rm SFR_p$ and $\rm Err_{sys}$ are the peak value of the Pop III star formation rate and the possible systematic errors due to the assumptions in Pop III population synthesis, respectively. $\rm Err_{sys}=1$ correspond to conventional parameters for Pop I stars. From the observation of the chirp signal of the coalescing Pop III BH-BHs, we can determine both the mass and the redshift of the binary for the cosmological parameters determined by Planck satellite. Our simulations suggest that the cumulative redshift distribution of the coalescing Pop III BH-BHs depends almost only on the cosmological parameters. We might be able to confirm the existence of Pop III massive stars of mass $\sim 30~\rm M_{\odot}$ by the detections of gravitational waves if the merger rate of the Pop III massive BH-BHs dominates that of Pop I BH-BHs.