Introduction
A very important qualitative feature of Hamiltonian systems is integrability [8]. According to the Liouville–Arnold theorem ([1], p. 270–272), if an autonomous Hamiltonian system with n degrees of freedom possesses n functionally independent first integrals in involution, then every regular, compact, and connected common level set of these integrals is diffeomorphic to an invariant n-dimensional torus Tn. Moreover, in a neighbourhood of such a torus, action–angle coordinates can be introduced.
Let us begin with a perturbed, autonomous, two-degree-of-freedom system given by a Hamiltonian of the general form. Given that the perturbed Hamiltonian is an analytic function in terms of the parameter ε, on an open interval around ε=0, it can consequently be expressed as a power series in ε ([2], p. 40),
H(J1,J2,w1,w2)=H0(J1,J2)+εH1(J1,J2,w1,w2)(1)
The perturbed Hamiltonian is written in curvilinear symplectic coordinates called action-angle variables in an open region of the phase space:
Ji∈S⊂Rn,wi mod2π∈Tn
More specifically, H1(J1,J2,w1,w2) is the perturbation, which is periodic in w1,w2 and can be expanded as a double Fourier series with period 2π. The function H0(J1,J2) is the Hamiltonian of the unperturbed problem. It is Liouville–Arnold integrable, with finite motions; that is, in the corresponding regular open region of phase space, the unperturbed system is locally foliated by invariant n-dimensional tori. These tori are invariant with respect to the phase flow of the unperturbed system. Furthermore, H0 satisfies the nondegeneracy condition ([2], p. 40):
det(∂Ji∂Jj∂2H0)=0,i,j=1,2
The motion of a phase point is defined by the solution of the system
J˙i=0,w˙i=ωi(J),ωi(J)=∂Ji∂H0,i=1,2(2)
That is,
Ji=Ji0,wi=ωit+wi0,i=1,2(3)
Each torus is determined by the value Ji=Ji0, while wimod2π are the angular coordinates .
In an abstract n-dimensional system, a torus Tn with frequencies
ω=(ω1,…,ωn)
is called non-resonant if k⋅ω=0 ∀k∈Zn∖{0},
and resonant if there exists a nonzero integer vector
k∈Zn∖{0}
such that
k⋅ω=0
Obviously, in a two-dimensional system, a torus T2 is said to be resonant if
∃(k1,k2)∈Z2∖{(0,0)}:k1ω1+k2ω2=0
For k1=0 and ω2=0, this gives
ω2ω1=−k1k2∈Q
Generally, the nondegenerate condition corresponds to the case in which the frequencies are functionally independent. In nondegenerate systems the non-resonant tori form an everywhere dense set of full measure. The resonant tori form a set of measure zero, which is also everywhere dense ([3], p. 273).
A phase trajectory on a non-resonant torus covers the torus densely. One can simply say that when the frequency ratio
ω2ω1∈Q,
the invariant tori contain periodic orbits that close after a finite time, whereas when the frequency ratio is irrational,
ω2ω1∈R∖Q,
the invariant tori have quasi-periodic orbits that cover the tori densely.
Poincaré's theorem on nonintegrability
Obviously,
H(J1,J2,w1,w2)
is an integral of motion. Let us suppose that H has a second integral of motion
Φ(J1,J2,w1,w2,ε)
in an open neighborhood of J1,J2, for w1,w2 from 0 to 2π, and for ε in a neighborhood of ε=0. Then Φ can be expanded in a convergent series around ε=0:
Φ=Φ0(J1,J2)+εΦ1(J1,J2,w1,w2)+ε2Φ2+O(ε3)(4)
If we take a non-resonant torus, note that Φ0 is independent of the angles. A detailed explanation can be found in ([3] p. 353).
Poincaré proves that Φ0 is an integral of H0, independent of it. Furthermore, up to order ε, since Φ is an integral of H, it holds that
That leads to
[Φ0,H0]=0,[Φ0,H1]+[Φ1,H0]=0(5)
The first relation in (5) states that Φ0 is an integral of H0. We shall examine whether Φ0 is independent of H0. The second relation in (5) yields the relation between the Fourier coefficients Bk1,k2 and Ak1,k2:
Bk1,k2(k1∂J1∂H0+k2∂J2∂H0)=Ak1,k2(k1∂J1∂Φ0+k2∂J2∂Φ0)(6)
Here,
H1=k1,k2∑Ak1,k2(J1,J2)exp(i(k1w1+k2w2))
and
Φ1=k1,k2∑Bk1,k2(J1,J2)exp(i(k1w1+k2w2))
are the Fourier series of H1 and Φ1 in w1 and w2.
Relation (6) is an identity on the phase space; that is, it holds for every Ji,wi in an open region d×Tn, at least ([2] p. 47).
According to equation (6), its left-hand side becomes zero on the dense set of resonant tori of the unperturbed system H0, since it is nondegenerate:
∂J2∂H0∂J1∂H0=ω2ω1=−k1k2(7)
If an infinite set of Fourier coefficients Ak1,k2=0 of H1 are different from zero on these dense resonant tori, then according to equation (6),
∂J1∂Φ0/∂J2∂Φ0=−k1k2(8)
If we combine (7) and (8), we get
D(J)=∂(J1,J2)∂(Φ0,H0)=0(9)
Therefore, Φ0 is functionally dependent on H0, which is a contradiction. We point out that (9) holds for the entire domain of the phase space. Since D is continuous, if it is zero on a dense set, then it is zero everywhere. Indeed, take any J∗ in the domain. Since the resonant set of tori is dense, there exists a sequence
J(n)→J∗
such that
D(J(n))=0.
Therefore,
D(J∗)=n→∞limD(J(n))=0,
so (9) holds throughout the entire domain. Thus, Poincaré’s theorem [7] states:
If an infinite set of Fourier coefficients Ak1,k2=0 of H1 are nonzero on a dense set of resonant tori of the unperturbed, integrable, nondegenerate system, then H does not possess a second integral of motion analytic in the real variables in an open neighborhood of J1,J2, for w1,w2 from 0 to 2π, and for ε in a neighborhood of ε=0.
Meletlidou–Ichtiaroglou nonintegrability criterion and equivalence to Poincaré's theorem on nonintegrability
A more detailed presentation is given in [4]. From the second relation in (5), we obtain
[Φ0,H1]=[H0,Φ1],
and therefore
dtdΦ1=[Φ0,H1](10)
Consider a member of the family of resonant tori for which
ω2ω1=mn,n,m∈Z∖{0}.
By integrating (10) over one period, we obtain the following. The right-hand side of this equation gives zero upon integration, since the corresponding term is the time derivative of Φ1 integrated over one period T on the resonant torus.
The period T is determined by
T=Ω2π,Ω=nω1=mω2
The left-hand side results in
(k2∂J1∂Φ0+k1∂J2∂Φ0)∂φ∂∫0TH1dt=0(11)
where
φ=k2w1+k1w2
We note that every orbit on a resonant torus is periodic and differs from the others only in terms of w10, w20. Thus, we obtain an entire family of periodic orbits on a specific resonant torus.
Along the unperturbed motion,
dtdφ=k2w˙1+k1w˙2=0.
Therefore, φ is constant on every periodic orbit. On a resonant torus of H0 there are infinitely many periodic orbits parametrized by φ.
The mean value of H1 is
⟨H1⟩(φ)=T1∫0TH1(t;φ)dt
It can either be the same for all periodic orbits of the torus, or it can vary as φ changes. In the second case, there exists a trajectory for which
∂φ∂∫0TH1dt=0
On a resonant torus, this leads to
∂J2∂Φ0∂J1∂Φ0=k2k1
This relation concerns the entire resonant torus. It means that Φ0 depends on H0 on this torus. This occurs for a dense set of resonant tori if
∂φ∂∫0TH1dt=0
on this dense set. But this is not possible according to Poincaré’s theorem, since Φ0 is independent of H0. Therefore, H is nonintegrable.
More specifically, the exact proposition of [4] is that if there exists a dense set of resonant tori of H0 such that, on each of them, the average value of H1 depends on the particular periodic orbit, then the perturbed Hamiltonian is nonintegrable for ε∈(−ε0,ε0)∖{0} in an open neighborhood of ε=0.
It is proven that the Meletlidou–Ichtiaroglou criterion and Poincaré’s theorem are equivalent.
A new theorem on nonintegrability in the case of two independent resonances
We are going to establish a new theorem on the nonintegrability of a perturbed Hamiltonian system with two resonances. We begin with the Hamiltonian
H=H0(J1,J2)+ε[A1(J1,J2)cos(k1w1−k2w2)+A2(J1,J2)cos(l1w1−l2w2)](12)
Here J1,J2,w1,w2 are the action-angle variables, H0 is the integrable, nondegenerate, unperturbed part in action-angle variables, and H1 is the perturbation with two resonances, k1,k2,l1,l2∈Z. and arbitrary. The vectors
k=(k1,−k2),l=(l1,−l2)
are independent.
We suppose, as in Poincaré’s theorem of nonintegrability, that there exists a second integral of motion independent of H, and we are going to use Poincaré’s lemma that Φ0 is independent of H0. To study the nonintegrability generated by the interaction of these two resonances, we apply a sequence of Lie-series symplectic transformations ([6]).
The Lie operator is defined by the Poisson bracket:
LWF=[H,F].
Our Lie generators Wmk,nk is a sequence of function in angle variables
Wmk,nk=cos(mkw1−nkw2),(13)
where
mk=pkk1+qkl1,nk=jkk2+skl2,
with
pk,qk,jk,sk,k1,k2,l1,l2∈Z−{0},k∈N−{0}i=1,2,
Then the transformed Hamiltonian is
H′=exp(εLWmk,nk)H=n=0∑∞n!εnLWmk,nknH==H0+εH1+ε[Wmk,nk,H0]+ε2[Wmk,nk,H1]+2ε2[Wmk,nk,[Wmk,nk,H0]]+O(ε3).(14)
According to Cauchy’s theorem for Lie series, the series for H′ converges, since it does not have any small divisors[5].
At each step we apply Wmk,nk to the initial Hamiltonian H. At each step we shall prove that eq. (16) is non-zero at the resonant torus with ω2ω1=k1+mkk2+nk. After applying this infinite sequence of Lie generators, we shall prove that (16) is non-zero at the resonant tori ω2ω1=k1+mkk2+nk where pk,qk,jk,sk take infinite values and ω2ω1 at these resonant tori form a dense set in action space.
We also apply the Poisson bracket with Wmk,nk to equation (5):
[Wmk,nk,[Φ0,H1]]+[Wmk,nk,[Φ1,H0]]=0
Using again Jacobi’s identity for
[Wmk,nk,[Φ0,H1]]=[Φ0,[H1,Wmk,nk]]+[[Wmk,nk,Φ0],H1]
and
[Wmk,nk,[Φ1,H0]]=[Φ1,[Wmk,nk,H0]]+[[Wmk,nk,Φ1],H0],
we obtain
[Φ0,[Wmk,nk,H1]]+[[Wmk,nk,Φ0],H1]+[Φ1,[Wmk,nk,H0]]+[[Wmk,nk,Φ1],H0]=0(15)
According to the criterion by E. Meletlidou and S. Ichtiaroglou, we integrate (16) on the periodic orbits of the unperturbed invariant resonant tori of H0, where the resonance condition holds:
k2+nkω1=k1+mkω2
And also
l2+nkω1=l1+mkω2
The results of the analysis using (α) also hold for (β). Let us consider the resonance condition (α). The period T is
T=Ω2π,Ω=k2(1+jk)+skl2ω1=k1(1+pk)+qkl1ω2
We integrate over one period:
∫0T([Wmk,nk,[Φ0,H1]]+[Wmk,nk,[Φ1,H0]])dt
Thus,
∫0T([Φ0,[Wmk,nk,H1]]+[[Wmk,nk,Φ0],H1]+[Φ1,[Wmk,nk,H0]]+[[Wmk,nk,Φ1],H0])dt=0(16)
Consider
w1=ω1t+w10,w2=ω2t+w20,
with frequencies of the form
ω1=[k2(1+jk)+skl2]Ω,
ω2=[k1(1+pk)+qkl1]Ω
For the first term of eq. (16) which is
∫0T[Φ0,[W,H1]]dt
and selecting the term cos((k1+m1)ω1−(k2+m2)ω2) as a general term, we obtain
21(∂J2∂A1nk−∂J1∂A1mk)(∂J2∂Φ0(k2+nk)−∂J1∂Φ0(k1+mk))⋅⋅∫0Tcos((k1+mk)ω1−(k2+nk)ω2)dt.
Provided that Φ0 is independent of H0, the term
(k2+nk)∂J2∂Φ0−(k1+mk)∂J1∂Φ0
is different from zero
Also the other coefficient is nonzero if
mknk=∂J2∂A1∂J1∂A1
as well as
The integral is nonzero since, on the above resonance (α),
∫0Tcos[(k1+mk)ω1−(k2+nk)ω2]dt=∫0Tcosϕkdt=Tcosϕk.(17)
where ϕk=(k1+mk)ω10−(k2+nk)ω20 since the resonance does not vanish the initial conditions ωi0.
Looking at the last term of (17), it can be proven that it is zero, since it is the time derivative of [Wmk,nk,Φ1]:
∫0T[[Wmk,nk,Φ1],H0]dt=∫0T(∂ω1∂[Wmk,nk,Φ1]∂J1∂H0+∂ω2∂[Wmk,nk,Φ1]∂J2∂H0)dt==∫0T[∂ω1∂[Wmk,nk,Φ1]ω˙1+∂ω2∂[Wmk,nk,Φ1]ω˙2]dt.
That is,
∫0Tdtd[Wmk,nk,Φ1]dt=0
The other terms are
∫0T([[Wmk,nk,Φ0],H1]+[Φ1,[Wmk,nk,H0]])dt(18)
The second term is
∫0T[[Wmk,nk,Φ0],H1]dt
and has the coefficient
21[∂J1∂A1mk(∂J2∂Φ0nk−∂J1∂Φ0mk)+∂J2∂A1nk(∂J1∂Φ0mk−∂J2∂Φ0nk)+A1(−∂J1∂J2∂2Φ0k1k2(pk+jk)+k1l2sk+k2l1qk+k1∂J12∂2Φ0mk+k2∂J22∂2Φ0nk)]
It can be proven that this it is nonzero when H0 does not have a linear integral in J1,J2, namely Φ0, that is continued in the perturbed system. Also, the integral (18) and the last term in (19) at the resonance is nonzero.
Also, the terms do not balance.
We shall prove in Appendix A. that the ratios
ω2ω1=k1+mkk2+nk
and
ω2ω1=l1+mkl2+nk
form a dense set in the relevant interval. Therefore, the resonant tori form a dense set in the phase space. After proving that the resonant tori are dense in the action space, the Hamiltonian given by (12) is nonintegrable
For H to possess a second independent first integral analytic in ε, equation (17) must vanish. Hence, the vanishing of (17) is a necessary condition for the existence of such an integral. Under the assumptions stated below, equation (17) is non-zero on a dense set of resonant tori, which yields a contradiction.
Finally,according to the criterion, the Hamiltonian system given by (12) is nonintegrable, except in the cases described above.
We formulate the following theorem:
Theorem: Consider the analytic autonomous Hamiltonian system
H(J,w;ε)=H0(J1,J2)+ε[A1(J1,J2)cos(k1w1−k2w2)+A2(J1,J2)cos(l1w1−l2w2)], where H0 is Liouville-integrable and non-degenerate, and
k=(k1,−k2),l=(l1,−l2) are linearly independent integer vectors.
Let
Wmk,nk=cos(mkw1−nkw2), with
mk=pkk1+qkl1,nk=jkk2+skl2, where
pk,qk,jk,sk∈Z−{0},k∈N−{0}. Assume that the corresponding resonance ratios
ω2ω1=k1+mkk2+nk(19)
or, analogously,
ω2ω1=l1+mkl2+nk
generate a dense set of resonant tori in the considered action domain. Assume further that, on this dense resonant set, the second-order equation (16) does not vanish. In particular, either for the relevant resonant contribution,
nk∂J2∂A1−mk∂J1∂A1=0,
with analogous non-vanishing condition for A2, or H0 does not possess a linear integral with respect to J1,J2 or that the remaining averaged Poisson-bracket contributions do not cancel each other we proved that for sufficiently small
ε=0,
the perturbed Hamiltonian H admits no second first integral which is functionally independent of H and analytic in ε in the considered open domain that is continued in the perturbed system. Consequently H is nonintegrable.
Remark: The theorem is inconclusive in cases:
Where all second-order resonant coefficient vanishes, as well as the coefficients,
nk∂J2∂Ai−mk∂J1∂Ai=0,i=1,2,k∈N−{0}
The two Fourier vectors are not independent, i.e.
k∥l;
then it is known that the system is integrable.
Failure of any of these hypotheses does not imply integrability; it only means that the present criterion is inconclusive.
Conclusions
In this work we have investigated a two-degree-of-freedom nearly integrable Hamiltonian system whose perturbation contains two linearly independent Fourier modes. By combining a Lie-series symplectic transformation with the non-integrability approach of Poincaré and the criterion of Meletlidou–Ichtiaroglou, we obtain a second-order resonant obstruction to the existence of an additional analytic first integral.
The Lie transformation generates combination harmonics at order O(ε2). When one of these harmonics becomes resonant with an unperturbed invariant torus, its average over the corresponding periodic orbit does not, in general, vanish and depends on the phase labelling the periodic orbit on that torus.
If this non-vanishing, orbit-dependent averaged contribution occurs on a dense set of resonant tori, and if the remaining second-order averaged terms do not cancel it, the necessary condition for the existence of a second first integral analytic in ε is violated. Consequently, for sufficiently small ε=0, no second first integral functionally independent of H and analytic in the stated class can exist in the considered domain.
Thus, the obstruction is not produced merely by the presence of two independent Fourier modes, but by their second-order interaction together with the density and non-vanishing conditions stated in the theorem.
The Lie-series formulation is particularly convenient for this analysis because it preserves the canonical structure and provides an explicit order-by-order description of the harmonics generated by the interaction of the perturbing terms.
Dense set of Resonant Tori
We prove that the resonance ratios
k1+mkk2+nk
form a dense set in the relevant positive frequency-ratio interval.
Recall that
mk=pkk1+qkl1,nk=jkk2+skl2,
where
pk,qk,jk,sk∈Z∖{0}.
Let
d1=gcd(k1,l1),d2=gcd(k2,l2).
By Bézout’s identity,
{pk1+ql1:p,q∈Z}=d1Z,
and similarly,
{jk2+sl2:j,s∈Z}=d2Z.
Therefore mk and nk may attain arbitrary integer multiples of d1 and d2, respectively.
Consider an arbitrary open interval
(α,β)⊂R+,0<α<β.
Since mk may take arbitrarily large values, we can choose mk such that
Dk:=k1+mk>0
and
Dkd2<β−α.(20)
Now choose nk such that
Nk:=k2+nk
which is a large enough multiple of d2 satisfying
Nk>αDk⇒DkNk>α
Therefore the set of integers
B={nk:k∈N, DkNk>α}
is nonempty. In the above set there is a minimum integer nk such that Sk=minnk, nk∈B for which it holds that Sk∈B.
We will show that for Sk it holds that
DkSk<β
Let us suppose the opposite
DkSk≥β
then,
DkSk−Dkd2≥β−Dkd2>β−(β−α)=α
which is a contradiction, so
DkSk<β
since Sk is the minimum of the set B, and Dk may attain arbitrary large values,
α<DkSk<β
The interval (α,β) is arbitrary, the set
{k1+mkk2+nk}
is dense in R+.
The same argument applies to the second family
{l1+mkl2+nk},
and therefore it is also dense in R+.
Finally, since the unperturbed Hamiltonian H0 is non-degenerate, the map
J↦ω(J)=∇H0(J)
is locally invertible. Hence the the above rational frequency ratios imply that the corresponding resonant invariant tori form a dense set in the considered regular action domain.