Dynamics of conduction blocks in a model of paced cardiac tissue

We study numerically the dynamics of conduction blocks using a detailed electrophysiological model. We find that this dynamics depends critically on the size of the paced region. Small pacing regions lead to stationary conduction blocks while larger pacing regions can lead to conduction blocks that travel periodically towards the pacing region. We show that this size-dependence dynamics can lead to a novel arrhythmogenic mechanism. Furthermore, we show that the essential phenomena can be captured in a much simpler coupled-map model.Comment: 8 pages 6 figure

A probabilistic peridynamic framework with an application to the study of the statistical size effect

This is the final version. Available from Elsevier via the DOI in this record. Data availability: No data was used for the research described in the article.The high computational expense of peridynamic models remains a major limitation, hindering ‘outer-loop’ applications that require a large number of simulations, for example, uncertainty quantification. This contribution presents a framework that makes such computations feasible. By employing a Multilevel Monte Carlo framework, where the majority of simulations are performed using a coarse mesh, and performing relatively few simulations using a fine mesh, a significant reduction in computational cost can be realised, and statistics of structural failure can be estimated. The maximum observed speed-up factor is 16 when compared to a standard Monte Carlo estimator, thus enabling the efficient forward propagation of uncertain parameters in a computationally expensive peridynamic model. Furthermore, the multilevel method provides an estimate of both the discretisation error and sampling error, thereby improving confidence in numerical predictions. The performance of the approach is demonstrated through an examination of the statistical size effect in quasi-brittle materials.Turing Fellowshi

Coexisting Pulses in a Model for Binary-Mixture Convection

We address the striking coexistence of localized waves (`pulses') of different lengths which was observed in recent experiments and full numerical simulations of binary-mixture convection. Using a set of extended Ginzburg-Landau equations, we show that this multiplicity finds a natural explanation in terms of the competition of two distinct, physical localization mechanisms; one arises from dispersion and the other from a concentration mode. This competition is absent in the standard Ginzburg-Landau equation. It may also be relevant in other waves coupled to a large-scale field.Comment: 5 pages revtex with 4 postscript figures (everything uuencoded

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

Identifying elastoplastic parameters with Bayes' theorem considering output error, input error and model uncertainty

We discuss Bayesian inference for the identification of elastoplastic material parameters. In addition to errors in the stress measurements, which are commonly considered, we furthermore consider errors in the strain measurements. Since a difference between the model and the experimental data may still be present if the data is not contaminated by noise, we also incorporate the possible error of the model itself. The three formulations to describe model uncertainty in this contribution are: (1) a random variable which is taken from a normal distribution with constant parameters, (2) a random variable which is taken from a normal distribution with an input-dependent mean, and (3) a Gaussian random process with a stationary covariance function. Our results show that incorporating model uncertainty often, but not always, improves the results. If the error in the strain is considered as well, the results improve even more

Phase synchronization and noise-induced resonance in systems of coupled oscillators

We study synchronization and noise-induced resonance phenomena in systems of globally coupled oscillators, each possessing finite inertia. The behavior of the order parameter, which measures collective synchronization of the system, is investigated as the noise level and the coupling strength are varied, and hysteretic behavior is manifested. The power spectrum of the phase velocity is also examined and the quality factor as well as the response function is obtained to reveal noise-induced resonance behavior.Comment: to be published in Phys. Rev.

Resonance between Noise and Delay

We propose here a stochastic binary element whose transition rate depends on its state at a fixed interval in the past. With this delayed stochastic transition this is one of the simplest dynamical models under the influence of ``noise'' and ``delay''. We demonstrate numerically and analytically that we can observe resonant phenomena between the oscillatory behavior due to noise and that due to delay.Comment: 4 pages, 5 figures, submitted to Phys.Rev.Lett Expanded and Added Reference

Structure And Dynamics Of Modulated Traveling Waves In Cellular Flames

We describe spatial and temporal patterns in cylindrical premixed flames in the cellular regime, $Le < 1$, where the Lewis number $Le$ is the ratio of thermal to mass diffusivity of a deficient component of the combustible mixture. A transition from stationary, axisymmetric flames to stationary cellular flames is predicted analytically if $Le$ is decreased below a critical value. We present the results of numerical computations to show that as $Le$ is further decreased traveling waves (TWs) along the flame front arise via an infinite-period bifurcation which breaks the reflection symmetry of the cellular array. Upon further decreasing $Le$ different kinds of periodically modulated traveling waves (MTWs) as well as a branch of quasiperiodically modulated traveling waves (QPMTWs) arise. These transitions are accompanied by the development of different spatial and temporal symmetries including period doublings and period halvings. We also observe the apparently chaotic temporal behavior of a disordered cellular pattern involving creation and annihilation of cells. We analytically describe the stability of the TW solution near its onset+ using suitable phase-amplitude equations. Within this framework one of the MTW's can be identified as a localized wave traveling through an underlying stationary, spatially periodic structure. We study the Eckhaus instability of the TW and find that in general they are unstable at onset in infinite systems. They can, however, become stable for larger amplitudes.Comment: to appear in Physica D 28 pages (LaTeX), 11 figures (2MB postscript file

Effect of External Noise Correlation in Optical Coherence Resonance

Coherence resonance occurring in semiconductor lasers with optical feedback is studied via the Lang-Kobayashi model with external non-white noise in the pumping current. The temporal correlation and the amplitude of the noise have a highly relevant influence in the system, leading to an optimal coherent response for suitable values of both the noise amplitude and correlation time. This phenomenon is quantitatively characterized by means of several statistical measures.Comment: RevTeX, 4 pages, 7 figure

Self-organized Vortex State in Two-dimensional Dictyostelium Dynamics

We present results of experiments on the dynamics of Dictyostelium discoideum in a novel set-up which constraints cell motion to a plane. After aggregation, the amoebae collect into round ''pancake" structures in which the cells rotate around the center of the pancake. This vortex state persists for many hours and we have explicitly verified that the motion is not due to rotating waves of cAMP. To provide an alternative mechanism for the self-organization of the Dictyostelium cells, we have developed a new model of the dynamics of self-propelled deformable objects. In this model, we show that cohesive energy between the cells, together with a coupling between the self-generated propulsive force and the cell's configuration produces a self-organized vortex state. The angular velocity profiles of the experiment and of the model are qualitatively similar. The mechanism for self-organization reported here can possibly explain similar vortex states in other biological systems.Comment: submitted to PRL; revised version dated 3/8/9