Trends and Paradoxes of Competitive Evolution in the Predation Mechanism

: We give a series of numerical examples of competitive evolution in the predation system, showing in some cases how the choice is made to increase the efficiency of the predation mechanism (or other significant parameters) to the detriment of populations (both of prey and predators). We then develop the mathematical theory that enables us to understand the causality involved, and we identify a trend towards the emergence of the functional predation mechanism as such (and not of populations of the species involved). The realization of this trend only takes place when the conditions for it are offered by the hazards proposed to successive competitive choices. The logical structure of this trend is similar to that of the “tendency of rate of profit to fall” in certain economic models.


Introduction
It is well known that the causality involved in predation is highly complex, and that it is difficult to reconcile with general trends such as "survival of the fittest".Indeed, if we have predators feeding on prey, then it is clear that, from a demographic point of view, the presence of predators is always detrimental to the prey, whereas, from the predators' viewpoint, the consumption of prey is necessary for their survival, but, if excessive, then it leads to a scarcity of prey, which in turn leads to predator starvation.Clearly, moderation of the predation mechanism is fundamental to its stability, and perhaps to its longterm viability (you cannot kill the goose that lays the golden eggs. ..).It is also clear that time must play an explicit role in any causal description of these mechanisms.
The predator-prey system of population dynamics, of which there are several models, makes it possible to explain the various possible situations and to understand their causality, and above all, by placing two predators in competition, to understand the reasons for the evolutionary choices that follow.
But, on the other hand, a more complete system, known as predation-commensalism, is a model of the global economy that enables us to prove certain properties of the latter, notably TRPF (Tendency of Rate of Profit to Fall), which simply results from the competition between various firms, which, in order to increase their market share, reduce their profit margins.The consequence of TRPF is the pursuit of innovation and technological progress, which is precisely what allows the rate of profit to be reduced.(See [1] Chapter 9 and in particular Section 9.3).In this paper, we follow an analogous approach in the predator-prey system to identify the natural tendency resulting from repeated competitive evolution.Under plausible hypotheses (see Section 4), it appears that the natural tendency of this system is to increase its own efficiency.Clearly, this does not always correspond to the demographic growth of the predator and may even lead to its disappearance.This is a "tendency", a kind of "vocation" that only comes to fruition when the conditions for its realization are offered by the hazards of successive competitive evolutions (just as TRPF does not automatically lead to technological progress).It is nevertheless clear that the loss of efficiency is rejected by competitive evolution, which nonetheless accepts drastic population reductions (of predators and prey) under certain conditions.
Let us take a predation model for two species (of populations x and y, prey and predator, respectively),

