Triggering synchronized oscillations through arbitrarily weak diversity in close-to-threshold excitable media

It is shown that arbitrarily weak (frozen) heterogeneity can induce global synchronized oscillations in excitable media close to threshold. The work is carried out on networks of coupled van der Pol-FitzHugh-Nagumo oscillators. The result is shown to be robust against the presence of internal dynamical noise.Comment: 4 pages (RevTeX 3 style), 5 EPS figures, submitted to Phys. Rev. E (16 aug 2001

Modeling rhythmic patterns in the hippocampus

We investigate different dynamical regimes of neuronal network in the CA3 area of the hippocampus. The proposed neuronal circuit includes two fast- and two slowly-spiking cells which are interconnected by means of dynamical synapses. On the individual level, each neuron is modeled by FitzHugh-Nagumo equations. Three basic rhythmic patterns are observed: gamma-rhythm in which the fast neurons are uniformly spiking, theta-rhythm in which the individual spikes are separated by quiet epochs, and theta/gamma rhythm with repeated patches of spikes. We analyze the influence of asymmetry of synaptic strengths on the synchronization in the network and demonstrate that strong asymmetry reduces the variety of available dynamical states. The model network exhibits multistability; this results in occurrence of hysteresis in dependence on the conductances of individual connections. We show that switching between different rhythmic patterns in the network depends on the degree of synchronization between the slow cells.Comment: 10 pages, 9 figure

Greenhouse gas balance over thaw-freeze cycles in discontinuous zone permafrost

Peat in the discontinuous permafrost zone contains a globally significant reservoir of carbon that has undergone multiple permafrost-thaw cycles since the end of the mid-Holocene (~3700 years before present). Periods of thaw increase C decomposition rates which leads to the release of CO2 and CH4 to the atmosphere creating potential climate feedback. To determine the magnitude and direction of such feedback, we measured CO2 and CH4 emissions and modeled C accumulation rates and radiative fluxes from measurements of two radioactive tracers with differing lifetimes to describe the C balance of the peatland over multiple permafrost-thaw cycles since the initiation of permafrost at the site. At thaw features, the balance between increased primary production and higher CH4 emission stimulated by warmer temperatures and wetter conditions favors C sequestration and enhanced peat accumulation. Flux measurements suggest that frozen plateaus may intermittently (order of years to decades) act as CO2 sources depending on temperature and net ecosystem respiration rates, but modeling results suggest that—despite brief periods of net C loss to the atmosphere at the initiation of thaw—integrated over millennia, these sites have acted as net C sinks via peat accumulation. In greenhouse gas terms, the transition from frozen permafrost to thawed wetland is accompanied by increasing CO2 uptake that is partially offset by increasing CH4 emissions. In the short-term (decadal time scale) the net effect of this transition is likely enhanced warming via increased radiative C emissions, while in the long-term (centuries) net C deposition provides a negative feedback to climate warming

Control spiral wave dynamics using feedback signals from line detectors

We numerically study trajectories of spiral-wave-cores in excitable systems modulated proportionally to the integral of the activity on the straight line, several or dozens of equi-spaced measuring points on the straight line, the double-line and the contour-line. We show the single-line feedback results in the drift of core center along a straight line being parallel to the detector. An interesting finding is that the drift location in $y$ is a piecewise linear-increasing function of both the feedback line location and time delay. Similar trajectory occurs when replacing the feedback line with several or dozens of equi-spaced measuring points on the straight line. This allows to move the spiral core to the desired location along a chosen direction by measuring several or dozens of points. Under the double-line feedback, the shape of the tip trajectory representing the competition between the first and second feedback lines is determined by the distance of two lines. Various drift attractors in spiral wave controlled by square-shaped contour-line feedback are also investigated. A brief explanation is presented.Comment: 6 pages and 7 figures; Accepted for publication in EPL; Figs.5 and 6 are in JPG forma

Limits and dynamics of stochastic neuronal networks with random heterogeneous delays

Realistic networks display heterogeneous transmission delays. We analyze here the limits of large stochastic multi-populations networks with stochastic coupling and random interconnection delays. We show that depending on the nature of the delays distributions, a quenched or averaged propagation of chaos takes place in these networks, and that the network equations converge towards a delayed McKean-Vlasov equation with distributed delays. Our approach is mostly fitted to neuroscience applications. We instantiate in particular a classical neuronal model, the Wilson and Cowan system, and show that the obtained limit equations have Gaussian solutions whose mean and standard deviation satisfy a closed set of coupled delay differential equations in which the distribution of delays and the noise levels appear as parameters. This allows to uncover precisely the effects of noise, delays and coupling on the dynamics of such heterogeneous networks, in particular their role in the emergence of synchronized oscillations. We show in several examples that not only the averaged delay, but also the dispersion, govern the dynamics of such networks.Comment: Corrected misprint (useless stopping time) in proof of Lemma 1 and clarified a regularity hypothesis (remark 1

Noise Induced Coherence in Neural Networks

We investigate numerically the dynamics of large networks of $N$ globally pulse-coupled integrate and fire neurons in a noise-induced synchronized state. The powerspectrum of an individual element within the network is shown to exhibit in the thermodynamic limit ($N\to \infty$) a broadband peak and an additional delta-function peak that is absent from the powerspectrum of an isolated element. The powerspectrum of the mean output signal only exhibits the delta-function peak. These results are explained analytically in an exactly soluble oscillator model with global phase coupling.Comment: 4 pages ReVTeX and 3 postscript figure

Magnetic Field-Induced Condensation of Triplons in Han Purple Pigment BaCuSi2_2O6_6

Besides being an ancient pigment, BaCuSi$_2$O$_6$ is a quasi-2D magnetic insulator with a gapped spin dimer ground state. The application of strong magnetic fields closes this gap creating a gas of bosonic spin triplet excitations called triplons. The topology of the spin lattice makes BaCuSi$_2$O$_6$ an ideal candidate for studying the Bose-Einstein condensation of triplons as a function of the external magnetic field, which acts as a chemical potential. In agreement with quantum Monte Carlo numerical simulations, we observe a distinct lambda-anomaly in the specific heat together with a maximum in the magnetic susceptibility upon cooling down to liquid Helium temperatures.Comment: published on August 20, 200

Dynamics of lattice spins as a model of arrhythmia

We consider evolution of initial disturbances in spatially extended systems with autonomous rhythmic activity, such as the heart. We consider the case when the activity is stable with respect to very smooth (changing little across the medium) disturbances and construct lattice models for description of not-so-smooth disturbances, in particular, topological defects; these models are modifications of the diffusive XY model. We find that when the activity on each lattice site is very rigid in maintaining its form, the topological defects - vortices or spirals - nucleate a transition to a disordered, turbulent state.Comment: 17 pages, revtex, 3 figure

Structural Stability and Renormalization Group for Propagating Fronts

A solution to a given equation is structurally stable if it suffers only an infinitesimal change when the equation (not the solution) is perturbed infinitesimally. We have found that structural stability can be used as a velocity selection principle for propagating fronts. We give examples, using numerical and renormalization group methods.Comment: 14 pages, uiucmac.tex, no figure