Emergent global oscillations in heterogeneous excitable media: The example of pancreatic beta cells

Using the standard van der Pol-FitzHugh-Nagumo excitable medium model I demonstrate a novel generic mechanism, diversity, that provokes the emergence of global oscillations from individually quiescent elements in heterogeneous excitable media. This mechanism may be operating in the mammalian pancreas, where excitable beta cells, quiescent when isolated, are found to oscillate when coupled despite the absence of a pacemaker region.Comment: See home page http://lec.ugr.es/~julya

Depinning of kinks in a Josephson-junction ratchet array

We have measured the depinning of trapped kinks in a ratchet potential using a fabricated circular array of Josephson junctions. Our ratchet system consists of a parallel array of junctions with alternating cell inductances and junctions areas. We have compared this ratchet array with other circular arrays. We find experimentally and numerically that the depinning current depends on the direction of the applied current in our ratchet ring. We also find other properties of the depinning current versus applied field, such as a long period and a lack of reflection symmetry, which we can explain analytically.Comment: to be published in PR

Dynamical phase diagram of the dc-driven underdamped Frenkel-Kontorova chain

Multistep dynamical phase transition from the locked to the running state of atoms in response to a dc external force is studied by MD simulations of the generalized Frenkel-Kontorova model in the underdamped limit. We show that the hierarchy of transition recently reported [Braun et al, Phys. Rev. Lett. 78, 1295 (1997)] strongly depends on the value of the friction constant. A simple phenomenological explanation for the friction dependence of the various critical forces separating intermediate regimes is given.Comment: 12 Revtex Pages, 4 EPS figure

Resonance phenomena of a solitonlike extended object in a bistable potential

We investigate the dynamics of a soliton that behaves as an extended particle. The soliton motion in an effective bistable potential can be chaotic in a similar way as the Duffing oscillator. We generalize the concept of geometrical resonance to spatiotemporal systems and apply it to design a nonfeedback mechanism of chaos control using localized perturbations.We show the existence of solitonic stochastic resonance.Comment: 3 postscript figure

A quasi-diagonal approach to the estimation of Lyapunov spectra for spatio-temporal systems from multivariate time series

We describe methods of estimating the entire Lyapunov spectrum of a spatially extended system from multivariate time-series observations. Provided that the coupling in the system is short range, the Jacobian has a banded structure and can be estimated using spatially localised reconstructions in low embedding dimensions. This circumvents the ``curse of dimensionality'' that prevents the accurate reconstruction of high-dimensional dynamics from observed time series. The technique is illustrated using coupled map lattices as prototype models for spatio-temporal chaos and is found to work even when the coupling is not strictly local but only exponentially decaying.Comment: 13 pages, LaTeX (RevTeX), 13 Postscript figs, to be submitted to Phys.Rev.

Array-induced collective transport in the Brownian motion of coupled nonlinear oscillator systems

Brownian motion of an array of harmonically coupled particles subject to a periodic substrate potential and driven by an external bias is investigated. In the linear response limit (small bias), the coupling between particles may enhance the diffusion process, depending on the competition between the harmonic chain and the substrate potential. An analytical formula of the diffusion rate for the single-particle case is also obtained. In the nonlinear response regime, the moving kink may become phase-locked to its radiated phonon waves, hence the mobility of the chain may decrease as one increases the external force.Comment: 4 figures, to appear in Phys. Rev.

Infrared spectroscopy of phytochrome and model pigments

Fourier-transform infrared difference spectra between the red-absorbing and far-red-absorbing forms of oat phytochrome have been measured in H2O and 2H2O. The difference spectra are compared with infrared spectra of model compounds, i.e. the (5Z,10Z,15Z)- and (5Z,10Z,15E)-isomers of 2,3,7,8,12,13,17,18-octaethyl-bilindion (Et8-bilindion), 2,3-dihydro-2,3,7,8,12,13,17,18-octaethyl-bilindion (H2Et8-bilindion), and protonated H2Et8-bilindion in various solvents. The spectra of the model compounds show that only for the protonated forms can clear differences between the two isomers be detected. Since considerable differences are present between the spectra of Et8-bilindion and H2Et8-bilindion, it is concluded that only the latter compound can serve as a model system of phytochrome. The 2H2O effect on the difference spectrum of phytochrome supports the view that the chromophore in red-absorbing phytochrome is protonated and suggests, in addition, that it is also protonated in far-red-absorbing phytochrome. The spectra show that protonated carboxyl groups are influenced. The small amplitudes in the difference spectra exclude major changes of protein secondary structure

