arxiv: v1 [math.ds] 17 Oct 2013

Transcription

1 DISSIPATIVE OUTER BILLIARDS: A CASE STUDY GIANLUIGI DEL MAGNO, JOSE PEDRO GAIVÃO, AND EUGENE GUTKIN This paper is dedicated to the memory of the third author, who unexpectedly passed away while the paper was being completed. arxiv: v1 [math.ds] 17 Oct 2013 Abstract. We study dissipative polygonal outer billiards, i.e. outer billiards about convex polygons with a contractive reflection law. We prove that dissipative outer billiards about any triangle and the square are asymptotically periodic, i.e. they have finitely many global attracting periodic orbits. A complete description of the bifurcations of the periodic orbits as the contraction rates vary is given. For the square billiard, we also show that the asymptotic periodic behavior is robust under small perturbations of the vertices and the contraction rates. Finally, we describe some numerical experiments suggesting that dissipative outer billiards about regular polygon are generically asymptotically periodic. 1. Introduction and motivation The outer billiard about a convex compact region Ω C is the map T : C \ Ω C \ Ω introduced by M. Day [7]. For z X := C \ Ω, let L z C be the ray through z intersecting tangentially Ω on the right with respect to the obvious orientation. Suppose first that Ω is strictly convex, and let t Ω be the intersection point of Ω and L z. Then T (z) L z is the point symmetric to z with respect to t. The transformation T : X X is an area preserving twist map [20]. Neumann asked if there are points z X whose ω-limit sets contain subsets of Ω or infinity [26]. Applying his famous KAM-type theorem, Moser showed that this is impossible if Ω is continuously differentiable 333 times and is strictly convex [23, 24]. In this case, in fact, X contains invariant circles arbitrarily close to Ω and arbitrarily far away from it. It follows that the orbit of every z X belongs to an annulus bounded by two invariant circles. Later on, R. Douady showed that Date: October 18, Mathematics Subject Classification. Primary: 37E15; Secondary: 37E99, 37D50. Key words and phrases. Outer billiards, Contracting Reflection Law, Piecewise contractions. The first two authors were partially supported by Fundacão para a Ciência e a Tecnologia through the Program POCI 2010 and the Project Randomness in Deterministic Dynamical Systems and Applications (PTDC- MAT ). 1

2 2 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN six continuous derivatives of Ω suffice for the existence of invariant circles [10]. However, the strict convexity of Ω is necessary [14]. The definition of the outer billiard map has to be modified if Ω is not strictly convex. Suppose for example that Ω is a convex k-gon P. For points z X such that the ray L z contains a side of P, the image T (z) is not uniquely defined. These points form the singular set S of T, which is a union of k half-lines. The outer billiard map is not defined on S, and has jump discontinuities across its lines. The complement X \ S is a disjoint union of k open cones A 1,..., A k. The apex of A k is a corner point w k P, and the restriction T Ak is the euclidean reflection about w k. Thus, the outer billiard about P is a planar piecewise isometry on k cones. The dynamics of some classes of piecewise isometries 1 have been studied in detail, mostly in two dimensions [2, 5, 12, 13, 17]. Polygonal outer billiards make a very special and beautiful class of piecewise isometries. Moser 2 pointed out that for these billiards, Neumann s question about the existence of orbits accumulating on the boundary of the table or at infinity was open [25]. In view of Moser s remark, the question has attracted much of attention, and several partial results in both directions were obtained [15, 21, 27, 32]. However, the full scope of the Neumann-Moser question remains open. For various aspects of polygonal outer billiards, we refer the reader to [3, 11, 16, 18, 29, 28]. We will now embed the outer billiard about the convex k-gon P into a one-parameter family of dynamical systems. As before, let A i be the open cones of C whose apex is the corner w i of P. Let λ (0, + ). Define T λ (z) = λz + (1 + λ)w i for z A i. (1.1) The map T λ is a piecewise affine dilation on the cones A 1,..., A k. For λ = 1, we recover the outer billiard map. For 0 < λ < 1 (resp. 1 < λ), the map is piecewise contracting (resp. piecewise expanding). If we denote by Q the polygon obtained by reflecting P about the origin, then it is easy to see that T 1 λ is conjugated to the map T 1/λ of the outer billiard about Q with contraction rate 1/λ. In particular, the ω- limit sets of T λ and T 1/λ coincide. Therefore, as far as the asymptotic properties of T λ for λ 1 are concerned, we will only consider the case 0 < λ < 1. We will call the map T λ with 0 < λ < 1 the dissipative outer billiard, and will call the outer billiard the map T λ with λ = 1 (see Fig. 1). The dynamics of the outer billiard when λ = 1 about the equilateral triangle, the square and the hexagon is particularly simple: every orbit 1 Interval exchanges transformation are one-dimensional piecewise isometries. This very special class of piecewise isometries makes a research subject on its own right, and it would take us too far afield to comment on it. 2 Moser presented this question as a caricature of the problem on the stability of motion in celestial mechanics [25].

3 OUTER BILLIARDS 3 T λ (z) z T 2 λ(z) T 3 λ(z) Figure 1. Example of a dissipative polygonal outer billiard. of these billiards is periodic [25]. In this paper, we investigate in detail the dissipative outer billiards about the equilateral triangle and the square. We show that both billiards are asymptotically periodic, i.e., their ω-limit sets consist of periodic orbits, and every orbit converges to one of them. We then perform a bifurcation analysis, by computing the number and the period of these orbits as λ varies between 0 and 1. Using a result by Catsigeras and Budelli on general piecewise contractions [6], we extend the results obtained for the square to the more general situation when the contraction rate λ may depend on the vertices of P, and to quadrilaterals sufficiently close to the square. By the affine equivariance of T λ, the results obtained in this paper are not limited to the equilateral triangle, the square and quadrilaterals close to the square, but extend also to the polygons obtained by applying to them affine transformations. Our motivation for studying dissipative polygonal outer billiards is threefold: (1) similarly to the inner billiard [11], little is known about periodic orbits of polygonal outer billiards [29]. By a recent result of Culter, every outer polygonal billiard has a periodic orbit [31]. A detailed knowledge of the periodic points of T λ for λ < 1 may help obtain information on periodic points when λ = 1 by studying the limit λ 1; (2) the dissipative polygonal inner billiard exhibits chaotic attractors with SRB measures [1, 8, 9, 22]. It would be instructive to compare these phenomena with the dynamics exhibited by dissipative polygonal outer billiards; (3) by a result of Bruin and Deane, almost every piecewise contraction is asymptotically periodic [4]. It would be desirable to obtain a similar result for dissipative polygonal outer billiards.

