Towards a Mathematical Description of Biodiversity Evolution

We outline in this work a mathematical description of biodiversity evolution based on a second-order differential equation (also known as the “inertial/Galilean view”). After discussing the motivations and explicit forms of the simplest “forces”, we are lead to an equation analogue to a harmonic oscillator. The known solutions for the homogeneous problem are then tentatively related to the biodiversity curves of Sepkoski and Alroy<i> et al.</i>, suggesting mostly an inertial behavior of the time evolution of the number of genera and a quadratic behavior in some long-term evolution after extinction events. We present the Green function for the dynamical system and apply it to the description of the recovery curve after the Permo-Triassic extinction, as recently analyzed by Burgess, Bowring and Shen. Even though the agreement is not satisfactory, we point out direct connections between observed drop times after massive extinctions and mathematical constants and discuss why the failure ensues, suggesting a more complex form of the second-order mathematical description.