Complexity

Volume 2017, Article ID 7209762, 6 pages

https://doi.org/10.1155/2017/7209762

## On the Periods of Parallel Dynamical Systems

^{1}Research Institute of Applied Mathematics in Science and Engineering, Ciudad Real, Spain^{2}Department of Mathematics, University of Castilla-La Mancha, Albacete, Spain

Correspondence should be addressed to Jose C. Valverde; se.mlcu@edrevlav.esoj

Received 17 June 2016; Revised 20 October 2016; Accepted 16 November 2016; Published 12 January 2017

Academic Editor: David Arroyo

Copyright © 2017 Juan A. Aledo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

In this work, we provide conditions to obtain fixed point theorems for parallel dynamical systems over graphs with (Boolean) maxterms and minterms as global evolution operators. In order to do that, we previously prove that periodic orbits of different periods cannot coexist, which implies that Sharkovsky’s order is not valid for this kind of dynamical systems.

#### 1. Introduction

It is beyond dispute that mathematical modeling of real phenomena is one of the main tasks of applied mathematics nowadays. Nevertheless, once a model fits a phenomenon, it is necessary to analyze such a model theoretically in order to apply the results found to any particular case under study.

In the last two decades,* graph dynamical systems* (GDS) have been revealed as an important tool for the mathematical modeling of computer processes. For the formulation of such models, these systems are decomposed into the lower unities of aggregation, called* entities*, where each entity has a (numerical) value which represents its state. When the possible states of each entity are activated or deactivated, the state value of an entity is formalized by considering . Indeed, the state values can belong to a more general Boolean algebra [1]. The relations among entities are represented by a graph called the* dependency graph* of the system, while the evolution or update of the system is implemented by* local (Boolean) functions* which together constitute a global (evolution) operator. That is, for the dynamic evolution of the state of any entity, the corresponding local function acts on such state and those ones corresponding to entities related to it.

When the states of the entities are updated in a synchronous manner, the system is called a* parallel dynamical system* (PDS) [2, 3], while if they are updated in an asynchronous way, the system is named* sequential dynamical system* (SDS) [4].

Other tools for mathematical modeling as cellular automata (CA) [5–11] or Boolean networks (BN) [12–14] are particular cases of PDS [15]. It shows the enormous versatility of this new paradigm, which has applications to different branches of sciences as biology [12, 16, 17], ecology [18–20], mathematics [21, 22], cryptography [23, 24], physics [25], and chemistry [26]. (The abbreviations GDS, PDS, SDS, CA, and BN will be written for the singular and plural forms of the corresponding terms, since it seems better from an aesthetic point of view.)

The main aim regarding the study of a dynamical system is to give a complete characterization of its orbit structure [27]. That is, to derive as much information as possible about the phase diagram, based on the initial states and the evolution operator of the system. Specifically, as pointed out in [28], it means to determine the length and number of coexisting limit cycles and which different initial states arrive in the same limit cycle.

Concerning these questions, in [3] PDS and SDS are studied, considering the simplest Boolean functions OR (resp., AND) and NOR (resp., NAND) as global evolution operators. These results are extended in [2] for PDS with any general Boolean maxterm (resp., minterm) as global functions; in particular, it is proved that for such PDS uniquely periodic or eventually periodic orbits of period lower than or equal to 2 are possible.

The results in [2] left open some questions among those set out by Derrida and Pomeau [28] regarding the phase diagram of such systems. One of them consists in studying the coexistence of periodic orbits of different periods for PDS with general maxterm (resp., minterm) functions as evolution operators in the same sense of Sharkovsky’s Theorem [29, 30]. In other words, it consists in guessing whether the existence of certain periods implies the appearance of other ones in a similar way to Sharkovsky’s order.

Concerning this question, in this work, we show that Sharkovsky’s Theorem and the order established with it do not apply in the context of PDS. On the contrary, we demonstrate that the coexistence of periodic orbits with different periods is not possible.

