Gaussian Mixture models for the local peculiar velocity distribution of neighbours of the Sun. Left: One- to ten-component models, together with Aikake and Bayesian information criteria values relative to the best fit. Middle: One- and two-component models, together with the decomposition of the latter. Right: Probability of each star belonging to the high-velocity component of the two-component model.