114 research outputs found
Invariant submanifold for series arrays of Josephson junctions
We study the nonlinear dynamics of series arrays of Josephson junctions in
the large-N limit, where N is the number of junctions in the array. The
junctions are assumed to be identical, overdamped, driven by a constant bias
current and globally coupled through a common load. Previous simulations of
such arrays revealed that their dynamics are remarkably simple, hinting at the
presence of some hidden symmetry or other structure. These observations were
later explained by the discovery of (N - 3) constants of motion, each choice of
which confines the resulting flow in phase space to a low-dimensional invariant
manifold. Here we show that the dimensionality can be reduced further by
restricting attention to a special family of states recently identified by Ott
and Antonsen. In geometric terms, the Ott-Antonsen ansatz corresponds to an
invariant submanifold of dimension one less than that found earlier. We derive
and analyze the flow on this submanifold for two special cases: an array with
purely resistive loading and another with resistive-inductive-capacitive
loading. Our results recover (and in some instances improve) earlier findings
based on linearization arguments.Comment: 10 pages, 6 figure
Chapman-Enskog method and synchronization of globally coupled oscillators
The Chapman-Enskog method of kinetic theory is applied to two problems of
synchronization of globally coupled phase oscillators. First, a modified
Kuramoto model is obtained in the limit of small inertia from a more general
model which includes ``inertial'' effects. Second, a modified Chapman-Enskog
method is used to derive the amplitude equation for an O(2) Takens-Bogdanov
bifurcation corresponding to the tricritical point of the Kuramoto model with a
bimodal distribution of oscillator natural frequencies. This latter calculation
shows that the Chapman-Enskog method is a convenient alternative to normal form
calculations.Comment: 7 pages, 2-column Revtex, no figures, minor change
Synchronization in populations of globally coupled oscillators with inertial effects
A model for synchronization of globally coupled phase oscillators including
``inertial'' effects is analyzed. In such a model, both oscillator frequencies
and phases evolve in time. Stationary solutions include incoherent
(unsynchronized) and synchronized states of the oscillator population. Assuming
a Lorentzian distribution of oscillator natural frequencies, , both
larger inertia or larger frequency spread stabilize the incoherent solution,
thereby making harder to synchronize the population. In the limiting case
, the critical coupling becomes independent of
inertia. A richer phenomenology is found for bimodal distributions. For
instance, inertial effects may destabilize incoherence, giving rise to
bifurcating synchronized standing wave states. Inertia tends to harden the
bifurcation from incoherence to synchronized states: at zero inertia, this
bifurcation is supercritical (soft), but it tends to become subcritical (hard)
as inertia increases. Nonlinear stability is investigated in the limit of high
natural frequencies.Comment: Revtex, 36 pages, submit to Phys. Rev.
Solvable model of a phase oscillator network on a circle with infinite-range Mexican-hat-type interaction
We describe a solvable model of a phase oscillator network on a circle with
infinite-range Mexican-hat-type interaction. We derive self-consistent
equations of the order parameters and obtain three non-trivial solutions
characterized by the rotation number. We also derive relevant characteristics
such as the location-dependent distributions of the resultant frequencies of
desynchronized oscillators. Simulation results closely agree with the
theoretical ones
Dynamics of fully coupled rotators with unimodal and bimodal frequency distribution
We analyze the synchronization transition of a globally coupled network of N
phase oscillators with inertia (rotators) whose natural frequencies are
unimodally or bimodally distributed. In the unimodal case, the system exhibits
a discontinuous hysteretic transition from an incoherent to a partially
synchronized (PS) state. For sufficiently large inertia, the system reveals the
coexistence of a PS state and of a standing wave (SW) solution. In the bimodal
case, the hysteretic synchronization transition involves several states.
Namely, the system becomes coherent passing through traveling waves (TWs), SWs
and finally arriving to a PS regime. The transition to the PS state from the SW
occurs always at the same coupling, independently of the system size, while its
value increases linearly with the inertia. On the other hand the critical
coupling required to observe TWs and SWs increases with N suggesting that in
the thermodynamic limit the transition from incoherence to PS will occur
without any intermediate states. Finally a linear stability analysis reveals
that the system is hysteretic not only at the level of macroscopic indicators,
but also microscopically as verified by measuring the maximal Lyapunov
exponent.Comment: 22 pages, 11 figures, contribution for the book: Control of
Self-Organizing Nonlinear Systems, Springer Series in Energetics, eds E.
Schoell, S.H.L. Klapp, P. Hoeve
How to suppress undesired synchronization
It is delightful to observe the emergence of synchronization in the blinking
of fireflies to attract partners and preys. Other charming examples of
synchronization can also be found in a wide range of phenomena such as, e.g.,
neurons firing, lasers cascades, chemical reactions, and opinion formation.
However, in many situations the formation of a coherent state is not pleasant
and should be mitigated. For example, the onset of synchronization can be the
root of epileptic seizures, traffic congestion in communication networks, and
the collapse of constructions. Here we propose the use of contrarians to
suppress undesired synchronization. We perform a comparative study of different
strategies, either requiring local or total knowledge of the system, and show
that the most efficient one solely requires local information. Our results also
reveal that, even when the distribution of neighboring interactions is narrow,
significant improvement in mitigation is observed when contrarians sit at the
highly connected elements. The same qualitative results are obtained for
artificially generated networks as well as two real ones, namely, the Routers
of the Internet and a neuronal network
Desynchronizing effect of high-frequency stimulation in a generic cortical network model
Transcranial Electrical Stimulation (TCES) and Deep Brain Stimulation (DBS)
are two different applications of electrical current to the brain used in
different areas of medicine. Both have a similar frequency dependence of their
efficiency, with the most pronounced effects around 100Hz. We apply
superthreshold electrical stimulation, specifically depolarizing DC current,
interrupted at different frequencies, to a simple model of a population of
cortical neurons which uses phenomenological descriptions of neurons by
Izhikevich and synaptic connections on a similar level of sophistication. With
this model, we are able to reproduce the optimal desynchronization around
100Hz, as well as to predict the full frequency dependence of the efficiency of
desynchronization, and thereby to give a possible explanation for the action
mechanism of TCES.Comment: 9 pages, figs included. Accepted for publication in Cognitive
Neurodynamic
Multiple dynamical time-scales in networks with hierarchically nested modular organization
Many natural and engineered complex networks have intricate mesoscopic
organization, e.g., the clustering of the constituent nodes into several
communities or modules. Often, such modularity is manifested at several
different hierarchical levels, where the clusters defined at one level appear
as elementary entities at the next higher level. Using a simple model of a
hierarchical modular network, we show that such a topological structure gives
rise to characteristic time-scale separation between dynamics occurring at
different levels of the hierarchy. This generalizes our earlier result for
simple modular networks, where fast intra-modular and slow inter-modular
processes were clearly distinguished. Investigating the process of
synchronization of oscillators in a hierarchical modular network, we show the
existence of as many distinct time-scales as there are hierarchical levels in
the system. This suggests a possible functional role of such mesoscopic
organization principle in natural systems, viz., in the dynamical separation of
events occurring at different spatial scales.Comment: 10 pages, 4 figure
Dynamical aspects of mean field plane rotators and the Kuramoto model
The Kuramoto model has been introduced in order to describe synchronization
phenomena observed in groups of cells, individuals, circuits, etc... We look at
the Kuramoto model with white noise forces: in mathematical terms it is a set
of N oscillators, each driven by an independent Brownian motion with a constant
drift, that is each oscillator has its own frequency, which, in general,
changes from one oscillator to another (these frequencies are usually taken to
be random and they may be viewed as a quenched disorder). The interactions
between oscillators are of long range type (mean field). We review some results
on the Kuramoto model from a statistical mechanics standpoint: we give in
particular necessary and sufficient conditions for reversibility and we point
out a formal analogy, in the N to infinity limit, with local mean field models
with conservative dynamics (an analogy that is exploited to identify in
particular a Lyapunov functional in the reversible set-up). We then focus on
the reversible Kuramoto model with sinusoidal interactions in the N to infinity
limit and analyze the stability of the non-trivial stationary profiles arising
when the interaction parameter K is larger than its critical value K_c. We
provide an analysis of the linear operator describing the time evolution in a
neighborhood of the synchronized profile: we exhibit a Hilbert space in which
this operator has a self-adjoint extension and we establish, as our main
result, a spectral gap inequality for every K>K_c.Comment: 18 pages, 1 figur
- …