4 4 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN The paper is structured as follows. In Section 2, we discuss the dissipative outer billiard about a segment, which is in a sense a onedimensional caricature of the dissipative polygonal outer billiard. Despite its somewhat artificial nature, this dynamical system shows in a nutshell some of the main features of our study. In Section 3, we introduce some notation and collect a few general results about the periodic orbits and the ω-limit set of dissipative outer billiards. The dissipative outer billiards about a triangle and about the square are studied in Sections 4 and 5, respectively. In Section 6, we address the problem of the robustness of the asymptotic periodicity under perturbations of the polygon and the contraction rates. Finally, in Section 7, we illustrate our numerical experiments and make a conjecture concerning dissipative outer billiards about general regular polygons. While writing this paper, we learned from R. Schwartz that Jeong obtained results similar to ours concerning the asymptotic periodicity of the square and the bifurcation analysis of its periodic orbits [19]. 2. One-dimensional billiards The outer billiards defined in the previous section are piecewise maps of C. In a similar fashion, one can define outer billiards on R. For these billiards, the choice of the region Ω is limited, since Ω can be only an interval I R. In this subsection, we will see that the dissipative outer billiard on R yields a family of one-dimensional transformations, exhibiting typical features of the dissipative polygonal outer billiard. We assume without loss of generality that I = [ 1, 1], and set X + = (1, + ) and X = (, 1). Then X = X + X and T λ : X ± X, where T λ (x) = λx (1 + λ) on X ±. In other words, every point of R \ I is reflected about the further endpoint of I. Since Tλ 2(X ±) X ±, by symmetry, it suffices to analyze Tλ 2 X +. For λ (0, 1), the map Tλ 2 X + is a contraction. Its unique fixed point is p λ = (1 + λ)/(1 λ). Furthermore, lim λ 0 p λ = 1 and lim λ 1 p λ = +. By the affine equivariance of our setting, the previous discussion yields the following. Proposition 2.1. Let I = [ a, a] R be any finite interval with a > 0, and let T λ : R \ I R \ I be the dissipative outer billiard about I. For 0 < λ < 1, the map T λ contracts R \ I with contraction coefficient λ 2 to the unique two-periodic orbit Moreover, we have ω λ = { a 1 + λ 1 λ, a1 + λ 1 λ }. lim ω λ = {a, a} and lim ω λ = {+, }. λ 0 λ 1 In a nutshell, the one-dimensional dissipative outer billiard is a smooth family of linear contractions. By iteration, each T λ contracts the whole

5 OUTER BILLIARDS 5 space to the unique periodic orbit. This orbit depends analytically on the dissipation parameter λ. As the parameter approaches the boundary {0, 1} of the range, the transformation T λ converges to the limit maps T 0 and T 1. The map T 0 is degenerate: T 0 (X) = { a, a}, whereas the map T 1 is the usual outer billiard about the interval [ a, a]. As the reader will see that the dissipative outer billiard about the equilateral triangle and the square retain these features, adding to them an infinite sequence of bifurcations. 3. General results We now extend the definition of the dissipative outer billiards, allowing it to have several contraction rates. Let P be a convex k-gon with vertices w = (w 1,..., w k ), and let λ = (λ 1,..., λ k ) with 0 < λ i < 1 for i = 1,..., k. The sup norm of R k is denoted by. We also denote by B r (z) C the closed ball of radius r > 0 centered at z C with respect to the norm. As before, let X = C \ P, and let A i be the cone with apex w i. The dissipative outer billiard map about P with contraction rates λ is given by T λ (z) = λ i z + (1 + λ i )w i for z A i. (3.1) As for the single contraction rate case, the singular set S of T λ consists of the half-lines containing the sides of P. For any integer n 1, let S n denote the set of points z X for which there exists an integer 0 i < n such that T i λ (z) belongs to S (so S 1 = S). We call S n the singular set of order n. Note that, S n is a finite union of half-lines, and {S n } is a filtration, i.e. S n S n+1 for every n 1. We say that a point z X is non-singular if z / S n for every n 1. For such a point z, we denote by ω(z) the ω-limit set of the T λ -orbit of z. We also denote by ω(t λ ) the union of the ω-limit sets of all nonsingular points of T λ. The following notion is central in this paper. Definition 3.1. We say that T λ is asymptotically periodic if ω(t λ ) is the union of finitely many periodic orbits, and the forward orbit of every non-singular point converges to one of them Forward invariant set and periodic orbits. Let a = λ and b = w. Denote by K λ the closed ball B r (0) of radius r = b(1 + a)/(1 a) 2. Lemma 3.1. Let z X be a non-singular point. Then there exists a sequence {i j } j N with i j {1,..., k} such that for every n N, n Tλ n (z) = ( 1)λ i1 λ in + ( 1) n j (1 + λ ij )λ ij+1 λ in w ij. (3.2) j=1 The proof of the lemma is straightforward, and so we omit it. Proposition 3.2. We have T λ (K λ ) K λ and ω(t λ ) K λ.

6 6 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN Proof. The claim follows easily from the inequality Tλ n (z) a n (1 + a) z + 1 a b for z X \ S n and n N, which is a direct consequence of (3.2). By the previous proposition, every periodic orbit of T λ is contained in K λ. Proposition 3.3. The map T λ has finitely many periodic orbits of period less than n N. Proof. Suppose that z is a periodic point of period n. Then by Lemma 3.1, there exists a sequence i 1,..., i n of integers contained in {1,..., k} such that n z = Tλ n (z) = ( 1) n λ i1 λ in z + ( 1) n j (1 + λ ij )λ ij+1 λ in v j. j=1 Since 0 < λ i1 λ in < 1, we have n j=1 z = ( 1)n j (1 + λ ij )λ λ ij+1 i n v j. (3.3) 1 + ( 1) n+1 λ i1 λ in Thus, z is determined by the sequence i 1,..., i n. Since there are at most k n of these sequences, we conclude that there are at most k n periodic points of period n. The simplest periodic orbit of the outer billiard about a convex k-gon has rotation number 3 1/k. The sequence of the vertices of P visited by this orbit is the infinite periodic sequence w 1, w 2,..., w k, w 1, w 2,.... This orbit is sometimes called Fagnano 4 orbit [30]. Corollary 3.4. There exists 0 < δ < 1 depending only on the vertices of P such that if λ < δ, then the Fagnano orbit exists, and is the global attractor for T λ. Proof. Suppose that v 1,..., v k is a Fagnano orbit of P. Relabeling the vectors v i, if necessary, and setting v k+1 = v 1, we can write w i = λ i v i + 1 v i+1 for i = 1,..., k. 1 + λ i 1 + λ i The previous equations define a linear map (v 1,..., v n ) (w 1,..., w n ) which is invertible provided λ 1 λ k ( 1) k. Since this condition is sufficient and necessary for the existence of the Fagnano orbit, the Fagnano orbit exists for dissipative outer billiards. Every periodic orbit of a dissipative outer billiard is locally attracting. Thus, to prove that the Fagnano orbit is a global attractor for 3 Given a periodic orbit γ of T λ, let Γ be the closed curve obtained by joining consecutive points of γ using straight lines. The rotation number of γ is defined to be the winding number of Γ divided by k.

