## Heterogeneous, Weakly Coupled Map Lattices

M Dolores Sotelo Herrera , Jes´s San Mart´ ın , Mason A. Porter a a u a b

- a Departamento de Matem´ atica Aplicada, E.T.S.I.D.I., Universidad Polit´cnica de Madrid, e Madrid, Spain
- b Oxford Centre for Industrial and Applied Mathematics, Mathematical Institute, University of Oxford, Oxford, United Kingdom

## Abstract

Coupled map lattices (CMLs) are often used to study emergent phenomena in nature. It is typically assumed (unrealistically) that each component is described by the same map, and it is important to relax this assumption. In this paper, we characterize periodic orbits and the laminar regime of type-I intermittency in heterogeneous weakly coupled map lattices (HWCMLs). We show that the period of a cycle in an HWCML is preserved for arbitrarily small coupling strengths even when an associated uncoupled oscillator would experience a period-doubling cascade. Our results characterize periodic orbits both near and far from saddle-node bifurcations, and we thereby provide a key step for examining the bifurcation structure of heterogeneous CMLs.

Keywords: heterogeneous CML, intermittency, period preservation, synchronization

## 1. Introduction

Numerous phenomena in nature - such as human waves in stadiums [1] and flocks of seagulls [2] - result from the interaction of many individual elements, and they can exhibit fascinating emergent dynamics that cannot arise in individual or even small numbers of components [3]. In practice, however, a key assumption in most such studies is that each component is described by the same dynamical system. However, systems with heterogeneous elements are much more common than homogeneous systems. For example, a set of interacting cars on a highway that treats all cars as the same ignores different types of cars (e.g., their manufacturer, their age, different levels of intoxication among the drivers, etc.), and a dynamical system that governs the behavior of different cars could include different parameter values or even different functional forms entirely for different cars. Additionally, one needs to use different functional forms to address phenomena such as interactions among cars, traffic

Email address: jsm@dfmf.uned.es (Jes´s San Mart´ ın) u

December 11, 2015

lights, and police officers. Unfortunately, because little is known about heterogeneous interacting systems [4, 5], the assumption of homogeneity is an important simplification that allows scholars to apply a plethora of analytical tools. Nevertheless, it is important to depart from the usual assumption of homogeneity and examine coupled dynamical systems with heterogeneous components.

The study of coupled map lattices (CMLs) [6, 7] is one important way to study the emergent phenomena (e.g., cooperation, synchronization, and more) that can occur in interacting systems. CMLs have been used to model systems in numerous fields that range from physics and chemistry to sociology, economics, and computer science [7, 8, 9, 10, 11]. In a CML, each component is a discrete dynamical system (i.e., a map). There are a wealth of both theoretical and computational studies of homogeneous CMLs [6, 7, 12, 13, 14, 15, 16, 17, 18], in which the interacting elements are each governed by the same map. Such investigations have yielded insights on a wide variety of phenomena. As we mentioned above, the assumption of homogeneity is a major simplification that often is not justifiable. Therefore, we focus on heterogeneous CMLs, in which the interacting elements are governed by different maps or by the same map with different parameter values. The temporal evolution of a heterogeneous coupled map lattice (CML) with p components is given by

/negationslash

$$X _ { i } ( n + 1 ) = f _ { R _ { i } } ( X _ { i } ( n ) ) + \varepsilon \sum _ { \substack { h = \frac { 1 } { h } \\ h \neq i } } ^ { p } f _ { R _ { h } } ( X _ { h } ( n ) ) \,, \quad i \in \{ 1, \dots, p \} \,, \quad ( 1 )$$

where X n i ( ) represents the state of the entity at instant n at position i of a lattice and ε &gt; 0 weights the coupling between these entities. We consider entities in the form of oscillators, where the i th oscillator evolves according to the map

$$X _ { i } ( n + 1 ) = f _ { R _ { i } } ( X _ { i } ( n ) ) \,, \quad i \in \{ 1, \dots, p \} \,, \quad \quad ( 2 )$$

where the f R i are, in general, different functions that depend on a parameter R i (where i ∈ { 1 , . . . , p } ). We assume that each f R i is a C 2 unimodal function that depends continuously on the parameter R i with a critical point C at R i . As usual, f m means that f is composed with itself m times. If an uncoupled oscillator X n i ( ) takes the value x i,n , then the evolution of this value under the map is x i,n +1 = f R i ( x i,n ).

In this paper, we examine heterogeneous, weakly coupled map lattices (HWCMLs). Weakly coupled systems can exhibit phenomena (e.g., phase separation because of additive noise [19]) that do not arise in strongly coupled systems, and one can even use weak coupling along with noise to fully synchronize nonidentical oscillators [20]. Thus, it is important to examine HWCMLs, which are amenable to perturbative approaches. In our paper, we characterize periodic orbits both far away from and near saddle-node (SN) bifurcations. Understanding periodic orbits is interesting by itself and is also crucial for achieving an understanding of more complicated dynamics (such as chaos) [21, 22]. We then characterize the laminar regime of type-I intermittency in our HWCMLs. Finally, we summarize our results and briefly comment on applications.

## 2. Theoretical Results

Before discussing our results, we need to define some notation. Let x i,n R | i denote the points in a periodic orbit of the i th uncoupled oscillator with control parameter R i . The parameter value r i is a bifurcation value of R i for the i th map, so x i,n r | i denotes the points in a periodic orbit at this parameter value.

Suppose that R i = r i + ε α , where ε is the same as in the coupling term of the CML (1) and α ∈ (0 , ∞ ) is a constant. We seek to derive results that are valid at size O ε ( ). We need to consider the following situations:

α &lt; 1 In this case, when we expand to size O ε ( ), the coupling term does not contribute at all. Therefore, the oscillators in (1) behave as if they were uncoupled at this order of the expansion.

- α &gt; 1 In this case, the coupling term controls the ε bifurcation terms. Thus, to size O ε ( ), we cannot study the behavior of the bifurcation.

α = 1 In this case, we are considering a perturbation of the same size as the coupling term, and we can simultaneously study the coupling and the bifurcation analytically.

/negationslash

To study orbits close to bifurcation points, we thus let R i = r i + , where ε ε is the same as in the coupling term of the CML (1). In our numerical simulations (see Section 3), we will also briefly indicate the effects of considering α = 1 (see Section 3.3).

## 2.1. Study of the CML Far from and Close to Saddle-Node Bifurcations

In this section, we examine heterogeneous CMLs in which the uncoupled oscillators have periodic orbits either far from or near SN bifurcations. As periodic orbits exhibit different dynamics from each other depending on whether they are near or far from SN bifurcations [23, 24], it is important to distinguish between these two situations.

A periodm SN orbit is a periodic orbit that is composed of m 'SN points' of the composite map f m r i . Each of these m SN points is a fixed point of f m r i at which f m r i undergoes an SN bifurcation. Periodm SN orbits play an important role in a map's bifurcation structure, because they occur at the beginning of periodic windows in bifurcation diagrams. Studying them is thus an important step towards examining the general bifurcation structure of a map.

When f m r i undergoes an SN bifurcation, the map f r i has two properties that we highlight. Let { x i, 1 | r i , x i, 2 | r i . . . , x i,m r | i } be an periodm SN orbit. It then follows that

1.

Consequently, orbits that are near an SN orbit satisfy

$$\frac { \partial f _ { r _ { i } } ^ { m } } { \partial x } ( x _ { i, j | r _ { i } } ) = 1 = \prod _ { k = j } ^ { j + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, k | r _ { i } } ) \\ \text{tlv. orbits that are near an SN orbit satisfy}$$

