Chemical self-replication of oligonucleotides and helical peptides show the so-called square root rate law. Based on this rate we extend our previous work on ideal replicators to include the square root rate and other possible nonlinearities, which we couple with an enzimatic sink. Although the nonlinearity is necessary for complex dynamics, the nature of the sink is the essential feature in the mechanism that allows temporal and spatial patterns. We obtain exact general relations for the Poincare-Adronov-Hopf and Turing bifurcations, and our generalized results include the Higgins, autocatalator, and templator models as specific cases.

The authors consider a minimal cross-catalytic self-replicating system of only two cross-catalytic templates that mimics the R3C ligase ribozyme system of Dong-Eu and Joyce [Chem. Biol. 11, 1505 (2004)]. This system displays considerably more complex dynamics than its self-replicating counterpart. In particular, the authors discuss the Poincare-Andronov-Hopf bifurcation, canard transitions, excitability, and hysteresis that yield birhythmicity between simple and complex oscillations.