We observe rare and recurrent large spiking events in a heterogeneous network of superconducting Josephson junctions (JJ) connected through a resistive load and driven by a radio-frequency (rf) current in addition to a constant bias. The intermittent large spiking events show characteristic features of extreme events (EE) since they are larger than a statistically defined significant height. Under the influence of repulsive interactions and an impact of heterogeneity of damping parameters, the network splits into three sub-groups of junctions, one in incoherent rotational, another incoherent librational motion and a third sub-group originating EE. We are able to scan the whole population of junctions with their distinctive individual dynamical features either in EE mode or non-EE mode in parameter space. EE migrates spatially from one to another sub-group of junctions depending upon the repulsive strength and the damping parameter. For weak repulsive coupling, all the junctions originate frequent large spiking events, in rotational motion when the average inter-spike-interval (ISI) is small, but it increases exponentially with repulsive interaction; it largely deviates from its exponential growth at a breakpoint where EE triggers in a sub-group of junctions. The probability density of inter-event-intervals (IEI) in the subgroup exhibits a Poisson distribution. EE originates via bubbling instability of in-phase synchronization.

To be announced