Hopf algebra structures for the backward error analysis of ergodic stochastic differential equations

技術要約：エルゴード的確率微分方程式の逆誤差解析のためのホップ代数構造

問題提起
決定論的な常微分方程式（ODE）における幾何学的数値積分において、基礎となるバックワード誤差解析（BEA）は、数値積分器を修正された微分方程式の厳密解として解釈することを可能にする。この枠組みは、特定のケース（例えば、2次までのミルスタイン・スキーム）における強近似および弱近似に対して、確率微分方程式（SDE）へと成功裏に拡張されてきたが、エルゴード力学（不変測度）のサンプリングには単純には一般化されない。

不変測度に関する既存のアプローチ（[33]や[3]によるもの）は、SDEの生成子と不変測度の密度を通じて定義される修正ベクトル場 $\tilde{f}$ の反復計算に依存している。これらの計算は複雑であり、統一的な代数的表現を欠いており、エキゾチックなブッチャー級数（Butcher series）を用いて3次を超えて拡張できることも証明されていない。さらに、決定論的なケースでは修正ベクトル場がブッチャー級数として表現されるのに対し、確率的な設定においては、修正ベクトル場を形式的な級数として明示的に表現する手法が存在しない。本論文は、弱収束および不変測度のサンプリングに関して、級数の合成と代入を扱うための厳密な代数的基盤を構築する必要性に対処するものである。

手法
著者らは、標準的なB級数およびアロマティックB級数を確率的な文脈へと一般化した**エキゾチック・アロマティックS級数（exotic aromatic S-series）**に基づく新しい代数的枠組みを開発している。この手法の核心は以下の通りである：

1. 組合せ構造： 本論文では、装飾されたアロマティック森林（decorated aromatic forests）、そして極めて重要なことに、**クランプされた森林（clumped forests）**を導入している。クランプされた森林は、アロマティック・ツリーの単項式を表し、アロマティック・ツリーの普遍包絡代数を形成する。標準的なアロマティック森林とクランプされた森林の間の主要な相違点は、「アロマ（aromas）」が根付き成分にどのように付着するかにある。
2. ホップ代数構造：
   • グロスマン・ラーソンソン・ホップ・アルボイド（Grossman-Larson Hopf Algebroid）： 著者らは、装飾されたアロマティック森林上のグロスマン・ラーソンソン・ホップ・アルボイドを構築している。この構造は、S級数の合成則に関連している。
   • CEM余積と代入： 代入則を扱うにあたり、「アロマ」の存在（具体的には、確率的ノイズやリアナ（lianas）の文脈において）によって複雑化するため、著者らはクランプされた森林に関する概念を導入する。彼らは、森林に対する余積 $\Delta_{CEM}$（Calaque–Ebrahimi-Fard–Manchon）を定義している。この余積により、ペアになった頂点（リアナ）やアロマの完全性を維持しながら、黒い頂点を他のツリーで置換することが可能になる。
   • D-代数： 本研究では、グラフト積（grafting products）と発散写像を定義するために、跡を持つ可換D-代数（tracial commutative D-algebras）の理論を利用し、Guin-Oudomのプロセスをこれらの構造へと拡張している。
3. 部分積分（IBP）： 本論文では、S級数上の同値関係 $\sim$ （部分積分）を用いて、原始元（森林）をツリーへと変換する手法を用いている。これは、不変測度の文脈において修正ベクトル場を定義するために不可欠である。

主な貢献と結果

1. 明示的な合成および代入則：

   • 定理 2.5（合成則）： 2つのエキゾチック・アロマティックS級数の合成が、ブッチャー・コンヌ・クレイマー余積（$\Delta_{BCK}$）に関する係数写像の畳み込み積に対応することを確立している。
   • 定理 2.8（代入則）： クランプされた森林上のCEM余積（$\Delta_{CEM}$）を用いた、エキゾチック・アロマティックS級数の代入則を提供する。これにより、アロマが存在する場合の級数の代入における困難さが解決される。

2. 不変測度のためのバックワード誤差解析：

   • 定理 3.10： S級数の係数写像がキャラクタクターである一貫した手法に対して、エキゾチックB級数として明示的に記述できる修正ベクトル場 $\tilde{f}$ が存在することを証明している。この級数の係数写像 $b$ は、森林をツリーへと変換する演算子 $A$（IBPによる）と代入則を用いた明示的な再帰公式によって与えられる。
   • この結果は、エルゴード的SDEに対する修正ベクトル場の、任意の次数における初の明示的な表現を提供するものであり、このベクトル場が直交等変性のような自然な幾何学的性質を満たすことを保証する。

3. サンプリングのための修正方程式：

   • 定理 3.14： 修正方程式の構成へとバックワード誤差解析を拡張している。これは、修正されたベクトル場（エキゾチックB級数として）が存在し、その修正された方程式に数値積分器を適用することで、不変測量に対して厳密になることを示している。これは、不変測度サンプリングのための数値解法の次数を向上させる経路を提供する。

4. 代数的統一：

   • 本論文は、確率的オーダー条件、積分器の合成、およびポストプロセッシングの研究を、エキゾチック・アロマティックS級数の形式の中で統一している。これらを行うための代数的構造（ホップ・アルボイドおよびホップ代数）は、決定論的な幾何学的積分で使用されるものに対する自然な拡張であることを示している。

意義と主張
本論文は、弱の意味での、および不変測量のための、確率的数値解析を体系的に研究するために必要なホップ代数の基礎を提供することを主張している。その主な意義は以下の通りである：

• 「煩雑な計算」問題の解決： 合成則および代入則に関する明示的な代数的公式を提供することで、修正ベクトル場の逐次的なオーダーごとの組合せ的導出という手間を排除している。
• 明示的な修正ベクトル場： エルゴード的SDEに対する、エキゾチックB級数としての修正ベクトル場の初の明示的な表現を提供しており、これは低次を超えてこれまで証明されていなかった結果である。
• 幾何学的整合性： この代数的枠組みは、導出された修正ベクトル場が、長期積分において極めて重要な直交等変性などの重要な幾何学的性質を満たすことを保証する。
• 一般化： 本フレームワークは、B級数およびアロマティックB級数の決定論的理論を、加法的ノイズ（リアナ）や多様体制約（ストロン/stolons）による特有の複雑さを扱う確率的な設定へと一般化している。

著者らは、理論がユークリッド空間（$\mathbb{T}^d$）については完全に開発されているものの、これらを一般的な多様体へ拡張するには、部分積分演算子および射影法の具体的な構造に関するさらなる研究が必要であると述べている。本論文は、これらの代数的ツールが、エルゴード的確率力学のための高次積分器の開発において不可欠であると結論付けている。