$$1 - \prod _ { k = j } ^ { j + m - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) = o ( 1 ) \,.$$

By contrast, if

$$\iota, \, \Pi & 1 - \prod _ { k = j } ^ { j + m - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) = O ( 1 ) \,. & ( 4 ) \\ t \, \text{an orbit is the far from} " a \, \text{SN orbit}$$

we say that an orbit is 'far from' a SN orbit.

2. Because f r i has a critical point at C, so does f m r i . Suppose that x i,n r | i is the point of the SN orbit that is closest to C. As ∂f r i ∂x (C) = 0, for sufficiently large periods, we can find SN orbits with arbitrarily small ∣ ∣ ∣ ∣ ∂f r i ∂x ( x i,n r | i ) ∣ ∣ ∣ ∣ (see Fig. 1), and in particular we can find examples where ∣ ∣ ∣ ∣ ∂f r i ∂x ( x i,n r | i ) ∣ ∣ ∣ ∣ &lt; ε . We use the term small-derivative SN orbits for such orbits. Additionally, a small-derivative SN orbit includes points that are not close to the critical point C, so that ∂f ∂x ( x i ) = O (1) in general, and the associated terms cannot be neglected.

The overall bifurcation pattern in a typical unimodal map of the interval is topologically equivalent to the bifurcation pattern in any other typical unimodal map of the interval [25], so it is sensible to focus on a particular such map. The standard choice for such a map is the logistic map. Orbits of any period occur in the logistic map, which contains infinitely many small-derivative SN orbits. In particular, such orbits include the periodq SN orbits from which supercycles with symbol sequences CRL q -2 originate. 1 Given this fact and the broad applicability of results for the logistic map, we note that our results are relevant in numerous situations.

## Lemma 1

Let | ε | &lt; 1 in the CML (1), and suppose that the map f q i R i has an SN bifurcation at R i = r i , such that the associated SN orbit of f r i is a smallderivative SN orbit. Additionally, suppose that R i = r i + ε for i ∈ { 1 , . . . , s } , but that R i for i ∈ { s +1 , . . . , p } are far away from r i . Consider the following initial conditions:

- · For i ∈ { 1 , . . . , s } , let X n i ( ) = x i,n r | i + εA i,n + O ε ( 2 ), where x i,n r | i is the point of the SN orbit closest to the critical point C of f R i at R i = r i .
- · For i ∈ { s +1 , . . . , p } , let X n i ( ) = x i,n R | i + εA i,n + O ε ( 2 ).

The temporal evolution of the CML (1) is then given by

1 Recall that a 'supercycle' is a periodic orbit that includes C; if its period is q (i.e., if f q R ( C ) = C for some parameter value R ), then ∂f q R ∂x ( x k ) = 0 for all points x k in the supercycle.

## 1. For i ∈ { 1 , . . . , s } ,

/negationslash

$$X _ { i } ( n + 1 ) = x _ { i, n + 1 | r _ { i } } + s \left ( \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n | r _ { i } } ) + \sum _ { \substack { h = 1 \\ h \neq i } } ^ { s } x _ { h, n + 1 | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + 1 | R _ { h } } \right ) + O ( \varepsilon ^ { 2 } ) \,, \quad m = 1 \,,$$

/negationslash

/negationslash

$$X _ { i } ( n + m ) = x _ { i, n + m | r _ { i } } + \left [ \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n + m - 1 | r _ { i } } ) + \sum _ { k = n } ^ { n + m - 2 } \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, k | r _ { i } } ) \prod _ { l = k + 1 } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \\ + \sum _ { k = n + 1 } ^ { n + m - 1 } \left ( \left ( \sum _ { \substack { s \\ h = 1 \\ h \neq i } } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) ^ { n + m - 1 } \prod _ { l = k } ^ { \partial f _ { r _ { i } } } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \right ) \\ + \sum _ { \substack { h = 1 \\ h \neq i } } ^ { s } x _ { h, n + m | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + m | R _ { h } } \right ] \varepsilon + O ( \varepsilon ) \,, \quad m \in \{ 2, \dots, q \}.$$

## 2. For i ∈ { s +1 , . . . , p } ,

/negationslash

$$X _ { i } ( n + 1 ) = x _ { i, n + 1 | R _ { i } } + \varepsilon \left ( \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, n | R _ { i } } ) A _ { i, n } + \sum _ { h = 1 } ^ { s } x _ { h, n + 1 | R _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i \\ ( 7 ) } } ^ { p } x _ { h, n + 1 | R _ { h } } \right ) + O ( \varepsilon ^ { 2 } ) \,, \quad m = 1 \,, \\ - \quad \dots \quad.$$

/negationslash

$$X _ { i } ( n + m ) = x _ { i, n + m | R _ { i } } + \left [ \prod _ { k = n } ^ { n + m - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) A _ { i, n } + \sum _ { k = n + 1 } ^ { n + m - 1 } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } \sum _ { h = s + 1 } ^ { p } + x _ { h, k | R _ { h } } \right ) \\ \prod _ { \begin{array} { c } n + m - 1 \\ \end{array} } ^ { n + m - 1 } \frac { \partial f _ { R _ { i } } } { \wedge } ( x _ { i, l | R _ { i } } ) + \stackrel { s } { \lambda } ^ { s } x _ { h \, n + m | r _ { i } } + \end{array} \right ] _ { \begin{array} { c } n } \cdot \\ \end{array} } \right ] _ { \varepsilon } + O ( \varepsilon ^ { 2 } ) \,. \quad m \in \{ 2. \dots. q \}.$$

/negationslash

$$\prod _ { l = k } ^ { n + m - 1 } \frac { \L L } { \partial x } ( x _ { i, l | R _ { i } } ) \Big ) + \sum _ { h = 1 } ^ { s } x _ { h, n + m | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, n + m | R _ { h } } \Big ] \varepsilon + O ( \varepsilon ^ { 2 } ) \,, \quad m \in \{ 2, \dots, q \}.$$

## Proof of Lemma 1

We proceed by induction. Substitute X n i ( ) = x i,n r | i + εA i,n + O ε ( 2 ) and X n i ( ) = x i,n R | i + εA i,n + O ε ( 2 ) into equation (1) and expand in powers of ε . Note that we need to consider i ∈ { 1 , . . . , s } and i ∈ { s +1 , . . . , p } separately.

- 1. We initiate the iteration at the point x i,n r | i of the SN orbit closest to the critical point C. Because we have a small-derivative SN orbit, ∂f r i ∂x ( x i,n r | i ) is arbitrarily small, although this is not true for other points in the SN orbit.

For i ∈ { 1 , . . . , s } , we have

$$X _ { i } ( n + 1 ) = f _ { \underset { \varepsilon } { r _ { i } } } ( x _ { i, n | r _ { i } } ) + \varepsilon \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, n | r _ { i } } ) A _ { i, n } + \varepsilon \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n | r _ { i } } )$$

/negationslash

$$+ \varepsilon \sum _ { \substack { h = 1 \\ h \neq i } } ^ { s } f _ { r _ { h } } ( x _ { h, n | r _ { h } } + O ( \varepsilon ) ) + \sum _ { h = s + 1 } ^ { p } f _ { R _ { h } } ( x _ { h, n | R _ { h } } + O ( \varepsilon ) ) + O ( \varepsilon ^ { 2 } ) \\ \cdot$$

/negationslash

$$& \quad ^ { n \ne i } \\ = & \, x _ { i, n + 1 | r _ { i } } + \varepsilon \left ( \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n | r _ { i } } ) + \sum _ { \substack { h = 1 \\ h \ne i } } ^ { s } x _ { h, n + 1 | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + 1 | R _ { h } } \right ) + O ( \varepsilon ^ { 2 } ) \,,$$

