This post is concerned with dynamical systems in the sense of systems of ordinary differential equations. It is well-known that as soon as the dimension of the system is at least three general dynamical systems can include very complicated behaviour (chaos, strange attractors etc.) which defy detailed understanding in the absence of extra restrictions. In many applications of dynamical systems to the real world (for instance to biology) we encounter the following situation. A dynamical system is postulated which involves a lot of unknowns (in particular more than three) and depends on many parameters. There is only limited information about the parameters which means that in reality we are confronted with a large class of systems. Moreover there is only a limited amount of information about which initial data lead to solutions relevant to the application of interest. This means that in effect a large class of solutions is involved when an attempt is made to study the problem. The door is wide open to the daunting complexity of general solutions of general systems. The standard reaction is to make some choice of parameters and initial data and solve for the evolution on the computer. Since only a finite number of cases can be computed there is no guarantee that what comes out is typical of what happens in general or relevant for the application.
It might be hoped that general properties of dynamical systems which arise in particular wide classes of applications might reduce some of these difficulties. An example of a class of this type is defined by the systems which describe the competition of species in ecology. In 1976 Stephen Smale (J. Math. Biol. 3, 5) showed in a three-page paper that the restriction to this type of system does not, on its own, lead to any improvement whatsoever. In more detail what he did was the following. He considered a system of ODE of the form where the
are smooth and
for all
. The inequality expresses the competitive property while the factor
on the right hand side ensures that if the population of a species is initially zero it stays zero. The region where all
are non-negative, which is the biologically interesting region, is invariant. Finally he assumes that the functions
are negative outside a large ball so that solutions cannot escape to infinity. Let
be the set where all
are non-negative and
. Smale shows that there is a dynamical system in the class under discussion here where
is an invariant subset, all solutions converge to
as
and the dynamical system on
is arbitrary. Thus any complication of a general dynamical system of dimension
can be built into a model of the dynamics of
competing species. There is no fundamental simplification if
is at least three. Soon afterwards an analogous result was proved for a class of systems in chemical physics describing continuous flow stirred tank reactors by Perelson and Wallwork (J. Chem. Phys. 66, 4390). A paper of Kaplan and Yorke (American Naturalist 111, 1030) mentions the relevance of Smale’s result to competitive exclusion, a topic mensioned in a previous post.
It is not my intention here to spread pessimism. I hope and believe that dynamical systems can be of immense value in many scientific applications and in biology and medicine in particular. I just want to emphasize that the act of writing down a dynamical system satisfying a few general assumptions is no guarantee that some interesting information about a scientific problem has been implemented. I would add that being able to calculate a few solutions on the computer is not a sufficient reason to be confident that something has been achieved. It should be a matter of priority to develop ways of ensuring that systems proposed do contain information.
May 2, 2009 at 11:41 am |
I discovered from the book ‘Monotone dynamical systems, an introduction to the theory of competitive and cooperative systems’ by H. Smith that there is a nice result for systems which are cooperative. These differ from competitive systems only by a sign. Reversing the sign of time converts one into the other and so it did not occur to me when I wrote the post that there could be a big difference. In the example of Smale the complicated dynamics occurs in an invariant manifold which describes the late-time behaviour of solutions. The early-time behaviour of cooperative systems is just as complicated. In contrast, the late-time behaviour of cooperative systems has better properties. In that case the
-limit set of a generic solution consists of stationary points.
June 11, 2009 at 6:32 pm |
Hello. Your last sentence reads “It should be a matter of priority to develop ways of ensuring that systems proposed do contain information.”
Along your blog entry, you never define or propose a definition for information content, a some sort of scale that measure how much information a systems of ODE’s is potentially capable to have.
Is there such a “concept” for ODEs, i. e. a concept of information independent of how well the ODEs solutions “fit” the experimental data?
June 12, 2009 at 9:29 am |
In this post I was using the word ‘information’ in an informal, intuitive way. I did not intend to suggest any technical use of the word and I am not aware of a relevant concept of this kind.
February 19, 2010 at 9:21 am |
It is implicit in the above remarks about the difference between competitive and cooperative systems that the variables
are assumed positive. Otherwise systems could be transformed from one class to the other just by changing the sign of all the
.
January 26, 2011 at 8:11 am |
[…] this post I will explore some connections between two themes I have previously written about, dynamical systems empty of information and chemical reaction network theory. In the first of these posts I mentioned a paper of Perelson […]