The numerical computation of multicomponent mixture critical points has been the subject of much study due to their theoretical and practical importance. Both deterministic and stochastic methods have been applied with varying degrees of reliability and robustness. In this work, we utilize numerical polynomial homotopy continuation (NPHC) to reliably identify all mixture critical points. This method is unique due to its robustness, initialization‐free nature and ease of parallelization. For a given system of equations, all complex solutions are found. Computational times are also found to be invariant to mixture composition. We validate this technique against previous work and extend the method to mixtures of up to eight components. NPHC is shown to be a modern and powerful technique which offers mathematical reliability at a moderate computational cost. © 2016 American Institute of Chemical Engineers AIChE J, 62: 4497–4507, 2016