In the last step, we have neglected terms that contain ε ∂f r i ∂x ( x i,n r | i ) because ∂f r i ∂x ( x i,n r | i ) is arbitrarily small.

- 2. For i ∈ { s +1 , . . . , p } , we have

/negationslash

$$X _ { i } ( n + 1 ) = x _ { i, n + 1 | R _ { i } } + \varepsilon \left ( \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, n | R _ { i } } ) A _ { i, n } + \sum _ { h = 1 } ^ { s } x _ { h, n + 1 | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, n + 1 | R _ { h } } \right ) + O ( \varepsilon ^ { 2 } ) \,.$$

When using the induction hypothesis, we need to distinguish the case i ∈ { 1 , . . . , s } from the case i ∈ { s +1 , . . . , p } . For the CML (1), equations (9,10) yield the following equations.

- 1. For i ∈ { 1 , . . . , s } , we write the induction hypothesis for m ≥ 2 as

/negationslash

/negationslash

$$1. & \text{ For } i \in \{ 1, \dots, s \}, \text{ we write the induction hypothesis for } m \geq 2 \text{ as} \\ X _ { i } ( n + m ) = x _ { i, n + m | r _ { i } } + \left [ \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n + m - 1 | r _ { i } } ) + \sum _ { k = n } ^ { n + m - 2 } \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, k | r _ { i } } ) \prod _ { l = k + 1 } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \\ & + \sum _ { k = n + 1 } ^ { n + m - 1 } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) \prod _ { l = k } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \right ) \\ & + \sum _ { h = 1 } ^ { s } x _ { h, n + m | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + m | R _ { h } } \right ] \varepsilon + O ( \varepsilon ^ { 2 } ), \quad ( 1 1 )$$

which implies that

/negationslash

/negationslash

$$\text{which implies that} \\ X _ { i } ( n + m + 1 ) & = f _ { R _ { i } } \left [ x _ { i, n + m | r _ { i } } + \left ( \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n + m - 1 | r _ { i } } ) + \sum _ { k = n } ^ { n + m - 2 } \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, k | r _ { i } } ) \prod _ { l = k + 1 } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \\ & + \sum _ { k = n + 1 } ^ { n + m - 1 } \left ( \left ( \sum _ { h = \overset { s } { \in i } } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) \prod _ { l = k } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \right ) \\ & + \sum _ { h = \overset { s } { \in i } } x _ { h, n + m | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + m | R _ { h } } \right ) \varepsilon \right ] + \varepsilon \sum _ { h = \overset { s } { \in i } } f _ { r _ { h } } ( x _ { h, n + m | r _ { h } } + O ( \varepsilon ) ) \\ & + \varepsilon \sum _ { h = s + 1 } ^ { p } f _ { R _ { h } } ( x _ { h, n + m | R _ { h } } + O ( \varepsilon ) ) \\ \text{We Taylor expand all occurrences of $f$ and its derivatives to obtain}$$

/negationslash

We Taylor expand all occurrences of f and its derivatives to obtain

/negationslash

/negationslash

$$\text{We Taylor expand all occurrences of $f$ and its derivatives to obtain} \\ X _ { i } ( n + m + 1 ) = x _ { i, n + m + 1 | r _ { i } } + \left [ \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n + m | r _ { i } } ) + \sum _ { k = n } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, k | r _ { i } } ) \prod _ { l = k + 1 } ^ { n + m } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \\ + \sum _ { k = n + 1 } ^ { n + m } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { i } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { i } } \right ) \prod _ { l = k } ^ { n + m } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \right ) \\ + \sum _ { h = 1 } ^ { s } x _ { h, n + m + 1 | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + m + 1 | R _ { h } } \right ] \varepsilon + O ( \varepsilon ^ { 2 } ) \,. \\ 2. \text{ For } i \in \{ s + 1, \dots, p \}, \text{ we write the induction hypothesis for $m \geq 2$ as}$$

2. For i ∈ { s +1 , . . . , p } , we write the induction hypothesis for m ≥ 2 as

X n i ( + m ) = x i,n + m R | i

/negationslash

$$X _ { i } ( n + m ) & = x _ { i, n + m | R _ { i } } \\ + & \left [ \prod _ { k = n } ^ { n + m - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) A _ { i, n } + \sum _ { k = n + 1 } ^ { n + m - 1 } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) \prod _ { l = k } ^ { n + m - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, l | R _ { i } } ) \right ) \\ & \quad + \sum _ { h = 1 } ^ { s } x _ { h, n + m | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, n + m | R _ { h } } \right ] \varepsilon + O ( \varepsilon ^ { 2 } ) \,, & ( 1 3 )$$

/negationslash

which implies that

X n

$$\text{which implies that} \\ X _ { i } ( n + m + 1 ) = \\ f _ { R _ { i } } \left ( x _ { i, n + m | R _ { i } } + \left [ \prod _ { k = n } ^ { n + m - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) A _ { i, n } + \sum _ { k = n + 1 } ^ { n + m - 1 } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) \\ \prod _ { l = k } ^ { n + m - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, l | R _ { i } } ) \right ) + \sum _ { h = 1 } ^ { s } x _ { h, n + m | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + m | R _ { h } } \right ] \varepsilon + O ( \varepsilon ^ { 2 } ) \right ) \\ + \varepsilon \sum _ { h = 1 } ^ { s } f _ { r _ { h } } ( x _ { h, n + m | r _ { h } } + O ( \varepsilon ) ) + \sum _ { h = s + 1 } ^ { p } f _ { R _ { h } } ( x _ { h, n + m | R _ { h } } + O ( \varepsilon ) ) \\ \text{We Taylor expand of all occurrences of $f$ and its derivatives to obtain} \\ \nu / \cdot \quad. \quad. \quad.$$

/negationslash

/negationslash

We Taylor expand of all occurrences of f and its derivatives to obtain

X n i ( + m +1) = x i,n + m +1 | R i

/negationslash

$$X _ { i } ( n + m + 1 ) & = x _ { i, n + m + 1 | R _ { i } } \\ & + \left [ \prod _ { k = n } ^ { n + m } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) A _ { i, n } + \sum _ { k = n + 1 } ^ { n + m } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) \prod _ { l = k } ^ { m } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, l | R _ { i } } ) \right ) \\ & + \sum _ { h = 1 } ^ { s } x _ { h, n + m + 1 | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, n + m + 1 | R _ { h } } \right ] \varepsilon + O ( \varepsilon ^ { 2 } ) \\ & & \Box$$

/negationslash

## Theorem 1

Let | ε | &lt; 1 in the CML (1), and suppose that the hypotheses of Lemma 1 are satisfied. That is, we assume that the map f q i R i has an SN bifurcation at R i = r i , such that the associated SN orbit of f r i is a small-derivative SN orbit, that R i = r i + ε for i ∈ { 1 , . . . , s } , and that R i for i ∈ { s + 1 , . . . , p } are far away from r i . Let { x i, 1 | r i , x i, 2 | r i , . . . , x i,q i | r i } be a periodq i orbit for the uncoupled oscillator X i for i ∈ { 1 , . . . , s } , and let { x i, 1 | R i , x i, 2 | R i , . . . x i,q i | R i } be a periodq i orbit for the uncoupled oscillator X i for i ∈ { s +1 , . . . , p } . Consider the following initial conditions:

- · For i ∈ { 1 , . . . , s } ,

$$\., s \}, & \text{let} \\ & X _ { i } ( n ) = x _ { i, n | r _ { i } } + \varepsilon A _ { i, n } + O ( \varepsilon ^ { 2 } ) \,,$$

