Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Stability analysis”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40Linked to original sources

Slow-fast dynamics in a model of population and resource growth.

"Models of the interaction of population, the economy, and the environment often contain nonlinear functional relationships and variables that move at different speeds. These properties foster apparent unpredictabilities in system behaviour. Using a simple deterministic model of demographic, economic and environmental interactions we illustrate the usefulness of geometric singular perturbation theory in environmental population economics. In contrast to local stability analysis, the theory of slow-fast dynamics helps to gain new insights into the global behaviour of the system. In particular, the knowledge of the basins of attraction of the stationary states enables one to determine the regions of sustainable future paths of resources and population." (SUMMARY IN FRE)

Conservation of Natural Resources↗

Dynamics of neural populations: stability and synchrony.

A population formulation of neuronal activity is employed to study an excitatory network of (spiking) neurons receiving external input as well as recurrent feedback. At relatively low levels of feedback, the network exhibits time stationary asynchronous behavior. A stability analysis of this time stationary state leads to an analytical criterion for the critical gain at which time asynchronous behavior becomes unstable. At instability the dynamics can undergo a supercritical Hopf bifurcation and the population passes to a synchronous state. Under different conditions it can pass to synchrony through a subcritical Hopf bifurcation. And at high gain a network can reach a runaway state, in finite time, after which the network no longer supports bounded solutions. The introduction of time delayed feedback leads to a rich range of phenomena. For example, for a given external input, increasing gain produces transition from asynchrony, to synchrony, to asynchrony and finally can lead to divergence. Time delay is also shown to strongly mollify the amplitude of synchronous oscillations. Perhaps, of general importance, is the result that synchronous behavior can exist only for a narrow range of time delays, which range is an order of magnitude smaller than periods of oscillation.

Action Potentials↗

Natural selection and Y-linked polymorphism.

Several population genetic models allowing natural selection to act on Y-linked polymorphism are examined. The first incorporates pleiotropic effects of a Y-linked locus, such that viability, segregation distortion, fecundity and sexual selection can all be determined by one locus. It is shown that no set of selection parameters can maintain a stable Y-linked polymorphism. Interaction with the X chromosome is allowed in the next model, with viabilities determined by both X- and Y-linked factors. This model allows four fixation equilibria, two equilibria with X polymorphism and a unique point with both X- and Y-linked polymorphism. Stability analysis shows that the complete polymorphism is never stable. When X- and Y-linked loci influence meiotic drive, it is possible to have all fixation equilibria simultaneously unstable, and yet there is no stable interior equilibrium. Only when viability and meiotic drive are jointly affected by both X- and Y-linked genes can a Y-linked polymorphism be maintained. Unusual dynamics, including stable limit cycles, are generated by this model. Numerical simulations show that only a very small portion of the parameter space admits Y polymorphism, a result that is relevant to the interpretation of levels of Y-DNA sequence variation in natural populations.

Biometry↗

The regulation of an age-structured population by a fatal disease.

A model describing the effect of a fatal disease on an age-structured population which would otherwise grow is presented and analysed. If the disease is capable of regulating host numbers, there is an endemic steady age distribution (SAD), for which an analytic expression is obtained under some simplifying assumptions. The ability of the disease to regulate the population depends on a parameter R(alpha), which is defined in terms of the given age-dependent birth and death rates, and where alpha is the age-dependent disease-induced death rate. If R(alpha) < 1 the endemic SAD is attained, while R(alpha) > 1 means the disease cannot control the population's size. The number R(0) is the expected number of offspring produced by each individual in the absence of the disease; for a growing population we require R(0) > 1. A stability analysis is also performed and it is conjectured that the endemic SAD is locally asymptotically stable whenever it is attained. This is demonstrated explicitly for a very simple example where all rates are taken as constant.

Age Factors↗

Interactions between a uniformly proliferating tumour and its surroundings: uniform material properties.

The stability of a planar tumour growing into neighbouring tissue is examined and, in particular, its dependence on the properties of the tumour and of the surrounding material studied. An abundant supply of nutrient is assumed, so the proliferation of cells is uninhibited (resulting in exponential growth). We consider two possible constitutive relations. Darcy's law and Stokes flow, in describing the deformation of the tissue and the resulting model takes the form of a coupled system comprising a nonlinear reaction-diffusion-convection equation for the tumour cell concentration and an elliptic system for the deformation and stress fields. Using a combination of linear-stability analysis, numerical methods and thin-film approximations. the evolution of the advancing tumour boundary is determined. It is shown that when the tumour and surrounding material properties are the same, a planar interface is always linearly unstable, with the Stokes flow problem being reducible to the Darcy one. We treat the subsequent (nonlinear) evolution and suggest possible extensions to this work.

