Predicting structural properties of dark matter haloes is one of the fundamental goals of modern cosmology. We use the suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. We find that in order to understand the structure of dark matter haloes and to make 1–2 per cent accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the ratio of the virial radius to the core radius in the Navarro–Frenk–White (NFW) profile. For massive haloes, the average density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter α in the Einasto approximation. We show that haloes progress through three stages of evolution. They start as rare density peaks and experience fast and nearly radial infall that brings mass closer to the centre, producing a highly concentrated halo. Here, the halo concentration increases with increasing halo mass and the concentration is defined by the α parameter with a nearly constant core radius. Later haloes slide into the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. At this stage, the concentration does not depend on halo mass. Once the rate of accretion and merging slows down, haloes move into the domain of declining concentration–mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. Accurate analytical fits are provided.

MultiDark simulations: the story of dark matter halo concentrations and density profiles