We determine the distribution of cosmic string loops directly from simulations, rather than determining the loop production function and inferring the loop distribution from that. For a wide range of loop lengths, the results agree well with a power law exponent -2.5 in the radiation era and -2 in the matter era, the universal result for any loop production function that does not diverge at small scales. Our results extend those of Ringeval, Sakellariadou, and Bouchet: we are able to run for 15 times longer in conformal time and simulate a volume 300-2400 times larger. At the times they reached, our simulation is in general agreement with the more negative exponents they found, -2.6 and -2.4. However, our simulations show that this was a transient regime; at later times the exponents decline to the values above. This provides further evidence against models with a rapid divergence of the loop density at small scales, such as ``model 3'' used to analyze LIGO data and predict LISA sensitivity.
Read full abstract