In this paper we deal with a model describing the evolution in time of the density of aneural population in a state space, where the state is given by Izhikevich’s two -dimensional single neuron model. The main goal is to mathematically describe theoccurrence of a significant phenomenon observed in neurons populations, thesynchronization. To this end, we are making the transition to phasedensity population, and use Malkin theorem to calculate the phase deviations of a weaklycoupled population model.