7 OUTER BILLIARDS 7 T λ, it suffices to show that if λ is sufficiently small, then every nonsingular point of K λ eventually enters a small neighborhood of the Fagnano orbit. This property is a direct consequence of the following fact: for every non-singular z K λ A i with i {1,..., k}, we have T λ (z) v i+1 = T λ (z) T λ (v i ) = λ i z v i λ diam(k λ ) Skew-product. In this subsection, we assume that P is a regular k-gon, and that all the components of λ are equal to some λ (0, 1). This situation is special, because the map T λ commutes with the rotation of C about the center of P by an angle 2π/k. This fact allows us to conjugate T λ to a skew product transformation. Since λ is fixed in this subsection, we will drop the label λ from our notation. Denote by R the clockwise rotation of C about the center of P by an angle 2π/k, and define π : i A i A 0 by π(z) = R i (z) if z A i. Let S = A 0 T 1 (S). Next, let f : A 0 \ S A 0 be the map defined by f(z) = π T (z) for z A 0 \ S. It is easy to see that A 0 \ S consists of finitely many connected open sets D 1,..., D [k/2]+1 with boundary contained in S, and that f is an affine transformation on each D i. Note that the sets D i are characterized by the property T (D i ) A i. Since f π = π T by definition, f is a factor of T. For every integer n > 0, the iterate f n is defined only on A 0 \ S n, where S n is the set of z A 0 such that f m (z) S for some integer 0 m < n. Let S = n 0 S n. Every iterate of f is defined on A 0 \ S. Define σ : A 0 \ S Z k by σ(z) = i whenever z D i. Then, the skew-product F : (A 0 \ S) Z k A 0 Z k is given by F (z, j) = (f(z), σ(z) + j mod k), (z, j) (A 0 \ S) Z k. If we set σ n (z) = σ(f n 1 z) + + σ(z) for n 1 and z A \ S n, then we have F n (z, j) = (f n (z), σ n (z) + j mod k). The transformations F and T are clearly conjugated. By abuse of notation, we also use the symbol π for the projection from A 0 Z k A 0 given by π(z, i) = z. Every periodic orbit of F (and so of T ) is mapped by π into a single periodic orbit of f. The converse relation is less obvious, and is clarified in the next lemma. Lemma 3.5. Suppose that z A 0 \ S is a periodic point of f, i.e., f m (z) = z for some integer m > 0. Let O f (z) be the f-orbit of z. Then the set π 1 (O f (z)) consists of gcd(k, σ m (z)) periodic orbits of F of period mk/ gcd(k, σ m (z)). Proof. Let x π 1 (O f (z)). Then x = (z, i) for some z O f (z) and some 0 i < k. For every integer j > 0, we have F jm (x) = (f jm (z ), jσ m (z ) + i mod k) = (z, jσ m (z) + i mod k).

8 8 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN E S λ 1 v C A v A 2 Figure 2. Dissipative outer billiard about the equilateral triangle It follows that x is a periodic point of F of period equal to p := m min{j > 0: jσ m (z) = 0 mod k} = mk/ gcd(k, σ m (z)). This shows that every element of π 1 (O f (z)) is a period point of F of period p. Since π 1 (O f (z)) consists of c := mk elements, we conclude that π 1 (O f (z)) is the union of c/p = gcd(k, σ m (z)) periodic orbits of F. 4. Triangular billiards In this section we study the dissipative outer billiard about a triangle. We consider the special case λ 1 = λ 2 = λ 3 = λ (0, 1). By affine equivariance, we do not lose generality by restricting our analysis to the equilateral triangle of unit side length. Let v = e 2π 3 i. The corners of are 0, 1, v. Let T λ be the outer billiard map on. Hereafter, we shall omit λ from our notation. Let A 0, A 1 and A 2 be the three cones for T. The apex of A 0 is the origin and its sides are the rays {a : a 0} and {bv : b 0}. The cones A 1 and A 2 are determined by symmetry. Recall from Subsection 3.2 that S is the set of points in A 0 that are mapped by T to the singular set S, i.e. S = A 0 T 1 (S). For the equilateral triangle, S is the ray {λ 1 v + a(1 + v) : a 0}. Note that, T (A 0 ) A 1 A 2. Let A = A 0 and let C = T 1 (T (A 0 ) A 1 ) and E = T 1 (T (A 0 ) A 2 ). A simple computation shows that C is the polygonal region in A with one finite side [0, λ 1 v] and two infinite sides: S and the positive real axis. Clearly, E = A \ C. Note that the sets C and E coincide with the sets D 1 and D 2 introduced in Subsection 3.2. These sets are depicted in Fig. 2. We now describe the map f. Rotating properly the sets T (C) and T (E) about the center of the triangle we can reduce T to a piecewise

9 OUTER BILLIARDS 9 affine map f : A \ S A whose defining regions are C and E. Let f C = f C and f E = f E. It is easy to see that f C (z) = 1 + λ (1 + v)z, f E (z) = v λ vz. Note that f(e) = {a + b(1 + v): a > 1 and b > 0} C. Indeed, since any z E can be written in the form z = λ 1 v + av + b(1 + v) for some a, b > 0, we have f E (z) = v λ v(λ 1 v + av + b(1 + v)) = 1 + λa(1 + v) + λb. Thus f(e) C. Also note that, C and E depend on λ while f(c) and f(e) do not. Let B be the set of points a+b(1+v) C such that 0 < b < λ 1 +λ 2. Also, let C 1 = B and C n, n 2 be the set of points a + b(1 + v) C such that λ λ 2n+2 < b < λ λ 2n. The sets B and C n are depicted in Fig. 3. Lemma 4.1. f 2 (C n ) C n 1 for every integer n 2. Proof. Let z C n for some integer n 2. Then z = a + b(1 + v) with a > λ 1 and b > λ 1 + λ 2. Thus f(z) = f C (z) = (λb 1)v + (1 + aλ)(1 + v). Since, λb 1 > λ 1 and 1+aλ > 0 we conclude that f(z) E. Iterating again we get, f 2 (z) = f E f C (z) = α + β(1 + v) S C 3 λ 1 v + (λ 2 + λ 3 + λ 4 )(u + v) E C 2 λ 1 v + λ 2 (u + v) λ 1 v C 1 = B 0 1 v Figure 3. Sets B and C n.

