Geometry of Almost-Conserved Quantities in Symplectic Maps. Part II: Recovery of approximate invariant

1. Problématique et Contexte

Ce travail, deuxième partie d'une série, s'inscrit dans le cadre de la dynamique des systèmes hamiltoniens et des applications en physique des accélérateurs de particules. Le problème central est la compréhension et la modélisation des invariants approximaux dans les applications symplectiques (systèmes discrets), en particulier au-delà du voisinage immédiat des points fixes.

• Limites des approches classiques : Les méthodes traditionnelles, comme le développement de Lie ou les formes normales, sont souvent limitées à des amplitudes faibles (régime infinitésimal). À mesure que l'ordre d'approximation augmente, le domaine valide de l'invariant tend à se contracter vers un point (théorème de KAM), rendant difficile l'analyse des mouvements de grande amplitude et des systèmes chaotiques.
• Le défi des résonances : Dans les systèmes chaotiques, les résonances (où le nombre de rotation est rationnel) créent des structures complexes (chaînes d'îlots) qui obscurcissent la dynamique globale. Les méthodes existantes peinent souvent à capturer fidèlement ces structures ou à fournir des invariants globaux robustes.
• Objectif : Développer une méthode perturbative capable de reconstruire des invariants approximaux qui respectent les symétries discrètes (réversibilité) et qui restent valables sur de larges domaines de l'espace des phases, y compris dans des régimes fortement non linéaires et chaotiques.

2. Méthodologie

Les auteurs proposent une méthode perturbative systématique pour construire des invariants approximaux $K^{(n)}[p, q]$ pour une application symplectique $T: \mathbb{R}^2 \to \mathbb{R}^2$.

• Ansatz polynomial : L'invariant est construit comme une série de polynômes homogènes d'ordre croissant :
  $$K^{(n)} = K_0 + \epsilon K_1 + \dots + \epsilon^n K_n$$
  où les coefficients sont déterminés en imposant que l'invariant soit conservé à un ordre donné $O(\epsilon^{n+1})$.
• Symétries discrètes : Contrairement aux méthodes purement analytiques, cette approche impose que l'invariant préserve les symétries discrètes induites par la réversibilité du système, créant un parallèle conceptuel avec le théorème de Noether (symétries continues $\to$ lois de conservation exactes).
• Procédure d'averaging (Moyennage) : C'est l'élément clé de la méthode. Le système d'équations pour déterminer les coefficients est sous-déterminé. Les auteurs utilisent une procédure de moyennage sur les coordonnées de phase polaires $(\rho, \psi)$ pour minimiser l'erreur résiduelle :
  $$\frac{d}{dC_k} \int_0^{2\pi} R_n^2[\rho, \psi] \, d\psi = 0$$
  Cette étape est cruciale pour résoudre l'indétermination, assurer la compatibilité avec les résonances et obtenir des résultats stables.
• Théorème de Danilov : Une fois l'invariant approximatif obtenu, le nombre de rotation non linéaire $\nu$ est calculé directement à partir de l'invariant via le théorème de Danilov, évitant ainsi le calcul coûteux et complexe de l'action $J$.

3. Contributions Clés et Résultats

Le papier valide et étend la méthode à travers plusieurs tests rigoureux :

A. Récupération des systèmes intégrables

Les auteurs testent la méthode sur des applications intégrables connues (formes de McMillan) :

• Cas polynomiaux (McMillan symétrique) : La méthode récupère exactement l'invariant polynomial global. Les coefficients d'ordre supérieur générés par la procédure de moyennage s'annulent correctement, confirmant la convergence exacte.
• Cas exponentiels et trigonométriques (Suris) : Pour des invariants infinis (séries de puissances), la méthode montre une convergence « presque partout ». Les coefficients calculés convergent vers les valeurs théoriques à mesure que l'ordre $n$ augmente, bien que la convergence ne soit pas finie (comportement similaire au phénomène de Gibbs).

B. Application aux systèmes chaotiques (Cartes de Hénon)

La méthode est appliquée aux cartes quadratiques et cubiques de Hénon, des modèles standards pour la dynamique non linéaire en physique des accélérateurs :

• Reconstruction des structures : L'invariant moyen permet de visualiser avec précision les chaînes d'îlots (islands) et les séparatrices, même pour des amplitudes importantes.
• Nombre de rotation : Les courbes de nombre de rotation dérivées de l'invariant moyen correspondent remarquablement bien aux résultats du suivi numérique d'orbites, y compris près des séparatrices où les méthodes classiques échouent souvent.
• Diagrammes isochrones : En traçant la région simplement connexe de stabilité en fonction des paramètres de la carte, les auteurs montrent que la méthode capture la structure fractale des résonances (séquences de Farey) avec une grande fidélité.

C. Comparaison avec la méthode de la Matrice Carrée (Square Matrix - SM)

Une comparaison directe est faite avec la méthode SM :

• Limites de SM : Bien que SM fonctionne localement près du point fixe, elle échoue à récupérer les invariants exacts dans les cas intégrables (les termes divergent) et manque de mécanismes intrinsèques pour traiter les résonances analytiquement.
• Avantage de la méthode PT avec moyennage : La méthode proposée offre une précision supérieure, surtout dans les régimes résonnants et fortement non linéaires. Elle capture mieux les frontières des régions stables et les structures d'îlots que la méthode SM, même à des ordres élevés.

4. Signification et Impact

• Cadre analytique unifié : La méthode fournit un cadre analytique compact pour organiser la dynamique non linéaire, reliant directement les coefficients de la carte aux portraits de phase, aux nombres de rotation et aux empreintes de tune.
• Applications en physique des accélérateurs : L'approche permet d'étendre le formalisme classique de Courant-Snyder (limité au régime linéaire) aux régimes non linéaires. Elle est directement applicable à l'optimisation de l'aperture dynamique, à la conception de réseaux d'accélérateurs et à la résonance d'extraction.
• Robustesse : La démonstration que l'invariant moyen peut capturer la géométrie des structures fractales et des chaînes d'îlots dans des systèmes chaotiques représente une avancée significative par rapport aux développements perturbatifs traditionnels.

En conclusion, ce travail établit que la construction d'invariants approximaux respectant les symétries discrètes, couplée à une procédure de moyennage rigoureuse, constitue un outil puissant et robuste pour l'analyse globale de la dynamique symplectique, comblant le fossé entre les théories intégrables locales et la complexité des systèmes chaotiques globaux.