The purpose of this post is to demonstrate some very beautiful (I think!) mathematics that arises form Darwinian evolutionary theory. It is a real shame that most courses and discussions dealing with evolution **never** introduce any type of mathematical formalism which is very strange, since at the most fundamental levels, evolution must also be governed by quantum mechanics and electromagnetism, from which chemistry and biochemistry arise via top-down and bottom-up causation. See this article by George Ellis for more on the role of top-down causation in the universe and the hierarchy of physical matter. Indeed, my personal belief is that if some biologists and evolutionary biologists like Dawkins, Coyne, and others took the time to explain evolution with some modicum of mathematical formalism to properly describe the underlying mechanics instead of using it as an opportunity to attack religious people, the world would be a much better place, and the dialogue between science and religion would be much more smooth and intelligible.

In this post today, I will describe some formalism behind the phenomena of *prebiotic evolution. *It turns out that there has been a very good book by Claudius Gros and understanding evolution as a complex dynamical system (dynamical systems theory is my main area of research), and the interested reader should check out his book for more details on what follows below.

We can for simplicity consider a quasispecies as a system of macromolecules that have the ability to carry information, and consider the dynamics of the concentrations of the constituent molecules as the following dynamical system:

,

where are the concentrations of molecules, is the autocatalytic self-replication rate, and are mutation rates.

From this, we can consider the following catalytic reaction equations:

,

,

are the concentrations, are the autocatalytic growth rates, and are the transmolecular catalytic rates. We choose such that

.

Clearly:

,

that is, this quick calculation shows that the total concentration remains constant.

Let us consider now the case of homogeneous interactions such that

, , ,

which leads to

,

which becomes

.

This is a one-dimensional ODE with the following invariant submanifolds:

,

.

With homogeneous interactions, the concentrations with the largest growth rates will dominate, so there exists a such that where

,

.

The quantities and are determined via normalization conditions that give us a system of equations:

,

.

**For large , we obtain the approximation**

**,**

**which is the number of surviving species.**

Clearly, this is non-zero for a finite catalytic rate . This shows the formation of a hypercycle of molecules/quasispecies.

These computations clearly should be taken with a grain of salt. As pointed out in several sources, hypercycles describe closed systems, but, life exists in an open system driven by an energy flux. But, the interesting thing is, despite this, the very last calculation shows that there is clear division between molecules which can be considered as a type of primordial life-form separated by these molecules belonging to the environment.