We perform cosmological hydrodynamics simulations with non-equilibrium primordial chemistry to obtain 59 minihalos that host first stars. The obtained minihalos are used as the initial conditions of local three-dimensional radiation hydrodynamics simulations to investigate the formation of the first stars. We find that two-thirds of the minihalos host multiple stars, while the other third has single stars. The mass of the stars found in our simulations are in the range of 1 M☉ ≲ M ≲ 300 M☉, peaking at severalx 10 M☉. Most of the very massive stars of ≳ 140 M☉ are born as single stars, although not all of the single stars are very massive. We also find a few stars of ≲ 1 M☉that are kicked by the gravitational three body interactions to the position distant from the center of mass. The frequency that a star forming minihalo contains a binary system is ∼50%. We also investigate the abundance pattern of the stellar remnants by summing up the contributions from the first stars in the simulations. Consequently, the pattern is compatible with that of the low metallicity damped Lyα systems or the extremely metal-poor (EMP) stars, if the mass spectrum obtained in our experiment is shifted to the low mass side by 0.2 dex. If we consider the case that an EMP star is born in the remnant of the individual minihalo without mixing with others, the chemical signature of the pair instability supernova is more prominent, because most of them are born as single stars.
early universe - stars: abundances - stars: formation - stars: Population III