our model reduces to a system of self-propelled particles with velocity alignment akin to the one proposed by Vicsek [@Vicsek1995] but in a continuous-time formulation [@Ratushnaya2007; @Farrell2012; @Chepizhko2013; @CAP; @CAPhard]. As such, our model should reproduce in this limit the main collective features of Vicsek-like models. At low values of $\text{g}$ the system sets into a disordered homogeneous state. For values above the flocking threshold, $\text{g}>\text{g}^c$, the system develops global polar order. Flocking is understood as a first-order-like phase transition, characterized by a coexistence region where high density bands coexist with a disordered low-density background. For higher values of the coupling strength, the system eventually leaves the coexistence region, rendering the traveling bands unstable and leading to the so-called fluctuating flocking state or Toner-Tu phase [@Toner1995; @Toner2005; @Ramaswamy2010; @Chate2008; @Mishra2010]. The aligning term, Eq. (\[theta\_dot\]), is reminiscent of the (noisy) Kuramoto model of synchronization phenomena and the planar XY model of magnetism. Indeed, our model can be thought of as a Kuramoto model on a time evolving network generated by the positions of the particles [@AlbertPRX]. The internal variable $\theta_i$, is equivalent to the phase of an oscillator which here dictates its self-propulsion direction. In the absence of motion, i.e. 