Cell Division↗

An epidemic model with a density-dependent death rate.

This paper deals with a mathematical model for a disease where the death rate depends on the number of people in the population. This sort of model would be suitable for diseases in developing countries or for diseases amongst animal and insect populations. We also assume that the population under consideration is regulated by the disease so that there is a mortality induced by the disease. We use a compartmental model and perform an equilibrium and stability analysis to find that there is a threshold condition. If the threshold is exceeded, then there is a unique equilibrium with disease present which is locally stable to small perturbations. We conclude by looking at several specific models and seeing how the results relate to previous work.

Animals↗

Specific phosphorylation of Torpedo 43K rapsyn by endogenous kinase(s) with thiamine triphosphate as the phosphate donor.

43K rapsyn is a peripheral protein specifically associated with the nicotinic acetylcholine receptor (nAChR) present in the postsynaptic membrane of the neuromuscular junction and of the electrocyte, and is essential for its clustering. Here, we demonstrate a novel specific phosphorylation of 43K rapsyn by endogenous protein kinase(s) present in Torpedo electrocyte nAChR-rich membranes and identify thiamine triphosphate (TTP) as the phosphate donor. In the presence of Mg(2+) and [gamma-(32)P]-TTP, 43K rapsyn is specifically phosphorylated with a (32)P-half-maximal incorporation at approximately 5-25 microM TTP. The presence of TTP in the cytosol and of 43K rapsyn at the cytoplasmic face of the postsynaptic membrane, together with TTP-dependent phosphorylation of 43K rapsyn without added exokinases, suggests that TTP-dependent-43K-rapsyn phosphorylation may occur in vivo. In addition, phosphoamino acid and chemical stability analysis suggests that the residues phosphorylated are predominantly histidines. Inhibition of phosphorylation by Zn(2+) suggests a possible control of 43K rapsyn phosphorylation state by its zinc finger domain. Endogenous kinase(s) present in rodent brain membranes can also use [gamma-(32)P]-TTP as a phosphodonor. The use of a phosphodonor (TTP) belonging to the thiamine family but not to the classical (ATP, GTP) purine triphosphate family represents a novel phosphorylation pathway possibly important for synaptic proteins.

Adenosine Triphosphate↗

Effects of fluorescein on tear breakup time and on tear thinning time.

A simple, noninvasive method of in vivo precorneal tear film stability analysis is described. The method involves observing the first catoptric (first Purkinje) image. The effect of fluorescein on the results and the relation with standard tear breakup time (BUT) are presented.

Adult↗

A model of neuronal bursting using three coupled first order differential equations.

We describe a modification to our recent model of the action potential which introduces two additional equilibrium points. By using stability analysis we show that one of these equilibrium points is a saddle point from which there are two separatrices which divide the phase plane into two regions. In one region all phase paths approach a limit cycle and in the other all phase paths approach a stable equilibrium point. A consequence of this is that a short depolarizing current pulse will change an initially silent model neuron into one that fires repetitively. Addition of a third equation limits this firing to either an isolated burst or a depolarizing afterpotential. When steady depolarizing current was applied to this model it resulted in periodic bursting. The equations, which were initially developed to explain isolated triggered bursts, therefore provide one of the simplest models of the more general phenomenon of oscillatory burst discharge.

Action Potentials↗

Dynamics of two feline retroviruses (FIV and FeLV) within one population of cats.

We present a deterministic model of the dynamics of two microparasites simultaneously infecting a single host population. Both microparasites are feline retroviruses, namely Feline Immunodeficiency Virus (FIV) and Feline Leukaemia Virus (FeLV). The host is the domestic cat Felis catus. The model has been tested with data generated by a long-term study of several natural cat populations. Stability analysis and simulations show that, once introduced in a population, FIV spreads and is maintained, while FeLV can either disappear or persist. Moreover, introduction of both viruses into the population induces an equilibrium state for individuals of each different pathological class. The viruses never induce the extinction of the population. Furthermore, whatever the outcome for the host population (persistence of FIV only, or of both viruses), the global population size at the equilibrium state is only slightly lower than it would have been in the absence of the infections (i.e. at the carrying capacity), indicating a low impact of the viruses on the population. Finally, the impact of the diseases examined simultaneously is higher than the sum of the impact of the two diseases examined separately. This seems to be due to a higher mortality rate when both viruses infect a single individual.

Animals↗

Chaos in a periodically forced chemostat with algal mortality.

