Constraining the Statistics of Population III Binaries

We perform a cosmological simulation in order to model the growth and evolution of Population III (Pop III) stellar systems in a range of host minihalo environments. A Pop III multiple system forms in each of the ten minihaloes, and the overall mass function is top-heavy compared to the currently observed initial mass function in the Milky Way. Using a sink particle to represent each growing protostar, we examine the binary characteristics of the multiple systems, resolving orbits on scales as small as 20 AU. We find a binary fraction of ~36%, with semi-major axes as large as 3000 AU. The distribution of orbital periods is slightly peaked at < 900 yr, while the distribution of mass ratios is relatively flat. Of all sink particles formed within the ten minihaloes, ~50% are lost to mergers with larger sinks, and ~50% of the remaining sinks are ejected from their star-forming disks. The large binary fraction may have important implications for Pop III evolution and nucleosynthesis, as well as the final fate of the first stars.