Entropy and equilibria in competitive systems

This paper investigates applicability of thermodynamic concepts and principles to competitive systems. We show that Tsallis entropies are suitable for characterisation of systems with transitive competition when mutations deviate from Gibbs mutations. Different types of equilibrium in competitive systems are considered and analysed. As competition rules become more and more intransitive, thermodynamic analogies are eroded and the behaviour of the system can become complex. This work analyses the phenomenon of punctuated evolution in the context of the competitive risk/benefit dilemma.


Introduction
The question of whether systems involving competition can be characterised by quantities resembling conventional thermodynamic parameters does not have a simple unambiguous answer. This problem was investigated in ref. [1] and it was found that such characterisation is possible under conditions of transitive competition but, as the system becomes more and more intransitive, the thermodynamic analogy weakens. The similarity with conventional thermodynamic principles is strongest when mutations present in the system belong to the class of Gibbs mutations. While deploying the conventional logarithmic definition of entropy, the analysis of ref. [1] misses an important point: when mutations deviate from the Gibbs mutations, the family of Tsallis entropies [2] represents a very convenient choice of entropy to treat these cases. This omission is rectified in the present work. We also note that Tsallis entropy has been recently used in modelling of biological replications [3].
Thermodynamics is strongly linked to the concept of equilibrium. Competitive systems allow the introduction of different types of equilibrium, possessing different degrees of similarity with the concept of equilibrium in conventional thermodynamics. The current work discusses possible cases of competitive equilibria and performs a detailed analysis based on Tsallis entropy of the equilibrium through a point of contact, which is more similar to conventional thermodynamics than the other cases.
From the thermodynamic perspective, the present work is only an example of using Tsallis entropy. We do not attempt to draw any general thermodynamic conclusions and the use of non-extensive entropy in other applications may well be different from our treatment of equilibria in completive systems. The problem of general consistency between physical equilibrium conditions and definitions of non-extensive entropies has been analysed by Abe [4,5]. Non-extensive statistical mechanics has been reviewed by Tsallis [2], while non-extensive entropies associated with this mechanics are discussed in refs. [2][3][4][5][6][7] and many other publications.
The last section deals with intransitive cases when the thermodynamic analogy weakens and the possibility of using entropy as a quantity that always tends to increase in time or remain constant is not assured. This section analyses the risk/benefit dilemma represented by a competitive system, whose evolution can be transitive or intransitive depending on the choice of the system parameters. In the intransitive case, the evolution of the system appears to be punctuated by sudden collapses and becomes cyclic. This punctuated evolution is similar to the concept of punctuated equilibrium in evolutionary biology [8], although in the context of thermodynamics the latter term might be misleading as the system is not in equilibrium and keeps evolving between the punctuations.