The main objective of this paper is to establish results for PDS with general maxterm (resp., minterm) functions as evolution operators in the sense of the fixed point theorem by Banach. In fact, as we show in Example 8, there exists PDS with several fixed points and so it naturally arises the problem of determining conditions which guarantee the uniqueness of a fixed point.

The results, which are initially given for state values of the entities in the basic Boolean algebra , are later generalized for state values belonging to any general Boolean algebra with elements, , .

The specific systems treated in this work appear when modeling real physical phenomena like the dynamics of genes. In particular, when studying gene regulatory networks, it is observed that the binding sites of genes regulated by more than one gene exhibit synergistic effects. This could be modeled by conjunctive Boolean networks, but they do not allow modeling inhibitory regulation, what needs from minterms evolution operators to formalize and analyze the corresponding dynamics (see [31]).

The paper is organized as follows. Section 2 gives some preliminaries concerning PDS on maxterm and minterm Boolean functions. In Section 3, we demonstrate that the coexistence of periodic orbits with different periods is not possible, so showing that Sharkovsky’s Theorem does not apply in this context. In Section 4, we establish a fixed point theorem for a PDS. In Section 5, we extend our results to PDS with state values belonging to any general Boolean algebra with elements. Finally, Section 6 provides some conclusions and future research directions.

#### 2. Preliminaries

In this section we revise some basic notions regarding* parallel dynamical system* (PDS) (see [15] and the references therein for the details).

As is well known, a PDS is totally determined by its* dependency graph* and its* evolution* or* updating function*, which we briefly describe below.

The dependency graph is given by a (undirected) graph , where is the vertex set and is the edge set. We will represent by the* state* of the vertex , meaning that the vertex is* activated* () or* deactivated* ().

We will assume that is connected; otherwise the reasonings and results which we will provide can be easily adapted by working on its connected components.

For every vertex , we will consider the set of vertices that are adjacent to ; namely, Then, given , we define

The evolution or update of the system is implemented by local functions which, for every , act only over .

*Definition 1. *Let be a graph on . Then a mapwhere is the updated state of the entity by applying a local function over the states of the entities in , constitutes a discrete dynamical system called* parallel dynamical system* (PDS) over .

Every is often the restriction of the global function acting only over the state of the entities in , as we are going to suppose in this document. Nevertheless, each could be seen as another function : such that .

In particular, the (global) evolution operator of a PDS can be given by means of a* Boolean function* of variables:where is obtained from using the logical AND (), the logical OR (), the logical NOT (), and the elements . Particular cases of Boolean functions are* maxterms* and* minterms*. Recall that a maxterm (resp., minterm) of variables is a Boolean function such as where or .

In this work, we deal with this special kind of PDS.

*Definition 2. *Let be a (undirected) graph with . Then the following (evolution) operatorwhere is the updated state of the entity by applying locally the Boolean function over the states of the vertices in , constitutes a PDS over with the Boolean function as evolution operator, which will be denoted by or -PDS when specifying that the dependency graph is not necessary.

In particular, we will represent by MAX-PDS (resp., MIN-PDS) a PDS whose evolution operator is the maxterm MAX (resp., the minterm MIN).

#### 3. Coexistence of Periodic Orbits in PDS over Graphs

In [2] it was shown that the periodic orbits of a MAX-PDS (resp., MIN-PDS) are fixed points or 2-periodic orbits. Moreover, when the maxterm MAX (resp., minterm MIN) has all the variables in its direct form, then only (eventually) fixed points can appear, while if it has all the variables in its complemented form, then only (eventually) 2-periodic orbits are possible.

However, some important questions remained open. One of them consists in studying the coexistence of periodic orbits with different periods in the same sense of Sharkovsky’s Theorem [29, 30] for PDS with general maxterm (resp., minterm) functions as evolution operators. In other words, it consists in guessing whether the existence of certain periods implies the appearance of other ones in a similar way to Sharkovsky’s order. Concerning this question, in this section, we prove that Sharkovsky’s Theorem and the order established with it do not apply in the context of PDS. In fact, we demonstrate that periodic orbits with different periods cannot coexist.

