We present a proposal for the full phase-space distribution of the Milky Way halo. The model is axially and reflection symmetric and its time evolution is self-similar. It describes the halo as a set of discrete dark matter flows with stated densities and velocity vectors everywhere. We first discuss the general conditions under which the time evolution of a cold collisionless self-gravitating fluid is self-similar, and show that symmetry is not necessary for self-similarity. When spherical symmetry is imposed, the model is the same as described by Fillmore and Goldreich, and by Bertschinger, twenty-three years ago. The spherically symmetric model depends on one dimensionless parameter $ϵ$ and two dimensionful parameters. We set $ϵ=0.3$, a value consistent with the slope of the power spectrum of density perturbations on galactic scales. The dimensionful parameters are determined by the galactic rotation velocity ($220\mathrm{km}/\mathrm{s}$) at the position of the Sun and by the age of the Galaxy (13.7 Gyr). The properties of the outer caustics are derived in the spherically symmetric model. The structure of the inner halo depends on the angular momentum distribution of the dark matter particles. We assume that distribution to be axial and reflection symmetric, and dominated by net overall rotation. The inner caustics are rings whose radii are determined in terms of a single additional parameter $j_{\max }$. We summarize the observational evidence in support of the model. The evidence is consistent with $j_{\max }=0.18$ in Concordance cosmology, equivalent to $j_{\max ,\mathrm{old}}=0.26$ in Einstein-de Sitter cosmology. We give formulas to estimate the flow densities and velocity vectors anywhere in the Milky Way halo. The properties of the first 40 flows at the location of the Earth are listed.

• Received 30 May 2008

DOI:https://doi.org/10.1103/PhysRevD.78.063508