Competitive systems
Competitive systems involve the process of competition in its most generic form. The elements of competitive systems compete with each other according to preset rules. The rules define the winners and losers for each competition round based on properties of the elements denoted here by y. The properties of the losers are lost while the winners utilise the resource vacated by the losers and duplicate their properties. The process of duplication is not perfect and involves random mutations, which are mostly detrimental for competitiveness of the elements. The expression A ≺ B (or equivalently y A ≺ y B ) indicates that element B with properties y B is the winner in competition with element A with properties y A . If two elements have equivalent strength y A y B , the winner is to be determined at random. In computer simulations, the elements are also called Pope particles and exchange of properties is called mixing by analogy with the conventions adopted in particle simulations of reacting flows. Two forms of mixing -conservative and competitive -can be distinguished. The former is predominantly used in the flow simulations while the latter is associated with competitive systems. The rest of this section introduces basic terms used in characterisation of competitive systems; further details can be found in refs. [1,9].
The competition rules are divided into two major categories: transitive and intransitive. In transitive competitions superiority of B over A and C over B inevitably demands superiority of C over A, that is As illustrated in figure 1a, transitive competitions enable introduction of an absolute ranking r(y), which is a numerical function that determines superior (stronger) and inferior (weaker) elements: The competitive transformations can be interpreted as reactions between the particles where B is different from B due to mutations. B is the winner in competition with A and, thus, B is entitled to occupy the resource (i.e. particle) previously occupied by A. Properties of A are lost and the properties of B are copied across into A. (Conservative properties, which are not considered here, would be transferred in the opposite direction from the loser to the winner.) The copying process is not perfect due to mutations, which are random alterations of properties of B during copying. If mutations are not present then B =B. Unlike random walks, mutations have a strong preference for the negative directions: it is likely (in case of non-positive mutations, it is certain) that B B. If rare positive mutations B B are present, the distribution of particles may escalate towards higher ranks when the leading particle (i.e the particle with the maximal absolute ranking in the group) is occasionally overtaken by a new leader. One of the main results of ref. [1] is linking absolute ranking to the entropy potential s y = s y (r) and, under some restrictions (e.g. Gibbs mutations), proving the associated competitive H-theorem.
The competition rules, however, are not necessarily transitive and the competition is deemed intransitive when at least one intransitive triplet exists (see figure 1b). Although intransitive competitions do not generally permit absolute ranking of elements, they can be characterised by a coranking function which, by definition, should be antisymmetric ρ(y A , y B ) = −ρ(y B , y A ). In the case of transitive competitions, the coranking function can be expressed in terms of absolute ranking by In addition to ρ(y A , y B ), it is useful to define sharp coranking In case of two-particle mixing, evolution of the system is determined by the sharp coranking R. The graded coranking can be useful in establishing relative ranks within each mixing group, when mixing of multiple particles is considered. The distributions of elements in the property space is characterised by the particle distribution function ϕ(y) = nf (y) where n is the total number of particles and f (y), which can be interpreted as the probability distribution function (pdf), satisfies the normalisation condition ∞ f (y)dy = 1 When a competitive system is divided into K subsystems I = 1, 2, ..., K and each subsystem I has the b I -th fraction of the particles, we may characterise each of these subsystems by its own normalised distribution φ I (y); that is ϕ I (y) = n I φ I (y), b I = n I n , The subsystems can be distinguished by having different domains D I or by other means. When subsystems are distinguished, it is useful to define the coranking of the distributions which indicate relative strength of subsystem distributions with respect to each other. We can say "the subsystem I is stronger than the subsystem J" and write Examples of systems using competitive mixing can be found in refs. [1,[9][10][11] 3. Competition and q-exponential distributions We consider transitive competition with elements possessing a scalar property y, which is selected to be aligned with ranking (that is r(y) is a monotonically increasing function and the absolute ranking is effectively specified by y). Hence, for any two elements A and B The problem is deemed to be uniform with respect to shifts along y. Assuming that mutations, which are originated at point y and distributed with the probability density function f m (y, y ), are uniform f m (y, y ) = f m (y − y ), the general competitive evolution equation [1] takes a more simple form given by The competitive evolution equation specifies a balance between mutations, given by the first term on the right hand side of this equation, and the losses due to competition, given by the second term. The function F is the cdf (cumulative distribution function) of the pdf f If mutations are non-positive then Equation (13) can be integrated to yield where F m is the cdf that corresponds to pdf f m . The Gibbs mutations [1] correspond to q = 1 implying that the distribution f m (y) is based in this case on the conventional exponent f m (y) = exp(y)H (−y) where is the Heaviside function. Note that there is no loss of generality in setting α = 1 in f m ∼ exp(αy) since the variable y can always be rescaled to eliminate α. In case of Gibbs mutations, the pdf f is given by where y * is the position of leading particle.
In this work we are interested in the case when the pdf f can be approximated by the q-exponential distribution where e y q = exp q (y) = (1 + (1 − q) y) 1 1−q is the so called q-exponent and ln q (y) = y 1−q − 1 1 − q is the corresponding q-logarithm. If q → 1 then the q-functions approach the conventional exp(y) and ln(y). The cdf, corresponding to pdf (20) is given by where Q = 1/(2 − q). The distribution (20) solves the governing equation (13) with an asymptotic precision of O((1 − q) 2 ) provided the mutations are distributed according to Figure 2 illustrates that, as expected, the cdf of simulated distributions are very close the corresponding q-exponents when q is close to unity. The q-exponential functions can also serve as very good approximations for distribution in competitive systems for a wide range of q. Consider the q-exponential distribution of mutations with the cdf F m (y), which is given by q-exponential functions similar to (21) with Q = 1/(2 − q ) The approximate solutions shown in figure 2 correspond to q-exponential (20) with The cdf shapes presented in Figure 2 indicate that, although q-exponential distributions are not necessarily exact for competitive systems, they are reasonably accurate and correspond very well to the physical nature of the problem when mutations deviate from Gibbs mutations. In the competitive system illustrated in figure 3, every location is taxed due to competition with superior elements and at the same time is supplied by mutations originated at superior elements. For Gibbs mutations, the competitive system schematically depicted in figure 3 is in the state of detailed balance: every location is taxed and supplied at the same rate by any given superior. In simple systems with constant a priori phase space A(y), the Gibbs mutations are distributed exponentially (q = 1). When mutations deviate from Gibbs mutations, the overall rates of taxing and supplying must negate each other under steady conditions but there is no detailed balance in relations with different groups of superiors. For long-tailed (superexponential) distributions with q > 1, weak particles are supplied more by the leaders and are taxed more by immediate superiors. For short-tailed (superexponential) distributions with q < 1, weak particles are supplied more by the immediate superiors and are taxed more by the leaders. Competitive systems are aimed at studying generic properties of systems with competition and mutations. Although we do not specifically intend to model distributions of biological mutations, these distributions are still of some interest here as real-world examples of complex competitive systems. Ohta [12] considered near-neutral genetic mutations and suggested that these mutations have exponential distributions. Modern works tend to use the Kimura distribution [13], which has a complicated mathematical form, deviates from pure exponents and, theoretically, corresponds to a genetic drift of neutral mutations. It seems that the reported distributions of genetic mutations tend to be slightly subexponential. Figure 4 illustrates that the experimental distribution of mutation A3243G of mitochondrial DNA in humans [13] is well approximated by q-exponential cdf with Q = 0.8. Since these mutations are known to be deleterious [14], they are shown as negative in the figure (in agreement with the notations adopted in the rest of the present work).

Tsallis entropy in competitive systems
Free Tsallis entropy in competitive systems is defined by with two likely choices of the exponent γ given by γ = 1 and γ = q. Here we denotef (y) = f (y)/A(y). The first term in the integral is the configurational entropy, which represents the randomising influence of mutations, while the second term involves the entropy potential s y (y), which reflects the influence of competition (i.e. s y increases with r reflecting higher likelihood of survival of more competitive elements). The termf (y) γ s y (y) in (26), which can be called the escort term, reflects nature's preference for elements with higher rankings (for example, in biological systems, ranking reflects fitness). Due to presence of the potential s y (y), the competitive entropy defined by (26) is analogous to free entropy of conventional thermodynamics, which is proportional to free energy (Gibbs or Helmholtz) taken with the negative sign. The distribution and entropy that correspond to s y = 0 (i.e. not affected by competition) can be termed "a priori" bearing some resemblance to prior probabilities in Bayesian inference. The definition of entropy in (26) is of Boltzmann type, i.e. implying validity of the Stosszahlansatz (stochastic independence of particles from each other). Variation of the distribution function results in where and χ y can be interpreted as the competitive potential of state y. Maximisation of S constrained by the normalisation and by the location of the leading element, that is f (y) = 0 for any y y * results in the following condition δS δf (y) where λ is the Lagrange multiplier, implying that the local competitive potential is the same everywhere in equilibrium. Consider a simple case of scalar y and constant a priori capacity A = const .
that with s y (y) = ky, A = c q q and λ = −ky * results in the pdf and cdf given by where Q = 1/q and q 2 = 2 − q. Since y * is arbitrary in this case, the distribution can be freely shifted along y. The location of y * is determined by (30).
2. The multiplicative case of γ = q. Equations (27) and (31) take the form with arbitrary value of Z q depending on λ. This value can be determined from the normalisation requiring that Z q = k −1 /(2 − q). The corresponding cdf F Q (y, y * ) is the same as (33) with the q-parameter given by Q = Z q k = 1/(2 − q).
In case of physical thermodynamics, Tsallis et. al. [6] recommend using γ = q in conjunction with the escort distribution for the energy constraints as the best option. Competitive thermodynamics , as considered here , does not have any energy constraints (assuming that the conservative properties are limited to the number of particles, we do not have any energy defined for the system) and selection of γ needs to be considered again. The choice of γ for competitive systems is determined by the physics of the problem and can be different for different processes. If infrequently positive mutations are present and the distribution with fixed number of particles escalates by gradually increasing y * in time, then γ = 1 is preferable. Indeed, while y * increases, the definition of entropy remains exactly the same and the escalation is seen as a natural process of increasing entropy in the system. If γ = q, the definition of entropy is dependent on the position of the leading particle. The choice of γ = q is more suitable for competitions between subsystems placed at fixed locations but with the numbers of particles that can be altered due to exchanges. Gibbs mutations correspond to q = 1 and the choices γ = 1 and γ = q coincide in this case. In the previous work [1], the Boltzmann-Gibbs entropy was used for non-Gibbs mutations by artificially making the phase volume dependent on the leading particle position A = A(y, y * ). Unlike the Tsallis entropy considered in the present work, the old treatment of the problem [1] did not allow for a unified definition of entropy valid for different y * (i.e. the Boltzmann-Gibbs entropy provides a unified, y * -independent definition of competitive entropy only for the Gibbs mutations).

Equilibria in competing systems.
A competitive system can be divided into subsystems and the question of equilibrium conditions between these subsystems appears. If the system is subdivided into K subsystems I = 1, 2, ..., K and subsystem I has the b I -th fraction of the particles, we may characterise each of these subsystems by its own normalised distribution φ I (y) as specified by equations (9) and (10). Assuming that equilibrium or steady-state conditions are achieved within each subsystem, the major equilibrium cases include: 1. Equilibria in isolated subsystems. Isolated subsystems do not exchange mutations and do not compete against each other. Equilibria are established in isolated subsystems independently of the other subsystems (see figure 5a) 2. Competing equilibria. Particles in these subsystems compete against each other but mutations do not cross the subsystem boundaries as illustrated in figure 5b. Competing equilibria tend to be less stable than connected equilibria considered below and, generally, are impossible in transitive competitions. Indeed, if y * I > y * J (i.e. the leading element of subsystem I is more competitive than the leading element of subsystem J), y * I cannot lose to any element of subsystem J while y * J will eventually lose to leading elements of I. There is no equilibrium in transitive competition depicted in figure 5d since the subsystem I = 1 is going to win all particle resources for itself. If y * I = y * J then the two leaders will eventually meet in competition and due to their equivalent strength, the winner of this round (which ultimately belongs to the winning subsystem) is to be selected at random. Competing equilibria are nevertheless possible in intransitive competitions. This, obviously, requires that for all I = 1, ..., K or, otherwise, n I would grow for R I > 0 and decrease for R I < 0. As discussed in Appendix of ref. [1], oscillations are to appear in competing equilibria between subsystems 1, ..., K, unless 3. Connected equilibria. In this case the subsystems I = 1, ..., K are connected by both competition and mutations. For Gibbs mutations (q = 1), the competitive H-theorem applies ensuring the detailed balance of the equilibrium state [1]. This implies that, in equilibrium, the connection between any two elements or groups of elements can be severed without any effect on the state of the system. Severing connection between two elements terminates both competition and mutations between these elements. This is illustrated in figure 5c, where the direct connection between points A and B is severed, although A and B remain connected through other elements as shown by the dashed line. The equilibrium conditions are given by the equivalence of all competitive potentials χ I = χ J for every I and J, where the formula for competitive potential [1,9] is obtained by differentiating entropy with respect to n I . The partition functions Z I are evaluated for each subsystem I as integral over the subsystem domain D I Equilibrium in competitive systems with Gibbs mutations resembles most the equilibria of conventional thermodynamics. In the case of general non-positive mutations (i.e. non-Gibbs mutations), the state of the system depends on the type of contact. Here, we distinguish two cases of interest: (a) Point of contact. Two subsystems I and J have a point of contact at y = y • when the elements from the vicinity of y = y • effectively belong to the both subsystems, while the other elements are isolated within their subsystems. Hence, at equilibrium the density of particles representing competing elements must be the same in both subsystems at the point of contact The phase volume associated with the distributions is likely to be the same on both sides A I (y • ) = A J (y • ). Existence of a single point of contact (or several points of contact that, as discussed below, do not form a loop while connecting several subsystems) changes n I but does not affect the distributions φ I (y). More than one point of contact between two systems with non-Gibbs mutations is likely to change not only n I but also the distributions φ I (y).
(b) Complete merger. The subsystems are merged into a single system with the overall stationary distribution f = f 0 (y). Unless mutations are limited to Gibbs mutations, the subsystems are likely to undergo complex adjustments changing their distributions. If the term equilibrium is used for this steady state, it should be remembered that, generally, there is no detailed balance in the system. The overall stationary distribution is inseparable: f 0 (y) may change if the contact between any two locations is severed. Note that, although unusual, inseparable systems exist in conventional thermodynamics: objects with negative heat capacity [15] may serve as an example.
Among different types of equilibrium in competitive systems, the equilibrium at a point of contact is most suitable for thermodynamic analysis, even when mutations substantially deviate from Gibbs mutations.

Entropy for equilibrium through a point of contact
Connections through a point of contact can be given different interpretations. Figure 6a, shows two subsystems with the same property y that are connected at location y = y • . Another interpretation, which is illustrated in figure 6b, is that y 1 and y 2 are internal properties of the subsystems, generally not related to each other, while the point of contact is an agreement that establishes correspondence of two locations y • 1 and y • 2 that are called open portals. Particles can freely move between these portals through the bridge connecting the portals. Note that subsystems can have more than one open portal (see figure  6d) as long as connections between these portals do not form a loop. Figure 6e illustrates such a loop that can make particle densities at two open portals that belong to a single subsystem inconsistent with each other. This would change the shapes of particle distributions φ I (y) within the subsystems.
The case that is most interesting from the thermodynamic perspective is shown in figure 6c: each subsystem has only one open portal -this ensures that the number of particles n I within each subsystem changes, while the subsystem distributions φ I (y) remain the same (presuming that each subsystem always converges to its internal steady state). Each portal can be connected to one or more of the portals that belong to the other subsystems. This connection is characterised by the subsystem particle numbers n I converging to their equilibrium values and by the detailed equilibrium between the subsystems (although the detailed balance is not necessarily achieved for the steady states within each subsystem). Assuming that the portal y • I of subsystem I is connected to the portal y • J of subsystem J, the equilibrium condition (40) is now rewritten as Let us consider how this equilibrium between K subsystems can be characterised by Tsallis entropy, which is defined as with constraints on the overall number of particles in the system: where A • I (y) is the effective phase volume in the subsystems defined by Here A I (y) is the true phase volume and a • I is the correcting coefficient, which depends on the location of the portal y • I . The value of a • I is determined later from the equilibrium conditions specified by (41). Note that the definition of entropy is extensive with respect to superposition of the subsystems but is generally non-extensive within each of the subsystems: As demonstrated by Abe [5], extensivity of entropy simplifies equilibrium analyses and removes the need to distinguish nominal and physical values of intensive properties. The general scaling laws that ensure extensivity of the entropy have been considered by Hanel and Thurner [7]. Here, extensivity with respect to superposition of the subsystems is simply enforced by the definition of the entropy in (45). Maximisation of S is conducted first over for the shape of ϕ I (y) under constraint that is where λ I are the Lagrangian multipliers associated with constraints (47) and the effective partition function Z • I is determined in terms of the true partition function Z I by the normalisation condition The substitution of ϕ I (y) into (46) results in the following problem of finding entropy extremum Maximisation of S in (51) yields where χ I is the competitive potential of I-th subsystem and λ • is the I-independent Lagrangian multiplier associated with fixing the overall number of particles to n in (51). The equilibrium distribution of particles between subsystems is then given by Note that equations (52) and (53) imply the detailed equilibrium χ I = χ J for any I and J. Consistency of (53) with (41) determines the correcting coefficients The constant in this equation is arbitrary (since it does not affect the equilibrium state) and can be set to unity without loss of generality. Assuming that A I = 1 and all Z I are the same, we obtain a • I = 1/φ I (y • I ) and the following expression for the overall entropy Thus, equilibrium through a point of contact results in defining the effective phase volumes of the subsystems, which are responsible for the equilibrium conditions. The competitors' perception of the phase volume of a subsystem depends on location of its portal. The perceived volume is larger and the subsystem has higher competitive potential χ I when the portal is located at lower ranks. Assuming that portal connections are consistent with true competitiveness of the elements, we conclude that more competitive subsystems tend to possess higher effective phase volume.

Intransitivity, transition to complexity and the risk/benefit dilemma
If competition becomes intransitive and intransitive triplets (4) exist, absolute ranking in not possible in such systems and there is no absolute entropy (since entropy potential is attached to the absolute ranking). Some intransitive systems may still retain local transitivity in smaller subdomains. In this case, the system may behave locally in the same way as transitive systems do and it is still possible to use local absolute ranking and local entropy. In this case, the analog of the zeroth law of thermodynamics becomes invalid allowing for intransitivity of competitive potentials, such as χ 1 ≺ χ 2 ≺ χ 3 ≺ χ 1 (consider the subsystems I = 1, 2, 3 shown in figure 5b assuming that these subsystems are connected), and for cyclic evolutions. The system shown in figure 1b is locally transitive and globally intransitive. Assuming that some positive mutations are present, this system evolves transitively in the vicinity of point A by escalating in the direction of increase of the local ranking, but the overall evolution appears to be cyclic moving from A to B then from B to C and finally from C back to A . When intransitivity becomes stronger (denser) and intransitive triplets can be found in vicinity of any point, even local evolution of the system may become inconsistent with the principles of competitive thermodynamics. In complex systems, this evolution may result in competitive degradation (a process accompanied by slow but noticeable gradual decrease of competitiveness) and in competitive cooperation (formation of structures with a reduced level of internal competition and violating the Stosszahlansatz). From the perspective of competitive thermodynamics, these processes are abnormal (see refs. [1,9] for further discussion).
In this section we consider a different example that involves punctuated evolution: for most of the time the system seems to behave transitively and escalate towards higher ranks and higher entropy. This escalation is nevertheless punctuated by occasional crisis events where the state of the system collapses to (or near to) the ground state. The system then repeats the slow growth / sudden collapse cycle. Note that only the cyclic component of evolution is considered here, while competitive evolutions may also involve a translational component (or components) and become spiral [9]. Cycles and collapses are common in real-world complex competitive systems of different kinds [17,18].
The present example of punctuated evolution is based on the risk/benefit dilemma (RBD): when comparing the available strategies, we would like to have low risk and high benefits; hence the problem two parameters: the risk is denoted by y (1) and the benefit denoted by y (2) . While high y (2) and low y (1) are most attractive, some compromises increasing risk to increase the benefit or lowering the benefit to lower the risk may be necessary. When comparing two strategies, y A and y B , the choice is performed according to the following coranking that is strategy A is preferred over strategy B when ρ (y A , y B ) > 0. We consider two choices of parameters One can easily see that choice RBD2 is transitive allowing for absolute ranking In the case RBD2 our assessment of the risk and the benefit is linear so that evolutions maximising the absolute ranking r are expected. In case RBD1, however, we tend to neglect small increases in risk and opt for higher benefits but a large increase in risk becomes the major concern that overweights even significant benefits. Choice RBD1 appears to be strongly (densely) intransitive: as illustrated in figure  7a, there are intransitive triplets (4) in vicinity of every point. Both cases RBD1 and RBD2 deploy the same mutations, which are predominately small but can be large on rare occasions. Mutations reaching the prohibited area are banned. Figure 7b shows the computational domain. The gray area y (2) > y (1) 1/3 is prohibited, reflecting the fact that one cannot have large benefits without being exposed to significant risks. The strategies superior with respect to A are in the small dark area causing the system to evolve to higher risks and higher benefits. In transitive case, the system grows to reach the equilibrium point maximising the absolute ranking r(y) and then remains in the this state of relatively high benefits and reasonable risks forever. In intransitive case, the system does not stay in equilibrium but collapses into a defensive strategy involving low risks and low benefits. The reason for this collapse is illustrated in figure 7b. the aggressive strategy A is preferred over defensive strategy C but as the system evolves even into a more aggressive strategy B, the risk associated with B becomes too high and at certain moment defensive strategy C becomes more attractive than B. This results in the collapse of the growth and rapid transition to defensive strategies.
For the transitive case, the entropy is defined by equation (26). The translational case γ = 1 with entropy potential depending linearly on ranking s y (y) = kr(y) is chosen. The parameters q = 1/Q = 1/1.2 and kq = 70 are selected to match the equilibrium distribution discussed below. The entropy definition takes the following form: In the case under consideration, the entropy is practically dominated by the ranking term and the difference between conventional logarithmic entropy and Tsallis entropy is not large. Figure 9 illustrates intransitive and transitive evolutions in the risk/benefit dilemma. The transitive branch is obtained by switching parameters from RBD1 to RBD2 at 410 time steps. The following intransitive and transitive evolutions seem to be very similar but only up to a point where maximal S is reached. The same definitions of ranking (57) and entropy (58) are used for both cases, transitive and intransitive. Then the evolutions diverge: the transitive branch remains in equilibrium state near the point of maximal entropy and maximal ranking while the intransitive branch falls down into the region of defensive strategies. Video files covering these evolutions between steps 1 and 590 is offered as an electronic supplement to this article (see Appendix for more details).
If the underlying competition rules and long-term history of the evolution are unknown, determining how a system is going to behave in the future by analysing the current distributions may be very difficult. Figure 8 illustrates this point. This figure shows the cdf of ranking r for transitive evolution (RBD2) and intransitive evolution (RBD1) at 590 time steps. Both distributions are very similar and can be approximated quite well by the q-exponential cdf (33) with Q = 1.2 and k/Q = 70.
The competitive mechanism represented by the risk/benefit dilemma, can be one of the forces enacting economic cycles in the real world. From the economic perspective, the strategies reflected by RBD2 are seen as rational behaviours of individual players (say, investment agents). The benefit is represented by returns on investments and ranking r is conventionally called utility in economics [16]. This utility weights different factors against each other and enforces transitivity of economic decisions. The intransitive strategies reflected by RBD1 would be seen by economists as semi-rational. Since risk and benefit do not represent directly comparable categories and evaluation of risk is always subject to greater uncertainty, overlooking small risks and being overly concerned with high risks is a plausible economic strategy for any individual or company. While switching from the linear RBD2 to non-linear RBD1 seems like a minor adjustment for an economic element, it has a major effect on functioning of the whole system: economic growth is interrupted by collapses and the system evolves cyclically. Competition forces the competing elements to take higher and higher risks until the risk becomes unsustainable.

Conclusions
In competitive systems with Gibbs mutations, the distributions tend to be exponential (assuming isotropy of the property space). This case is described by the strongest similarity to conventional thermodynamics and the existence of detailed balance in the system. When the distribution of mutations deviates from that of Gibbs mutations, the q-exponents become very good approximations characterising the existence of long or short tails in the distributions caused by biases in taxing and supplying. In competitive thermodynamics, this corresponds to replacing conventional Boltzmann-Gibbs entropy by Tsallis entropy.
Unlike in conventional thermodynamics, competitive systems allow different types of equilibria possessing different degrees of similarity with the conventional thermodynamic equilibrium. Competition between subsystems without exchange of mutations tends to be less stable than the connected equilibria where subsystems exchange particles through both competition and mutations. Among connected equilibria, the case of Gibbs mutations bears the highest resemblance to conventional thermodynamics. When mutations are not of the Gibbs type, the point of contact equilibrium preserves this resemblance more than the other cases. The point of contact equilibrium has been analysed using Tsallis entropy. This analysis results in equilibrium conditions determined by equivalence of competitive potentials. These potentials are linked to the introduced effective phase volumes of subsystems that depend on location of the point of contact.
The thermodynamic analogy requires transitivity of competition rules. In case of intransitive competition rules, the system may behave anomalously when considered from the perspective of competitive thermodynamics. This involves formation of structures, competitive degradations and cycles. The present work uses the example of the competitive risk/benefits dilemma and analyses the case of punctuated evolutions. For most of the time, the evolution of an intransitive competitive system, which represents the dilemma, closely resembles evolutions of transitive systems, which increase ranking and the associated entropy. At some moments, however, this evolution is punctuated and results in an abrupt collapse, which decreases ranking and the associated competitive entropy -this cannot possibly happen when the competition is transitive. Then the system starts to grow and repeats the cycle again. While consideration of competitive processes in this work is generic, similar behaviours can be found in biological, economic and other systems.
Competition is intransitive in RBD1 and transitive in RBD2. The cases are branched apart at 410 time steps with the same distribution of particles. The format of the videos is explained in figure 10. In the intransitive simulations of the risk/benefits dilemma, competition forces competitors to undertake more and more aggressive strategies, while the distribution moves from D to A. This leads to unsustainably high risk and punctuation of continuous evolution by a sudden collapse of the system by elements seeking refuge in defensive strategies near D. While the evolution is punctuated in the intransitive case, the transitive version of the simulations safely reaches equilibrium and remains there forever. In spite of principal differences, the ascending fragments of both simulations are very similar. absolute ranking     [13] distribution of the A3243G mtDNA mutation (solid line) compared to the q-exponential cdf exp Q (ky/Q) with Q = 0.8 and k/Q = 6 (dashed line). Since these mutations are deleterious, their extent is shown as negative.  Figure 6. Point of contact equilibrium when the subsystems have a) the same ranking and a direct connection, b) autonomous ranking and agreed connection through portals, c) multiple connections through portals that are unique for each subsystem, d) multiple connections through multiple portals without forming a loop, and e) multiple connections through multiple portals with a loop. Note that the last case can be inconsistent with the point of contact equilibrium.  Figure 7. Intransitivity in the risk/benefit dilemma: a) strong intransitivity of intransitive triplets A, B and C densely present in the domain; b) intransitivity of aggressive strategy B winning over A but losing to defensive strategy C, which is considered to be inferior to A.  Figure 8. Cdf for the rank distribution at 590 time steps (the same simulation as in figure  9, r max ≈ 0.65 ). Solid curve -intransitive (RBD1) simulation; dashed curve -transitive (RBD2) simulation; solid curve with dots -approximation by the q-exponent with Q = 1.2 and k/Q = 70.