We study the possibility of chaotic dynamics in the externally driven Droop model. This model describes a phytoplankton population in a chemostat under periodic nutrient supply. Previously, it has been proven under very general assumptions, that such systems are not able to exhibit chaotic dynamics. We show that the simple introduction of algal mortality may lead to chaotic oscillations of algal density in the forced chemostat. Our numerical simulations show that the existence of chaos is intimately related to plankton overshooting in the unforced model. We provide a simple measure, based on stability analysis, for estimating the amount of overshooting. These findings are not restricted to the Droop model but also hold for other chemostat models with mortality. Our results suggest periodically driven chemostats as a simple model system for the experimental verification of chaos in ecology.

Eukaryota↗

Hydrodynamics of granular gases of viscoelastic particles.

Our study examines the long-time behaviour of a force-free granular gas of viscoelastic particles, for which the coefficient of restitution depends on the impact velocity, as it follows from the solution of the impact problem for viscoelastic spheres. Starting from the Boltzmann equation, we derived the hydrodynamic equations and obtained microscopic expressions for the transport coefficients in terms of the elastic and dissipative parameters of the particle material. We performed the stability analysis of the linearized set of equations and found that any inhomogeneities and vortices vanish after a long time and the system approaches the flow-free stage of homogeneous density. This behaviour is in contrast to that of a gas consisting of particles which interact via a (non-realistic) constant coefficient of restitution, for which inhomogeneities (clusters) and vortex patterns have been proven to arise and to continuously develop.

Colloids↗

Stability of Bödewadt flow.

The stability of the flow produced over an infinite stationary plane in a fluid rotating with uniform angular velocity at an infinite distance from the plane is considered. The basic flow is an exact solution of the Navier-Stokes equations making it amenable to theoretical study. An asymptotic investigation is presented in the limit of large Reynolds number. It is shown that the stationary spiral instabilities observed experimentally can be described by a linear inviscid stability analysis. The prediction obtained for the wave angle of the disturbances is found to agree well with the available experimental and numerical results.

Journal Article↗

Geometric visual hallucinations, Euclidean symmetry and the functional architecture of striate cortex.

This paper is concerned with a striking visual experience: that of seeing geometric visual hallucinations. Hallucinatory images were classified by Klüver into four groups called form constants comprising (i) gratings, lattices, fretworks, filigrees, honeycombs and chequer-boards, (ii) cobwebs, (iii) tunnels, funnels, alleys, cones and vessels, and (iv) spirals. This paper describes a mathematical investigation of their origin based on the assumption that the patterns of connection between retina and striate cortex (henceforth referred to as V1)-the retinocortical map-and of neuronal circuits in V1, both local and lateral, determine their geometry. In the first part of the paper we show that form constants, when viewed in V1 coordinates, essentially correspond to combinations of plane waves, the wavelengths of which are integral multiples of the width of a human Hubel-Wiesel hypercolumn, ca. 1.33-2 mm. We next introduce a mathematical description of the large-scale dynamics of V1 in terms of the continuum limit of a lattice of interconnected hypercolumns, each of which itself comprises a number of interconnected iso-orientation columns. We then show that the patterns of interconnection in V1 exhibit a very interesting symmetry, i.e. they are invariant under the action of the planar Euclidean group E(2)-the group of rigid motions in the plane-rotations, reflections and translations. What is novel is that the lateral connectivity of V1 is such that a new group action is needed to represent its properties: by virtue of its anisotropy it is invariant with respect to certain shifts and twists of the plane. It is this shift-twist invariance that generates new representations of E(2). Assuming that the strength of lateral connections is weak compared with that of local connections, we next calculate the eigenvalues and eigenfunctions of the cortical dynamics, using Rayleigh-Schrödinger perturbation theory. The result is that in the absence of lateral connections, the eigenfunctions are degenerate, comprising both even and odd combinations of sinusoids in straight phi, the cortical label for orientation preference, and plane waves in r, the cortical position coordinate. 'Switching-on' the lateral interactions breaks the degeneracy and either even or else odd eigenfunctions are selected. These results can be shown to follow directly from the Euclidean symmetry we have imposed. In the second part of the paper we study the nature of various even and odd combinations of eigenfunctions or planforms, the symmetries of which are such that they remain invariant under the particular action of E(2) we have imposed. These symmetries correspond to certain subgroups of E(2), the so-called axial subgroups. Axial subgroups are important in that the equivariant branching lemma indicates that when a symmetrical dynamical system becomes unstable, new solutions emerge which have symmetries corresponding to the axial subgroups of the underlying symmetry group. This is precisely the case studied in this paper. Thus we study the various planforms that emerge when our model V1 dynamics become unstable under the presumed action of hallucinogens or flickering lights. We show that the planforms correspond to the axial subgroups of E(2), under the shift-twist action. We then compute what such planforms would look like in the visual field, given an extension of the retinocortical map to include its action on local edges and contours. What is most interesting is that, given our interpretation of the correspondence between V1 planforms and perceived patterns, the set of planforms generates representatives of all the form constants. It is also noteworthy that the planforms derived from our continuum model naturally divide V1 into what are called linear regions, in which the pattern has a near constant orientation, reminiscent of the iso-orientation patches constructed via optical imaging. The boundaries of such regions form fractures whose points of intersection correspond to the well-known 'pinwheels'. To complete the study we then investigate the stability of the planforms, using methods of nonlinear stability analysis, including Liapunov-Schmidt reduction and Poincaré-Lindstedt perturbation theory. We find a close correspondence between stable planforms and form constants. The results are sensitive to the detailed specification of the lateral connectivity and suggest an interesting possibility, that the cortical mechanisms by which geometric visual hallucinations are generated, if sited mainly in V1, are closely related to those involved in the processing of edges and contours.