10 10 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN S λ 1 v B 1 B 2 B 3 0 Figure 4. The sets B n. where α = 1 + λ + aλ 2 and β = λ 2 b λ 1. Taking into account that a > λ 1 and λ λ 2n+2 < b < λ λ 2n we get that α > 1 and λ λ 2n+4 < β < λ λ 2n+2. Thus, if z C n then f 2 (z) C n 1. Now we study the images of B under iteration of f. On can easily show that f(b) C (a consequence of Lemma 4.2 below). Let B n, n 1 be the points in B that are mapped to C n after one iteration of f, i.e. B n = f 1 (B) C n, n 1. Lemma 4.2. For every integer n 1 the set B n consists of points a + b(1 + v) B such that λ λ 2n+1 < a + b < λ λ 2n 1. Proof. Let z = a + b(1 + v) B. Then f C (z) = α + β(1 + v) where α = 1 λb and β = λ(a + b). From a > λ 1, a + b > 0 and 0 < b < λ 1 + λ 2 we conclude that f C (B) = { α + β(1 + v): λ 1 < α < 1, β > 0 and α + β > 0 }. Thus, f(b) C and f(b) C n for every n 1. Now f(z) C n if and only if λ λ 2n+2 < β < λ λ 2n. Taking into account that β = λ(a + b) we get the desired inequality. The sets B n are depicted in Fig. 4. Let B = n 1 B n. By previous considerations, we obtain that every point in B comes back to B after an odd number of iterations of f. In fact, for every n 1 we have B n = {z B : f 2n 1 (z) B and f k (z) / B if 1 k < 2n 1}. On B we define the Poincaré map φ: B B, i.e. the first return map to B. Let φ n denote the restriction of φ to B n. Lemma 4.3. For every integer n 1 φ n (z) = λ 2n 2 v(1 + + λ 2n 3 ) + λ 2n 1 (1 + v)z. Proof. By Lemma 4.1, φ n = (f E f C ) n 1 f C. The proof goes by induction on n and is left to the reader.

11 OUTER BILLIARDS 11 In the following lemma we prove that φ has some monotonicity property. Lemma 4.4. For every integer n 1 the following holds: (1) φ(b n ) B 1 B n+1, (2) φ(φ(b n ) B ) B 1 B n. Proof. We start by proving item (1). Let z = a + b(1 + v) B n. By Lemma 4.3, we have φ n (z) = α + β(1 + v) where α = λ 2n 2 bλ 2n 1, β = λ 2n 1 (a + b) 1 λ 2n 3. Taking into account the definition of B n+1 we only need to show that α+β < λ 2 + +λ 2n 3. But since, b > 0 and a+b < λ 2 + +λ 2n 1 we conclude that α + β = λ 2n 2 + λ 2n 1 (a + b) bλ 2n 1 < λ 2 + λ λ 2n 2 < λ λ 2n 3 for every λ (0, 1), and the first claim follows. To prove item (2), let z = a + b(1 + v) B n φ 1 (B ). We want to show that φ 2 (z) B 1 B n. By item (1), we only need to consider two cases: (a) z B n φ 1 (B n ). Thus φ 2 (z) = φ 2 n(z) = α + β (1 + v) where α + β = λ 2n 2 (1 + λ + + λ 2n 1 bλ 2n ). Since b > 0 we get α + β < λ λ 2n 1 for every λ (0, 1). Thus (2) follows. (b) z B n φ 1 (B n+1 ). In this case φ 2 (z) = φ n+1 φ n (z). As before we have φ n+1 φ n (z) = α + β (1 + v) where α + β = λ 2n (1 + λ + + λ 2n 1 bλ 2n ). Taking into account that b > 0 we get α + β < λ λ 2n 1 for every λ (0, 1). Thus (2) follows again in this case. Lemma 4.5. For every λ (0, 1) there exists an integer N = N(λ) > 0 such that φ(b ) B 1 B N. Proof. The claim follows from the fact that, Re(φ(z)) < λ 2n < 1 1 λ z B n.

12 12 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN Let S be the set of points in B whose orbit eventually intersect a discontinuity point of φ. Let ψ n denote the affine extension of φ n to the complex plane and denote by z n the fixed point of the extension ψ n. Also denote by w n the unique fixed point of the map ψ n+1 ψ n. Note that u n := ψ n (w n ) is the unique fixed point of the map ψ n ψ n+1. Finally, let ψ = {ψ n } and denote by Per(ψ) the union of all z n, w n and u n. Proposition 4.6. For every z B \ S there exists an integer n 1 and z Per(ψ) such that lim k φ 2k (z) = z B n. Proof. Given z B \ S, by Lemma 4.4, there exist integers k 0 0 and n 1 such that φ 2k (z) B n for every k k 0. Let w = φ 2k 0 (z). Thus φ 2k (w) B n for every k 0. Again, by Lemma 4.4, φ 2k+1 (w) B n 1 B n B n+1 for every k 0. Thus we have three cases: (1) There exists k 1 0 such that φ 2k1+1 (w) B n 1. Then, by Lemma 4.4 we have that φ 2k+1 (w) B n 1 for every k k 1. Thus φ 2k (w) = (φ n 1 φ n ) k (w) for every k k 1, which implies that lim k φ 2k (z) = u n 1. (2) There exists k 1 0 such that φ 2k1+1 (w) B n and φ 2k+1 (w) / B n 1 for every k 0. Then, again by Lemma 4.4 we have that φ 2k+1 (w) B n for every k k 1. Thus φ 2k (w) = φ 2k n (w) for every k k 1, which implies that lim k φ 2k (z) = z n. (3) φ 2k+1 (w) B n+1 for every k 0. Then, φ 2k (w) = (φ n+1 φ n ) k (w) for every k 0, which implies that lim k φ 2k (z) = w n. Using the terminology introduced in Section 1, we conclude that φ is asymptotically periodic. Let ω(φ) denote the ω-limit set of φ. Corollary 4.7. For every 0 < λ < 1, let N = N(λ) be the integer in Lemma 4.5. Then ω(φ) {z 1, w 1, u 1,..., z N, w N, u N }. Proof. The claim follows from Lemma 4.5 and Proposition 4.6. Next, we perform a bifurcation analysis of the ω-limit set of φ. Lemma 4.8. Let z = a + b(1 + v) such that z = ψ n (z). Then z B n if and only if λ λ 2n+1 < a + b < λ λ 2n 1. Proof. By Lemma 4.3, ψ n (z) = a+b(1+v) where b = λ 2n 1 (a+b) 1 λ 2n 3. Thus 0 < b < λ 1 + λ 2 if and only if λ λ 2n+1 < a + b < λ λ 2n 1. The claim follows. Proposition 4.9. z n is a fixed point of φ n if and only if λ (λ 2n 1, 1) where the sequence {λ m } is defined in the Appendix.