where x i,n r | i is the point of the SN orbit closest to the critical point of f R i at R i = r i , and A i,n is an arbitrary O (1) value.

/negationslash

- · For i ∈ { s +1 , . . . , p } , let

$$X _ { i } ( n ) = x _ { i, n | R _ { i } } + \varepsilon A _ { i, n } + O ( \varepsilon ^ { 2 } ) \,,$$

where x i,n R | i is a point of the orbit, and

/negationslash

/negationslash

$$& \text{where $x_{i,n|R_{i}}$ is a point of the orbit, and} \\ A _ { i, n } = \left [ \sum _ { k = n + 1 } ^ { n + q - 1 } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) \prod _ { l = k } ^ { n + q - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, l | R _ { i } } ) \right ) \\ & \quad + \sum _ { h = 1 } ^ { s } x _ { h, n + q | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, n + q | R _ { h } } \right ] \left ( \frac { 1 } { 1 - \prod _ { k = n } ^ { n + q - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) } \right ). \\ \\ & \text{The CML} \left ( 1 \right ) \text{has the solution}$$

The CML (1) has the solution

$$X _ { i } ( n + m ) = \begin{cases} x _ { i, n + m | r _ { i } } + \varepsilon A _ { i, n + m } + O ( \varepsilon ^ { 2 } ) \,, \quad i \in \{ 1, \dots, s \} \,, m \in \{ 1, \dots, q \} \\ x _ { i, n + m | R _ { i } } + \varepsilon A _ { i, n + m } + O ( \varepsilon ^ { 2 } ) \,, \quad i \in \{ s + 1, \dots, p \} \,, m \in \{ 1, \dots, q \} \,, \end{cases}$$

where the coefficients A i,n + m are periodic with period q = lcm( q 1 , q 2 , . . . , q p ) and satisfy the following formulas:

$$1. \text{ For } i \in \{ 1, \dots, s \},$$

/negationslash

$$A _ { i, n + 1 } = \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n | r _ { i } } ) + \sum _ { \substack { h = 1 \\ h \neq i } } ^ { s } x _ { h, n + 1 | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + 1 | R _ { h } } \,, \quad m = 1 \quad ( 1 7 )$$

(18)

/negationslash

$$A _ { i, n + m } = \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n + m - 1 | r _ { i } } ) + \sum _ { k = n } ^ { n + m - 2 } \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, k | r _ { i } } ) \prod _ { l = k + 1 } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \\ + \sum _ { k = n + 1 } ^ { n + m - 1 } \left ( \left ( \sum _ { h \neq i } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) \prod _ { l = k } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \right ) \\ + \sum _ { h = 1 } ^ { s } x _ { h, n + m | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + m | R _ { h } } \,, \quad m \in \{ 2, \dots, q \} \,.$$

/negationslash

2. For i ∈ { s +1 , . . . , p } ,

/negationslash

/negationslash

$$2. \text{ For } i \in \{ s + 1, \dots, p \}, \\ A _ { i, n + m } & = \left [ \sum _ { k = n + m + 1 } ^ { n + m + q - 1 } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, k | R _ { h } } \right ) ^ { n + m + q - 1 } \prod _ { l = k } ^ { 0 } \frac { 0 f R _ { i } } { \partial x } ( x _ { i, l | R _ { i } } ) \right ) \\ & + \sum _ { h = 1 } ^ { s } x _ { h, n + m + q | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, n + m + q | R _ { h } } \right ] \left ( \frac { 1 } { 1 - \prod _ { k = n + m } ^ { n + m + q - 1 } \frac { 0 f R _ { i } } { \partial x } ( x _ { i, k | R _ { i } } ) } \right ), \\ & \quad m \in \{ 1, \dots, q \}.$$

Remark. Although the initial conditions given in the statement of Theorem 1 may seem restrictive, our numerical computations demonstrate that - independently of the type of the orbit (i.e., either close to or far away from the SN) - it is sufficient to take as an initial condition any point of the unperturbed orbit plus a perturbation of size O ε ( ).

## Proof of Theorem 1.

We need to consider i ∈ { 1 , . . . , s } and i ∈ { s +1 , . . . , p } separately.

- 1. Using Lemma 1, it follows from X n i ( ) = x i,n r | i + εA i,n + O ε ( 2 ) that

/negationslash

$$X _ { i } ( n + 1 ) = x _ { i, n + 1 | r _ { i } } + \varepsilon \left ( \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n | r _ { i } } ) + \sum _ { h = 1 } ^ { s } x _ { h, n + 1 | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + 1 | R _ { h } } \right ) + O ( \varepsilon ^ { 2 } )$$

and

/negationslash

$$\text{and} \\ X _ { i } ( n + q + 1 ) = x _ { i, n + q + 1 | r _ { i } } + \left [ \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n + q | r _ { i } } ) + \sum _ { k = n } ^ { n + q - 1 } \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, k | r _ { i } } ) \prod _ { l = k + 1 } ^ { n + q } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \\ & + \sum _ { k = n + 1 } ^ { n + q } \left ( \left ( \sum _ { \hbar { = } i } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) \prod _ { l = k } ^ { n + q } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \right ) \\ & + \sum _ { \hbar { = } i } ^ { s } x _ { h, n + q + 1 | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + q + 1 | R _ { h } } \right ] \varepsilon + O ( \varepsilon ^ { 2 } ) \,.$$

/negationslash

Because n + q ∏ l = k ∂f r i ∂x ( x i,l r | i ) includes the arbitrarily small term ∣ ∣ ∣ ∣ ∂f r i ( x i,n r | i ) ∂x ∣ ∣ ∣ ∣ , it follows from (22) that

/negationslash

With q = lcm( q , q 1 2 , . . . , q p ), we have

$$X _ { i } ( n + q + 1 ) = x _ { i, n + q + 1 | r _ { i } } + \left ( \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n + q | r _ { i } } ) + \sum _ { \substack { h = 1 \\ h \neq i } } ^ { s } x _ { h, n + q + 1 | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + q + 1 | R _ { h } } \right ) \varepsilon + O ( \varepsilon ^ { 2 } ) \,. \\ \text{With } q = \text{lcm} ( q _ { 1 }, q _ { 2 }, \dots, q _ { p } ), \text{ we have}$$

$$\psi$$

$$x _ { i, n + q + 1 | r _ { i } } & = x _ { i, n + 1 | r _ { i } } \,, \\ x _ { i, n + q + 1 | R _ { i } } & = x _ { i, n + 1 | R _ { i } } \,,$$

/negationslash

$$\sum _ { \substack { h = 1 \\ h \neq i } } ^ { s } x _ { h, n + q + 1 | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, n + q + 1 | R _ { h } } = \sum _ { \substack { h = 1 \\ h \neq i } } ^ { s } x _ { h, n + 1 | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, n + 1 | R _ { h } } \,,$$

/negationslash

/negationslash

/negationslash because x i,j is a point of a q i -period orbit. Consequently, equations (21) and (22) become the same equation. From equations (9) and (21), we can write equation (6) in Lemma 1 as X n i ( + m ) = x i,n + m r | i + εA i,n + m + O ε ( 2 ) to obtain

/negationslash

/negationslash

$$\text{to obtain} \quad & \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \ \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \, \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd. \quad \end{array} \quad \text{to obtain} \quad & \quad \partial f _ { r _ { i } } ( x _ { i, n + m - 1 } ) + \sum _ { k = n } ^ { p } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \prod _ { l = k + 1 } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \right ) \\ & + \sum _ { \substack { h = 1 \\ h = 1 \\ h = 1 } } ^ { s } x _ { h, n + m | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + m | R _ { h } } \,, \quad \quad m \in \{ 1, \dots, q \} \,. \quad \quad \quad ( 2 3 ) \\ 2. \ \text{With} \ X _ { i } ( n + m ) = x _ { i, n + m | R _ { i } } + \varepsilon A _ { i, n + m } + O ( \varepsilon ^ { 2 } ), \ \text{Lemma} \ 1 \text{ implies that}$$