Animals↗

DevR-DevS is a bona fide two-component system of Mycobacterium tuberculosis that is hypoxia-responsive in the absence of the DNA-binding domain of DevR.

Two-component systems play a central role in the adaptation of pathogenic bacteria to the environment prevailing within host tissues. The genes encoding the response regulator DevR (Rv3133c/DosR) and the cytoplasmic portion (DevS(201)) of the histidine kinase DevS (Rv3132c/DosS), a putative two-component system of Mycobacterium tuberculosis, were cloned and the protein products were overexpressed, purified and refolded as N-terminally His(6)-tagged proteins from Escherichia coli. DevS(201) underwent autophosphorylation and participated in rapid phosphotransfer to DevR in a Mg(2+)-dependent manner. Chemical stability analysis and site-directed mutagenesis implicated the highly conserved residues His(395) and Asp(54) as the sites of phosphorylation in DevS and DevR, respectively. Mutations in Asp(8) and Asp(9) residues, postulated to form the acidic Mg(2+)-binding pocket, and the invariant Lys(104) of DevR, abrogated phosphoryl transfer from DevS(201) to DevR. DevR-DevS was thus established as a typical two-component regulatory system based on His-to-Asp phosphoryl transfer. Expression of the Rv3134c-devR-devS operon was induced at the RNA level in hypoxic cultures of M. tuberculosis H37Rv and was associated with an increase in the level of DevR protein. However, in a devR mutant strain expressing the N-terminal domain of DevR, induction was observed at the level of RNA expression but not at that of protein. DevS was translated independently of DevR and induction of devS transcripts was not associated with an increase in protein level in either wild-type or mutant strains, reflecting differential regulation of this locus during hypoxia.

Bacterial Proteins↗

Analytical approach for the Floquet theory of delay differential equations.

We present an analytical approach to deal with nonlinear delay differential equations close to instabilities of time periodic reference states. To this end we start with approximately determining such reference states by extending the Poincaré-Lindstedt and the Shohat expansions, which were originally developed for ordinary differential equations. Then we systematically elaborate a linear stability analysis around a time periodic reference state. This allows us to approximately calculate the Floquet eigenvalues and their corresponding eigensolutions by using matrix valued continued fractions.

Journal Article↗

Magnetic effect for electrochemically driven cellular convection.

Hydrodynamic instability analogous to Rayleigh-Bénard convection is observed in an electrolytic solution between two parallel copper wire electrodes. The laser interferometric technique can reveal the dissipation structure created by the motion of the fluid, which is controlled electrochemically. It is shown that under the presence of horizontal magnetic field the roll cells move horizontally along the electrodes. The electrochemically driven convection is simply controlled and monitored by setting and measuring the electrochemical parameters and forms many kinds of spatiotemporal patterns, especially under the magnetic field. The phenomenon is modeled by considering a Boussinesq fluid under a concentration gradient. The stability of the resulting equations is studied by linear stability analysis. The time dependent nonlinear system is investigated numerically and the main features of the experimental response are reproduced.

Journal Article↗

Colored-noise-induced chaotic array synchronization.

The effect of a time-correlated Gaussian noise on one-dimensional arrays consisting of diffusively coupled chaotic cells is analyzed. A resonance effect between the time scale of the chaotic attractor and the colored Gaussian noise has been found. As well, depending on the number of cells, coupling, and noise strength, an improvement of the synchronization or a poor synchronization between cells within the array can occur for some values of the time correlation. These nonlinear cooperative effects are studied in terms of a linear stability analysis around the uniform synchronized behavior.

Journal Article↗