13 OUTER BILLIARDS 13 Proof. We want to show that z n B n if and only if λ (λ 2n 1, 1). By Lemma 4.3, z n = a n + b n (1 + v) where a n + b n = By Lemma 4.8, z n B n if and only if λ 2n 2 λ 4n 3 (1 λ) (1 λ 2n 1 + λ 4n 2 ). λ λ 2n+1 < a n + b n < λ λ 2n 1. Substituting the expression for a n + b n into the previous inequalities we get λ 2 (1 λ 2n 2 ) < λ4n 1 (1 λ 2n 1 ) 1 λ 2n 1 + λ 4n 2 < 1 λ2n. We study each inequality separately. The right-hand side inequality is equivalent to 1 λ 2n 1 λ 2n + λ 4n 2 > 0, which holds for every λ (0, 1). Indeed, 1 λ 2n 1 λ 2n + λ 4n 2 > (1 λ 2n 1 ) 2 > 0. The left-hand side inequality is equivalent to p 2n 1 (λ) = 1 λ 2n 2 λ 2n 1 + λ 4n 2 < 0. By Lemma A.1, we conclude that p 2n 1 (λ) < 0 for every λ (λ 2n 1, 1). This concludes the proof of the proposition. Proposition w n is a fixed point of φ n+1 φ n if and only if λ (γ 2n 1, 1) where the sequence {γ m } is defined in the Appendix. Proof. We want to show that w n B n and u n B n+1 if and only if λ (γ 2n 1, 1). By Lemma 4.3, w n = a n + b n (1 + v) and u n = ψ n (w n ) = α n + β n (1 + v) where a n + b n = λ 2n (1 λ 6n 1 ) (1 λ) (1 + λ 4n + λ 8n ), α n + β n = λ 2n 1 a n + λ 2n 2. By Lemma 4.8, w n B n and u n B n+1 if and only if λ λ 2n+1 < a n + b n < λ λ 2n 1, λ λ 2n 1 < α n + β n < λ λ 2n 3. Substituting the expressions for a n + b n and α n + β n into the previous inequalities we get λ 2 (1 λ 2n 2 ) < λ4n+1 (1 λ 6n 1 ) 1 + λ 4n + λ 8n < 1 λ 2n, λ 2 (1 λ 2n ) < λ4n+1 (1 λ 6n+1 ) 1 + λ 4n + λ 8n < 1 λ 2n+2.

14 14 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN It turns out that, among these four inequalities, there are two which are equivalent. In fact, we can reduce the above chain of inequalities to the following system: 1 λ 2n 2 λ 4n 1 + λ 4n λ 6n + λ 8n < 0, 1 λ 2n + λ 4n λ 4n+1 λ 6n + λ 8n > 0, 1 λ 2n λ 4n 1 + λ 4n λ 6n + λ 8n < 0. We claim that this system of inequalities holds for every λ (γ 2n 1, 1). Let I 1, I 2 and I 3 denote the first, second and third polynomials in the previous system. Note that I 1 (λ) < I 3 (λ) and I 2 (λ) > 0 for every λ (0, 1). Indeed, I 3 (λ) I 1 (λ) = λ 2n 2 (1 λ 2 ) > 0. Moreover, since λ 2n + λ 6n > 2λ 4n and λ 4n > λ 4n+1 for every λ (0, 1) we get I 2 (λ) > 1 4λ 2n + 6λ 4n 4λ 6n + λ 8n = (1 λ 2n ) 8 > 0. Finally, by Lemma A.2, I 3 (λ) < 0 for every λ (γ 2n 1, 1) and the proof of the proposition is complete. Remark Taking into account the return time to B, we conclude that the fixed point of φ n corresponds to a period orbit of period 2n 1 of f and the fixed point of φ n+1 φ n corresponds to a periodic orbit of period 4n of f. For every λ (0, 1) let m 1 = m 1 (λ) > 0 and m 2 = m 2 (λ) > 0 be the following integers, m 1 = max{n N: λ 2n 1 < λ λ 2n+1 }, m 2 = max{n N: γ 2n 1 < λ γ 2n+1 }. The following result completely describes the bifurcations of the ω-limit set of φ. Theorem For every λ (0, 1) we have, ω(φ) = {z 1,..., z m1, w 1, u 1,..., w m2, u m2 }. Proof. When λ / {λ 2n 1 } {γ 2n 1 }, the claim follows from Proposition 4.6, 4.9 and So, suppose that λ = λ 2n 1 for some n 2. We want to show that z n / ω(φ). It follows from the proof of Proposition 4.9, that z n = a n + b n (1 + v) where a n = λ λ 2n+1 and b n = 0. Thus, z n is located at the bottom left corner of B n. Since, by Lemma 4.3, ψ n is a composition of an anticlockwise rotation of angle π/3 and a contraction towards z n, we conclude that under the iteration of φ n every point close to z n will eventually enter B n 1. Thus, moving away from z n. This proves that z n / ω(φ). Now suppose that λ = γ 2n 1 for some n 2. As before we want to show that w n, u n / ω(φ). It is sufficient to show that u n / ω(φ). It follows from the proof of Proposition 4.10, that u n = α n +β n (1+v) where

15 OUTER BILLIARDS 15 Figure 5. Dissipative outer billiard about the equilateral triangle for λ = The billiard map has three periodic orbits Γ 1,1, Γ 1,2 and Γ 2,1 of periods 3, 9 and 12, respectively. The basin of attraction of each periodic orbit is the union of all the regions of the same color. α n +β n = λ 2 + +λ 2n 1. Hence u n belongs to the common boundary of B n and B n+1. Since, by Lemma 4.3, ψ n ψ n+1 is composition of an anticlockwise rotation of angle 2π/3 and a contraction towards u n, we conclude that under the iteration of φ n φ n+1 every point in B n+1 close to u n will eventually enter B n. Thus, by Lemma 4.4, every point close to u n will eventually enter B 1 B n and never leave this set. This proves that u n / ω(φ) and concludes the proof of the theorem. Denote by dist the Euclidean distance on C. The following result is a consequence of Theorem 4.12 and Lemma 3.5. Corollary For every λ (0, 1), the map T has exactly m 1 = m 1 (λ) distinct periodic orbits Γ 1,i, i = 1,..., m 1, each one having period 3(2i 1) and exactly m 2 = m 2 (λ) distinct periodic orbits Γ 2,i, i = 1,..., m 2, each one having period 12i. Moreover, for every nonsingular z X, there exist i {1, 2} and 1 j m i such that dist(t n (z), Γ i,j ) 0 as n +. (see Fig. 5). Proof. Let 0 < λ < 1. Since T and F (see Subsection 3.2) are conjugated, we prove the corollary with T replaced by F. By Theorem 4.12 and Lemma 3.5, it follows that O f (z 1 ),..., O f (z m1 ), O f (u 1 ),..., O f (u m2 ) are the only periodic orbits of f. Since z i is a fixed point of φ i, it follows from the definition of φ i that the period of O f (z i ) is equal to 2i 1, and that f 2i 1 (z i ) = (f E f C ) i 1 f C (z i ). An easy computation then