2. With X n i ( + m ) = x i,n + m R | i + εA i,n + m + O ε ( 2 ), Lemma 1 implies that

$$X _ { i } ( n + m + q ) = x _ { i, n + m + q | R _ { i } } + \prod _ { k = n + m } ^ { n + m + q - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) A _ { i, n + m } \,.$$

/negationslash

/negationslash

$$X _ { i } ( n + m + q ) & = x _ { i, n + m + q | R _ { i } } + \left [ \prod _ { k = n + m } ^ { n + m + q - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) A _ { i, n + m } \\ & \quad + \sum _ { k = n + m + 1 } ^ { n + m + q - 1 } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) ^ { n + m + q - 1 } \prod _ { l = k } ^ { \partial f _ { R _ { i } } } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, l | R _ { i } } ) \right ) \\ & \quad + \sum _ { h = 1 } ^ { s } x _ { h, n + m + q | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + m + q | R _ { h } } \right ] \varepsilon + O ( \varepsilon ^ { 2 } ) \,. \\ & 1 1$$

By taking q = lcm( q , q 1 2 , . . . , q p ), we obtain x i,n + m r | i = x i,n + m q r + | i and x i,n + m R | i = x i,n + m q R + | i because x i,j is a point of a periodic orbit. Consequently, X n i ( + m ) -X n i ( + m + ) = q O ε ( ) whenever

/negationslash

/negationslash

$$\text{query}, \, \Lambda _ { i } ( n + m ) - \Lambda _ { i } ( n + m + q ) = \mathcal { U } _ { \infty } ( \varepsilon ) \text{ whenever} \\ A _ { i, n + m } & = \prod _ { k = n + m } ^ { n + m + q - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) A _ { i, n + m } + \sum _ { k = n + m + 1 } ^ { n + m + q - 1 } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, k | R _ { h } } \right ) \\ & \prod _ { l = k } ^ { n + m + q - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, l | R _ { i } } ) \right ) + \sum _ { h = 1 } ^ { s } x _ { h, n + m + q | r _ { h } } + \sum _ { \substack { h = s + 1 \\ h \neq i } } ^ { p } x _ { h, n + m + q | R _ { h } } \cdot \, ( 2 4 ) \\ \text{Furthermore} \quad A \colon \dots \text{ is neriodic}$$

Furthermore, A i,n + m is periodic. Equation (24) now implies that

/negationslash

/negationslash

$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{	\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{ $\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{#\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{  $\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text}$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text#\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$
$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\test$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text'$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text<$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text.$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text`$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$$

It follows that A i,n + m has period q because it is given by sums and products of q -periodic functions evaluated at points of a q -periodic orbit.

Observe that the formula for A i,j for i ∈ { 1 , . . . , s } in equation (19) does not contain the term   1 -j + q -1 ∏ k = j ∂f r i ∂x ( x i,k r | i )   in the denominator [see equation

(3)]. Otherwise, A i,j would be of size O (1 /ε ), and the expansion that we used to prove Theorem 1 would not be valid. By contrast, the formula for A i,j for

$$& \text{to prove $1$ theorem $1$ would not be valid. By contrast, the formula for $A_{i,j}$ for} \\ & i \in \{ s + 1, \dots, p \} \text{ in equation (20) includes the term } \left [ 1 - \prod _ { k = j } ^ { j + q - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) \right ] \\ & \text{in the denominator because the oscillators are far from SN bifurcations for} \\ & i \in \{ s + 1, \dots, p \}. \text{ Therefore.}$$

in the denominator because the oscillators are far from SN bifurcations for i ∈ { s +1 , . . . , p } . Therefore,

$$1 - \prod _ { k = j } ^ { j + q - 1 } \frac { \partial f _ { R _ { i } } } { \partial x } ( x _ { i, k | R _ { i } } ) = O ( 1 ) \,, \\. \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cap \, \cdot$$

and it follows that A i,j also has size O (1).

,

## 2.2. Type-I Intermittency Near Saddle-Node Bifurcations

Theorem 1 concerns the behavior of the CML (1) with a mixture of periodic oscillators that are near the SN bifurcation with others that are far from the SN bifurcation. If an SN orbit takes place at R i = r i , then the oscillators with R i = r i + ε are the ones that are close to the SN orbit.

We now want to study the behavior of the CML (1) when an uncoupled oscillator has type-I intermittency [26] at R i = r i -ε (i.e., just to the left of where it undergoes an SN bifurcation). Type-I intermittency is characterized by the alternation of an apparently periodic regime (a so-called 'laminar phase'), whose mean duration follows the power law 〈 〉 l ∝ ε -1 2 (so the laminar region becomes longer as ε becomes smaller), and chaotic bursts. As R i = r i -ε , we expand f r i -ε in powers of ε to obtain

$$f _ { r _ { i } - \varepsilon } ^ { q _ { i } } ( x _ { j | r _ { i } } ) = f _ { r _ { i } } ^ { q _ { i } } ( x _ { j | r _ { i } } ) - \varepsilon \frac { \partial f _ { r _ { i } } ^ { q _ { i } } } { \partial r } ( x _ { j | r _ { i } } ) + O ( \varepsilon ^ { 2 } ) \,,$$

where x j a point of a periodq i SN orbit. Therefore the laminar phase is driven by the periodq i SN orbit associated with the SN bifurcation. Thus, as ε becomes smaller, the orbit spends more iterations in the laminar regime, and it thus more closely resembles the periodq i SN orbit. In particular, ∣ ∣ ∣ x j r | i -f q i r i -ε ( x j r | i ) ∣ ∣ ∣ = O ε ( ) .

To approximate the temporal evolution of the laminar regime using the periodq i SN orbit, we proceed in the same way as in Theorem 1, except that we replace R i = r i + ε by R i = r i -ε . We thus write

/negationslash

$$X _ { i } ( n + 1 ) = x _ { i, n + 1 | r _ { i } } + \left [ - \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n | r _ { i } } ) + \sum _ { \substack { h = 1 \\ h \neq i } } ^ { s } x _ { h, n + 1 | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + 1 | R _ { h } } \right ] \varepsilon + O ( \varepsilon ^ { 2 } ) \,, \quad m = 1 \,,$$

/negationslash

$$X _ { i } ( n + m ) = x _ { i, n + m | r _ { i } } + \left [ - \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n + m - 1 | r _ { i } } ) - \sum _ { k = n } ^ { n + m - 2 } \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, k | r _ { i } } ) \prod _ { l = k + 1 } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \\ + \sum _ { k = n + 1 } ^ { n + m - 1 } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) \prod _ { l = k } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \right ) \\ \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat{ \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat$$

$$( n + m ) = x _ { i, n + m | r _ { i } } + \left [ - \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, n + m - 1 | r _ { i } } ) - \sum _ { k = n } ^ { n + m - 2 } \frac { \partial f _ { r _ { i } } } { \partial r } ( x _ { i, k | r _ { i } } ) \prod _ { l = k + 1 } ^ { n + m - 1 } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \\ + \sum _ { k = n + 1 } ^ { n + m - 1 } \left ( \left ( \sum _ { h = 1 } ^ { s } x _ { h, k | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, k | R _ { h } } \right ) ^ { n + m - 1 } \prod _ { l = k } ^ { \partial f _ { r _ { i } } } \frac { \partial f _ { r _ { i } } } { \partial x } ( x _ { i, l | r _ { i } } ) \right ) \\ + \sum _ { h = 1 } ^ { s } x _ { h, n + m | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + m | R _ { h } } \right ] \varepsilon + O ( \varepsilon ^ { 2 } ) \,, \quad m \in \{ 2, \dots, q \} \\ \dots \,.$$