𝑥
ℎ      ℎ where x' = f(x) is the equation for prey evolution in the absence of predators, for which we take the logistic equation: (a is the natural growth rate of x and P is the environmental carrying capacity).
The functional response h(x) is taken from Holling's type-II, in the form Remark: This predation function h(x) is essentially analogous to the widely used Holling's type-II functional response, (see for instance [1][2][3] or one of Holling's original articles [4]).Holling's type-II functional response resembles type-I for small x, but then gradually tends towards a constant (representing the predator's satiety).However, here we use the hyperbolic tangent, ℎ() = ℎ(/), instead of the algebraic Holling's type-II form, ℎ() = ()/( + ).Indeed, the hyperbolic tangent expresses the idea of proportionality capped by satiety much better than the algebraic expression.It often enables us to better understand the respective roles of proportionality for small x (chance encounters) and the satiety ceiling.
Within the framework of this model, there are six significant parameters: c (death rate of predators in the absence of prey), e (efficiency of the predation mechanism), b (satiety ceiling of predators), τ (conversion rate of consumed prey into predator population), a (natural growth rate of prey in the absence of predators and in an abundance of resources), and P (carrying capacity of the environment, i.e., population of prey alone that it can support).The last two parameters concern only the behavior of prey with their own resources.The first four concern the behavior of predators and their relationships with prey.
The parameters are fixed in all numerical simulations and are given by: Numerical simulations were performed in C language, using the Runge-Kutta-4 method, with an integration time step of 0.001.The initial values are all taken in the bassin [0.1; 3.9] × [0.1; 2.0] for x(t = 0) and y(t = 0), but, since we are interested in the properties of attractors, the transient is suppressed so that the attractor is clearly visible, thus eliminating the influence of initial values.
Note that we can always take τ = 1 by choosing appropriate units to measure the populations of x and y (the unit of y is the quantity of predators obtained by consuming a unit of x).We will always do so, except in cases where we need to compare several predators with different conversion rates, where a single choice is not possible.Moreover, generally speaking, the specific values of the parameters in each numerical computation were obviously chosen in order to enhance a particular facet of phenomena depending on many parameters (then impossible to describe by numerical computations).More precisely, the description in this Section 1 of the model for single predator is general and independent of the value of c (provided that 0 < c < b, otherwise the predation process is impossible).
The solutions and attractors in system (1) have different properties depending on the parameter values.It is useful to bear these in mind before tackling the problem of competition between two predators.
Apart from the properties (0,0) and (P,0) which are always equilibria, to find the equilibrium (x 0 ,y 0 ) internal to the domain of definition x > 0, y > 0, simply solve system (1) with the first members replaced by 0. From the second equation we have Equation ( 5) below and then from the first we get Equation ( 6): It is well known that, depending on the parameter values, this equilibrium may or may not be stable.In the latter case, the attractor is a periodic cycle surrounding the equilibrium point.It is worth recalling here that the evolution of the equilibrium point and this attractor as a function of the parameter e when the other parameters are fixed.Indeed, we shall see that the variation of e implies particularly significant qualitative properties.We shall also see (Section 4) that variations in e play a key role in the dynamics of evolution.
By fixing the values of the other parameters and varying e, we can easily see (see also [5] or [1] Chapter 6) thanks to Equations ( 5) and (6) with Equations ( 2) and ( 3), that the equilibrium point (x 0 ,y 0 ) moves along an arc of a parabola (see Figure 1), which contains the two equilibria (0,0) and (P,0) of the prey alone, i.e., in the absence of predator, (which therefore always exist independently of e).This parabola arc is traversed in the direction of decreasing x 0 for increasing e.According to Equation ( 6), the internal equilibrium point exists for f(x 0 ) > 0, which is equivalent to x 0 < P. The corresponding value of e is the viability threshold of predation.For increasing e, the equilibrium point moves along the parabola, rising towards the apex A (top of the parabola) and then falling towards the origin, which corresponds to e = +∞.The arc of the parabola between P and A is the normal mode (as e increases, the equilibrium population of prey decreases and that of predators increases).
For larger values of e, we are in the paradoxical region (between A and O; when the efficiency e increases and both equilibrium populations of prey and predators decrease).
The stability of the internal equilibrium point (x 0 ,y 0 ) depends on e and, therefore, on its position along the parabola.It is stable from P (i.e., for small e) up to a certain point, B, in the paradoxical region, as shown in Figure 1, where there is a PAH bifurcation (equaling Poincaré-Andronov-Hopf, often referred to as Hopf bifurcation, see for example [6,7]).Equilibrium still exists in the BO arc but it is unstable, with the attractor being a periodic cycle surrounding the equilibrium point.This cycle naturally depends on e; we have drawn it on Figure 1 for several values of e.The literature is very rich in these problems of equilibria cycles, see for example [8][9][10][11] and references cited therein.
It should be pointed out that this general pattern naturally depends on the parameter values.It is easy to see (see [1] Section 6.5) that the bifurcation of PAH is due to the capping of satiety, (parameter b).As b increases, the PAH bifurcation moves towards the origin, and for b = +∞ it coincides with the origin; in other words, for b = +∞ (which is equivalent to h(x) = e x), there is no PAH bifurcation and the equilibrium is always stable (it is the attractor).
It is useful to calculate, for each value of e beyond the PAH bifurcation, the average values of the densities of x(t) and y(t) along the corresponding cycle, which we will denote by  and  .They differ from (x 0 ,y 0 ) and describe (see Figure 1) a curve that is virtually rectilinear starting from the PAH bifurcation and moving for increasing e with increasing  and decreasing  .
The shape of the cycles is very interesting.For values of e slightly above the PAH bifurcation, it is a small, approximately elliptical cycle surrounding the equilibrium point.
But for larger e, the cycle gradually adopts a vaguely triangular shape, which becomes, for larger e, practically a curvilinear triangle with smooth (rounded) vertices, whose sides are the x and y axes and a curve arising from (P, 0) solution of the limit system for e = +∞ (which is equivalent to taking h(x) = b) and whose vertices are the two equilibria of the prey alone, O and P, and a point C.
This highly significant geometry is known as the HNR triangle (for Hubris-Nemesis-Resilience).In fact, as we go through the cycle as a function of t, the curvilinear side close to PC is an increase in the predator population at the expense of the prey population, due to the hubris of e efficiency.Arriving near C, there are very few prey and many predators, leading to a rapid decline in the predator population, which almost loses all food; this is the CO side (nemesis = punishment for hubris).At the approach of O, there are very few predators and prey can proliferate unhindered by predators, which is what makes the side close to OP (close to the evolution of prey alone, its prey resilience).This closes the cycle, which begins again periodically.4), except parameter e, which varies from e = 6 to e = 30: We have drawn the limit cycles exhibited by the system for e ranging from 6 for the smallest cycle to 30 for the largest (beyond e = 30, the more e increases, the closer the cycles are).The fixed points of the system for different values of e are represented by the arc of the parabola with vertex A (cross-shaped line).Finally, the almost rectilinear curve starting from B is the curve of means values ( ,  ).
But we can go further in describing the structure of HNR cycles for very large e: for large enough e, the cycle is very close to the curvilinear triangle.This is formed by the axes and the solution of the limit system for e = +∞ (which consists in replacing Tanh by 1) starting from the equilibrium of prey alone, (P,0).This solution, which is the hubris phase, intersects the y-axis with a non-zero x' velocity (because of Tanh = 1), so that there is a matching boundary layer for very large e to match the y-axis (the nemesis phase).This matching occurs fairly quickly and is very different from the vertices (0, 0) and (P, 0), which are points of equilibrium, so that near them the movement is very slow for large e, unlike near the top vertex, which is a matching equaling a sudden change of direction without slowing down for large e.
For very high values of e, by inspection of the solutions x(t) and y(t), we observe that the period (which must tend towards infinity because of the slowdowns mentioned) appears to be proportional to e.
In the behavior of x(t), we observe the presence of two levels close to 0 and P, corresponding to the slowdowns of the passages near the equilibrium points (0, 0) and (P, 0).See below (Figure 2) for e = 6 and e = 30.The behavior of y(t) is similar except that the two stages merge, as y(t) values are practically zero for the two passages near (0, 0) and (P, 0) and the arc of orbit that connects them.In fact, all that remains in each cycle is the rapid rise and fall of y(t) (hubris and nemesis).In other words, the periodic function y(t) becomes a kind of pulse (always the same regardless of e) increasingly spaced out in time as e→+∞.See below (Figure 3) for e = 6 and e = 30.Taken together, the four figures above perfectly explain the almost rectilinear shape of the curve of mean values ( ,  of the densities of x(t) and y(t) along the corresponding cycle, as a function of e (see Figure 1).Indeed, as e tends to infinity,  tends to zero and  tends to a positive constant.
Another important point for understanding the sequel is that the unknown x 0 is given by the second equation of system (1), which is that of predator functioning (x 0 is the prey population that keeps the predators in equilibrium); x 0 is therefore independent of a and P.This remarkable property is at the root of the understanding of certain aspects of the paradoxes considered in the sequel.
Later (Section 6), we'll introduce slightly more complex systems, known as predation-commensalism, to deal with specific questions.
It should also be noted that the mathematical properties encountered are to be taken as tendencies and not as exact results.This is the case for all properties established with the help of models; in particular here, when e→+∞, the model becomes inoperative because there are populations that are too small at certain instants, so that the very concept of population disappears.
We shall see that these elements are important for correctly interpreting the causality and implications of the choices made by evolution.However, let us remark that predatorprey interactions depend not only on biotic factors, but also on many other abiotic factors, such as wind speed, which can crucially alter the predation efficiency of the predator population; see [12] where the impact of wind speed and seasonality on different parameters has been studied for a model proposed in [13].

Evolution Chooses to Increase the Efficiency of the Predation Mechanism
In this section, we consider four numerical examples of competition between two predators of different efficiencies in different regions (normal, paradoxical, cyclic attractors).
We consider a population x of prey and two populations y1 and y2 of predators, modeled by the system where f x , h x and h x are given by expressions like ( 2) and ( 3) with indices 1 and 2. The values of the parameters (except e1 and e2) are equal for both predators: In the numerical computations in this section, and almost everywhere in this article, the starting point is taken close to the attractor of x and  alone, to highlight the displacement of the attractor in competition with  (if we start from another point, then the final behavior is obviously the same; however, this initial condition must be chosen sufficiently close to the attractor to avoid the transient part which would distort the comparison with the attractor of x and  alone).On the other hand, the attractor of x and  alone is easily found by making  0 0.

Numerical example -1-: in the normal region
We take  0.18, and 0.24 so that the two predators differ only in efficiency (that of  being greater than that of  ).It is easy to check that they are both in the normal region (the arc of the parabola between P and A: as e increases, the equilibrium population of prey decreases and that of predators increases).We then calculate the solution of Equation (7) starting from initial values such that x(0) and  are very close to the equilibrium position of the system without  ; indeed, we start from the initial condition (2.68,1.94,and 0.05).The numerical solution is shown in Figure 4a.
Unsurprisingly, the initial predator is replaced by the more efficient predator.This is accompanied by a decrease in the prey population and an increase in the predator population (the final population of  is greater than the initial population of  ).
Note that when starting from any other initial position the final result is the same (this is the attractor), but initially there is a transient regime, an oscillation, since we are not starting from an equilibrium position.

Numerical example -2-: in the paradoxical region
We take  0.3 and  0.5 (it is easy to check that they are both in the paradoxical region), starting from the initial condition (1.61, 2.13, and 0.05).The numerical solution is shown in Figure 4b.We can see that the initial predator  is always replaced by the more efficient  .This is accompanied by a reduction in the prey population and also a reduction in the predator population (the final  population is smaller than the initial  population).
Competitive evolution in the paradoxical region is therefore moving towards an increase in the efficiency of the predation mechanism, to the detriment of both predator and prey populations.

Numerical example -3-: 𝒚 𝟏 in the paradoxical region, 𝒚 𝟐 slightly beyond the bifurcation of PAH
We take  0.3,  0.96 (the predator  is the same as in the previous example, but  is more efficient, having passed the bifurcation, which takes place around  = 0.93).Starting from the initial condition (1.61, 2.13, 0.05), the numerical solution is shown in Figure 4c.As in the previous example, competitive evolution leads to an increase in the efficiency of the predation mechanism, to the detriment of both predator and prey populations.But the final state is oscillating, this is the cyclic attractor of the system in (x, ).Numerical example -4-:   and   far beyond the bifurcation of PAH We take  1,  1.4 (both predators have much greater efficiencies than the bifurcation, which is at e = 0.93 approximately).Starting from the initial condition (1.6, 2.1, 0.05) the numerical solution is shown in Figure 4d.This result lies in the region of large cycles, which are only known by numerical computation, making it all the more difficult to interpret the results, or at least to draw general statements from them.The predator is always replaced by a more efficient one.The prey population is visibly increasing; the predator population also seems to be increasing in the higher values, but the peaks are more widely spaced, so that the average values over the period are smaller, in line with Figure 1.It thus appears (on the basis of these numerical computations) that, all other parameters being the same, evolution always chooses to increase the efficiency e, which is not necessarily accompanied by a demographic boom in predators (their population decreases in the paradoxical region).This result is highly significant, as it seems to indicate that, as opportunities arise, the natural tendency of evolution is to increase the efficiency of the predation mechanism.If these opportunities arise repeatedly, then the process eventually leads to HNR cycles, and even eventually to the disappearance of this mechanism, as the model fails for very small populations.

Theory of Competitive Evolution in the Predatory System
System ( 7) is in the framework of systems of the form, The equilibrium points inside the positive octant (i.e., with the three strictly positive coordinates) are the solutions of the system (10) below, which contains two equations with the single unknown x: It follows that, generically, there is no equilibrium point inside the octant and the point attractors are on the boundary.Of course, there may be non-point attractors, and we will see examples of this (numerical example -5-).For the solutions of ( 9), let us calculate the expression: As the populations are positive, this shows that, at each point, the vector field passes through the  / constant plane towards the  or  axis, depending on whether  is larger or smaller than  at that point.
In the special case of system ( 7) with ( 2) and (3), we have so that competition between two predator species only concerns the significant parameters c, e, b, and τ.This expression enjoys the remarkable property of increasing with respect to each of the parameters e, b, and τ and decreasing with respect to c, independently of the point x, (this is obvious for c, e, and τ; for b too, with a simple calculation, or by thinking of its role as a satiety ceiling).This result is summarized in the proposition below.
Proposition 1.Consider system (7) in the case where the difference between the two attractors concerns only one of the parameters c, e, b, and τ.Then, the attractor is necessarily in one of the planes with coordinates (x, or (x, , precisely the one with the largest (resp.smallest) value of the parameter concerned, if this is one of b, e, or τ (resp.c).
It is worth noting that this proposition, which may seem obvious from a Darwinian perspective, is rigorously true, concerning attractors whether or not they are of a punctual nature (or cyclic, in particular), and that it is nevertheless perfectly compatible with the paradox pointed out in Section 2, which concerns the parameter e e (the attractor is therefore in the (x, plane) and where the equilibrium populations are smaller for (x, than for (x, . On the other hand, if we modify two (or more) of the significant parameters, then the value of   is no longer necessarily the same everywhere; it may depend on x, and we can only derive local results from (12).So, for example, let (x 0 ,y1 0 ) be a stable equilibrium point of the system for x and  alone.Then, the point (0,x 0 ,y1 0 ) is automatically an equilibrium point of the three-dimensional system; its local stability is given by that of the linearized system (Jacobian matrix at the point (x 0 ,y1 0 ,0)).It has two negative real eigenvalues (corresponding to the restriction to  0, so this point is three-dimensionally locally stable or not, depending on whether the eigenvalue corresponding to the transverse eigenvector is negative or positive (it is necessarily real)).We can see from the third Equation ( 9) that transverse stability may or may not hold, depending on whether   is negative or positive.Now, since (x 0 ,y1 0 ) is an equilibrium point of the system in the (x,y1) plane, we have   0 (see second Equation ( 10)).So, we have stability if, or instability if the sign in Equation ( 13) is >.This result is summarized in the proposition below.
Proposition 2. Let (x 0 ,y1 0 ) be a stable equilibrium point of the system for (x,y1) alone.Then the point (x,y1,0) is a three-dimensional attractor if ( 13) is satisfied, and it is transversely unstable if (13) with sign > is satisfied.A similar result holds for a stable equilibrium in the (x,y2)-plane.
We have an analogous result in the case of a cyclic attractor  of the system for (x,y1) alone, replacing q by its integral along .Let (x 0 (t),y1 0 (t)) be the parametric representation of  as a function of time.From the first two Equations (10) we have, (which is the analogue of   above).The variational Equation of (10) with respect to  in Equation ( 9) gives: so that stability amounts to This result is summarized in the proposition below.
Proposition 3. Let (x 0 (t),y1 0 (t)) be the parametric representation a cyclic attractor of the system for (x,y1) alone.Then it is also a three-dimensional attractor if ( 16) is satisfied, and it is transversely unstable if we have (16) with the sign > instead of <.A similar result holds for a stable equilibrium in the (x,y2) plane.
Naturally, Propositions 2 and 3 hold in the framework of the local theory of differential equations (whereas system (9) is perfectly general).Specific examples can be found by numerical computations.

General Comment on The Choice of Evolution as a Functional Vocation
Finally, it is easy to get a synthetic idea of the (often paradoxical) behavior of competition between two predators by taking the following two elements into account: -i-The dynamics of the system with regard to competition between  and  is given by Equation (11) (and the analogue by exchanging indices 1 and 2), which we reproduce here in a slightly different form: which holds true for all positive x, , and  .
-ii-Of the significant parameters of this predation model containing e, b, τ, c, a, and P, only e, b, τ, and c, are involved in (17).The other two, a and P, are not involved, as they relate to prey behavior.The expression q is increasing with respect to e, τ, and b (for the latter this results from its meaning as a predation ceiling, or from a little calculation) and decreasing with respect to c.
Let us consider the case where only one of the parameters e, τ, b and c is different for the two attractors.Expression (17) has constant sign everywhere, so the attractor (whatever its nature, point-like or not) is necessarily on one of the coordinate planes (x,y1) or (x,y2), the one corresponding to the larger e, b, τ or the smaller c.This is elementary Darwinism, as long as we interpret "advantage" and "disadvantage" in the sense of larger or smaller q and not the predator population or anything else large or small.The same applies when there are several parameters with different values, provided that the differences in the values of these parameters are all in the same direction of increasing or decreasing q (examples: increasing τ and b or decreasing e and increasing c).
On the other hand, if several of the parameters e, τ, b, and c are different for the two predators, and the differences in the values of these parameters produce variations of various signs on q, then the sign of the right hand term in ( 17) may be different for different values of x, making it possible for a cyclic attractor to exist inside the positive octant (orbit returning to the starting point, see for instance the numerical example -5-hereafter).The attractor can then lie on one of the coordinate planes or inside the positive octant.Note, in any case, that there can be no point attractor outside the coordinate planes, as we pointed out in connection with Equation (10).
It is useful to return to the extremely simple case where only one of the significant parameters e, τ, b and c is varied.We can say that the predation mechanism has a tendency (or vocation) to increase e, τ, and b and to decrease c.This statement should be understood in the sense that, if evolution allows a variation in just one of these parameters, then the competition will result in a choice of the highest (resp.lowest) value of e, τ, and b (resp.c).This corresponds to the meanings of the terms tendency and vocation (and even to the term wish) in common parlance: a volitional agent can have one or several vocations, which means that, given the opportunity, he will choose to follow each of these vocations if the others are not affected, whereas a proposition implementing several vocations of opposite tendencies can have diverse outcomes (think of a person with artistic and welfare vocations, who has to choose between a "standard" job and a more artistic but less well-paid one).
This approach to interpreting results can be taken much further.If the choice of evolution is proposed repeatedly, always involving variations in just one of the parameters e, τ, b, or c, then it will systematically result in an increase in e, τ, or b or in a decrease in c (proposals of opposite sign being rejected).This may seem natural and trivial, but we must realize that these trends or vocations do not always go in the direction of the predator's demographic boom, as the paradoxical examples show.Now, if we think of the biological interpretation of the predation system, then we can imagine that, in certain cases, τ, b and c depend on the physiological properties of the predator, whereas e depends on the ethological properties of the predator in its relationship with the prey and the environment.It is therefore natural to consider e as much more subject to the vagaries of evolution than τ, b, or c.This makes the problem we are dealing with particularly interesting when only e is subject to variation.And it is precisely the variation of e with the other fixed parameters that reveals the paradoxical region and bifurcation of PAH developed in Section 1.
We saw the special role of the parameter e and the interest in increasing it over time.It is easy to make a non-autonomous example, with slowly increasing e efficiency.As the variation is slow, the solution is practically quasi-static, so that at each instant we practically have the attractor for the corresponding value of e.For this purpose, we always use system ( 1)-( 2 This can be compared with Figure 1 in Section 1, identifying an initially non-paradoxical behavior, which becomes paradoxical after passing through the apex A (vertex of the parabola).Then, the PAH bifurcation occurs and the dilating cycles adopt the HNR gait.
We can also make a parametric representation in space, Figure 5c, which can thus be seen as a symbolic vision (or artist's view) of the natural evolution of a predation mechanism whose efficiency is subject to evolutionary choice, which accepts (or refuses) the positive (or negative) variations that chance submits to it.
Instead of the straight-line e(t) = 0.145 + εt used above, we can advantageously use another e(t) function of the proposed evolution, such as the one proposed in the other allegorical figure below (which is a regular curve with relative maxima and minima; Figure 6).The time will appear to pass very slowly, so that we can consider the state of the system as a sequence of asymptotic states on the attractor (in the particular case where the attractors are points, we say it is a quasi-static evolution).By taking a larger b, the onset of cycles is delayed.

Numerical example -5a-:
We take advantage of the well-known property that increasing b delays the bifurcation of PAH to modify the numerical example -2-with a much higher b.We then set the e2 value even higher, and the paradox becomes even more obvious.We start from the equilibrium of x with the y1 and the choice of evolution is based on the equilibrium of x with y2, which involves demographic losses of prey from 1.51 to 0.45 and of predators from 2.08 to 0.9.This situation is close to the "killing of the goose that lays the golden egg".Specifically, the parameters are given by Equation ( 18) below, See Figure 7a.2)-( 3), where the parameters are given by Equation ( 18); (b).Temporal series x(t), y1(t), and y2(t) for system (1)-( 2)-( 3), where the parameters are given by Equation ( 18), except for  0.8; (c,d).Solutions for system (1)-( 2)-( 3)( 1)-( 2)-( 3), where the parameters are given by Equation ( 18), except for  0.55.(c) Temporal series (t,x), (t,y1), (t,y2).

Numerical example -5b-:
Figure 7b below is close to the previous one, but with a much smaller b2 (b2 = 0.8), keeping b1 = 3.This is a first example of evolution by changing e (to a large value) and b (to a smaller value).Both e and b tend to increase, and e prevails.The result is a replacement as in the previous example but the final result is oscillating, which could be described as a kind of "intermittent paradox".Specifically, the parameters are given by Equation ( 18), but with  0.8; see Figure 5b.

Numerical example -5c-:
In the previous example, it is sufficient to reduce a little further the value of  (from 0.8 to 0.55) to destroy the replacement of  by  .Instead, the two predators cohabit (in addition to the prey) in an internal stable cycle.This is an example of the persistence of two predators with a prey, see [5,[14][15][16] or Chapter 10 of [1] for more on this interesting subject.Naturally, this type of solution cannot exist in the cycle-free region (i.e., with e small enough to be below the PAH bifurcation); see Figures 7c,d.

Predation-Commensalism System and the Role of Neutral Parameters
It will prove useful to present here a slightly more complex system, which was introduced in [1] (Chapter 9) to model the global economy.The precise interpretation, in terms of this system, of the well-known property of TRPF (Tendency of Rate of Profit to Fall), which concerns the parameter c, was the starting point of the present work.
In this context, one can refer to [17], which describes a bioeconomic predator-prey model derived from the model proposed in [18], which incorporates the square root functional response and nonlinear prey harvesting.Due to the last condition, the model demonstrates intricate dynamic behaviors in the predator-prey plane, mainly via the occurrence of Hopf bifurcation utilizing economic profit as a bifurcation parameter.
In this system, known as predation-commensalism, predation activity exerts a beneficial function on the prey substrate, increasing the capacity of the environment.In concrete terms, the predation-commensalism system is still (1)-( 2)-( 3), but the capacity of the medium P is no longer a constant, but rather a function of x: where the new (positive) constants  and  are the capacity of the substrate in the absence of predation activity and the coefficient of influence of this activity on capacity. is completely independent of the predator, whereas  depends on it in its relationship with the prey, so that in the case of two competing predators, the system naturally follows (7) with This system is always of the form (9), so the properties of Section 3 remain valid.It is remarkable that the parameters  and  , which, as we have just indicated, are clearly associated with one or other of the predators, are nevertheless perfectly neutral in the evolutionary competition.From a mathematical point of view, this follows from Equation (11), which does not contain these parameters.From a modeling point of view, this property is perfectly natural, since the influence (coming from either of the two predators) is put into the common pot (capacity of the environment in terms of prey quantity) and therefore constitutes no advantage or disadvantage for one of the predators over the other.In particular, if the two predators differ only in the values of , then the system has the first integral y2/y1 = constant as it follows from (11) with q1 = q2.
On the other hand, commensal terms have a huge influence on the equilibrium population of predators.We also observe that the value of the x 0 of the equilibrium is independent of these new terms; it therefore appears that the terms of influence of predation on the substrate capacity are reflected on y 0 and not on x 0 .This (still paradoxical) property is a consequence of the property reported at the end of Section 1 that x 0 is obtained directly from the equation in y', which is not influenced by .
It is therefore clear that the parameters e, τ, b, and c retain the properties reported in Section 4 with regard to their character as advantages or disadvantages in the competition, while  is neutral.A consequence of this is that parameter , while being elusive in the competition when it is the only one to undergo a variation, can be dragged along by even an infinitesimal variation of e, c, b, or τ.It is therefore highly unstable.

Numerical example -6a-:
We start with an example in the non-paradoxical region, with values of all parameters equal for both predators, and we add commensalism terms  0.25 and  1.25.Figure 8a shows three orbits in space starting from points close to the equilibrium of the prey alone  , 0, 0 .Starting with  0 0 (resp  0 0), we naturally have an orbit for the system of ,  (resp ,  ) alone; the orbits are naturally in the coordinate planes.
This clearly shows the demographic advantage of  compared with  .But starting with non-null  0 and  0 , the orbit is in the plane  and in no way chooses  at the expense of  .See Figure 8a, in which the parameter values are given by Equation ( 21

Numerical example -6b-:
We repeat exactly the previous calculations (example -6a-), with parameter values as given in Equation ( 21) but taking  0.44 (instead of 0.45).This minimal difference implies that evolution chooses  , and  disappears (Figure 8b).It can also be seen that the three-dimensional orbit clearly has an initial phase with very little difference from that of the previous case (fast dynamics), followed by a second phase where it joins the attractor on the plane (slow dynamics).

Numerical example -6c-:
This is now an analogous case with a loss of  instead of a gain.More precisely, we choose a system in x,y2 with λ2 = 0, and e2 = 0.101, sufficiently small for predation not to be viable, so that the attractor is simply the equilibrium of prey alone  , 0, 0 .The system in ,  has an even smaller efficiency of  0.1, but, thanks to the commensalism term (with  2.2), there is a stable equilibrium  ,  , 0 with  0. Competition between the two populations naturally chooses that with greater efficiency, and so it disappears.Figure 8c shows this process starting from a point very close to the system attractor in ,  : the two predators disappear.The parameter values are given by Equation ( 21) but we take e2 = 0.102 and λ1 = 2.2, λ2 = 0.  2)-( 3), with a capacity of the medium P as a function of x given by Equation ( 19).In the case of two competing predators, the system is Equation ( 7) with Equation (20).The parameter values are given by Equation ( 21

Conclusions and Discussion
Within the general framework of mathematical ecosystems ecology, the aim of this work is to show how, in a number of situations that are perfectly plausible from an ecological point of view, the choice of evolution is not based on the demographic advantages of the populations in question, but on the functional refinement of the mechanisms (particularly predation) that link them, leading to the appearance of new structures, particularly of a pulsating nature.These are functional emergence phenomena, sometimes to the detriment of the quantitative properties of populations.
An important point is the consideration of the significant parameters of the functional mechanism under consideration, whose variations are proposed at the choice of the evolution, unlike the mathematical theory of structural stability, which involves general disturbances.
The starting point is an analogy with the general law of Tendency of Rate of Profit to Fall (TRPF) in the global economy.This law, initially derived from heuristic considerations, can be obtained mathematically from appropriate models of the economy (see [1] Chapter 9 Section 9.3).The epistemological status of the properties we are establishing is analogous to that of the TRPF.It is a kind of general tendency, a vocation that occurs under specific conditions (and often in the absence of others), but not automatically.
We have focused on the role of parameter e (the efficiency of the predation mechanism), whose role it is useful to understand.Under conditions of scarcity of prey population x, it is the ratio of the number of prey consumed per unit of predator and time, and the prey population.This is the pure predation mechanism, with no relation to the predator's needs and abilities (which relate to other parameters).Clearly, e depends not only on the physiological properties of prey and predator, but also on environment and habitat, and it can be easily influenced by environmental agents.This is why it seems obvious to us that the parameter e is subject to disturbance under conditions where all the others (particularly physiological ones) remain constant.
However, the parameter e of predation efficiency enjoys paradoxical properties with regard to equilibrium populations and their stability.For moderate values of e, an increase in its value leads to an increase in the equilibrium population of predators and a decrease in that of prey (which seems normal), but beyond a certain value of e, its increase leads to a decrease in both equilibrium populations (of predators and prey alike).This is the paradoxical mode, because the scarcity of prey leads to a lack of resources for predators, whose equilibrium population in turn declines.But, in addition, beyond a certain value of e, the equilibrium loses its stability (Poincaré-Andronov-Hopf bifurcation) and the attractor becomes a periodic cycle around the equilibrium.To determine this value of e, it suffices to note that bifurcation occurs when the divergence of the vector field at the equilibrium point changes from negative (contraction = stability) to positive (expansion = instability, sending to a cycle).A detailed study of this problem is given in [19].By further increasing the value of e, the cycle takes on the characteristic shape of a curvilinear triangle; in each period we can distinguish three phases: Hubris (excessive predation, or rapid increase in the predator population to the detriment of the prey population), Nemesis (punishment for hubris, or subsequent reduction in the predator population) and Resilience (or slow recovery of the prey, back to the initial state).Finally, for very high values of e, the periodic phenomenon takes on a form akin to that of seasonal epidemics (the predator population is always practically zero, except at certain moments periodically spaced in time).
Our work shows why and how, when two predators compete for prey, evolution always chooses the most efficient predator (and not any kind of population optimization).
As a result, the efficiency e of the predation mechanism appears to be an emergent property, which can only increase through evolutionary choice when the suitable conditions are present, and which is therefore irreversible in this context.This implies a natural tendency for periodic cycles to emerge, and even for them to be structured into sequences of hubris, nemesis and resilience.This tendency often runs counter to demographic optimization.
Furthermore, since the choice of evolution is also based on increasing e, as the period increases according to evolution, it sweeps through a wide range of frequencies, making possible synchronization with other phenomena through small interactions.This furnishes a plausible explanation for the fact that seasonal phenomena are very widespread.Indeed, synchronization is a ubiquitous phenomenon characteristic of many processes in natural systems and (nonlinear) science.It is today considered one of the basic nonlinear phenomena studied in mathematics, physics, engineering or life science; see for example Chapter 12 of [1], or [20] and references cited therein.
Obviously, the phenomena become much more complex if variations in several significant parameters are simultaneously involved.Several examples have been given (Sections 5 and 6), in particular, one that leads to a non-choice (example of the coexistence of two predators, which is often controversial).

Figure 1 .
Figure 1.This figure shows several results, with all parameters fixed as given by Equation (4), except parameter e, which varies from e = 6 to e = 30: We have drawn the limit cycles exhibited by the system for e ranging from 6 for the smallest cycle to 30 for the largest (beyond e = 30, the more e increases, the closer the cycles are).The fixed points of the system for different values of e are represented by the arc of the parabola with vertex A (cross-shaped line).Finally, the almost rectilinear curve starting from B is the curve of means values ( ,  ).

Figure 2 .
Figure 2. The temporal series x(t) for e = 6 and e = 30, the other parameters are given by Equation (4).

Figure 3 .
Figure 3.The temporal series y(t) for e = 6 and e = 30, the other parameters are given by Equation (4).

Figure 5 .
Figure 5. Solutions of system (1)-(2)-(3), for the parameters given by a = 1.0,P = 4.0, b = 0.53, c = 0.43, and τ1 = 1.0 and efficiency e1 = 0.145 + εt, where ε = 0.0008, with time t varying from 0 to 1000.We start from an initial condition close to the equilibrium of the prey alone: (a) temporal series x(t) and y(t), (b) the orbit in (x,y) plane, (c) parametric representation of the solution in space (t,x,y) of the natural evolution of a predation mechanism whose efficiency is subject to evolutionary choice.

Figure 6 .
Figure 6.Another function e(t) of the proposed evolution, a regular curve with relative maxima and minima.