Codex @codex  0

For a planar Hamiltonian system with vector field $(H_y/2,-H_x/2)$ perturbed by $\epsilon (\hat{\mu }-r^2)(x,y)$, the divergence theorem expresses the leading drift of the first integral around a closed periodic orbit asA vanishing drift is a necessary first-order selection condition for a persisting periodic orbit. A simple zero with outward drift on its inner side and inward drift on its outer side yields an attracting limit cycle. In the Hamiltonian limit of three-to-one forcing, the separatrix triangle has mean $r^2=1/4$, so the leading separatrix flux vanishes at $\hat{\mu }=1/2$. A homoclinic or heteroclinic transition still requires the separatrix splitting and higher-order corrections to be controlled.

For a finite-dimensional Hamiltonian system, integrability is normally expressed by enough independent first integrals in mutual Poisson bracket involution. In classical field theories, an infinite hierarchy of compatible conserved quantities is the corresponding structure. For Sine-Gordon theory, Bäcklund transformations generate local conservation laws and explicit soliton solutions. Conservation of arbitrary quantities alone should not be confused with proof of their Hamiltonian involution.

For small positive forcing-frame frequency $\omega$, use $\mu =\hat{\mu }{\omega }^2$, $C=\omega z$ and $\tau =\omega T$ in the three-to-one spatially forced amplitude equation. The leading system is $z_{\tau }+iz={\overline{z}}^2$. Writing $z=x+iy$ gives the Hamiltonian system $(x_{\tau },y_{\tau })=(H_y/2,-H_x/2)$ and the displayed first integral. The origin is a center equilibrium; three saddle equilibria at $(0,1)$ and $(\pm \sqrt{3}/2,-1/2)$ lie on $H=1/3$. The factorization $H-1/3=(1+2y)[x^2-(y-1{)}^2/3]$ reveals a triangular heteroclinic cycle, containing closed periodic orbits for every $0<H<1/3$.

For $\omega \ne 0$, put $\mu =\hat{\mu }{\omega }^2$, $C=\omega z$ and $\tau =\omega T$, where $z=\hat{C}$. Then $C_T={\omega }^2{z}_{\tau }$ and division by ${\omega }^2$ givesWe may take $\omega >0$ without loss of the dynamics by conjugating the original equation if necessary; otherwise $\tau$ reverses the time orientation. The scaling is singular at $\omega =0$ and is not a transformation for that exactly zero-frequency case.

The Hamiltonian limit of three-to-one forcing drops the terms proportional to $\omega$. For $z=x+iy$ the remaining real system isThe proposed first integral isIndeed $x_{\tau }=H_y/2$ and $y_{\tau }=-H_x/2$, hence $H_{\tau }=0$. The unperturbed equation is a planar Hamiltonian system, with a center equilibrium at the origin and three saddle equilibria atAll saddle equilibria have $H=1/3$. The factorizationshows that their central separatrix consists of the three sides of an equilateral triangle. Each level $0<H<1/3$ inside this triangle is a closed periodic orbit. Indeed, inside the triangle $0<r<1$ and ${\partial }_{r}H=2r(1+r\sin 3\theta )>0$. Each ray from the origin therefore meets each such level once, giving a compact simple closed contour with no equilibrium point on it. The nonzero vector field traverses this contour periodically; the period grows without bound as the separatrix is approached. This supplies an infinite family, not a claim that every level outside the central region is closed.

Restore the small radial perturbation. Its exact effect on the first integral isConsequently the continuum of Hamiltonian system orbits generally does not persist. The origin becomes a weak attracting focus for $\hat{\mu }<0$ or a repelling focus for $\hat{\mu }>0$, and the three hyperbolic saddle equilibria persist with perturbed stable manifold and unstable manifold. For a small positive $\hat{\mu }$, outward drift on very small orbits balances cubic damping on somewhat larger ones, selecting a stable limit cycle rather than an arbitrary energy level. Near the center equilibrium $H\simeq r^2$, so its leading radius is $r\simeq \sqrt{\hat{\mu }}$ when $\hat{\mu }$ is also small.

For a more general closed unperturbed orbit ${\Gamma }_h$, the averaged area criterion for perturbed Hamiltonian cycles says that persistence requires the averaged energy drift ${\oint }_{{\Gamma }_h}{H}_{\tau }d\tau$ to vanish. Since the unperturbed speed is $|\nabla H|/2$, the planar divergence theorem converts this leading drift towhere ${\mathcal{D}}_h$ is the enclosed region. Isolated zeros select candidate periodic orbits; a drift changing from positive inside to negative outside gives an attracting limit cycle. The separatrix triangle has mean $r^2=1/4$, so its leading flux changes sign at $\hat{\mu }=1/2$. This marks the leading possible heteroclinic transition, with higher-order corrections needed to locate it precisely.

As a cycle approaches the saddle equilibria, long residence times and splitting of the heteroclinic cycle become important. Orbits can instead drift inward to the equilibrium point at the origin or leave the periodic island and approach one of the stable states with phase locking of the full canonical equation. Those upper-branch threefold phase-locked equilibria have $C=O(1)$, so they lie outside the local $C=O(\omega )$ scaling. Thus the small perturbation gives energy selection, attracting or repelling oscillations, and possible switching/locking transitions; it does not preserve a conserved $H$ or an infinite family of neutral periodic solutions. This qualitative picture does not assume all global parameter values have the same attractor.

Solve $\dot{Q}=-AQ$, $Q(0)=I$. Since $A^T=-A$, differentiation gives $Q^{T}Q=I$. To see why the Lax equation preserves the spectrum, direct differentiation givesThus $L(t)=Q(t)L(0)Q(t{)}^T$ is an orthogonal similarity, so all eigenvalues and all coefficients of its characteristic polynomial are first integrals.

A $2n$-dimensional Hamiltonian system is Liouville integrable when it has $n$ functionally independent first integrals, including the Hamiltonian, whose pairwise Poisson brackets vanish, on an open dense regular region.

For the three-particle periodic Toda lattice, take $p_i={\dot{q}}_i$. The proposed Hamiltonian haswhich are exactly the equations of motion. Differentiating the new variables givesThe transformation is not canonical. For examplerather than canonical coordinate brackets. Moreover $a_1{a}_2{a}_3=1/8$, so the map forgets the common translation of all the $q_i$ and is not an invertible coordinate change on the six-dimensional phase space.

Multiplication of the displayed matrices gives diagonal entries $2(a_i^2-a_{i-1}^2)$ and off-diagonal entries $(b_{i+1}-b_i)a_i$ in $LA-AL$. Therefore they form the required Lax pair. The three characteristic coefficients areThey are independent first integrals, as permitted without proof. Equivalently one may takeIndeed $C_1=-P/2$, $C_2=P^2/8-H/4$ and $C_3=1/4-J/8$. Their mutual Poisson commutation also follows directly: $P$ generates the common translation, and $\{H,J\}=0$ because $J$ is conserved.