arxiv: v1 [math.ds] 8 Nov 2017

Transcription

1 Normal forms for three-parametric families of area-preserving maps near an elliptic fixed point Natalia Gelfreikh arxiv: v1 [math.ds] 8 Nov 017 Department of Higher Mathematics and Mathematical Physics, Saint-Petersburg State University, Russia st00454@spbu.ru Abstract We study dynamics of area-preserving maps in a neighbourhood of an elliptic fixed point. We describe simplified normal forms for a fixed point of co-dimension 3. We also construct normal forms for a generic three-parameter family which contains such degeneracy and use the normal form theory to describe generic bifurcations of periodic orbits in these families. 1 Introduction In a Hamiltonian system small oscillations around a periodic orbit are often described using the normal form theory [1, ] which provide an important tool for the study of local dynamics (see e.g. [1, 3, 4]). In the case of two degrees of freedom the Poincaré section is used to reduce the problem to studying a family of area-preserving maps in a neighbourhood of a fixed point. The Poincaré map depends on the energy level and possibly on other parameters involved in the problem. A sequence of coordinate changes is used to transform the map to a normal form. Our approach to the normal form of a map is similar to [5]. In the absence of resonances the normal form is a rotation of the plane, and the angle of the rotation depends on the amplitude. In a generic oneparameter family of area-preserving maps, the normal form provides a description for a chain of islands which is born from the origin when the multiplier of the fixed point crosses a resonant value [1,, 6, 7, 8]. In [9] unique normal forms for two-parametric families were constructed and used to analyse bifurcations of n-periodic orbits. 1

2 In the paper we study three-parametric families of APM with fixed point of elliptic type. In such families two possibilities of degenerations are possible: (1) degenerations in the twist terms; () degeneration in the leading resonance term. Normal form in the case of (1) was constructed in [9] (for arbitrary number of parameters). Normal forms for () are constructed in the paper first for an individual map (Section ) and then for families (Section 3). In Sections 4 and 5 the normal forms are used to inverstigate biffurcations. 1.1 Individual maps Let F 0 : R R be an area-preserving map (APM) which also preserves orientation. Let the origin be a fixed point: F 0 (0) = 0. Since F 0 is areapreserving det DF 0 (0) = 1. Therefore the two eigenvalues of the Jacobian matrix DF 0 (0) are λ 0 and λ 1 0. We will consider an elliptic fixed point, i.e. the case of non real λ 0. As the map is real λ 1 0 = λ 0. Consequently the multipliers of an elliptic fixed point belong to the unit circle: λ 0 = 1, i.e. λ 0 = e iα 0. There is a linear area-preserving change of variables such that the Jacobian of F 0 takes the form of a rotation: ( ) cos α0 sin α DF 0 (0) = R α0, where R α0 = 0. (1) sin α 0 cos α 0 It is well known that APM with elliptic fixed point can be represented in Birkhoff normal form [1], i.e. there is an area-preserving change of coordinates which transforms F 0 into the resonant normal form N 0 that commutes with the rotation: N 0 R α0 = R α0 N 0. The change of coordinates and the map N 0 are formal series. The linear part of a normal form N 0 is R α0. Following the method suggested in [5] by Takens (see e.g. [10, 8]), we consider a formal series H 0 such that N 0 = R α0 Φ 1 H 0, () where Φ t H 0 is a flow generated by the Hamiltonian H 0. The Hamiltonian has the Takens normal form, i.e. it is invariant with respect to the rotation: H 0 R α0 = H 0. Our goal is to transform formal series H 0 to the most simple form. We use changes of variables which commutate with the rotation R α0. Then in

