  We consider the inertial Kuramoto model of $N$ globally coupled oscillators
characterized by both their phase and angular velocity, in which there is a
time delay in the interaction between the oscillators. Besides the academic
interest, we show that the model can be related to a network of phase-locked
loops widely used in electronic circuits for generating a stable frequency at
multiples of an input frequency. We study the model for a generic choice of the
natural frequency distribution of the oscillators, to elucidate how a
synchronized phase bifurcates from an incoherent phase as the coupling constant
between the oscillators is tuned. We show that in contrast to the case with no
delay, here the system in the stationary state may exhibit either a subcritical
or a supercritical bifurcation between a synchronized and an incoherent phase,
which is dictated by the value of the delay present in the interaction and the
precise value of inertia of the oscillators. Our theoretical analysis,
performed in the limit $N \to \infty$, is based on an unstable manifold
expansion in the vicinity of the bifurcation, which we apply to the kinetic
equation satisfied by the single-oscillator distribution function. We check our
results by performing direct numerical integration of the dynamics for large
$N$, and highlight the subtleties arising from having a finite number of
oscillators.