/negationslash

$$+ \sum _ { \substack { h = 1 \\ h \neq i } } ^ { s } x _ { h, n + m | r _ { h } } + \sum _ { h = s + 1 } ^ { p } x _ { h, n + m | R _ { h } } \Big { ] } \, \varepsilon + O ( \varepsilon ^ { 2 } ) \,, \quad m \in \{ 2, \dots, q \}$$

which determines the temporal evolution of the CML (1) in the laminar regime.

## 3. Numerical Computations

Theorem 1 proves the existence of an approximately periodic orbit. In principle, one can deduce the existence of a periodic orbit by using the Implicit Function Theorem (IFT). However, the IFT fails at the SN bifurcation (i.e., at R i = r i ) for free oscillators and consequently fails near an SN bifurcation (i.e., for R i = r i + ) of the HWCML (1), because the Jacobian determinant vanishes. ε

Had we expanded all terms in Theorem 1, we would have obtained terms of size O ε ( 2 ) that depend on the coefficients of the terms of size O ε ( ) (i.e., as functions of the A i,n + m terms in Theorem 1), so terms of size O ε ( 2 ) would have the same period as the A i,n + m . We could then obtain terms of size O ε ( 3 ) as a functions of the coefficients of lower-order terms. These terms would also have the same period as A i,n + m , and the same is true for all higher-order terms if we continued the expanding in powers of ε . This reasoning suggests the existence of a periodic orbit of period q = lcm( q , . . . , q 1 p ) (not just an approximate one), and our numerical simulations successfully illustrate the existence of such periodic orbits.

For simplicity, we consider a pair of coupled oscillators,

$$X ( n + 1 ) & = f ( X ( n ) ) + \varepsilon g ( Y ( n ) ) \,, \\ Y ( n + 1 ) & = g ( Y ( n ) ) + \varepsilon f ( X ( n ) ) \,,$$

where f ( x ) = R x 1 (1 -x ) and g y ( ) = cos( R y 2 ). We initially fix the coupling to be ε = 0 0001, though we will later consider 2 . ε , 3 , and so on. ε The uncoupled oscillator Y ( n ) has a fixed period of 4 and is far away from a SN bifurcation for R 2 = 1 9. . We use values of R 1 such that the uncoupled oscillator X n ( ) is near an SN bifurcation, and we consider SN orbits with different periods.

## 3.1. Uncoupled Oscillator X n ( ) with a Period-3 Orbit

For the oscillator X n ( ), we fix R 1 = r 1 +2 , where ε r 1 ≈ 3 828427 is an SN . bifurcation point of f . When there is no coupling, the free oscillator X n ( ) has a period-3 SN orbit, and the free oscillator Y ( n ) has a period-4 orbit. When coupled, both X n ( ) and Y ( n ) have a periodic orbit with period q = lcm(3 4) = , 12 (see Fig. 2).

At R 1 = r 1 + , the HWCML (29) exhibits type-I intermittency associated ε with the SN bifurcation (see Fig. 3). However, for larger R 1 (e.g., r 1 + 2 , ε r 1 +3 , ε . . . , r 1 +7 ), the periods of the uncoupled oscillators ε X n ( ) and Y ( n ) are preserved because we are farther away from the bifurcation point. We observe type-I intermittency for R 1 = r 1 +0 , ε R 1 = r 1 -ε , R 1 = r 1 -2 . ε

Remark. When R 1 = r 1 +2 , we calculate 1 ε -∏ j + m -1 k = j ∂f r i ∂x ( x i,k r | i ) ≈ 0 24 for . ε = 0 0001. (For . R 1 = r 1 + , we obtain a smaller value for the second quantity). ε Recall the quantifications of 'far from' and 'near' in section 2.1. Although ε can be very small, the periodic windows that are born with an SN orbit can be even smaller than ε . Thus, from the dynamical standpoint, a very small

value of the coupling parameter ε can nevertheless be large as a variation on a bifurcation parameter.

In Section 2.2, we determined the temporal evolution of the oscillators in the laminar regime of type-I intermittency up to size O ε ( ). By comparing Fig. 2 (which depicts the dynamics for a parameter value slightly larger than the SN bifurcation point) and Fig. 4 (which depicts the dynamics right before the bifurcation), we observe periodic behavior just after the bifurcation and laminar behavior just before it.

## 3.2. Uncoupled Oscillator X n ( ) with a Period-5 Orbit

We proceed as in Section 3.1 and obtain similar results.

For the oscillator X n ( ), we fix R 1 = r 1 +2 , where ε r 1 ≈ 3 738173 is an SN . bifurcation point of f . When there is no coupling, the free oscillator X n ( ) has a period-5 SN orbit, and the free oscillator Y ( n ) has a period-4 orbit. When coupled, both X n ( ) and Y ( n ) have a periodic orbit with period q = lcm(5 4) = , 20 (see Fig. 5).

At R 1 = r 1 + , the HWCML (29) exhibits type-I intermittency associated ε with the SN bifurcation (see Fig. 6). However, for larger R 1 (e.g., r 1 + 2 , ε r 1 +3 , ε r 1 +4 ε . . . ), the periods of the uncoupled oscillators X n ( ) and Y ( n ) are preserved because we are farther away from the bifurcation point. We observe type-I intermittency for R 1 = r 1 +0 , ε R 1 = r 1 -ε , R 1 = r 1 -2 . ε

## 3.3. Summary of HWCML Dynamics

Our results allow us to deduce the dynamics of the HWCML (29) when R i = r i + ε α . We worked with a coupling strength of ε = 0 0001 and a control . parameter of R i = r i + kε . In our numerical computations, we observed the following behavior:

- (a) intermittency for R i ≤ r i + ; ε
- (b) periodic behavior for R i ≥ r i +2 . ε

Therefore, the following occurs.

- (i) If we choose R i = r i + ε α with α &gt; 1, then R &lt; r i i + , and the HWCML ε exhibits intermittent behavior according to (a).
- (ii) If we choose R i = r i + ε α with 0 &lt; α &lt; 1, then R &gt; r i i +2 ; this holds ε even for α close to 1, as long as ε α &gt; 2 ε (e.g., 0 &lt; α /lessorapproxeql 0 92 for . ε = 0 0001). . Therefore, the HWCML exhibits periodic behavior according to (b).

Based on our numerical computations, we can thus establish the following statement: 'Under the hypotheses of Theorem 1, the oscillators in the CML (1) have periodic orbits that persist with the same period as in Theorem 1 for perturbations of size O ε ( ). That is, higher-order terms do not change the period, as we heuristically stated at the beginning of Section 3.'

We now discuss the consequences of all oscillators in an HWCML having the same period q = lcm( q , q 1 2 , . . . , q p ), where q , . . . , q 1 p are the periods of the free oscillators. One can adjust the parameters to obtain periods q 1 . . . q p so

/negationslash that q = lcm( q 1 , q 2 , . . . , q p ) remains constant. For example, if q 1 = 3 and q 2 = 2 , then k q = lcm(3 2 ) = 3 , k × 2 k (for integers k &gt; 0). If the first oscillator undergoes a period-doubling cascade, then its period is 3, 3 × 2, 3 × 2 , and so 2 on. However, the period m of the HWCMLs is q = lcm(3 2 ) = lcm(3 , k × 2 2 ) = , k · · · = lcm(3 × 2 k , 2 ) = 3 k × 2 , k so it does not change even after an arbitrary number of period-doubling bifurcations. That is, for arbitrarily small ε = 0, the HWCML period remains the same even amidst a period-doubling cascade.