3 new variables the map is still in Birkhoff normal form and the corresponding Hamiltonian remains in the Takens normal form. It is convinient to use complex variables defined by z = x + iy and z = x iy. (3) A fixed point is called resonant if there exists n N such that λ n 0 = 1. The least positive n is called the order of the resonance [1]. The rotation R α0 takes the form (z, z) (e iα 0 z, e iα0 z). As the map N 0 commutes with R α0 it contains only resonant terms: N 0 (z, z) = λ 0 z + f kl z k z l. k+l, k l=1 (mod n) Corresponding Hamiltonian has the Takens normal form [8]: H 0 (z, z) = h kl z k z l, h kl = h lk. (4) k+l 3, k=l (mod n) It was established in [8], [9] that if h n0 0 then Hamiltonian H 0 can be transformed to a normal form H 0 (z, z) = k a k z k z k + (z n + z n ) k 0 b k z k z k. In the paper we consider the case of h n0 = 0 but h 0 and h n+1,1 4h h n,0 0. In order to investigate Hamiltonian it is convenient to use the symplectic polar coordinates (I, ϕ) given by { x = I cos ϕ, y = I sin ϕ or { z = Ie iϕ, z = Ie iϕ. The Hamiltonian in Takens normal form (4) takes the form : H 0 (I, ϕ) = I a k I k + a jk I jn/+k cos(jnϕ + β jk ), (6) k 0 j 1,k 0 where a k, a jk and β jk are real coefficients. Proposition 5 of Section implies the following theorem. 3 (5)

4 Theorem 1 Let N 0 be BNF for some APM F 0 with resonance of order n at the origin. Let for corresponding Hamiltonian in TNF (6) a 10 = 0, a 0 0 and a 11e iβ 11 8a 0 a 0 e iβ 0 0. Then there is a formal canonical change of coordinates such that in new variables Ñ0 = R α0 Φ 1 H0 and H 0 (I, ϕ) = I k 0 a k I k + I n/ k 1 b k cos(nϕ + ψ k )I k + I n cos nϕ k 0 c k I k. The coefficients in the form (7) are not unique: the coefficients b k are replaced by b k after rotation by π/n. There is also an alternative normal form (only a 0 0 is required), which contains fewer terms of low orders: Theorem Let N 0 be BNF for some APM F 0 with resonance of order n at the origin. Let for corresponding Hamiltonian in TNF (6) a 0 0. Then there is a formal canonical change of coordinates such that in new variables Ñ 0 = R α0 Φ 1 H0 and H 0 (I, ϕ) = I k 0 a k I k + k 1 c k I nk/ cos(knϕ + ψ k ). (7) Corresponding Proposition 6 is proved in Section. 1. Families Let us consider a three-parametric family F µ of APM with fixed point of an elliptic type at the origin with λ µ = e iαµ, µ = (µ 1, µ, µ 3 ) R 3. We assume that for µ = (0, 0, 0) the function F 0 has resonance of order n (i.e. λ 0 = e iα 0, λ n 0 = 1). After a linear change of coordinates the map F µ takes the form F µ = R αµ Φ, where Φ is a tangent-to-identity APM. In complex variables (3) F µ = (f µ, f µ ), f µ (z, z) = λ µ z + f kl (µ)z k z l. It is natural to use the value k+l ε = α 0 α µ (8) 4

5 as one of parameters. Let µ 3 be expressed in terms of (ε, µ 1, µ ). Then f µ (z, z) can be presented as a series in three variables (z, z, ε) with coefficients depending on (µ 1, µ ): f µ (z, z) = λ 0 z + f klm (µ 1, µ )z k z l ε m. k+l+m After an appropriate change of coordinates the map f µ can be written in Birkhoff normal form which contains only resonant terms: f µ (z, z) = λ 0 z + f klm (µ 1, µ )z k z l ε m. (9) k+l+m k l=1 (mod n) Interpolation theorem for families (see for example [8]) gives F µ = R α0 Φ 1 H µ, where Hamiltonian H µ in complex variables has the form H µ (z, z) = h klm (µ 1, µ )z k z l ε m, (10) k+l+m 3 k=l (mod n) where h 111 (µ 1, µ ) = 1 as ε is determined by (8). coordinates (5) Or, in the symplectic H µ (I, ϕ) = εi + I k 0 a km (µ 1, µ )I k ε m m 0 + j 1 a jkm (µ 1, µ )I k+j n ε m cos(jnϕ + ψ jkm (µ 1, µ )). k,m 0 (11) The form (11) is not unique. Rotations and time-one shifts Φ 1 χ with resonant Hamiltonian χ preserve its structure. Our goal is to derive the most simple form of H µ, i.e. to eliminate as many terms of H µ as possible. If a 100 (µ 1, µ ) = h n00 (µ 1, µ ) = 0 for all values of (µ 1, µ ) then Hamiltonian can be transformed to the form [?]: H µ (z, z) = εz z+ a km (µ 1, µ )z k z k ε m +(z n + z n ) b km (µ 1, µ )z k z k ε m k,m 0 k,m 0 k,m 0 or, in coordinates (I, ϕ): H µ (I, ϕ) = εi + a km (µ 1, µ )I k ε m + I n/ cos nϕ b km (µ 1, µ )I k ε m. 5 k,m 0 (1)

6 In this paper the main attention is paid to the case of Hamiltonian (10) when the following conditions are valid: h n00 (0, 0) = 0, h 0 (0, 0) 0, h n+1,1,0(0, 0) 4h 0 (0, 0)h n,0,0 (0, 0) 0. (13) Note that these conditions correspond to a family which is a three-parametric unfolding of F 0 from Theorem 1. We show that there is such canonical change of variables that Hamiltonian has the form (11) with the last sum containing only two harmonics, i.e. terms with j = 1 and j =. Moreover a 1,l+1,m = 0, a,l+1,m = 0 and ψ,l,m = 0 for all l and m. For the sake of receiving ψ 00 (µ 1, µ ) = 0, as it will be shown ( in Section 3, it is necessary to make a rotation by an angle ϕ = 1 arg h n n,0,0 h n+1,1,0 4h 0 ). After the rotation the small coefficient of z n in (10) is h n00 (µ 1, µ ) exp ( i ( )) arg h n,0,0 (µ 1, µ ) h n+1,1,0(µ 1, µ ) = ν(µ 1, µ ), 4h 0 (µ 1, µ ) (14) ν(µ 1, µ ) = ν 1 (µ 1, µ ) + iν (µ 1, µ ) = γ(µ 1, µ )e iβ(µ 1,µ ). Let µ 1 and µ be expressed by ν 1 and ν : µ 1 = µ 1 (ν 1, ν ), µ = µ (ν 1, ν ). Then we can consider (ν 1, ν ) as new parameters instead of (µ 1, µ ). Below we use the following notations: µ = (µ 1, µ, µ 3 ), Υ = (ε, ν 1, ν ), m = (m 1, m, m 3 ), m = m 1 + m + m 3, Υ m = ε m 1 ν m 1 ν m 3 The following theorem gives a simplification of the Hamiltonian to a normal form. Theorem 3 Let F µ be a smooth (C or analytic) family of area preserving maps with fixed point of elliptic type at the origin such that (1) F 0 has a resonant of order n at the origin: λ 0 = e iα 0, λ n = 1; () coefficients of Hamiltonian H µ in the TNF (10) satisfy conditions (13); (3) parameters µ = (µ 1, µ, µ 3 ) can be expressed by (ε, ν 1, ν ) defining by (8) and (14). Then there is a formal Hamiltonian H Υ and formal canonical change of variables which conjugates F µ with R α0 Φ 1 HΥ. Moreover, H Υ in coordinates. 6

7 (I, ϕ) has the following form: H Υ (I, ϕ) = εi + γi n/ cos(nϕ + β) + I +I n/ k, m 0, k+m 1 1 +I n k, m 0 a km I k Υ m b km I k Υ m cos(nϕ + ψ km ) k, m 0 c km I k Υ m cos(nϕ). The theorem follows from Proposition 8 of Section 3. From Proposition 9 follows an alternative variant, useful for study of bifurcations: Theorem 4 If F µ is a smooth (C or analytic) family of area preserving maps such that F 0 has a resonant elliptic fixed point at the origin and in (11) a 00 (0, 0) 0, then there is a formal Hamiltonian H Υ and formal canonical change of variables which conjugates F µ with R α0 Φ 1 HΥ. Moreover, H Υ has the following form: H Υ (I, ϕ) = εi + γi n/ cos(nϕ + β) + I k, m 0 a kmi k Υ m +I n/ k+m 1 1, c km I kn/ Υ m cos((k + )nϕ + ψ km ). m,m 3 0 (15) Then we use the normal forms (1) and (15) to study biffurcations. We keep lower order terms and skip those, which do not change the picture qualitatively. The detailed discussion of typical level sets and their bifurcations for the case of degeneracy in the main resonant term are presented in Section 4. In Section 5 there is a brief description of some possible bifurcations for three-parametric families when in (1) a 0 (0, 0) = a 30 (0, 0) = 0, b 00 (µ 1, µ ) 0. This case is different from two-parametric families considered in [9] only in tiny domain in the space of parameters. Simplification of a formal Hamiltonian with h n,0 = 0 for individual map In this section we construct two degenerate resonant normal forms (17) and (7). The first one provides particularly simple form of Hamiltonian in the 7

8 symplectic polar coordinates: it contains only two lowest harmonics of the angle variable, namely nϕ and nϕ. An alternative normal form (7) contains fewer terms of low orders. Proposition 5 If H(z, z) = k+l 4, k,l 0, k=l (mod n) h kl z k z l (16) is a formal series such that h kl = h lk, h n0 = 0, h 0 and h n+1,1 4h h n,0 0 then there exists a formal canonical change of variables which transforms the Hamiltonian H into H(z, z) = k 0 a k (z z) k+ + k 1 (b k z n + b k z n )(z z) k + k 0 c k (z z) k (z n + z n ), where a k, c k R, b k C. Moreover, a 0 = h and c 0 = h n,0 h n+1,1. 4h (17) Proof. The key idea behind the proof is based on studying terms of formal power series in the order of their δ-degree. For a resonant monomial (z k z l, k = l mod n) we define its δ-degree by δ(z k z l ) = k l n + min{ k, l } = 1 n (k + l) k l. (18) n Grouping in (16) terms of the same δ-degree we get H(z, z) = m h m (z, z), where h m (z, z) is a homogeneous resonant polynomial of δ-degree m. The term of the lowest δ-degree is h (z, z) = h z z + h n+1,1 z n+1 z + h 1,n+1 z z n+1 + h n,0 z n + h 0,n z n. Let χ be a resonant polynomial. After the substitution (z, z) Φ 1 χ(z, z) the Hamiltonian takes the form H = H + L χ H + k 1 k! Lk χh, (19) 8

9 where L χ H = i{h, χ}. The δ-degrees of monomials in Poisson brackets satisfy the following relation: Indeed, {z k1 z l 1, z k z l } = (k 1 l k l 1 )z k 1+k 1 z l 1+l 1 δ(z k 1+k 1 z l 1+l 1 ) δ(z k1 z l 1 ) + δ(z k z l ) 1. (0) δ(z k 1+k 1 z l 1+l 1 ) = 1(k 1 + k + l 1 + l ) n k n 1 l 1 + k l 1(k 1 + l 1 ) n k n 1 l 1 + 1(k + l ) n k n l 1. Let χ = αz n + α z n whith α = h n+1,1 inh. It is not difficult to see that after the substitution (z, z) Φ 1 χ(z, z) the coefficients of Hamiltonian H = k+l 4 h kl z k z l of δ-degree are h = h, h n+1,1 = h n+1,1 + inαh = 0, h n,0 = h n,0 + inαh n+1,1 n α h = h n,0 h n+1,1 4h. After rotation z e i arg h n,0 n z the coefficient h n,0 c 0 = h n,0. Then in new variables h (z, z) = a 0 z z + c 0 (z n + z n ), (1) where a 0 = h, c 0 = h n,0 h n+1,1. 4h So the terms of δ-degree has the declared form. We proceed by induction. Let h m (z, z) has the declared form for m < p, namely p 3 h j+ (z, z) = a j (z z) j+ + c j (z z) j (z n + z n ), for 0 j,() p h j+1 (z, z) = a j 1 (z z) j+1 + (z z) j (b j z n + b j z n ), for 1 j.(3) 9

10 For the sake of receiving such form for the term of δ-degree p we use several changes of variables consequently. Any homogeneous polynomial χ of δ-degree p 1 generates the change of variables (z, z) Φ 1 χ(z, z). Formulae (19) and (0) imply h m = h m for m p 1 and for m = p h p = h p + L(χ). Here L(χ) is the homological operator: L(χ) = [L χ h ] p, (4) where h is determined by (1) and [ ] p denotes terms of δ-degree p. Let p be δ-degree of resonant monomial z k z l and j = k l n. Let and for 1 j p Q p0 (z, z) = z p z p Q pj (z, z) = z p+j(n 1) z p j, Q p, j (z, z) = z p j z p+j(n 1). The homological operator (4) acts on monomial Q p 1,j (z, z) of δ-degree p 1 with j > 0 by and for j = 0: L(Q p 1,j ) = ia 0 njq pj inc 0 (p 1 j)q p,j+ L(Q p 1,0 ) = inc 0 (p 1)(Q p, Q p, ). We denote coefficients of resonant polynomials h p (z, z) and h p (z, z) by g j and g j respectively: h p = g 0 Q p0 + (g j Q pj + gj Q p, j ), hp = g 0 Q p0 + ( g j Q pj + g j Q p, j ). 1 j p 1 j p The change of variables generated by resonant polynomial χ k of δ-degree p 1 χ k = α k Q p 1,k + α k Q p 1, k, 1 k p 1, χ 0 = α 0 Q p 1,0 (5) 10

11 transforms coefficients g j by the following way: g j = g j, j / {k, k + }, g k = g k + ia 0 nkα k, g k+ = g k+ inc 0 (p 1 k)α k. Below we describe the order of this transformations wich provides the declared form of the terms of δ-degree p. Let k 0 = 1 for even p and k 0 = for odd p. Change of variables generated by χ k (5) with α k = g k eliminates g ia 0 kn k, changes g k+ and does not affect any other terms of δ-degree p: g j = g j, j {k, k + }, g k = 0, g k+ = g k+ + c 0(p 1 k) ka 0 g k. Starting with k = k 0, then k = k 0 + etc. proceeding up to k = p 1 one get corresponding terms g k = 0 for all k k 0 for which k k 0 is even. To eliminate terms with odd k k 0 we use changes of variables generated by χ k (5) with g α k = k+, 1 k p, inc 0 (p k+1) α 0 = Ig. nc 0 (p 1) Each of these substitutions eliminates g k+ (for k = 0 only image part of g eliminates), changes g k and does not affect any other terms of δ-degree p: g j = g j, j {k, k + }, g k+ = 0, (for k = 0 I g = 0), g k = g k + ka 0g k+ c 0 (p k+1). Starting from k = p, then k = p 4 etc. down to k = k 0 1 we get corresponding g k+ = 0 and g R. Remark about uniqueness. The kernel of the homological operator is empty for odd p and one-dimensional for even p. So only [h k ] k ker L which implies uniqueness of the coefficients a k, c k, b k and arg b k (mod π) in formula (17). But rotation by π/n change arg b k arg b k + π. An alternative normal form is derived in the next proposition. 11

12 Proposition 6 If H(z, z) = a 0 z z + k+l 5, k,l 0, k=l (mod n) h kl z k z l (6) is a formal series such that h kl = h lk and a 0 0, then there exists a formal tangent-to-identity canonical change of variables which transforms the Hamiltonian H into H(z, z) = (z z) k 0 a k (z z) k + k 1(c k z nk + c k z nk ), (7) where a k R, c k C. Proof. Let χ(z, z) = αz k z l + α z l z k, k > l. After the substitution (z, z) Φ 1 χ(z, z) the Hamiltonian takes the form (19), where L χ H = i{h, χ} = ia 0 α(k l)z k+1 z l+1 ia 0 α (k l)z l+1 z k+1 + O k+l+3, where O k+l+3 denotes terms of degree k + l + 3 and higher. Then h st = h st for 4 s + t k + l + 1 and for s + t = k + l + if (s, t) (k + 1, l + 1). And h k+1,l+1 = h k+1,l+1 + ia 0 α(k l). h k+1,l+1 Let α = ia 0 (k l). Then h k+1,l+1 = 0. Repeating this substitutions one gets by induction h kl = 0 with the exception of h kk and h nk,0 = h 0,nk. 3 Simplification of a formal Hamiltonian for families with small h 11 and h n0 Now we consider the three-parametric unfolding of Hamiltonian (16): H(z, z; ε, µ 1, µ ) = εz z + h klm (µ 1, µ )z k z l ε m, (8) k+l 3,m 0 k=l (mod n) 1

13 where h klm (µ 1, µ ) = h lkm (µ 1, µ ), h n00 (0, 0) = 0, h 0 (µ 1, µ ) 0, h n+1,1,0(µ 1, µ ) 4h 0 (µ 1, µ )h n,0,0 (µ 1, µ ) 0. The purpose of the section is to prove that after some canonical formal transformations Hamiltonian (8) takes the form (34). An alternative form (37) is also derived. As in Section δ-degree is defined by (18). Then H(z, z; ε, µ 1, µ ) = h (p)m (z, z; µ 1, µ )ε m, (9) p 1,m 0 where h (p)m is homogeneous resonant polynomial on (z, z) of δ-degree p with coefficients depending on (µ 1, µ ). The terms of δ-degree 1 are h (1)1 (z, z; µ 1, µ ) = z z + h n01 (µ 1, µ )z n + h n01(µ 1, µ ) z n, h (1)m (z, z; µ 1, µ ) = h n0m (µ 1, µ )z n + h n0m(µ 1, µ ) z n, m 1. (30) The all transformations described below preserve the form of h (1)m (i.e. h 110 = 1 and h 11m = 0 for m 1) although coefficients h n0m are changed. As for individual map, normalisation of terms of δ-degree p = is different from p 3. So it is convenient to consider the case of p = separately. Lemma 7 Let H(z, z; ε, µ 1, µ ) be a formal series (9) with h (1)m (z, z; µ 1, µ ) as in (30) and h ()m (z, z; µ 1, µ ) = h m (µ 1, µ )z z + h n+1,1,m (µ 1, µ )z n+1 z +h n+1,1,m(µ 1, µ )z z n+1 + h n,0,m (µ 1, µ )z n + h n,0,m(µ 1, µ ) z n. If h 0 (µ 1, µ ) 0 and h n+1,1,0(µ 1, µ ) 4h 0 (µ 1, µ )h n,0,0 (µ 1, µ ) 0 then there exists a formal canonical change of variables which transforms the Hamiltonian H into H = h p 1,m 0 (p)m ε m with h n00 (µ 1, µ ) = h n00 (µ 1, µ ) exp ( i ( arg h n,0,0 (µ 1, µ ) h n+1,1,0(µ 1, µ ) 4h 0 (µ 1, µ ) (31) h 11m (µ 1, µ ) = h 11m (µ 1, µ ) = 0 for m 1, h 111 (µ 1, µ ) = h 111 (µ 1, µ ) = 1, )), h ()m (z, z; µ 1, µ ) = h m (µ 1, µ )z z + h n,0,m (µ 1, µ )(z n + z n ), (3) where h n,0,m (µ 1, µ ) R for all m 0, h n,0,0 (µ 1, µ ) = h n,0,0 (µ 1, µ ) h n+1,1,0 (µ 1,µ ). 4h 0 (µ 1,µ ) 13

14 Proof. The lemma is proved by induction. For m = 0 we use change of variables (z, z) Φ 1 χ(z, z) with χ = βz n + β z n, β = h n+1,1,0 inh 0. From (19) we obtain coefficients of terms of δ-degree 1: and of δ-degree, m = 0 : h n0m = h n0m for m = 0 and m, h n01 = h n01 + inβ = h n01 h n+1,1,0 h 0 h 0 = h 0 for n 4, h 0 = h 0 in β h 300 for n = 3, h n+1,1,0 = h n+1,1,0 + inβh 0 = 0, h n,0,0 = h n,0,0 + inβh n+1,1,0 n β h 0 = h n,0,0 h n+1,1,0 4h 0 0. Note that h 0 is changed only for the case of resonance of the order n = 3 but h 0 (0, 0) = h 0 (0, 0) as h( 300 (0, 0) = 0. ) After rotation z z exp the coefficient h n,0,0 h n,0,0. i arg h n,0,0 n Thus h ()0 has the form (3). Let h ()m has the form (3) for all m M 1. In order to normalise the term h ()M we use successively two substitutions of the type (z, z) Φ εm χ (z, z). Then H = H + ε M L χ H + ε km k! Lk χh. k First we assume χ = βz n + β z n and then χ = αz z. Choosing appropriate value of β we eliminate h n+1,1,m and we get Ih n,0,m = 0 by choosing α. After the first substitution with β = h n+1,1,m inh 0 the coefficients of terms of δ-degree 1 are: h n0m = h n0m for m M + 1, h n,0,m+1 = h n,0,m+1 + inβ = h n,0,m+1 h n+1,1,m h 0 and for terms of δ-degree for m M 1 the coefficients are not changed and for m = M: h M = h M for n 4, h M = h M in β h 300 for n = 3, h n+1,1,m = h n+1,1,m + inβh 0 = 0, h n,0,m = h n,0,m + inβh n+1,1,0 = h n,0,m h n+1,1,m h n+1,1,0 h 0. 14

15 The change (z, z) Φ εm χ (z, z) with χ = αz z leads to h klm = h klm for m M 1 and h klm = h klm iα(k l)h kl0 for k > l. So, h n,0,m = h n,0,m inαh n00 and h M = h M, h n+1,1,m = h n+1,1,m = 0, h n,0,m = h n,0,m inαh n,0,0. Choosing α = Ih n,0,m nh n,0,0 we get h n,0,m R. Now we introduce new parameters (ν 1, ν ) instead of (µ 1, µ ). Let h n00 (µ 1, µ ) = ν = ν 1 + iν, where h n00 (µ 1, µ ) is determined by the formula (31). Let µ 1 and µ be expressed in terms of (ν 1, ν ) as power series. Then we get Hamiltonian H in the form of the series in five variables (z, z; ε, ν 1, ν ). Let Υ = (ε, ν 1, ν ), m = (m 1, m, m 3 ), m = m 1 + m + m 3, Υ m = ε m 1 ν m 1 ν m 3. After collecting terms of the same δ-degree Hamiltonian H takes the form: H = εz z + νz n + ν z n + (b 0m z n + b 0m z n )Υ m m 1 1 m,m az z + c(z n + z n ) + (a 0m z z + c 0m (z n + z n ))Υ m (33) + m 1 p 3, m 0 h (p)m (z, z)υ m. Here the terms of δ-degree 1 are in the form (30) and the terms of δ-degree are already in the form (3), a = h 0 (0, 0), c = h n,0,0 (0, 0) h n+1,1,0 (0,0) 4h 0 (0,0). The next proposition completes transformation of Hamiltonian to the normal form. 15

16 Proposition 8 Let (33) be a formal series where h (p)m are real-valued resonant polynomials on (z, z) of δ-degree p and ac 0. There exists a formal canonical change of variables which transforms the Hamiltonian H into H(z, z; Υ) = εz z+νz n +ν z n +a(z z) +c(z n + z n )+(z z) a km (z z) k Υ m + k+m 1 1 m,m 3 0 (z z) k (b km z n + b km z n )Υ m + (z n + z n ) where a km, c km R, b km C. k+ m 1 k+ m 1 (34) c km (z z) k Υ m, Proof. Let δ-degree of monomial z k z l is defined by (18) as before. Now we introduce a new δ Υ -degree for monomial on five variables: δ Υ (z k z l Υ m ) = δ(z k z l ) + m. (35) Then H(z, z; Υ) = s h s (z, z; Υ), where h s is a homogeneous polynomial of δ Υ -degree s. The main term of δ Υ -degree s = is already in the normal form: h (z, z; Υ) = a(z z) + c(z n + z n ). We proceed by induction. Let S 3 and h s (z, z; Υ) has the declared form for s S 1. The term of δ Υ -degree S has the form: h S (z, z; Υ) = h (S M )M(z, z)υ M, (36) 0 M S 1 where h (S M )M (z, z) is a polynomial of δ-degree S M. The terms of δ- degree 1 and are already in the normal form. So we consider S M 3, i.e. 0 M S 3. Let χ be a homogeneous polynomial on (z, z) of δ-degree P 1 = S M 1. After the substitution (z, z) Φ ΥM χ (z, z) the Hamiltonian takes the form H = H + Υ M L χ H + Υ Mj L j j! χh. j 16

17 Then and h s (z, z; Υ) = h s (z, z; Υ) for s S 1 h S (z, z; Υ) = h S (z, z; Υ) + Υ M L(χ), where L is the homological operator (4) from Section. So h (P )m = h (P )m for m M and h (P )M = h (P )M + L(χ). It was been shown in Section that for odd P there exists such χ that and for even P h (P )M = a P M z P z P + (z z) P 1 (b P M z n + b P M z n ) h (P )M = a P M z P z P + c P M (z z) P (z n + z n ), where a P M, c P M R, b P M C. So, taking sequentially all sets of (m 1, m, m 3 ) such that 0 m 1 + m + m 3 S 3 we transform all terms of δ Υ-degree S to the declared form. The following proposition establishes the normal form of Hamiltonian useful for study of bifurcations. Proposition 9 Let (33) be a formal series where a 0 and h (p)m are realvalued resonant polynomials of δ-degree p. There exists a formal canonical change of variables which transforms the Hamiltonian H into H(z, z; Υ) = εz z + a(z z) + (z z) k+ m 1 a km(z z) k Υ m +νz n + ν z n + k 1,k+m 1, m 0 (c kmz kn + c km zkn )Υ m (37) where a km R, c km C, c 0 = c. Proof. We use δ Υ -degree introduced by δ Υ (z k z l Υ m ) = k + l + m and normalized order by order as in Proposition 8 but now we choose the compliment to the image of the homological operator as in the proof of Proposition 6. 17

18 4 Bifurcations for the case of degeneracy in the leading resonant term We discuss typical level sets of Hamiltonian (37) for 0 < ε + ν ɛ 0 with sufficiently small ɛ 0. In order to investigate biffurcations we keep lower order terms and skip those, which do not change the picture qualitatively: H(z, z; ε, ν) = εz z + a(z z) + νz n + ν z n + c(z n + z n ), where ν = ν + εc ε k c 1k00, k = n. Applying the symplectic polar coordinates (5) and assuming ν = 1 γeiβ we get the model Hamiltonian in the form H(I, ϕ; ε, γ, β) = εi + I + γi n/ cos(nϕ + β) + I n cos nϕ. (38) Here the coefficients a and c are normalized to unity with the help of a scaling applied to the variable I, the parameters ε and γ, and the Hamiltonian function H. The case of n = 3 is different from the case of n Typical level sets of the model Hamiltonian for n 4 Since I H = ε + I + O(γI n/ 1 ) + O(I n 1 ), (39) when ε > 0 for n 4 Hamiltonian H has not any critical points near the origin and level sets are closed curves. Typical level sets for ε < 0 are depending on value of parameters ε and γ. We consider separatly the case of γ ε n/ and γ = O( ε n/ ). If ε n/ = o(γ) then the last term in (38) can be omitted. Then critical points of Hamiltonian are located near I = ε, cos(nϕ + β) = ±1. There is a chain of n islands on the distance of order ε 1/ from the origin. The size of islands is of order ( γ 1/ ε n/4. For γ = O ε n/) the last term in (38) is also essential. In order to study bifurcations we introduce the scaling: ( γ = 4 ε ) n/b, ε I = ( + ε ) n/j, ( ψ = nϕ, H = ε n H. (40) ) 18

19 Figure 1: Bifurcation diagram on the complex plane of ν = γe iβ for ε < 0. The domains D and D are separated by the vertical line segment. After this scaling the Hamiltonian takes the form (scipping the constant term) H(J, ψ) = H ( ( ) ) ε n/ 1 0 (J, ψ) + O, where H 0 (J, ψ) = J + 4b cos(ψ + β) + cos ψ. (41) Equilibrium points of H 0 (J, ψ) locate at the points for which J = 0 and b sin(ψ + β) + sin ψ cos ψ = 0. (4) To figure out how many solutions has this equation we denote sin ψ = x and cos ψ = y. Then equation (4) is equvalent to the system { xb cos β + yb sin β + xy = 0 x + y = 1 First equation corresponds to hyperbole (or straight lines if cos β = 0, or sin β = 0, or b = 0). One branch of the hyperbole (or one of stright lines) passes through the origin. So at least two solutions exist for arbitrary b and β. There are two more solutions if second branch of hyperbole (or second straight line) crosses the unit circle. Tangency condition for the unit circle and the second brunch of the hyperbole is (b cos β) /3 + (b sin β) /3 = 1. 19

20 (a) (b) (c) (d) Figure : Critical level sets of the model Hamiltonian for ε < 0 depending on value of ν indicated on bifurcation diagram on Fig. 1: (a) ν D, (b) ν D, β = 0, (c) ν is in the boundary between D and D (β = ±π/), (d) ν is in the boundary between D and D 1. On the complex plane of γe iβ = ν = ν 1 + iν corresponding line is given by astroid: ν 1 /3 + ν /3 = 4/3 ( ε ) n/3 (43) and is presented on Fig. 1. Typical critical level sets of Hamiltonian corresponding different values of parameter ν (provided ε < 0) are shown on Fig.. Fragments on the figure are repeated n times. On the boundary between D and D (vertical line segment) there is a chain of n islands (Fig. (c)). Outside a domain bounded by the astroid there is a chain of n islands. 4. Typical level sets of model Hamiltonian for n = 3 For ε < 0 the case of n = 3 essentially differs from the case of n 4 as the third term in (39) cannot be omitted. The first three terms in (39) are of the same order if I and ε are of the order of γ. So we use the following scaling in (38) and get H = H 0 + O(γ ), where ε = aγ, I = γ J, H = γ 4 H (44) H 0 = aj + J 3/ cos(3ϕ + β) + J. (45) 0

21 Figure 3: Bifurcation diagram on the plane (γ, ε) for n = (a) (b) (c) (d) Figure 4: Critical level sets of the model Hamiltonian depending on value of (γ, ε): (a) (γ, ε) D 1, (b) the boudary between D 1 and D 1 (ε = 0), (c) (γ, ε) D 1, (d) the boudary between D 1 and D 0 (ε = 9 3 γ ). The critical points of H 0 are located at the points for which cos(3ϕ+β) = σ ϕ = ±1, and a + 3σ ϕ J 1/ + J = 0. The last equation can be solved explicitly. There are no real solutions for a > a 0, where a 0 = 9. If a (0; a 3 0) there are two solutions for σ ϕ = 1 and no one for σ ϕ = 1. And if a < 0 there is one solution for σ ϕ = 1 and one for σ ϕ = 1. The line ε = 9 3 γ. on the plane (γ, ε) bounded the domain in which level sets of the Hamiltonian (45) are closed curves (the domain D 0 on Fig. 3). The critical level sets of the Hamiltonian H 0 are illustrated on Figure 4. 1

22 For ε < 0 the case of n = 3 is not different from n 4 (see the analysis in the previous subsection). In particular, the boundary between D 1 and D is given by (43). On the Fig. 3 we assume ν 1 = ν = γ, ν = 0 which correspond to β = 0 (ν = γe iβ ). Then the boundary between D 1 and D is given by ε = 1/3 γ /3. For (γ, ε) D level sets of the model Hamiltonian are presented on Fig.. 5 Families with twist degeneracy The theorem about normal forms for families with twist degeneracy can be found in [9]. For three-parametric families it takes the following form. Proposition 10 Let Υ = (ε 1, ε, ε 3 ), m = (m 1, m, m 3 ), Υ m = ε m 1 1 ε m ε m 3 3. If H(z, z; Υ) = ε 1 z z + ε (z z) + ε 3 (z z) 3 + h klm z k z l Υ m (46) k+l>, k,l, m 0, k=l (mod n) is a formal series such that h klm = h lkm, h m = h 33m = 0 and h 440 h n00 0, then there exists a formal tangent-to-identity canonical change of variables which transforms the Hamiltonian H into H(z, z; Υ) = ε 1 z z + ε (z z) + ε 3 (z z) 3 + (z z) 4 a km (z z) k Υ m where a km, b km R, b 00 = h n00. Moreother if 3 n 8 k, m 0 +(z n + z n ) k, m 0 b km (z z) k Υ m, a km = 0 for k = n 5 (mod n) and b km = 0 for k = n 1 (mod n) if n 8 b km = 0 for k = 3 (mod 4)

23 (a) (b) Figure 5: Bifurcation diagram for stable case on the unit sphere in the space of parameters (ε 1, ε, ε 3 ). The coefficients a km and b km are defined uniquely. In the symplectic polar coordinates (5) the model Hamiltonian is H(I, ϕ; Υ) = ε 1 I + ε I + ε 3 I 3 + I 4 + b 0 I n/ cos nϕ. (47) Its typical level sets depend on which of two last terms is the main. If the term I 4 is smaller then the last term (n 7 or n = 8 and b 0 > 1) then for Υ = 0 the origin is not stable. If n 9 or n = 8 and b 0 < 1 then the origin is stable for Υ = 0. Bifurcations in stable and unstable cases are briefly considered below. 5.1 Stable case Let n 9. For 0 I 1, the Hamiltonian (47) can be considered as a small perturbation of H 0 = ε 1 I + ε I + ε 3 I 3 + I 4. The level sets of H 0 are circles for all values of the parameters. The Hamiltonian (47) does not posses this symmetry. The implicit function theorem implies that the last term can affect on the shape of level sets of H only near zeroes of I H 0. The equation I H 0 = ε 1 + ε I + 3ε 3 I + 4I 3 = 0 (48) has from 0 to 3 solutions depending on (ε 1, ε, ε 3 ). Let D k be a domain in the parameter space such that equation (48) has k roots, i.e. H 0 has kn stationary points. These domains on the unit sphere in the space (ε 1, ε, ε 3 ) are shown on the Figure 5. 3

24 (a) (b) (c) (d) Figure 6: The typical critical level sets for the Hamiltonian (47) in the stable case if parameters (ε 1, ε, ε 3 ) D 3. (a) (b) Figure 7: Bifurcation diagram for unstable case on the unit sphere in the space of parameters (ε 1, ε, ε 3 ). In D 0 all level sets of H are closed curves around the origin. In D 1 there is one chain of islands. In D typical level sets are similar to ones in two-parametric families [9]. In neighbourhood of the point ε 1 = 0, ε = 0, ε 3 = 1 there is a tiny domain D 3. When parameters are in D 3 the equation (48) have 3 roots. They correspond to three sets of hyperbolic and three sets of elliptic stationary points of the model Hamiltonian (47). Some of the possible corresponding critical level sets are shown on Figure 6. The fragment shown on the figure is repeated n times around the origin. 4

25 (a) (b) (c) Figure 8: Some critical level sets for n = 8, b 0 =. 5. Unstable case For n 7 or n = 8 and b 0 > 1 the bifurcation diagram on the unit sphere in the space of parameters (ε 1, ε, ε 3 ) is shown in Figure 7. In D 1 and D 1 typical level sets are the same as in one-parametric families. In D and D typical level sets and their bifurcations are essentially the same as in two-parametric families. In a neighbourhood of the points ε 1 = ε = 0, ε 3 = ±1 there are tiny domains with additional sets of stationary points. Some possible critical level sets for this case are shown on Figure 8. References [1] Arnold, V.I., Kozlov, V.V., Neishtadt, A.I., Mathematical aspects of classical and celestial mechanics. Dynamical systems. III. Third edition. Encyclopaedia of Mathematical Sciences, 3. Springer-Verlag, Berlin, 006. [] Siegel, C., Moser, J. Lectures on celestial mechanics, Springer, 1971, xii+90 pp. [3] Bryuno, A.D. The normal form of a Hamiltonian system. Russian Math. Surveys 43 (1988), no. 1, [4] Kuznetsov, Y.A., Elements of applied bifurcation theory. Third edition. Applied Mathematical Sciences, 11. Springer-Verlag, New York, pp. 5

26 [5] Takens, F., Forced Oscillations and Bifurcations, Global Analysis of Dynamical Systems, Festschrift dedicated to Floris Takens for his 60th birthday, Bristol: Inst. Physics, 001, pp (a reprint of Communications of the Mathematical Institute Rijksuniversiteit Utrecht, 3 (1974)). [6] Meyer, K.R. Generic bifurcation of periodic points. Trans. Amer. Math. Soc. 149 (1970) [7] Simó, C., Vieiro, A. Resonant zones, inner and outer splittings in generic and low order resonances of area preserving maps, Nonlinearity no.5 (009) [8] Gelfreich, V., Gelfreikh, N., Unique resonant normal forms for area preserving maps at an elliptic fixed point, Nonlinearity no.4 (009) [9] Gelfreich, V., Gelfreikh, N., Unique normal forms near a degenerate elliptic fixed point in two-parametric families of area-preserving maps. Nonlinearity 7 no.7 (014) [10] Bazzani, A., Giovannozzi, M., Servizi, G., Todesco, E., Turchetti, G. Resonant normal forms, interpolating Hamiltonians and stability analysis of area preserving maps. Phys. D 64 (1993), no. 1-3,