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.