Altai State Technical University, Barnaul Sobolev Institute of Mathematics, Novosibirsk

1. On Nonuniqueness of Cycles inDissipative Dynamical Systems of Chemical KineticsA.A.Akinshin,V.P.Golubyatnikov Altai State Technical University, Barnaul Sobolev Institute of Mathematics, Novosibirsk 8 June 2012, Novosibirsk. “Solitons , Collapses and Turbulence”

2. Our aim is to give a mathematical explanations and predictions of numerical experiments with periodic trajectoriesof nonlinear dynamical systems of chemical kinetics considered as models of gene networks regulated by simplecombinationsof negativeand positive feedbacks. OUTLINE: 1.Existence;2. Stability; 3. Non-Uniquenessof cycles. 2/57

3. Contemporary mathematics, v. 553, 2011. Russian Journal Numerical Analysis and Math.Modeling, v.26, N 4, 2011. Novosibirsk State University Herald, v.10, NN 1, 3, 2010; v.12, N 2, 2012. A.N.Kolmogorov, I.G.Petrovskii, N.S.PiskunovMoscow University Herald, 1937.

4. 4/57 1.Some simple gene networks models We study odd-dimensional dynamical systems Functions are smooth and monotonically decreasing. This corresponds to the negative feedbacks in the gene networks. Lemma 1. Each systemof this typehas exactly one stationary point in the positive octant: Lemma 2. is an invariant domain of the system. [(2k+1). +(2k)]

5. 5/57 We describe below the case (2k+1) Or, in vector form: dX/dt = F(X) – X. We Sincediv(F(X) – X) ≡ – (2k + 1), (2k+1)-dimensional volume of any bounded domain W in Qdecreases exponentially: Vol (W(t))=exp(–t(2k+1))·Vol(W(0)). All these domains collapse, but their attractors are not points!In most interesting cases these attractors contain cycles.

6. Non-convex 3D invariant domain incomposed by six triangle prisms Q 6/57 0

7. Theorem1.If the stationary point is hyperbolic then the system (2k+1) has at least one periodic trajectory in the invariant domain Q. The following diagram (D) shows the discrete scheme of some of the trajectories of the system(2k+1). 7/57 We reduce this invariant domain Q to the union of 4k+2 triangle prisms in order to localize the position of this cycle (“first’’cycle). Potential level=1.

8. 8/57 A trajectory and a limit cycle. ! !

9. 9/57 Below we demonstrate projections of trajectories of symmetric 5-D system onto 2-D and 3-D planes. Theorem 1′.If the dynamical system (2k+1)in the Th.1 is symmetric with respect to the cyclic permutation of the variables then the system has a cycle with corresponding symmetry.

10. 10 Projections of trajectories of 5-D systemonto the 2-D plane corresponding to

11. Projections of trajectories of the 5-D system onto the 3-D plane corresponding to the eigenvalues with positive real parts and the negative eigenvalue of the linearization matrix. 11 The blue spot shows the position of projection of the stationary point.

12. 12 The characteristic polynomial of the linearization of the system (2k+1) at the stationary point has the form Here is the product of all derivatives at the point . We arrange the eigenvalues of this linearization according to the values of their real parts: The eigenvalue is real and negative. So, If the point is hyperbolic then none of these real parts vanishes. If k=2 then

13. 13/57 symmetric dynamical system Eigenvalues of one 9-D

14. 14/57 Trajectories of 9-D symmetric system projected onto 3D-planes corresponding to different eigenvectors of the linearization of this system near the stationary point. The trajectories are contained in Q.

15. Similar results can be obtained for the systems of the types, etc. 15/57 M.Hirsch (1987).

16. 16/57 2. Stability questions (VM)=(2k+1) η >0. The eigenvalues of A can be expressed explicitly:

17. The transfer matrix Let 17/57 be the Jacobi matrix of , and let i = 1,2,… be its norm.

18. 18 Russel Smithhas shown that if then the system(VM)has a stable cycle(1987). Actually, he notes that this is not a sharp estimate!! Theorem 2.If the system (2k+1) satisfies the conditions of the theorem 1 and for some positive η then the invariant domain Q contains a stable cycle of this system. Potential level equals 1.

19. 19/57 3.Nonuniqueness of cycles in the system(2k+1) According to the Grobman-Hartmann theorem, each nonlinear dynamical system can be linearized in some neighborhood W of its hyperbolic point. Consider in W 2-D planes corresponding to pairs of the eigenvalues with positive real parts. These planes are composed by unwinding trajectories of the dynamical system (2k+1). Conjecture 1: Outside of W different 2-D planes generatedifferentcycles.

