We aim to locate previously unknown stellar clusters using the VISTA variables in the Via Lactea Survey (VVV) catalogue data. The method fits a mixture model of Gaussian densities and background noise and uses the expectation maximization algorithm to pre-filtered near-infrared survey stellar catalogue data; it was developed by the authors for the UKIDSS Galactic Plane Survey (GPS). The search located 88 previously unknown candidates, most of which are embedded stellar cluster candidates, and 39 previously unknown sites of star formation in the 562deg^2^ covered by VVV in the Galactic bulge and the southern disk.