Velocity correlation function for a sample of “stars” extracted from the stellar halos in the Aquarius simulations (Lowing et al. 2015). The stars are located in a sphere of 2.8 kpc and have magnitudes in the range 8 ≤ G ≤ 12.5. The sample size is the same as our observed halo dataset. A comparison to Fig. 3 shows a similar excess of pairs in the simulations as in the data, although there is some variance in the simulations.