16 16 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN shows that σ 2i 1 (z i ) = 3i 2 mod 3. Since 3i 2 and 3 are coprime, Lemma 3.5 implies that each set Γ 1,i := π 1 (O f (z i )) is a periodic orbit of F of period 3(2i 1). A similar argument shows that Γ 2,i := π 1 (O f (u i )) is a periodic orbit of F of period 12i. Since f is a factor of F, all periodic orbits of F are contained in the union of the sets Γ i,j, i = 1, 2 and j = 1,..., m i. Now, let z be a non-singular point of F. Theorem 4.12 implies that there exists z Per(ψ) such that dist(f n (z), O f (z )) 0 as n. By previous considerations, π 1 (O f (z )) = Γ i,j for some i = 1, 2 and j = 1,..., m i. Finally, note that f is a factor of F and that dist(w, π 1 (D)) = dist(π(w), D) for every w X and every set D A. Combining all the previous observations together, we conclude that lim dist(f n (z), Γ i,j ) = lim dist(f n(π(z)), O f (z )) = 0. n + n + 5. Square billiard We now study the dissipative outer billiard about the square. As for the triangular billiard, we consider the special case λ 1 = = λ 4 = λ (0, 1). In our analysis, we follow closely the approach used for the dissipative outer billiard about the equilateral triangle. In particular, we will construct a first return map φ on a strip B and obtain a detail description of the dynamics of the dissipative outer billiard map T. For the square, the geometry of the action of the map φ is simpler, because the singular set of the billiard map consists of horizontal and vertical lines. In this section, we will formulate a result without proof, every time that the proof can be easily recovered from the corresponding ones for the equilateral triangle. For convenience, we assume that the corners of the square are given by w 0 = 0, w 1 = i, w 2 = 1 i and w 3 = 1. Let A i be the open conical region corresponding to w i (see Fig. 6). So A 0 = {Re(z) > 0 and Im(z) > 0} is the positive quadrant of C. Let T λ be the outer billiard map on the square with λ (0, 1). From now on, we will drop the label λ from our notation. We have T (A 0 ) A 1 A 2, and rotating clockwise A 1 T (A 0 ) and A 2 T (A 0 ) about the center of the square by an angle of 90 and 180 degrees, respectively, we bring these regions back to A 0 (see Fig. 6). We now compute explicitly the map f for the square. Let B = T 1 (A 1 T (A 0 )), and let C = T 1 (A 2 T (A 0 )). These sets were denoted by D 1 and D 2 in Subsection 3.2. It is easy to check that B = {Re(z) > 0, Im(z) < 1/λ} and C = {Re(z) > 0, Im(z) > 1/λ} (see Fig. 6).

17 OUTER BILLIARDS 17 Figure 6. Outer billiard map about the unit square and reduction to the first quadrant. Next, define f B = f B and f C = f C. Note that f A and f B are precisely the maps f 1 and f 2 defined in Subsection 3.2. A simple computation shows that f B (z) = iλz + 1, f C (z) = λz + 1 i. (5.1) Let C 1 = B, and let C n+1 = f 1 C (C n) C for n 1. For n 2, the set C n consists of elements of C hitting B exactly after n iterations of f. In the next lemma, we prove that every non-singular point of C visits B eventually. Lemma 5.1. C \ S n 2 C n. Proof. We argue by contradiction. Suppose that there exists z C such that z k := fc k(z) C for all k 1. Thus lim inf k + Im(z k ) 1/λ. However, by (5.1), we have z k (1 i)/(1 λ) as k, yielding a contradiction. Denote by B n the subset of B formed by elements whose first return time to B is equal to n 1. From the definition of C n, it follows that B n = f 1 B (C n) B for n 1. Let B = n B n, and denote by φ: B B the first return map to B induced by f. Also, let φ n = φ Bn. Lemma 5.2. We have (1) φ(b n ) = φ n (B n ) n i=1 B i for n 1, (2) N = N(λ) such that φ(b ) N i=1 B i. Proof. Similar to the proof of Lemma 4.4. See also Fig. 7. Let ψ n be the affine extension of φ n to the whole C. Since ψ n is a strict contraction, it has a unique fixed point z n. Note that z n is a fixed

18 18 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN Figure 7. Rectangles B n. point of φ n (and so of φ) if and only if z n B n. It is also clear that if z n is a fixed point of φ n, then it is also a periodic point of f and T. The following proposition is an immediate consequence of Lemma 5.2. Proposition 5.3. Let N(λ) be the constant in Proposition 5.2. For every z B \ S, there exists an integer 0 < n N(λ) such that lim k + φ k (z) = z n B n. Using the terminology introduced in Section 1, we can rephrase the conclusion of the previous proposition by saying that φ is asymptotically periodic. Let ω(φ) denote the ω-limit set of φ. The next corollary follows directly from Proposition 5.3. Corollary 5.4. We have ω(φ) {z 1,..., z N(λ) } for every λ (0, 1). A bifurcation analysis of the fixed points {z n } n 1 can be performed similarly to the corresponding analysis for the equilateral triangle. From this analysis, we can derive a precise description of ω(φ) and ω(t ). The following result can be proved by adapting the arguments used in Section 4, in particular the proofs of Proposition 4.9 and Theorem Let {λ i } i 1 be the sequence in Lemma A.1. Theorem 5.5. Let m(λ) > 0 be the integer defined by λ m(λ) < λ λ m(λ)+1. Then ω(φ) = {z 1,..., z m(λ) } for every λ (0, 1). The following is a corollary of Theorem 5.5 and Lemma 5.1. Its proof is similar to the one of Corollary Corollary 5.6. For every λ (0, 1), the map T has exactly m = m(λ) distinct periodic orbits Γ 1,..., Γ m, and the period of Γ i is equal to 4i. Moreover, for every non-singular z X, there exists 1 i m such that dist(t n (z), Γ i ) 0 as n +. (see Fig. 8).

