This paper aims to identify the optimal number of flux modulation pole (FMP) in Vernier permanent magnet synchronous machines (VPMSMs) with concentrated tooth-coil windings. First, the influence of number of FMP on machine performance is revealed based on the unified permanent magnet magnetomotive force permeance analytical model. It is found that for a fixed stator slot number, there always exists an optimal number of FMP to maximize the torque, thanks to the enhanced contribution by the fundamental permeance harmonic. Then, the influence of critical dimensional parameters, including the machine diameter, the permanent magnet thickness, and the width of FMP, on the electromagnetic performance of VPMSMs with different numbers of FMP are parametrically analyzed by both analytical and finite element analyses. Results show that the VPMSMs with the optimal number of FMP always have advantages of a higher torque density together with a reduced PM volume, particularly for those with a relatively large diameter. Finally, two VPMSM prototypes with different numbers of FMP are manufactured and tested to validate the analyses.