Stabilizing function of skeletal muscles: an analytical investigation.

Stability is the ability of a system to return to its original state after a disturbance. Taking vertical oscillations of the centre of mass of a human bending his legs as an example we prove that the intrinsic mechanical properties of musculature can stabilize the oscillatory movement (preflex) without reflexive changes in activation. The human is represented by a model consisting of a massless two-segment linkage system (knee) topped by a point mass. Conditions for stability are calculated analytically based on the theory of Ljapunov and the results are illustrated by numerical examples. In order to guarantee a self-stabilizing ability of the muscle-skeletal system, the muscle properties such as force-length relationship, force-velocity relationship and the muscle geometry must be tuned to the geometric properties of the linkage system.     

基于比率的具有反馈控制的非自治捕食系统的动力学行为

讨论了一类基于比率的具有反馈控制的非自治捕食系统,所有的参数都是时滞的.先研究了该系统的一致持久性和全局渐近稳定性,并通过构造适当的Lyapunov函数,得到了系统存在惟一渐近稳定的正概周期解的充分性条件.最后,通过一个例子说明了结论的可行性.     

A new necessary condition for Turing instabilities

Reactivity (a.k.a initial growth) is necessary for diffusion driven instability (Turing instability). Using a notion of common Lyapunov function we show that this necessary condition is a special case of a more powerful (i.e. tighter) necessary condition. Specifically, we show that if the linearised reaction matrix and the diffusion matrix share a common Lyapunov function, then Turing instability is not possible. The existence of common Lyapunov functions is readily checked using semi-definite programming. We apply this result to the Gierer-Meinhardt system modelling regenerative properties of Hydra, the Oregonator, to a host-parasite-hyperparasite system with diffusion and to a reaction-diffusion-chemotaxis model for a multi-species host-parasitoid community.     

The PSO family: deduction, stochastic analysis and comparison

The PSO algorithm can be physically interpreted as a stochastic damped mass-spring system: the so-called PSO continuous model. Furthermore, PSO corresponds to a particular discretization of the PSO continuous model. In this paper, we introduce a delayed version of the PSO continuous model, where the center of attraction might be delayed with respect to the particle trajectories. Based on this mechanical analogy, we derive a family of PSO versions. For each member of this family, we deduce the first and second order stability regions and the corresponding spectral radius. As expected, the PSO-family algorithms approach the PSO continuous model (damped-mass-spring system) as the size of the time step goes to zero. All the family members are linearly isomorphic when they are derived using the same delay parameter. If the delay parameter is different, the algorithms have corresponding stability regions of any order, but they differ in their respective force terms. All the PSO versions perform fairly well in a very broad area of the parameter space (inertia weight and local and global accelerations) that is close to the border of the second order stability regions and also to the median lines of the first order stability regions where no temporal correlation between trajectories exists. In addition, these regions are fairly similar for different benchmark functions. When the number of parameters of the cost function increases, these regions move towards higher inertia weight values (w=1) and lower total mean accelerations where the temporal covariance between trajectories is positive. Finally, the comparison between different PSO versions results in the conclusion that the centered version (CC-PSO) and PSO have the best convergence rates. Conversely, the centered-progressive version (CP-PSO) has the greatest exploratory capabilities. These features have to do with the way these algorithms update the velocities and positions of particles in the swarm. Knowledge of their respective dynamics can be used to propose a family of simple and stable algorithms, including hybrid versions.     

Hopf and torus bifurcations,torus destruction and chaos in population biology

