  The synchronization of different $\gamma$-rhythms arising in different brain
areas has been implicated in various cognitive functions. Here, we focus on the
effect of the ubiquitous neuronal heterogeneity on the synchronization of PING
(pyramidal-interneuronal network gamma) and ING (interneuronal network gamma)
rhythms. The synchronization properties of rhythms depends on the response of
their collective phase to external input. We therefore determined the
macroscopic phase-response curve for finite-amplitude perturbations (fmPRC),
using numerical simulation of all-to-all coupled networks of integrate-and-fire
(IF) neurons exhibiting either PING or ING rhythms. We show that the intrinsic
neuronal heterogeneity can qualitatively modify the fmPRC. While the
phase-response curve for the individual IF-neurons is strictly positive (type
I), the fmPRC can be biphasic and exhibit both signs (type II). Thus, for PING
rhythms, an external excitation to the excitatory cells can, in fact, delay the
collective oscillation of the network, even though the same excitation would
lead to an advance when applied to uncoupled neurons. This paradoxical delay
arises when the external excitation modifies the internal dynamics of the
network by causing additional spikes of inhibitory neurons, whose delaying
within-network inhibition outweighs the immediate advance caused by the
external excitation. These results explain how intrinsic heterogeneity allows
the PING rhythm to become synchronized with a periodic forcing or another PING
rhythm for a wider range in the mismatch of their frequencies. We demonstrate a
similar mechanism for the synchronization of ING rhythms. Our results identify
a potential function of neuronal heterogeneity in the synchronization of
coupled $\gamma$-rhythms, which may play a role in neural information transfer
via communication through coherence.