We illustrate the above phenomenon with a simple computation. Consider the HWCML (29) and suppose that R 1 = 3 83 and . R 2 = 1 9. When . ε = 0 (i.e., when there is no coupling), the free oscillator X n ( ) has a period-3 orbit and the free oscillator Y ( n ) has a period-4 orbit. However, when ε = 0 001, both . X n ( ) and Y ( n ) have a periodic orbit with period q = lcm(3 4) = 12. As we show in , Table 1, the free oscillator X n ( ) undergoes period-doubling bifurcations, but the HWCML exhibits synchronization and still has period-12 orbits for ε = 0 001. .

## 4. Conclusions and Discussion

We have examined heterogeneous weakly coupled map lattices (HWCMLs) and have given results to describe periodic orbits both near and far from saddlenode orbits and to describe the temporal evolution of the laminar regime in typeI intermittency. All periodic windows of the bifurcation diagram of unimodal maps originate from SN bifurcations, so it is important to explore the dynamics near such bifurcation points.

An important implication of our results is that HWCMLs of oscillators need not behave approximately like their associated free-oscillator counterparts. In particular, they can have periodic-orbit solutions with completely different periods even for arbitrarily small coupling strengths ε = 0.

/negationslash

Our numerical calculations illustrate an important result about period preservation when oscillator parameters change. Even when one varies the parameters R i of the functions f R i such that the uncoupled oscillator X i undergoes a perioddoubling cascade, the periods of each of the coupled oscillators are preserved as long as the least common multiple of the periods remains constant. That is, the oscillation period is resilient to changes.

Period preservation is a rather generic phenomenon in CMLs. Suppose, for example, that one oscillator has period of q × 2 n , which can originate either from period doubling or from an SN bifurcation [27]. One can then change parameters so that different individual oscillators (if uncoupled) would undergo a period-doubling cascade, whereas the least common multiple of the periods of those oscillators will remain constant until one oscillator (if uncoupled) has period p × 2 n +1 . In a CML, a very large number of oscillators can each undergo a period-doubling cascade, so the period of a CML can be very resilient even in situations when other conditions - in particular, the values of the parameters in the CML - are changing a lot. Moreover, one can adjust the parameters to obtain oscillations of arbitrary periods q 1 . . . q p with q = lcm( q 1 , q 2 , . . . , q p ) = constant. Consequently, period preservation is a very common phenomenon: it

is not limited to the aforementioned period-doubling cascade but rather appears throughout a bifurcation diagram.

Periodic orbits anticipated by Theorem 1 and confirmed in Section 3 correspond to traveling waves in a one-dimensional HWCML and to periodic patterns in a multidimensional HWCML. Such patterns have been studied in homogeneous CMLs [13, 29], and our results can help to describe such dynamics in heterogeneous CMLs both near and far from bifurcations. Our observation about period resilience implies that there will be many different patterns with the same period. Small changes in an HWCML can change the specific pattern, but the period itself is rather robust.

Our results also have implications in applications. A toy macroscopic traffic flow model, governed by the logistic map, was proposed in [30]. The derivation of the model is based on very general assumptions involving speed and density. When these assumptions are satisfied, one can use the model to help examine the evolution of flows of pedestrians, flows in a factory, and so on. When such flows interact weakly, then equations of the form that we discussed in Section 2.1 can be useful for such applications. For example, one could do a simple examination of the temporal evolution of two groups of football fans around a stadium (or of sheep around an obstacle [31]). The two groups have different properties, so suppose that they are governed by an HWCML. From our results, if each group is regularly entering the stadium on its own (i.e., their behavior is periodic), then both groups considered together would continue to enter regularly at the same rate, provided that the interaction between the two groups is weak. This suggests that it would be interesting to explore a security strategy that models erecting a light fence to ensure that the interaction between the two groups remains weak.

The model in Ref. [30] also admits chaotic traffic patterns. One can construe the intermittent traffic flow in a traffic jam as being formed by regular motions (i.e., a laminar regime) and a series of acceleration and braking (i.e., chaotic bursts). Our results give the temporal evolution of such a laminar regime in a chaotic intermittent flow if the interaction between entities is weak (i.e., when the laminar regime is long, as we discussed in Section 2.2). Indeed, as has been demonstrated experimentally for the flow of sheep around an obstacle [31], it is possible to preserve laminar behavior for a longer time through the addition of an obstacle.

## Acknowledgements

We are grateful to the anonymous referee for his/her enlightening and detailed suggestions. We also thank Daniel Rodr´ ıguez P´rez for his help in the e preparation of this manuscript.

## References

## References

- [1] Farkas I, Helbing D, Vicsek T. Human waves in stadiums. Physica A 2003; 330:18-24
- [2] Sumpter DJT. Collective Animal Behavior. Princeton: Princeton University Press; 2010.
- [3] Boccara N. Modelling Complex Systems. Second Edition. Berlin: SpringerVerlag; 2010.
- [4] Coca D, Billings SA. Analysis and reconstruction of stochastic coupled map lattice models. Phys. Lett. A 2003; 315:61-75
- [5] Pavlov EA, Osipov GV, Chan CK, Suykens JAK. Map-based model of the cardiac action potential. Phys. Lett. A 2011; 375:2894-2902.
- [6] Kaneko K. Theory and Applications of Coupled Map Lattices. New York: John Wiley&amp; Sons; 1993.
- [7] Kaneko K. From globally coupled maps to complex-systems biology. Chaos 2015; 25:097608
- [8] Special issue. Chaos 1992; 2(3).
- [9] Special issue. Physica D 1997; 103.
- [10] Tang Y, Wang Z, Fang J. Image encryption using chaotic coupled map lattices with time-varying delays. Commun. Nonlinear Sci. Numer. Simulat. 2010; 15:2456-2468.
- [11] Wang S, Hu G, Zhou H. A one-way coupled chaotic map lattice based selfsynchronizing stream cipher. Commun. Nonlinear Sci. Numer. Simulat. 2014; 19:905-913.
- [12] Kaneko K. Spatiotemporal chaos in one- and two-dimensional coupled map lattices. Physica D 1989; 37:60-82.
- [13] He G. Travelling waves in one-dimensional coupled map lattices. Commun. Nonlinear Sci. Numer. Simulat. 1996; 1:16-20.
- [14] Franceschini V, Giberti C, Vernia C. On quasiperiodic travelling waves in coupled map lattices. Physica D 2002; 164:28-44.
- [15] Cherati ZR, Motlagh MRJ. Control of spatiotemporal chaos in coupled map lattice by discrete-time variable structure control. Phys. Lett. A 2007; 370:302-305.
- [16] Jakobsen A. Symmetry breaking bifurcation in a circular chain of N coupled logistic maps. Physica D 2008; 237:3382-3390.