One of the simplest population biological models displaying a Hopf bifurcation is the Rosenzweig–MacArthur model with Holling type II response function as essential ingredient. In seasonally forced versions the fixed point on one side of the Hopf bifurcation becomes a limit cycle and the Hopf limit cycle on the other hand becomes a torus, hence the Hopf bifurcation becomes a torus bifurcation, and via torus destruction by further increasing relevant parameters can follow deterministic chaos. We investigate this route to chaos also in view of stochastic versions, since in real world systems only such stochastic processes would be observed.However, the Holling type II response function is not directly related to a transition from one to another population class which would allow a stochastic version straight away. Instead, a time scale separation argument leads from a more complex model to the simple 2 dimensional Rosenzweig–MacArthur model, via additional classes of food handling and predators searching for prey. This extended model allows a stochastic generalization with the stochastic version of a Hopf bifurcation, and ultimately also with additional seasonality allowing a torus bifurcation under stochasticity.Our study shows that the torus destruction into chaos with positive Lyapunov exponents can occur in parameter regions where also the time scale separation and hence stochastic versions of the model are possible. The chaotic motion is observed inside Arnol’d tongues of rational ratio of the forcing frequency and the eigenfrequency of the unforced Hopf limit cycle.Such torus bifurcations and torus destruction into chaos are also observed in other population biological systems, and were for example found in extended multi-strain epidemiological models on dengue fever. To understand such dynamical scenarios better also under noise the present low dimensional system can serve as a good study case.     

A Lyapunov function for piecewise-independent differential equations: stability of the ideal free distribution in two patch environments

In this article we construct Lyapunov functions for models described by piecewise-continuous and independent differential equations. Because these models are described by discontinuous differential equations, the theory of Lyapunov functions for smooth dynamical systems is not applicable. Instead, we use a geometrical approach to construct a Lyapunov function. Then we apply the general approach to analyze population dynamics describing exploitative competition of two species in a two-patch environment. We prove that for any biologically meaningful parameter combination the model has a globally stable equilibrium and we analyze this equilibrium with respect to parameters.      

Mathematical analysis of a metapopulation model with space-limited recruitment

We investigate a mathematical aspect of a multi-species' sessile metapopulation model with space-limited recruitment proposed by Iwasa et al. in 1986. We define some basic reproduction numbers to show the threshold condition for the stability of trivial steady state and the existence of coexistent steady state. We show the existence of steady state where all species exist when some reproduction numbers are greater than one by the fixed point theorem. And we construct the Lyapunov function to show the global stability of trivial steady state when some basic reproduction numbers are not greater than one.     

Effects of deterministic and random refuge in a prey-predator model with parasite infection

Most natural ecosystem populations suffer from various infectious diseases and the resulting host-pathogen dynamics is dependent on host's characteristics. On the other hand, empirical evidences show that for most host pathogen systems, a part of the host population always forms a refuge. To study the role of refuge on the host-pathogen interaction, we study a predator-prey-pathogen model where the susceptible and the infected prey can undergo refugia of constant size to evade predator attack. The stability aspects of the model system is investigated from a local and global perspective. The study reveals that the refuge sizes for the susceptible and the infected prey are the key parameters that control possible predator extinction as well as species co-existence. Next we perform a global study of the model system using Lyapunov functions and show the existence of a global attractor. Finally we perform a stochastic extension of the basic model to study the phenomenon of random refuge arising from various intrinsic, habitat-related and environmental factors. The stochastic model is analyzed for exponential mean square stability. Numerical study of the stochastic model shows that increasing the refuge rates has a stabilizing effect on the stochastic dynamics.     

Dynamics of coupled thalamocortical modules

We develop a model of thalamocortical dynamics using a shared population of thalamic neurons to couple distant cortical regions. Behavior of the model is determined as a function of the connection strengths with shared and unshared populations in the thalamus, either within a relay nucleus or the reticular nucleus. When the coupling is via the reticular nucleus, we locate solutions of the model where distant cortical regions maintain the same activity level, and regions where one region maintains an elevated activity level, suppressing activity in the other. We locate and investigate a region where both types of solutions exist and are stable, yielding a mechanism for spontaneous changes in global activity patterns. Power spectra and coherence are computed, and marked differences in the coherence are found between the two kinds of modes. When, on the other hand, the coupling is via a shared relay nuclei, the features seen with the reticular coupling are absent. These considerations suggest a role for the reticular nucleus in modulating long distance cortical communication.     

Nonlinear dynamics and bifurcations of a coupled oscillator model for calling behavior of Japanese tree frogs (Hyla japonica)