Resonant-Cavity-Induced Phase Locking and Voltage Steps in a Josephson Array

We describe a simple dynamical model for an underdamped Josephson junction array coupled to a resonant cavity. From numerical solutions of the model in one dimension, we find that (i) current-voltage characteristics of the array have self-induced resonant steps (SIRS), (ii) at fixed disorder and coupling strength, the array locks into a coherent, periodic state above a critical number of active Josephson junctions, and (iii) when $N_a$ active junctions are synchronized on an SIRS, the energy emitted into the resonant cavity is quadratic with $N_a$. All three features are in agreement with a recent experiment [Barbara {\it et al}, Phys. Rev. Lett. {\bf 82}, 1963 (1999)]}.Comment: 4 pages, 3 eps figures included. Submitted to PRB Rapid Com

Exploring hypotheses of the actions of TGF-beta 1 in epidermal wound healing using a 3D computational multiscale model of the human epidermis

In vivo and in vitro studies give a paradoxical picture of the actions of the key regulatory factor TGF-beta 1 in epidermal wound healing with it stimulating migration of keratinocytes but also inhibiting their proliferation. To try to reconcile these into an easily visualized 3D model of wound healing amenable for experimentation by cell biologists, a multiscale model of the formation of a 3D skin epithelium was established with TGF-beta 1 literature-derived rule sets and equations embedded within it. At the cellular level, an agent-based bottom-up model that focuses on individual interacting units ( keratinocytes) was used. This was based on literature-derived rules governing keratinocyte behavior and keratinocyte/ECM interactions. The selection of these rule sets is described in detail in this paper. The agent-based model was then linked with a subcellular model of TGF-beta 1 production and its action on keratinocytes simulated with a complex pathway simulator. This multiscale model can be run at a cellular level only or at a combined cellular/subcellular level. It was then initially challenged ( by wounding) to investigate the behavior of keratinocytes in wound healing at the cellular level. To investigate the possible actions of TGF-beta 1, several hypotheses were then explored by deliberately manipulating some of these rule sets at subcellular levels. This exercise readily eliminated some hypotheses and identified a sequence of spatial-temporal actions of TGF-beta 1 for normal successful wound healing in an easy-to-follow 3D model. We suggest this multiscale model offers a valuable, easy-to-visualize aid to our understanding of the actions of this key regulator in wound healing, and provides a model that can now be used to explore pathologies of wound healing

Dynamics of a Josephson Array in a Resonant Cavity

We derive dynamical equations for a Josephson array coupled to a resonant cavity by applying the Heisenberg equations of motion to a model Hamiltonian described by us earlier [Phys. Rev. B {\bf 63}, 144522 (2001); Phys. Rev. B {\bf 64}, 179902 (E)]. By means of a canonical transformation, we also show that, in the absence of an applied current and dissipation, our model reduces to one described by Shnirman {\it et al} [Phys. Rev. Lett. {\bf 79}, 2371 (1997)] for coupled qubits, and that it corresponds to a capacitive coupling between the array and the cavity mode. From extensive numerical solutions of the model in one dimension, we find that the array locks into a coherent, periodic state above a critical number of active junctions, that the current-voltage characteristics of the array have self-induced resonant steps (SIRS's), that when $N_a$ active junctions are synchronized on a SIRS, the energy emitted into the resonant cavity is quadratic in $N_a$, and that when a fixed number of junctions is biased on a SIRS, the energy is linear in the input power. All these results are in agreement with recent experiments. By choosing the initial conditions carefully, we can drive the array into any of a variety of different integer SIRS's. We tentatively identify terms in the equations of motion which give rise to both the SIRS's and the coherence threshold. We also find higher-order integer SIRS's and fractional SIRS's in some simulations. We conclude that a resonant cavity can produce threshold behavior and SIRS's even in a one-dimensional array with appropriate experimental parameters, and that the experimental data, including the coherent emission, can be understood from classical equations of motion.Comment: 15 pages, 10 eps figures, submitted to Phys. Rev.