Cold initial conditions with errors

First I wanted to demonstrate that I have sorted out the normalization problem I had, the density no longer depends on the binning as shown by this plot, which has 10 particles per bin and 300 particles per bin plotted on top of one another:


Next, the cold initial condition runs are no longer contaminated by the warm runs, and they have also been run to 700 ATUs instead of 400. The plots are below. The color coding used here is the same as I used initially (in this post).


A few observations: While it follows the same scaling relation in the density, the red run seems to have many fewer particles in the interior core of the halo than the others (this can be seen because the innermost data point containing the innermost 10 particles is significantly farther out than in the other runs). The red and the yellow, which are the scenarios that form via merger, seem to have a more rounded relation than the blue and green, which form monolithically. If the particles are binned more coarsely, however, this trend seems to be less apparent. Here are the x>0 portions of the 4 runs plotted together, with coarser binning.


You will see that run 4 has now had a chance to virialize and looks much like the other runs. As yet I have not established the extent to which these quantitatively follow an $r^{-1/2}$ scaling relation, but by eye (with coarse enough binning) they seem to.