Possible Indirect Confirmation of the Existence of Pop III Massive Stars by Gravitational Wave

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.