19 OUTER BILLIARDS 19 Figure 8. Dissipative outer square billiard for λ = The billiard map has three periodic orbits of period 4, 8, 12. The basin of attraction of each orbits is the union of all the regions of the same color. Remark 5.7. It is worth observing that the bifurcation value λ 2i 1 is the same for the periodic orbits of period 4(2i 1) of the dissipative outer billiard about the square and the periodic orbits of period 3(2i 1) of the dissipative outer billiard about the equilateral triangle. 6. Persistency In this section, we prove that the dissipative square billiard remains asymptotically periodic under small perturbations of its vertices and the contraction rates. The precise statement is given in Theorem 6.1. Let P k be the collection of all pairs (P, λ) consisting of a convex k- gon P with no more than two vertices lying on the same side and one vertex at the origin of C. Let λ (0, 1) k. We can naturally identify a k-gon P with the vector w R 2k 2 whose components are the vertices w 1,..., w k 1 of P chosen so that if P is positively oriented, then w 1 follows the vertex at the origin, and w i+1 follows w i for i = 1,..., k 1. In virtue of this identification, P k is an open subset of R 3k 2 and a locally compact metric space with the sup norm of R 3k 2. Now, let P be the square, and let Λ = {λ 1, λ 2,...} be the set of all bifurcation values for the dissipative square billiards (see Theorem 5.5). In this section, we will write T P,λ instead of T λ to emphasize the dependence of the billiard map on the polygon P. Theorem 6.1. Let P be the square, and suppose that λ i = λ (0, 1)\Λ for i = 1,..., 4. Then for every ɛ > 0, there exists δ > 0 such that if (P, λ ) P 4 and (P, λ ) (P, λ) < δ, then (1) T P,λ is asymptotically periodic,

20 20 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN (2) card(ω(t P,λ )) = card(ω(t P,λ)), (3) the Hausdorff distance of ω(t P,λ ) and ω(t P,λ) is less than ɛ, (4) corresponding periodic orbits of T P,λ and T P,λ have the same period. Proof. The theorem follows from Propositions 6.4 and 6.5. The next corollary is an immediate consequence of the previous theorem. Corollary 6.2. The dissipative outer billiard about any quadrilateral sufficiently close to the square with contraction coefficients sufficiently close to some real in (0, 1) \ Λ is asymptotically periodic. We now recall the notions of piecewise locally contracting maps and persistency introduced in [6]. We present these notions only for twodimensional maps, which is the situation of interest here. Let K be a compact subset of C with the metric inherited from C. We call 5 Φ = {Φ i } a piecewise contraction on K if there exist a constant 0 < ρ < 1 and a finite compact cover C = {C 1,..., C m } of K such that C i C j i C i C j for i j, and Φ i : C i Φ(C i ) K is a bijective strict contraction with Lipschitz constant L Φi ρ for every i. Let L Φ = max i L Φi. The map Φ is defined by Φ(A) = i Φ i(a C i ) for every A K. The singular set of Φ is denoted by S Φ = i j C i C j. Note that Φ is multivalued transformation on S Φ. We will write (Φ, C) to emphasize the role of the cover C in the definition of Φ. The space of all piecewise contractions on K is endowed with a topology defined by the following base. Given a piecewise contraction (Φ, C) on K and a number ɛ > 0, the element B(Φ, ɛ) of the base is given by all piecewise contractions (Ψ, D) on K such that (1) the covers C and D have the same cardinality, (2) sup z Ci D i Φ(z) Ψ(z) < ɛ for every i, (3) L Φ L Ψ < ɛ, (4) the Haussdorf distance of C i and D i is less than ɛ for every i. Suppose that Φ is piecewise contraction that is asymptotically periodic. We say that the ω-limit set of Φ is persistent provided that Φ is asymptotically periodic, and that for every ɛ > 0, there exists δ > 0 such that if Ψ B(Φ, δ), then (1) Ψ is asymptotically periodic, (2) card(ω(ψ)) = card(ω(φ)), (3) the Hausdorff distance of ω(ψ) and ω(φ) is less than ɛ, (4) corresponding periodic orbits of Ψ and Φ have the same period. The following proposition, proved by Catsigeras and Budelli, gives a sufficient condition for the persistency of ω(φ) [6, Lemma 3.3 and Remark 3.4]. 5 These transformations are called piecewise continuous locally contracting maps [6].

21 OUTER BILLIARDS 21 Proposition 6.3. If Φ n (K) S Φ = for some integer n > 0, then Φ is asymptotically periodic, and ω(φ) is persistent. We now reformulate Proposition 6.3 in the context of billiards maps. Given (P, λ) P k for some k, let K P,λ be the forward invariant compact set for the map T P,λ as in Proposition 3.2, and let S TP,λ be the singular set of T P,λ. Proposition 6.4. Let (P, λ) P k, and suppose that T n P,λ (K P,λ) S TP,λ = for some integer n > 0. Then T P,λ is asymptotically periodic, and ω(t P,λ ) is persistent with respect to the topology of P k. Proof. Let O be a compact neighborhood of P k containing (P, λ). Define a = sup λ and b = sup w(p ), (P,λ) O (P,λ) O where w(p ) denotes the vector formed the vertices of P as explained at the beginning of this section. Since O is compact, we have 0 < a < 1 and 0 < b <. It follows easily from Proposition 3.2 that there exists a compact set W C forward invariant for every T P,λ with (P, λ ) O. Thus T P,λ is a piecewise contraction6 on W in the sense of Catsigeras and Budelli for every (P, λ ) O. Now, denote by Π the transformation associating a pair (P, λ ) O to the piecewise contraction T P,λ on W. If O is endowed with the Euclidean metric of R 3k 2, then it follows that Π is continuous. From the definition of W, we see immediately that K P,λ is contained in W. Then, Proposition 3.2 implies that TP,λ m (W ) K P,λ. Using TP,λ n (K P,λ) S TP,λ =, we obtain T n+m P,λ (W ) S T P,λ =. The wanted conclusion now follows from Proposition 6.3 and the continuity of Π. Next, we show that the hypothesis of Proposition 6.4 is satisfied for the square with λ 1 = = λ 4 = λ (0, 1) \ Λ. Proposition 6.5. Let P be the square, and suppose that λ i = λ (0, 1) \ Λ for i = 1,..., 4. Then there exists an integer n > 0 such that TP,λ n (K P,λ) S TP,λ =. Proof. First, we extend T P,λ from A i to A i. This is clearly possible, because T P,λ is linear on each A i. Then, we redefine T P,λ (E) = k i=1 T P,λ(E A i ) for every set E X. Note that T P,λ is now a multivalued map on S TP,λ with T P,λ (z) consists of at most of as many points as the number of singular lines meeting at z. The ω-limit set ω(z) has to be understood now as the set of accumulations points of the sequence of sets {T n P,λ (z)} n 0. Of course, we still have T P,λ (K P,λ ) K P,λ. 6 The map T P,λ is not defined on P. This can be easily fixed, by choosing u P, and defining T P,λ (z) = min i λ i (z u) for z P. We also need to extend the map T P,λ on each domain of continuity up to its boundary.