Synchronization has been observed in various systems, including living beings. In a previous study, we reported a new phenomenon with antisynchronization in calling behavior of two interacting Japanese tree frogs. In this paper, we theoretically analyse nonlinear dynamics in a system of three coupled oscillators, which models three interacting frogs, where the oscillators of each pair have the property of antisynchronization; in particular, we perform bifurcation analysis and Lyapunov function analysis.     

Nonlinear dynamics and bifurcations of a coupled oscillator model for calling behavior of Japanese tree frogs (Hyla japonica)

Synchronization has been observed in various systems, including living beings. In a previous study, we reported a new phenomenon with antisynchronization in calling behavior of two interacting Japanese tree frogs. In this paper, we theoretically analyse nonlinear dynamics in a system of three coupled oscillators, which models three interacting frogs, where the oscillators of each pair have the property of antisynchronization; in particular, we perform bifurcation analysis and Lyapunov function analysis.     

Dendritic synchrony and transient dynamics in a coupled oscillator model of the dopaminergic neuron

Transient increases in spontaneous firing rate of mesencephalic dopaminergic neurons have been suggested to act as a reward prediction error signal. A mechanism previously proposed involves subthreshold calcium-dependent oscillations in all parts of the neuron. In that mechanism, the natural frequency of oscillation varies with diameter of cell processes, so there is a wide variation of natural frequencies on the cell, but strong voltage coupling enforces a single frequency of oscillation under resting conditions. In previous work, mathematical analysis of a simpler system of oscillators showed that the chain of oscillators could produce transient dynamics in which the frequency of the coupled system increased temporarily, as seen in a biophysical model of the dopaminergic neuron. The transient dynamics was shown to be consequence of a slow drift along an invariant subset of phase space, with rate of drift given by a Lyapunov function. In this paper, we show that the same mathematical structure exists for the full biophysical model, giving physiological meaning to the slow drift and the Lyapunov function, which is shown to describe differences in intracellular calcium concentration in different parts of the cell. The duration of transients was long, being comparable to the time constant of calcium disposition. These results indicate that brief changes in input to the dopaminergic neuron can produce long lasting firing rate transients whose form is determined by intrinsic cell properties.     

Extinction of species in age-structured, discrete noncooperative systems

Using a multi-species model of Ebenman for the dynamics of a discretely reproducing population that consists of noncooperation between juveniles and adults, we obtain exclusion principles by Lyapunov function methods. In the very general age-structured model, we show that, if there is an adult whose growth function is always larger than that of all the other species in the system, then it dominates the system by driving all the others to extinction. This result confirms a general folklore. We also develop a quasi-dominance concept and show that it implies the extinction of all the quasi-dominated species. The quasi-dominance concept applies even if there is no species whose adults always grow faster than all the others. In addition, a notion of weak dominance is developed. We show with specific examples that weak dominance does not necessarily imply extinction of species. If all the growth functions are exponential functions, then weak dominance is equivalent to quasi-dominance.     

Self-sustained pH oscillations in a compartmentalized enzyme reactor system

This work represents our continued effort toward fulfilling the need to discover a model system for experimental investigations of temporal oscillations in an enzyme-membrane system. In this paper, the regions in the parameter space where self-sustained pH oscillations can be induced for a compartmentalized enzyme reactor system, which consists of a well-stirred reactor, a reservoir and a membrane containing no enzyme, were determined via numerical simulation with two proteolytic enzymes: papain (EC 3.4.22.2) and alpha-chymotrypsin (EC 3.4.21.1). The sizes of the regions were qualitatively compared with those associated with enzymic membrane system. As a result, we found that the possibility of experimentally observing self-sustained oscillations in the compartmentalized papain reactor system, as well as in the papain-membrane system, is high. However, self-sustained pH oscillations are less likely in the compartmentalized alpha-chymotrypsin reactor system than in the alpha-chymotrypsin-membrane system.     

Pattern selection in a cooperative biochemical in vitro amplification system: the role of parasites

