No abstract

We simulate two-dimensional lattices of pulse-coupled oscillators with random natural frequencies, resembling pacemaker cells in the heart. As coupling increases, this system seems to undergo two phase transitions in the thermodynamic limit. First, the largest cluster of frequency entrained oscillators becomes macroscopic. Second, all oscillators frequency entrain, except possibly some isolated ones. Between the two transitions, the system has features indicating self-organized criticality. To our knowledge, the first transition and the intermediate phase have never been observed before. It remains to be seen if they are generic for large lattices of limit cycle oscillators.

A long chain of pulse-coupled oscillators was studied. The oscillators interacted via a phase response curve similar to those obtained from pacemaker cells in the heart. The natural frequencies were random numbers from a distribution with finite bandwidth. Stable frequency-entrained states were shown always to exist above a critical coupling strength. Below the critical coupling, the probability to have such states was shown to be zero if the number of oscillators is infinite. This discontinuity establishes the existence of a phase transition in the thermodynamic limit. For weak coupling, clusters of frequency-entrained oscillators emerged. The cluster sizes were exponentially distributed, even when the critical coupling was approached. At this coupling, the mean cluster size diverged to infinity according to a power law. The standard deviation of the distribution of mean frequencies in the chain converged to zero, also according to a power law.

We study oscillator chains of the form phi; (k) = omega(k) +K[Gamma( phi(k-1) - phi(k) )+Gamma( phi(k+1) - phi(k) )], where phi(k) epsilon[0,2pi) is the phase of oscillator k. In the thermodynamic limit where the number of oscillators goes to infinity, for suitable choices of Gamma(x), we prove that there is a critical coupling strength K(c), above which a stable frequency-entrained state exists, but below which the probability is zero to have such a state. It is assumed that the natural frequencies are random with finite bandwidth. A crucial condition on Gamma(x) is that it is nonodd, i.e., mid R: Gamma(x)+Gamma(-x) mid R: not equal 0. The interest in the results comes from the fact that any chain of limit-cycle oscillators can be described by equations of the above form in the limits of weak coupling and narrow distribution of natural frequencies.

A real-space renormalization transformation is constructed for lattices of nonidentical oscillators with dynamics of the general form dvarphi_{k}/dt=omega_{k}+g summation operator_{l}f_{lk}(varphi_{l},varphi_{k}) . The transformation acts on ensembles of such lattices. Critical properties corresponding to a second-order phase transition toward macroscopic synchronization are deduced. The analysis is potentially exact but relies in part on unproven assumptions. Numerically, second-order phase transitions with the predicted properties are observed as g increases in two structurally different two-dimensional oscillator models. One model has smooth coupling f_{lk}(varphi_{l},varphi_{k})=phi(varphi_{l}-varphi_{k}) , where phi(x) is nonodd. The other model is pulse coupled, with f_{lk}(varphi_{l},varphi_{k})=delta(varphi_{l})phi(varphi_{k}) . Lower bounds for the critical dimensions for different types of coupling are obtained. For nonodd coupling, macroscopic synchronization cannot be ruled out for any dimension D> or =1 , whereas in the case of odd coupling, the well-known result that it can be ruled out for D<3 is regained.

No abstract

No abstract

We study and classify firing waves in two-dimensional oscillator lattices. To do so, we simulate a pulse-coupled oscillator model aimed to resemble a group of pacemaker cells in the heart. The oscillators are assigned random natural frequencies, and we focus on frequency entrained states. Depending on the initial condition, three types of wave landscapes are seen asymptotically. A concentric landscape contains concentric waves with one or more foci. Spiral landscapes contain one or more spiral waves. A mixed landscape contains both concentric and spiral waves. Mixed landscapes are only seen for moderate coupling strengths g, since for higher g, spiral waves have higher frequency than concentric waves, so that they cannot mix in frequency entrained states. If the initial condition is random, the probability to get a concentric landscape increases with increasing coupling strength g, but decreases with increasing lattice size. The g dependence of the probability enables hysteresis, where the system jumps between the two landscape types as g is continuously changed. For moderate g, spiral tips rotate around a suppressed oscillator that never fires. We call such an oscillator an oscillator defect. A spiral may also rotate around a point defect situated between the oscillators. In that case all oscillators fire at the entrained frequency. For larger g, a spiral tip either moves around a row of suppressed oscillators, a row defect, or around an open curve situated between the oscillators, which may be called a line defect. The length of a row or line defect increases with g. Our results may help understand sinus node reentry, where the natural pacemaker of the heart suddenly shifts to a higher frequency. Some of the observed phenomena seem generic, based on simulations of other models.