22 22 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN For the extended map T P,λ, Corollary 5.6 can be reformulated as follows: let m(λ) and Γ 1,..., Γ m(λ) be the positive integer and the periodic orbits, respectively, as in Corollary 5.6. If λ / Λ, then ω(z) m(λ) i=1 Γ i for every z X. This statement can be easily deduced from the proof of the corollary, which is similar to that of Corollary Note that the statement holds for every z X, and not just nonsingular z X as in Corollary 5.6. If λ Λ, then the statement fails, because there exists a point z S P,λ such that z TP,λ n (z) for some n > 0, and so ω(z) S P,λ. The statement could have been easily modified to cover the case λ Λ as well, but that is not required for the purpose of this proof. Since each periodic orbit Γ i is asymptotically stable, the previous reformulation of Corollary 5.6 implies that for every z K λ, there exists a neighborhood U X of z and an integer n > 0 such that T i (U) S P,λ = for all i n. By a standard compact argument, we conclude that TP,λ n (K P,λ) S P,λ = for some integer n > 0. This completes the proof. Remark 6.6. Following Jeong [19], one can show that the dissipative outer billiard about the square, the equilateral triangle and the hexagon is asymptotically periodic whenever λ < 1. The argument is quite simple, and we will illustrate it in the rest of this section. Let T λ be the dissipative outer billiard about the unit square. When λ 1 = = λ k = 1, the billiard has an explicit integral of motion I(z) defined as follows. Given z A i, let x i = x i (z) and y i = y i (z) be the unique positive numbers such that z = w i +x i (w i 1 w i )+y i (w i w i+1 ) for i = 0,..., 3, where w 4 = w 0. Define h(x, y) = [x] + [y], where [a] is the integer part of a R. A simple computation shows that I(T λ (z)) = I(z) for z X \ S, i.e., I(z) := h(x i (z), y i (z)) is an integral of motion. However when λ < 1, we have I(T λ (z)) I(z), and I becomes a Lyapunov function for T λ. Now, since the sequence of the vertices visited by an orbit contained in a set {I = α} is completely determined by α N, it is not difficult to show that the map T λ is asymptotically periodic for every λ (0, 1) 4. The proof that the dissipative outer billiards about the equilateral triangle and the regular hexagon are both asymptotically periodic for λ < 1 is analogous. The Lyapunov function for both billiards is defined as for the square, but with h = h for the equilateral triangle, and with h = h + for the hexagon, where [ x ] [ [ ] y x ± y + 1 h ± (x, y) := ] 2

23 OUTER BILLIARDS 23 To complete this remark, we observe that while the previous argument proves the asymptotic periodicity for the dissipative outer billiards considered, it does not give a complete description of the bifurcations of the ω-limit set, for which the analysis described in Sections 4 and 5 or an alternative one is still required. 7. Concluding remarks According to Bruin and Deane [4], almost every piecewise contraction is asymptotically periodic. This result, however, does not apply to polygonal dissipative outer billiards. The reason is that different polygons generate different partitions of the domain of definition of the corresponding outer billiard map, whereas in Bruin and Deane s setting, the partition is fixed a priori for the entire family of piecewise contractions. Not every dissipative outer polygonal billiard is asymptotically periodic. There exist quadrilaterals whose dissipative outer billiard has a Cantor set lying on the singular set, which attracts nearby points [19]. Numerical experiments we performed on dissipative outer billiards about regular polygons suggest the following conjecture. Conjecture 7.1. For almost every λ (0, 1), the dissipative outer billiard map T λ of any regular polygon is asymptotically periodic. Let K λ be the closed ball associated to T λ defined in Subsection 3.1. We say that the singular set of T λ stabilizes if there exists an integer n 1 such that S n+1 K λ = S n K λ. Proposition 7.2. If the singular set of T λ stabilizes, then ω(t λ ) is finite. Proof. If the singular set of T λ stabilizes then there exists an integer n 1 such that the map Tλ k is continuous on each connected component of K λ \ S n for every k > 0. Since K λ \ S n has finitely many such components, say Y 1,..., Y m, and the map T λ contracts uniformly every Y i into some Y j, we conclude that ω(t λ ) is finite. We cannot claim that every element of ω(t λ ) in the previous proposition is a periodic orbit, because some of these elements may belong to the singular set of T λ. However, such an element does generate a periodic orbit for a proper extension of T λ up to S (degenerate periodic orbit). We believe that degenerate periodic orbits exist only in very specific circumstances, for example at the bifurcation of a true periodic orbit. Even though the stabilization of the singular set is not strictly speaking a sufficient condition for the asymptotic stability, it still provides a useful criterion for investigating asymptotic periodicity that can be effectively implemented numerically. We wrote a script in Mathematica that generates the singular sets of any order of a polygonal outer billiard, and checks whether they

24 24 G. DEL MAGNO, J. P. GAIVÃO, AND E. GUTKIN (a) (b) Figure 9. (A) Pentagon with λ = 0.95 and (B) magnification of (A). stabilize. When the singular set stabilizes at some order n, the script attributes a color to the connected components of K λ \ S n. The union of all regions of the same color gives the basin of attraction in K λ of a single attracting periodic orbit. To test numerically Conjecture 7.1, we run our code with different regular k-gons with 5 k 12, and different contraction rates λ i (0, 1). As the number of sides of the polygon increases or the contraction rates get closer to 1, the code execution time increases drastically. All experiments showed that the singular set stabilizes, in agreement with our conjecture. Figs 9-11 show some snapshots of our numerical experiments. In these figures, the polygon P is located at the center in white color. The contraction rates λ i are all equal to λ (0, 1). The set K λ surrounds P, and the lines forming the stabilized singular set S n are drawn in black color. The domains of continuity of Tλ n are the polygonal regions cut by S n. The different colors indicate the basins of attraction of different periodic orbits. The basins of attractions of the periodic orbits in the dissipative outer billiard about the pentagon can be seen in Fig. 9: the four colors correspond to the basins of attractions of four distinct attracting periodic orbits. In Fig. 9(B), the blue regions form the basin of attraction of the Fagnano orbit. In the figure, there is also another periodic orbit of period five, but winding around the pentagon twice. This orbit is located on the largest polygonal regions in green color. The basin of attraction of a period orbit of period 10 is plotted in brown color. Finally, the tiny regions in pink color form the basin of attraction of a period orbit of period 35.