Beyond classical theories: An integrative mathematical model of mating dynamics and parental care.

Classical theories, such as Bateman's principle and Trivers' parental investment theory, attempted to explain the coevolution of sexual selection and parental care through simple verbal arguments. Since then, quantitative models have demonstrated that it is rarely that simple because many non-intuitive structures and non-linear relationships are actually at play. In this study, we propose a new standard for models of mating dynamics and parental care, emphasizing the clarity and use of mathematical and probabilistic arguments, the meaning of consistency conditions, and the key role of spatial densities and the law of mass action. We used adaptive dynamics to calculate the evolutionary trajectory of the total care duration. Our results clearly show how the outcomes of parental care evolution can be diverse, depending on the quantitative balance between a set of dynamical forces arising from relevant differences and conditions in the male and female populations. The intensity of sexual selection, synergy of care, care quality, and relative mortality rates during mating interactions and caring activities act as forces driving evolutionary transitions between uniparental and biparental care. Sexual selection reduces the care duration of the selected sex, uniparental care evolves in the sex that offers the higher care quality, higher mortality during mating interactions of one sex leads to more care by that sex, and higher mortality during caring activities of one sex favours the evolution of uniparental care in the other sex. Both synergy and higher overall mortality during mating interactions can stabilize biparental care when sexual selection reduces the care duration of the selected sex. We discuss how the interaction between these forces influences the evolution of care patterns, and how sex ratios can vary and be interpreted in these contexts. We also propose new directions for future developments of our integrative model, creating new comparable analyses that share the same underlying assumptions and dynamical frameworks.