In order to do that, next, we describe which structure a PDS must have in order to admit (eventually) fixed points. We will outline the reasonings for the case of a MAX-PDS, although all of them can be dually rewritten for a MIN-PDS.

Theorem 3. *Let MAX-PDS be the PDS over a dependency graph associated with the maxterm MAX. All the periodic orbits of this system are fixed points if, and only if, every entity whose state variable in MAX is in complemented form is adjacent to an entity whose state variable in MAX is in directed form.*

*Proof. *Let us consider the set (resp., ) of vertices in such that the corresponding variables in MAX appear in direct (resp., complemented) form.

First, assume that all the periodic orbits of this system are fixed points.

In the case , the system has two fixed points: the one with all the variables activated and the one with all of them deactivated (see [2, 3]).

Thus, assume that and take a fixed point, where represents the (fixed) value of the vertex . Note that, for all , it must be . Otherwise (i.e., if ), it would change to after the following iteration.

Suppose that there exists such that . In such a case and so for every it is . But this is not possible, since in the following iteration the value of would change to 0, which is a contradiction.

To prove the converse implication, let us suppose that for all it is . We will write to indicate the state value of the entity after iterations of the evolution operator MAX. Thus, let us consider an arbitrary initial value for the variables . Since the dependency graph is finite (and so is the state space), note that after a certain number of iterations, let us say , the states of all the vertices become fixed. Let us take and let us prove that it is . In fact, let us suppose that and take . Then it would be (since we are assuming that the state value of is fixed from the iteration ). But then, since , it must be , which is a contradiction. Thus, it is for all and these state values do not change as can be easily inferred from the reasoning above.

Therefore, all the variables of the system become fixed after iterations and the proof finishes.

Dually, we have the following.

Theorem 4. *Let MIN-PDS be the PDS over a dependency graph associated with the minterm MIN. All the periodic orbits of this system are fixed points if, and only if, every entity whose state variable in MIN is in complemented form is adjacent to an entity whose state variable in MIN is in directed form.*

Analogously, dual results can be obtained for all the following results. For the sake of conciseness, we will not state them.

Since all the periodic orbits of the system are fixed points or 2-periodic orbits (see [2]), as a consequence of Theorem 3 we have the following.

Theorem 5. *Let MAX-PDS be the PDS over a dependency graph associated with the maxterm MAX. All the periodic orbits of this system are 2-periodic orbits if, and only if, there exists an entity whose state variable in MAX is in complemented form which is not adjacent to any entity whose state variable in MAX is in directed form.*

*Proof. *First, assume that all the periodic orbits of this system are 2-periodic orbits. If the theses were not true, that is, if for all , it is , then from Theorem 3 we have that all the periodic orbits of this system are fixed points, which is a contradiction.

Conversely, let us suppose that there exists such that . If the system has a fixed point, reasoning as in the proof of Theorem 3, we get that every entity whose state variable in MAX is in complemented form is adjacent to an entity whose state variable in MAX is in directed form, which is a contradiction. Hence, all the periodic orbits of the system must be 2-periodic orbits.

*Remark 6. *Note that the cases when the evolution operator of the system is AND, OR, NAND, and NOR studied in [3] can be immediately obtained as particular cases of Theorems 3 and 5.

As a direct consequence of Theorems 3 and 5, we get the main result of this section.

Theorem 7. *Let MAX-PDS be the PDS over a dependency graph associated with the maxterm MAX. Then (eventually) fixed points and (eventually) 2-periodic orbits cannot coexist.*

#### 4. A Fixed Point Theorem for PDS over Graphs

As commented in the introduction, our main objective is to obtain a fixed point theorem for PDS. Observe that, although fixed points and 2-periodic orbits cannot coexist, there are PDS whose state spaces contain more than one fixed point as shown in the following example.

