We analyze a large system of globally coupled phase oscillators whose natural frequencies are bimodally distributed. The dynamics of this system has been the subject of long-standing interest. In 1984 Kuramoto proposed several conjectures about its behavior; ten years later, Crawford obtained the first analytical results by means of a local center manifold calculation. Nevertheless, many questions have remained open, especially about the possibility of global bifurcations. Here we derive the system's stability diagram for the special case where the bimodal distribution consists of two equally weighted Lorentzians. Using an ansatz recently discovered by Ott and Antonsen, we show that in this case the infinite-dimensional problem reduces exactly to a flow in four dimensions. Depending on the parameters and initial conditions, the long-term dynamics evolves to one of three states: incoherence, where all the oscillators are desynchronized; partial synchrony, where a macroscopic group of phase-locked oscillators coexists with a sea of desynchronized ones; and a standing wave state, where two counter-rotating groups of phase-locked oscillators emerge. Analytical results are presented for the bifurcation boundaries between these states. Similar results are also obtained for the case in which the bimodal distribution is given by the sum of two Gaussians.

A method to characterize dynamical interdependence among nonlinear systems is derived based on mutual nonlinear prediction. Systems with nonlinear correlation will show mutual nonlinear prediction when standard analysis with linear cross correlation might fail. Mutual nonlinear prediction also provides information on the directionality of the coupling between systems. Furthermore, the existence of bidirectional mutual nonlinear prediction in unidirectionally coupled systems implies generalized synchrony. Numerical examples studied include three classes of unidirectionally coupled systems: systems with identical parameters, nonidentical parameters, and stochastic driving of a nonlinear system. This technique is then applied to the activity of motoneurons within a spinal cord motoneuron pool. The interrelationships examined include single neuron unit firing, the total number of neurons discharging at one time as measured by the integrated monosynaptic reflex, and intracellular measurements of integrated excitatory postsynaptic potentials ͑EPSP's͒. Dynamical interdependence, perhaps generalized synchrony, was identified in this neuronal network between simultaneous single unit firings, between units and the population, and between units and intracellular EPSP's. ͓S1063-651X͑96͒04012-3͔

We design and analyze the dynamics of a large network of theta neurons, which are idealized type I neurons. The network is heterogeneous in that it includes both inherently spiking and excitable neurons. The coupling is global, via pulselike synapses of adjustable sharpness. Using recently developed analytical methods, we identify all possible asymptotic states that can be exhibited by a mean field variable that captures the network's macroscopic state. These consist of two equilibrium states that reflect partial synchronization in the network and a limit cycle state in which the degree of network synchronization oscillates in time. Our approach also permits a complete bifurcation analysis, which we carry out with respect to parameters that capture the degree of excitability of the neurons, the heterogeneity in the population, and the coupling strength (which can be excitatory or inhibitory). We find that the network typically tends toward the two macroscopic equilibrium states when the neuron's intrinsic dynamics and the network interactions reinforce one another. In contrast, the limit cycle state, bifurcations, and multistability tend to occur when there is competition among these network features. Finally, we show that our results are exhibited by finite network realizations of reasonable size.