Abstract. A generally applicable, automatic method for the efficient computation of a database of global dynamics of a multiparameter dynamical system is introduced. An outer approximation of the dynamics for each subset of the parameter range is computed using rigorous numerical methods and is represented by means of a directed graph. The dynamics is then decomposed into the recurrent and gradient-like parts by fast combinatorial algorithms and is classified via Morse decompositions. These Morse decompositions are compared at adjacent parameter sets via continuation to detect possible changes in the dynamics. The Conley index is used to study the structure of isolated invariant sets associated with the computed Morse decompositions and to detect the existence of certain types of dynamics. The power of the developed method is illustrated with an application to the two-dimensional, density-dependent, Leslie population model. An interactive visualization of the results of computations discussed in the paper can be accessed at the website http://chomp.rutgers.edu/database/, and the source code of the software used to obtain these results has also been made freely available.

No abstract

In this note an idea of quasi-homogeneous normal form theory using new grading functions is introduced, the definition of N th order normal form is given and some sufficient conditions for the uniqueness of normal forms are derived. A special case of the unsolved problem in a paper of Baider and Sanders for the unique normal form of Bogdanov Takens singularities is solved. AcademicPress, Inc.

No abstract

No abstract

In this paper we investigate equilibria of continuous differential equation models of network dynamics. The motivation comes from gene regulatory networks where each directed edge represents either down-or upregulation, and is modeled by a sigmoidal nonlinear function. We show that the existence and stability of equilibria of a sigmoidal system is determined by a combinatorial analysis of the limiting switching system with piece-wise constant non-linearities. In addition, we describe a local decomposition of a switching system into a product of simpler cyclic feedback systems, where the cycles in each decomposition correspond to a particular subset of network loops.

Abstract. We develop a rigorous computational method for estimating the Lyapunov exponents in uniformly expanding regions of the phase space for onedimensional maps. Our method uses rigorous numerics and graph algorithms to provide results that are mathematically meaningful and can be achieved in an efficient way.

The purpose of this and forthcoming papers is to obtain a better understanding of complicated bifurcations for multiple homoclinic orbits. We shall take one particular type of codimension two homoclinic orbits called orbit-flip and study bifurcations to multiple homoclinic orbits appearing in a tubular neighborhood of the original orbit-flip. The main interest of the present paper lies in the occurrence of successive homoclinic doubling bifurcations under an appropriate condition, which is a part of the entire bifurcation for multiple homoclinic orbits. Since this is a totally global bifurcation, we need the aid of numerical experiments for which we must choose a concrete set of ordinary differential equations that exhibits the desired bifurcation. In this paper we employ a family of continuous piecewise-linear vector fields for such a model equation. In order to explain the cascade of homoclinic doubling bifurcations theoretically, we also derive a two-parameter family of unimodal maps as a singular limit of the Poincaré maps along homoclinic orbits. We locate bifurcation curves for this family of unimodal maps in the two-dimensional parameter space, which basically agree with those for the piecewise-linear vector fields. In particular, we show, using a standard technique from the theory of unimodal maps, that there exists an infinite sequence of doubling bifurcations which corresponds to the sequence of homoclinic doubling bifurcations for the piecewise-linear vector fields described above. Since our unimodal map has a singularity at a boundary point of its domain of definition, the doubling bifurcation is slightly different from that for standard quadratic unimodal maps, for instance the Feigenbaum constant associated with the accumulation of the doubling bifurcations is different from the standard value 4.6692.…