Post-processed frozen-flow methods for the long time sampling of ergodic dynamics on Riemannian manifolds

Technische Samenvatting: Post-processed frozen-flow methodes voor de langetermijn-sampling van ergodische dynamica op Riemanniaanse variëteiten

Probleemstelling
Het artikel behandelt de uitdaging van het ontwerpen van hogere-orde numerieke methoden voor het samplen van de invariante maat van ergodische stochastische differentiaalvergelijkingen (SDE's) op Riemanniaanse variëteiten. Specifiek beschouwen de auteurs SDE's van de vorm:
$$dX(t) = F(X(t))dt + \sqrt{2} \sum_{d=1}^D E_d(X(t)) \circ dW_d(t)$$
waarbij $M$ een gladde, volledige, verbonden Riemanniaanse variëteit is, $E_d$ een globale orthonormale frame vormt, en $F$ een vectorveld is (dat in het geval van Riemanniaanse Langevin-dynamica uit een potentiaal $V$ kan voortvloeien).

Bestaande benaderingen kampen met aanzienlijke beperkingen:

1. Extrinsieke methoden (bijv. projectiemethoden) embedden de variëteit in een hoger-dimensionale Euclidische ruimte. Hoewel deze handig zijn, leiden ze vaak tot strikte tijdstapbeperkingen en stabiliteitsproblemen.
2. Bestaande intrinsieke methoden richten zich doorgaans op sterke convergentie of zwakke orde één, zijn beperkt tot specifieke variëteiten (zoals Lie-groepen), of geven geen prioriteit aan de langetermijn-samplingefficiëntie die vereist is voor ergodische maten.

Het doel is om intrinsieke methoden te ontwikkelen (die uitsluitend vertrouwen op geometrische operaties zoals geodesen en parallel transport, zonder embeddings of lokale coördinaten) die een hoge-orde nauwkeurigheid bereiken voor de invariante maat, waarbij de nadruk ligt op de efficiëntie van langetermijn-sampling boven de nauwkeurigheid in eindige tijd.

Methodologie
De auteurs stellen een nieuw kader voor op basis van post-processed frozen-flow methoden. Deze benadering combineert twee kernconcepten:

1. Frozen-flow methoden: Een stochastische uitbreiding van de Crouch-Grossman en commutator-vrije Lie-groep methoden. Deze methoden evolueren op de variëteit met behulp van elementaire geometrische operaties (geodesen en parallel transport) geassocieerd met een specifieke krommingsvrije verbinding, de Weitzenböck-verbinding, in plaats van de Levi-Civita-verbinding.
2. Post-processing: Een willekeurige perturbatie van de identiteitsafbeelding die slechts één keer aan het einde van het traject wordt toegepast. Deze techniek, geïnspireerd door de Leimkuhler-Matthews methode in de Euclidische ruimte, maakt een hogere-orde nauwkeurigheid voor de invariante maat mogelijk zonder de computationele kosten per stap aanzienlijk te verhogen.

Om de ordevoorwaarden voor deze methoden af te leiden, maken de auteurs gebruik van een algebraïsch formalisme gebaseerd op exotische Lie-Butcher-reeksen. Zij introduceren:

• Planaire exotische bossen (planar exotic forests): Een generalisatie van standaard Butcher-bomen naar variëteiten, waarbij "exotische" decoraties worden toegevoegd om de Laplace-Beltrami-operator te representeren.
• Integratie door delen (IBP) op bossen: Een nieuwe algebraïsche operatie afgeleid van de Riemanniaanse integratie door delen formule tegen de Gibbs-maat. Deze operatie stelt de auteurs in staat om de adjunct van de numerieke operatoren ($A_j^*$) uit te drukken in termen van irreducibele bossen.
• Gemodificeerde Shuffle Algebra: De auteurs definiëren een karakter op een gemodificeerde shuffle algebra van exotische bossen om de ordevoorwaarden voor de invariante maat te karakteriseren.

Belangrijkste bijdragen

1. Intrinsieke hogere-orde karakterisering: Het artikel biedt een rigoureus criterium voor hogere-orde nauwkeurigheid van de invariante maat op Riemanniaanse variëteiten. Het generaliseert eerdere Euclidische resultaten (bijv. [2, 24, 65]) naar de setting van variëteiten, en bewijst dat als specifieke algebraïsche voorwaarden op de coëfficiënten van de methode worden voldaan, de fout in de invariante maat vervalt als $O(h^p)$.
2. Nieuwe numerieke schema's: De auteurs leiden drie nieuwe methoden af en presenteren deze:
   • Methode 1: Een post-processed frozen-flow methode die de Leimkuhler-Matthews methode generaliseert naar Riemanniaanse variëteiten. Deze vereist slechts één evaluatie van het vectorveld $F$ per stap en drie exponentiële afbeeldingen. Het is van zwakke orde één maar tweede orde voor de invariante maat.
   • Methode 2 & 3: Alternatieve schema's (Stochastische frozen-flow Heun en 2-staps frozen-flow Runge-Kutta) die eveneens tweede orde zijn voor de invariante maat zonder post-processing.
3. Algebraïsch kader: Het artikel introduceert de integratie door delen operatie op planaire exotische bossen en definieert een reductie-operator (RED) die algemene bossen naar irreducibele bossen brengt. Dit maakt de systematische afleiding van ordevoorwaarden mogelijk, waarbij wordt aangetoond dat het aantal vereiste voorwaarden voor de invariante maat aanzienlijk lager is dan voor de zwakke fout (bijv. 5 voorwaarden voor orde 2 versus 8 voor de zwakke fout).
4. Kostenefficiëntie: De voorgestelde methoden zijn ontworpen om computationeel goedkoper te zijn dan bestaande extrinsieke projectiemethoden of andere intrinsieke hogere-orde schema's, omdat ze meerdere evaluaties van het vectorveld en complexe coördinatentransformaties vermijden.

Resultaten
De theoretische claims worden gevalideerd door numerieke experimenten op twee klassieke variëteiten:

1. $SO(p)$ (Speciale Orthogonale Groep): De auteurs testten de methoden op een compacte Lie-groep met kwadratische en sextische potentialen. De resultaten bevestigden dat de nieuwe methoden (1, 2 en 3) de verwachte tweede-orde convergentie voor de invariante maat bereiken. Met name Methode 1 vertoonde een lagere computationele kost vergeleken met de SFF2-methode (een referentie intrinsische methode van zwakke orde 2) en bereikte snel de Monte-Carlo-drempel in het kwadratische geval.
2. $S^2$ (Sfeer): Experimenten werden uitgevoerd met Von-Mises Fisher-dynamica. De foutcurves vertoonden opnieuw de verwachte tweede-orde convergentie voor de invariante maat. Methode 1 toonde opnieuw superieure prestaties, waarbij de Monte-Carlo-drempel vrijwel onmiddellijk werd bereikt.

In beide gevallen presteerden de methoden beter dan voorgaande extrinsieke en intrinsieke benaderingen in termen van de kosten die nodig zijn om een bepaalde nauwkeurigheid te bereiken.

Betekenis en Claims
Het artikel beweert een eenvoudige generalisatie van de Leimkuhler-Matthews methode naar de Riemanniaanse setting te bieden, waarbij tweede-orde nauwkeurigheid voor de invariante maat wordt bereikt met minimale computationele overhead (één evaluatie van het vectorveld).

De auteurs benadrukken dat hun benadering:

• Volledig intrinsiek is, waarbij embeddings en lokale coördinaten worden vermeden.
• Prioriteit geeft aan de efficiëntie van langetermijn-sampling, wat cruciaal is voor ergodische dynamica, in plaats van de nauwkeurigheid van trajecten in eindige tijd.
• Een systematische algebraïsche methodologie (via exotische bossen en post-Lie algebra's) biedt voor het ontwerpen van integratoren van elke gewenste hoge orde voor de invariante maat.

Het artikel concludeert dat hoewel de huidige methoden steunen op de Weitzenböck-verbinding (die krommingsvrij is maar verschilt van de Levi-Civita-verbinding), toekomstig werk ernaar streeft deze frozen-flow methoden te laten rusten op louter Riemanniaanse operaties (Levi-Civita-verbinding) en de analyse uit te breiden naar SDE's met multiplicatieve ruis. Het ontwikkelde algebraïsche formalisme wordt ook aangewezen als zijnde potentieel toepasbaar buiten de numerieke analyse, in velden zoals rough paths en renormalisatietheorie.