Previously, numerical simulations have shown that evolving systems can be stabilized against emerging parasites by pattern formation in spatially extended flow reactors. Hence, it can be argued that pattern formation is a prerequisite for any experimental investigation of the biochemical evolution of cooperative function. Here, we study a model of an experimental biochemical system for the cooperative in vitro amplification of DNA strands and show that emerging parasites can induce a complex pattern formation even when no pattern formation occurs without parasites. In an adiabatic approximation where the cooperative amplification reaction is assumed to adapt fast to slowly emerging parasites, the parasite concentration itself acts as a Steuer parameter for the selection of various complex patterns. Without such an adiabatic approximation only transient patterns emerge. As any species can grow for very low concentrations, the parasite is able to infect the entire reactor and the system is finally diluted out. In the experimental biochemical system, however, the species are individual molecules and the growth of spatially separated, non-infected regions becomes feasible. Hence a cutoff threshold for the minimal concentration is applied. In these simulations the otherwise lethal infection by parasites induces the formation of spatiotemporal spirals, and this spatial structure help the host and parasitoid species to survive together. These theoretical results describe an inherent property of cooperative reactions and have an important impact on experimental investigations on the molecular evolution and complex function in spatially extended reactors. Since the formation of the complex pattern is restricted either to a rather large cutoff value or a special choice of the kinetic parameters, we, however, conclude that the persistence of evolving cooperative amplification is not possible in a simple reaction-diffusion reactor. Experimental set-ups with patchy environments, e.g. biomolecular amplification in coupled microstructured flow chambers or in microemulsion, are eligible candidates for the observation of such a self-organized pattern selection.     

A comparison of two predator-prey models with Holling's type I functional response

In this paper, we analyze a laissez-faire predator-prey model and a Leslie-type predator-prey model with type I functional responses. We study the stability of the equilibrium where the predator and prey coexist by both performing a linearized stability analysis and by constructing a Lyapunov function. For the Leslie-type model, we use a generalized Jacobian to determine how eigenvalues jump at the corner of the functional response. We show, numerically, that our two models can both possess two limit cycles that surround a stable equilibrium and that these cycles arise through global cyclic-fold bifurcations. The Leslie-type model may also exhibit super-critical and discontinuous Hopf bifurcations. We then present and analyze a new functional response, built around the arctangent, that smoothes the sharp corner in a type I functional response. For this new functional response, both models undergo Hopf, cyclic-fold, and Bautin bifurcations. We use our analyses to characterize predator-prey systems that may exhibit bistability.     

A mathematical study of a two-regional population growth model

The paper provides a mathematical study of a model of urban dynamics, adjusting to an ecological model proposed by Lotka and Volterra. The model is a system of two first-order non-linear ordinary differential equations. The study proposed here completes the original proof by using the main tools such as a Lyapunov function.     

State-dependent impulsive models of integrated pest management (IPM) strategies and their dynamic consequences

A state-dependent impulsive model is proposed for integrated pest management (IPM). IPM involves combining biological, mechanical, and chemical tactics to reduce pest numbers to tolerable levels after a pest population has reached its economic threshold (ET). The complete expression of an orbitally asymptotically stable periodic solution to the model with a maximum value no larger than the given ET is presented, the existence of which implies that pests can be controlled at or below their ET levels. We also prove that there is no periodic solution with order larger than or equal to three, except for one special case, by using the properties of the LambertW function and Poincare map. Moreover, we show that the existence of an order two periodic solution implies the existence of an order one periodic solution. Various positive invariant sets and attractors of this impulsive semi-dynamical system are described and discussed. In particular, several horseshoe-like attractors, whose interiors can simultaneously contain stable order 1 periodic solutions and order 2 periodic solutions, are found and the interior structure of the horseshoe-like attractors is discussed. Finally, the largest invariant set and the sufficient conditions which guarantee the global orbital and asymptotic stability of the order 1 periodic solution in the meaningful domain for the system are given using the Lyapunov function. Our results show that, in theory, a pest can be controlled such that its population size is no larger than its ET by applying effects impulsively once, twice, or at most, a finite number of times, or according to a periodic regime. Moreover, our theoretical work suggests how IPM strategies could be used to alter the levels of the ET in the farmers' favour.