20. Let (2k+1) = p·q Then the phase portrait of the symmetric system contains p-dimensional and q-dimensional invariant planes. If the stationary point of this system is hyperbolic for its restrictions to these planes, then each of these planes contain a cycle. They are not contained in Q! (Potential level > 1.) If the conditions of Theorem 2 are satisfied, then these two cycles are stable within corresponding planes, not in the ambient phase portrait! If p ≠q, then this system has at least 3 different cycles. 20/57

21. 21/57 symmetric dynamical system(what happens in 5-D?) Eigenvalues of one 9-D

22. 22/57 Two cycles and some chaos in 5-D (!!) symmetric system

23. “two” cycles in 5-D symmetric system

24. 24/57 Projections of two different cycles of 9-D symmetric system onto 3-D planeThe second cycle is not contained in Q.Its potential level equals 3. The stationary point is at the top of the picture.

25. Projections of two different cycles of 11-D symmetric system onto two different 3-D planes left;right. 25/57

26. 26/57 Trajectory of 9-D system attracted by the first cycle Chaotic oscillations

27. 27/57 Projections of 3 cycles of 15-dimensional system onto the plane

28. 26/57 3 cycles in symmetric 15-D system Potential level = 5 P. level = 3 1

29. 29/57 Trajectory of 15-D system attracted by the first cycle first cycle chaotic oscillations

30. 30/57 Similar phenomenon chaotic oscillations→ ↓ first cycle

31. 31/57 4 cycles in 15-D symmetric system

32. 32/57 4. Model of 3-D gene network regulated by a simple combination ofnegative andpositive feedbacks system(ffΛ): ; ; Smooth monotonically decreasing for or more general unimodal function.

33. 33/57 Letbe the maximal valueof andis the inverse function to Lemma 3. Let and either or for Then the system (ffΛ) has exactly one stationary point in the positive octant. (see Appendix A.) Let be defined by

34. 33/57 Linearizationof system (1) at this point is described by the matrix with one negative eigenvalue. Its other eigenvaluesare complex. Consider the case (+)‏ Theorem 3.If the condition (+) is satisfied thenthe system (ffΛ):has at least one periodic trajectory.

35. The proof is based on existence of an invariant domain of the system (ffΛ). This is the parallelepiped 35/57 Actually, one can construct essentially smaller invariant domain (see below). Now, existence of periodic trajectories follows from the Brower fixed point theorem, as usual. Recall that

36. Trajectories of the system (ffΛ)right:left: 36/57

37. 37/57 5.More complicated gene networks modelsregulated bycombinationsof positive and negative feedbacks system(fΛΛ): ; for or other unimodal functions. ;

38. 38/57 Stationary points of the system (fΛΛ). I II III IV V

40. 40/57 Stationary points I andIIof the system(fΛΛ): I II

41. 41/57 Analogs of the theorems 1 and 2 about existence of a cycle and existence of a stable cycle hold in the neighborhoods of the stationary pointsI andIII. The stationary pointVis stable. The stationary pointsII and IV have topological index +1.

42. 42 Stationary points and cycles of the system(fΛΛ): рублей.

43. 43/57 Stationary points and cycles of the system (fΛΛ)(same parameters, other trajectories).

44. 6. Glass-Mackey-typesystems Ricker function. 44/57 Glass-Mackey function. (GM)‏ Logistic function.(L)‏ (w)=rw(–w)

45. 45/57 (Λ)‏ Each of the systems listed here has exactly 7 stationary pointsin some invariant domain If the parameters of the system are sufficiently large. The origin is also stationary point of each of these systems, but it does not seem to be so interesting.

46. 46/57 - + - + - + - Positions of the stationary points of the Glass-Mackey system (GM).

47. Topological indices of the stationary points 47/57 Topological indices of the points marked by “+” equal +1, their Conley indices are : . Topological indices of the points marked by “-” equal -1, their Conley indices are: . . The cycles of the Glass-Mackey type systemsdo appear near the stationary points with negative indices.

48. 48/57 For the stationary point marked by green minus, we have proved analogs of the theorems 1 and 2 about existence of a (stable) cycle. Numerical experiments show existence of cycles near the 1-st, 3-d and the 7-th stationary points marked by blue minuses.

49. Cyclesof the system (Λ). =10, m=5, q=0. 49/57 Similar pictures were observed in the systems (GM), (L).

50. Cyclesof the system (GM). α=4.3, γ=17.25. Projections onto the planeZ=0. 50/57 - + - + - + -