*Example 8 (PDS with more than a fixed point). *Let us consider the graph with and , and let us take the evolution operator The fixed points of this PDS are , , and .

Thus, it would be desirable to find conditions to assure that the system has a unique fixed point. In order to do that, assume that the MAX-PDS has at least a fixed point; that is, for every vertex it is (see Theorem 3).

Let be the connected components which result from when we remove all the vertices in and the edges adjacent to those vertices. We will say that a vertex in and a connected component are adjacent if is adjacent to any vertex of . Let us consider the subsets and of given by

Let be a fixed point of the system. As we have showed in the proof of Theorem 3, for all , it is .

Furthermore, if , , is such that , then for all it is . In order to see that, note that in the fixed point all the vertices in must be either activated or deactivated. Let us take . Then, if it is for all and since it is for all , then the state of will change to 0 after the next iteration, which is a contradiction.

Therefore, if for every , , it is , then the system has a unique fixed point: the one with all the variables activated.

On the other hand, if for a , , it is , then in a fixed point two situations are possible: either all the vertices in are activated or all of them are deactivated. Regarding this last comment, we must point out that given , if , are the connected components adjacent to , then not all these components can be deactivated simultaneously in the fixed point; otherwise the value of would change from 1 to 0 in the following iteration.

In view of the explanations above, we have proved the following result.

Theorem 9. *Let MAX-PDS be the PDS over a dependency graph associated with the maxterm MAX and assume that for all it is . This PDS has a unique fixed point if, and only if, for every , , it is . In this situation, the only fixed point is the one with all the vertices activated.*

#### 5. Fixed Point Theorem for PDS over Graphs with General Boolean States

In the recent work [1], a generalization of PDS over graphs is introduced by considering that the states of the entities can take values in an arbitrary Boolean algebra with elements, , . We represent such a PDS as (resp., ) to indicate the dependency graph and the maxterm MAX (resp., minterm MIN) taken as evolution operator. Note that the previous study corresponds to the particular case where ; that is, .

In such a work, it was proved that the periodic orbits of these more general PDS are also fixed points or 2-periodic orbits. The key idea to prove this result consists in using the Stone Representation Theorem for Boolean Algebras in a suitable way to* decompose* this general dynamical system into PDS where the entities take values in (see [1, Section ] for the details).

Then, from our study above, a straightforward reasoning allows checking that all the results in Sections 3 and 4 are also true for PDS where the states of the entities take values in an arbitrary Boolean algebra. That is, (eventually) fixed points and (eventually) 2-periodic orbits cannot coexist.

Theorem 10. *Let be the PDS over a dependency graph associated with the maxterm MAX. Then (eventually) fixed points and (eventually) 2-periodic orbits cannot coexist.*

In particular, we have a fixed point theorem in this more general context as stated below

Theorem 11. *Let be the PDS over a dependency graph associated with the maxterm and assume that for all it is . This PDS has a unique fixed point if, and only if, for every , , it is . In this situation, the only fixed point is the one with the state of all the vertices equal to .*

*Proof. *Since this general dynamical system can be decomposed into PDS where the entities take values in and all of these PDS are defined on the same dependency graph , the system has a unique fixed point, resembling the one of , just identifying the value with and with .

#### 6. Conclusions and Future Research Directions

This paper provides a significant advance in the description of the phase diagram of PDS over graphs on (global) Boolean operators, by analyzing the veracity of two classical theorems in dynamical systems theory in this framework.

In particular, we demonstrate that the coexistence of periodic orbits with different period is impossible in such systems. This means that Sharkovsky’s Theorem for discrete dynamical systems associated with continuous functions is not true in the context of PDS. We also obtain a fixed point theorem for this kind of dynamical systems, providing conditions which assure the uniqueness of a fixed point for the PDS. In other words, we characterize the PDS where any initial state evolves to an equilibrium which, for real models, is a meaningful piece of information.