| [17]   | Sotelo Herrera D, San Mart´n ı J. Analytical solutions of weakly coupled map lattices using recurrence relations. Phys. Lett. A 2009; 373:2704-2709.                                                             |
|--------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [18]   | Xu L, Zhang G, Han B, Zhang L, Li MF, Han YT. Turing instability for a two-dimensional Logistic coupled map lattice. Phys. Lett. A 2010; 374:3447-3450.                                                          |
| [19]   | Angelini L, Pellicoro M, Stramaglia S. Phase ordering in chaotic map lattices with additive noise. Phys. Lett. A 2001; 285:293-300.                                                                              |
| [20]   | Lai YM, Porter MA. Noise-induced synchronization, desynchronization, and clustering in globally coupled nonidentical oscillators. Phys. Rev. E 2013; 88:012905                                                   |
| [21]   | Moehlis J, Josi´ c K, Shea-Brown ET. Periodic orbit. Scholarpedia 2006; 1(7):1358                                                                                                                                |
| [22]   | Cvitanovi´ c P, Artuso R, Mainieri R, Tanner G, Vattay G, Whelan N, Wirzba A. Chaos: Classical and Quantum. Version 14. (available at http://chaosbook.org ) (2012).                                             |
| [23]   | Sotelo Herrera D, San Mart´n ı J. Travelling waves associated with saddle- node bifurcation in weakly coupled CML. Phys. Lett. A 2010; 374:3292- 3296.                                                           |
| [24]   | Sotelo Herrera D, San Mart´n ı J, Cerrada L. Saddle-node bifurcation cas- cades and associated travelling waves in weakly coupling CML. Int. J. Bif. Chaos. 2012; 22:1250172.                                    |
| [25]   | Guckenheimer J, Holmes P. Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields Berlin: Springer-Verlag; 1983.                                                                            |
| [26]   | Pomeau Y, Manneville P. Intermittent Transition to Turbulence in Dissi- pative Dynamical Systems. Commun. Math. Phys. 1980; 74:189-197.                                                                          |
| [27]   | San Mart´n ı J. Intermittency cascade. Chaos Solitons & Fractals 2007; 32:816-831.                                                                                                                               |
| [28]   | Milnor J. Self-similarity and hairiness in the Mandelbrot set. Lect. Notes Pure App. Math. 1989; 114:211-257.                                                                                                    |
| [29]   | dos S. Silva FA, Viana RL, de L. Prado T, Lopes SR. Characterization of spatial patterns produced by a Turing instability in coupled dynamical systems. Commun. Nonlinear Sci. Numer. Simul. 2014; 19:1055-1071. |
| [30]   | Lo S-C, Cho H-J. Chaos and control of discrete dynamic traffic model. J. Franklin Inst. 2005; 342:839-851.                                                                                                       |
| [31]   | Zurigel I et al. Clogging transition of many-particle systems flowing through bottlenecks. Sci. Rep. 2014; 4:7324.                                                                                               |

## Figure Legends

Figure 1: Left: The maps f 3 and f 4 (where f is the logistic map) and the fixed points at which SN bifurcations occur . Observe that there are SN points far away from the critical point C. Because 4 &gt; 3, the extremum of f 4 near the critical point C is narrower than the extremum of f 3 near C. Right: Magnification of the extrema near the critical point C . We show the distances d between the SN point and the critical point for both f 3 and f 4 . Observe that the distance between this pair of points decreases as the period increases.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0303v4-figure-3451aff7-9fa2-4615-a4d5-2b180d831fb7.jpg)

Figure 2: Temporal evolution of the HWCML (29) for R 1 = r 1 + 2 ε , where r 1 ≈ 3 828427 . (which is an SN bifurcation point) and R 2 = 1 9. . The uncoupled oscillators have (a) period 3 and (b) period 4. When ε = 0 0001 (i.e., weak coupling), . the oscillators X ε ( n ) and Y ε ( n ) both have period lcm(3 , 4) = 12. In panels (c) and (d), we plot X ε ( n ) -X 0 ( n ) and Y ε ( n ) -Y 0 ( n ) (i.e., the solution in the coupled case minus the solution in the ε = 0 case) to better observe the period-12 dynamics.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0303v4-figure-bccfe865-ab3c-4712-aba0-a78994a50ffa.jpg)

Figure 3: Temporal evolution of the HWCML (29) for R 1 = r 1 + ε , where r 1 ≈ 3 828427 . (which is an SN bifurcation point) and R 2 = 1 9 . . The uncoupled oscillators have (a) period 3 and (b) period 4. When ε = 0 0001 (i.e., weak coupling), . the oscillators X ε ( n ) and Y ε ( n ) exhibit type-I intermittency. In panels (c) and (d), we plot X ε ( n ) -X 0 ( n ) and Y ε ( n ) -Y 0 ( n ) (i.e., the solution in the coupled case minus the solution in the ε = 0 case) to better observe the intermittency dynamics.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0303v4-figure-8beb1bae-2960-4d92-a9af-9c8859b691ac.jpg)

Figure 4: Magnification of the laminar regime of type-I intermittency from Fig. 3d . We can clearly see the resemblance with the temporal evolution of the oscillator in the periodic regime (see Fig. 2d).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0303v4-figure-046340fb-3267-46b0-a81f-03dcfeba3766.jpg)

Figure 5: Temporal evolution of the HWCML (29) for R 1 = r 1 + 2 ε , where r 1 ≈ 3 738173 . (which is an SN bifurcation point) and R 2 = 1 9 . . The uncoupled oscillators have (a) period 5 and (b) period 4. When ε = 0 0001 (i.e., weak coupling), . the oscillators X ε ( n ) and Y ε ( n ) both have period lcm(5 , 4) = 20. In panels (c) and (d), we plot X ε ( n ) -X 0 ( n ) and Y ε ( n ) -Y 0 ( n ) (i.e., the solution in the coupled case minus the solution in the ε = 0 case) to better observe the period-20 dynamics.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0303v4-figure-2c67c631-c292-4356-96ac-1f85dd172849.jpg)

Figure 6: Temporal evolution of the HWCML (29) for R 1 = r 1 + ε , where r 1 is an SN bifurcation point and R 2 = 1 9. . The uncoupled oscillators have (a) period 5 and (b) period 4. When ε = 0 0001 (i.e., . weak coupling), the oscillators X ε ( n ) and Y ε ( n ) show intermittency. In panels (c) and (d), we plot X ε ( n ) -X 0 ( n ) and Y ε ( n ) -Y 0 ( n ) (i.e., the solution in the coupled case minus the solution in the ε = 0 case) to better observe the intermittent dynamics.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0303v4-figure-1f4027b9-cf14-4d0f-9eae-cd969cf0d9e1.jpg)

## Tables

Table 1: Period of the CML (29) for r 2 = 1 9 . and ε = 0 001. The parameter . r 1 indicates when the logistic map, which is satisfied by the free oscillator X n ( ), exhibits orbits of various periods during a period-doubling cascade in the window of period-3 orbits in the bifurcation diagram. Although the period of X n ( ) changes, the period of the CML remains the same.

| r 1        | Period of X ( n )   | Period of the CML     |
|------------|---------------------|-----------------------|
| 3 . 831874 | 3                   | lcm(3 , 4) = 12       |
| 3 . 844568 | 3 × 2               | lcm(3 × 2 , 4) = 12   |
| 3 . 848344 | 3 × 2 2             | lcm(3 × 2 2 , 4) = 12 |
Dolores Sotelo Herrera Jesús San Martín Mason A. Porter
Heterogeneous, Weakly Coupled Map Lattices
Coupled map lattices (CMLs) are often used to study emergent phenomena in
nature. It is typically assumed (unrealistically) that each component is
described by the same map, and it is important to relax this assumption. In
this paper, we characterize periodic orbits and the laminar regime of type-I
intermittency in heterogeneous weakly coupled map lattices (HWCMLs). We show
that the period of a cycle in an HWCML is preserved for arbitrarily small
coupling strengths even when an associated uncoupled oscillator would
experience a period-doubling cascade. Our results characterize periodic orbits
both near and far from saddle--node bifurcations, and we thereby provide a key
step for examining the bifurcation structure of heterogeneous CMLs.