Finally, by using the techniques provided in this work, the results here obtained invite studying the same problems in other types of GDS, as those on directed graphs [32] and/or those where the evolution operator is given by means of local independent Boolean functions [33].

Other tools for mathematical modeling as cellular automata (CA) or Boolean networks (BN) are particular cases of PDS when considering the same evolution operator. That is, CA and BN on maxterm (resp., minterm) Boolean functions over undirected graphs are particular cases of PDS with maxterm (resp., minterm) that we treat in this work. It shows the applicability to different branches of sciences and engineering that this new paradigm of PDS can have.

#### Competing Interests

The authors declare that they have no competing interests.

#### Acknowledgments

Juan A. Aledo was supported by Junta de Comunidades de Castilla-La Mancha (Grant PEII-2014-049-P). Luis G. Diaz, Silvia Martinez, and Jose C. Valverde were supported by FEDER OP2014-2020 of Castilla-La Mancha (Spain) under the Grant GI20163581. Silvia Martinez and Jose C. Valverde were also supported by the Ministry of Economy and Competitiveness of Spain (Grant MTM2014-51891-P).

#### References

- J. A. Aledo, S. Martinez, and J. C. Valverde, “Graph dynamical systems with general Boolean states,”
*Applied Mathematics & Information Sciences*, vol. 9, no. 4, pp. 1803–1808, 2015. View at Publisher · View at Google Scholar · View at MathSciNet - J. A. Aledo, S. Martinez, F. L. Pelayo, and J. C. Valverde, “Parallel dynamical systems on maxterm and minterm Boolean functions,”
*Mathematical and Computer Modelling*, vol. 35, pp. 666–671, 2012. View at Google Scholar - C. L. Barrett, W. Y. Chen, and M. J. Zheng, “Discrete dynamical systems on graphs and Boolean functions,”
*Mathematics and Computers in Simulation*, vol. 66, no. 6, pp. 487–497, 2004. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - H. S. Mortveit and C. M. Reidys,
*An Introduction to Sequential Dynamical Systems*, Springer, New York, NY, USA, 2007. View at MathSciNet - A. Ilachinski,
*Cellular Automata. A Discrete Universe*, World Scientific, Singapore, 2001. - J. Kari, “Theory of cellular automata: a survey,”
*Theoretical Computer Science*, vol. 334, no. 1–3, pp. 3–33, 2005. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - J. Maynard-Smith,
*Evolution and the Theory of Games*, Cambridge University Press, Cambridge, UK, 1982. - P. Sarkar, “A brief history of cellular automata,”
*ACM Computing Surveys*, vol. 32, no. 1, pp. 80–107, 2000. View at Publisher · View at Google Scholar · View at Scopus - J. L. Schiff,
*Cellular Automata: A Discrete View of the World*, John Wiley & Sons, New York, NY, USA, 2008. View at MathSciNet - S. Wolfram, “Statistical mechanics of cellular automata,”
*Reviews of Modern Physics*, vol. 55, no. 3, pp. 601–644, 1983. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - S. Wolfram,
*Cellular Automata and Complexity*, Addison-Wesley, New York, NY, USA, 1994. - S. A. Kauffman, “Metabolic stability and epigenesis in randomly constructed genetic nets,”
*Journal of Theoretical Biology*, vol. 22, no. 3, pp. 437–467, 1969. View at Publisher · View at Google Scholar · View at Scopus - S. A. Kauffman,
*Origins of Order: Self-Organization and Selection in Evolution*, Oxford University Press, Oxford, UK, 1993. - I. Shmulevich and S. A. Kauffman, “Activities and sensitivities in Boolean network models,”
*Physical Review Letters*, vol. 93, no. 4, Article ID 048701, 4 pages, 2004. View at Publisher · View at Google Scholar · View at Scopus - J. A. Aledo, S. Martinez, and J. C. Valverde, “Parallel dynamical systems over graphs and related topics: a survey,”
*Journal of Applied Mathematics*, vol. 2015, Article ID 594294, 14 pages, 2015. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - A. Deutsch and S. Dormann,
*Cellular Automaton Modelling of Biological Pattern Formation*, Birkhäuser Boston, Cambridge, Mass, USA, 2004. - Z. Toroczkai and H. Guclu, “Proximity networks and epidemics,”
*Physica A: Statistical Mechanics and Its Applications*, vol. 378, no. 1, pp. 68–75, 2007. View at Publisher · View at Google Scholar · View at Scopus - P. Hogeweg, “Cellular automata as a paradigm for ecological modeling,”
*Applied Mathematics and Computation*, vol. 27, no. 1, pp. 81–100, 1988. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - U. Dieckman, R. Law, and J. A. J. Metz,
*The Geometry of Ecological Interactions. Simplifying Spatial Complexity*, Cambridge University Press, Cambridge, UK, 2000. - J. Hofbauer and K. Sigmund,
*Evolutionary Games and Population Dynamics*, Cambridge University Press, Cambridge, UK, 2003. - G. Chiaselotti, T. Gentile, and P. A. Oliverio, “Parallel and sequential dynamics of two discrete models of signed integer partitions,”
*Applied Mathematics and Computation*, vol. 232, pp. 1249–1261, 2014. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - G. Chiaselotti, G. Marino, P. A. Oliverio, and D. Petrassi, “A discrete dynamical model of signed partitions,”
*Journal of Applied Mathematics*, vol. 2013, Article ID 973501, 10 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - A. Fúster-Sabater and P. Caballero-Gil, “On the use of cellular automata in symmetric cryptography,”
*Acta Applicandae Mathematicae*, vol. 93, no. 1–3, pp. 215–236, 2006. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - S. Nandi, B. K. Kar, and P. Pal Chaudhuri, “Theory and applications of cellular automata in cryptography,”
*IEEE Transactions on Computers*, vol. 43, no. 12, pp. 1346–1357, 1994. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - B. Chopard and M. Droz,
*Cellular Automata for Modeling Physics*, Cambridge University Press, Cambridge, UK, 1998. - L. B. Kier, P. G. Seybold, and C.-K. Cheng,
*Modeling Chemical Systems Using Cellular Automata*, Springer, New York, NY, USA, 2005. - S. Wiggins,
*Introduction to Applied Nonlinear Dynamical Systems and Chaos*, Springer, New York, NY, USA, 1990. View at Publisher · View at Google Scholar · View at MathSciNet - B. Derrida and Y. Pomeau, “Random networks of automata: a simple annealed approximation,”
*Europhysics Letters*, vol. 1, no. 2, pp. 45–49, 1986. View at Publisher · View at Google Scholar - A. N. Sharkovsky, “Co-existence of cycles of a continuous mapping of a line onto itself,”
*Ukrainian Mathematical Journal*, vol. 16, pp. 61–71, 1964. View at Google Scholar - P. Štefan, “A theorem of Šarkovskii on the existence of periodic orbits of continuous endomorphisms of the real line,”
*Communications in Mathematical Physics*, vol. 54, no. 3, pp. 237–248, 1977. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - I. Park, K. H. Lee, and D. Lee, “Inference of combinatorial Boolean rules of synergistic gene sets from cancer microarray datasets,”
*Bioinformatics*, vol. 26, no. 12, pp. 1506–1512, 2010. View at Publisher · View at Google Scholar · View at Scopus - J. A. Aledo, S. Martinez, and J. C. Valverde, “Parallel dynamical systems over directed dependency graphs,”
*Applied Mathematics and Computation*, vol. 219, no. 3, pp. 1114–1119, 2012. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - J. A. Aledo, S. Martinez, and J. C. Valverde, “Parallel discrete dynamical systems on independent local functions,”
*Journal of Computational and Applied Mathematics*, vol. 237, no. 1, pp. 335–339, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus