[論文レビュー] Approximating linear response by nonintrusive shadowing algorithms

Angxiu Ni|arXiv (Cornell University)|Mar 22, 2020

Theoretical and Computational Physics参考文献 40被引用数 4

ひとこと要約

この論文は、非侵襲的シャドウイングアルゴリズムがカオス的系における線形応答の 'シャドウイング寄与' のみを計算することを確立している。全微分ではなく、残りの '不安定寄与' がシャドウイング法の系統的誤差であることを証明し、不安定次元が低い場合にはこの誤差が小さいことを示しており、流体力学のような高次元系におけるシャドウイングの成功を裏付けている。

ABSTRACT

Nonintrusive shadowing algorithms efficiently compute $v$, the difference between shadowing trajectories, then use $v$ to compute derivatives of averaged objectives of chaos with respect to parameters of the dynamical system. However, previous proofs of shadowing methods wrongly assume that shadowing trajectories are representative. In contrast, the linear response formula is proved rigorously, but is more difficult to compute. We prove that $v$ gives only a part, called the shadowing contribution, of the linear response; hence, the other part, the unstable contribution, is the systematic error of shadowing methods. For systems with a small ratio of unstable dimensions, with some further statistical assumptions, we show that the unstable contribution is small. We also briefly describe an algorithm for the unstable contribution, which is simpler to derive but less efficient than the fast linear response algorithm. Moreover, we prove the convergence of the nonintrusive shadowing algorithm, the fastest shadowing algorithm, to $v$ and to the shadowing contribution.

研究の動機と目的

• カオス的系における線形応答を近似する非侵襲的シャドウイングアルゴリズムの誤差を厳密に分析すること。
• シャドウイング法の系統的誤差を '不安定寄与' と特定すること——これは全線形応答の欠落部分である。
• 非侵襲的シャドウイングアルゴリズムが、全線形応答ではなくシャドウイング寄与に収束することを証明すること。
• 不安定寄与が小さい条件を提供し、高次元系におけるシャドウイングの経験的成功を正当化すること。
• シャドウイング寄与と不安定寄与の分解を通じて、高速な線形応答アルゴリズムの理論的基盤を構築すること。

提案手法

• 線形応答を二つの部分に分解する：シャドウイング寄与（シャドウイングによって計算される）と不安定寄与（誤差）。
• 非侵襲的シャドウイングアルゴリズムを適用し、$ m $ 個の接続方程式の解を用いてシャドウイング方向 $ v $ を計算する。ここで $ m $ は不安定次元である。
• シャドウイング寄与の誤差をサンプリングベースで分析し、相関の崩壊仮定を用いてサンプリング誤差を評価する。
• 仮定3として、$ \Phi_u v $ の相関が指数関数的に減少することを導入し、サンプリング誤差 $ \|\tilde{e}^S\| $ の定量的評価を可能にする。
• ホルダー連続性と一様双曲的系における相関の崩壊を用いて、非侵襲的シャドウイングアルゴリズムの収束限界を導出する。
• 後続の研究で、2階接続方程式を用いて不安定寄与を効率的に計算する高速な線形応答アルゴリズムを提案する。

実験結果

リサーチクエスチョン

• RQ1 非侵襲的シャドウイングアルゴリズムは、線形応答のどの部分を実際に計算しているのか？
• RQ2 シャドウイング法は、既知の理論的制限があるにもかかわらず、流体力学のような高次元カオス的系でなぜ成功するのか？
• RQ3 シャドウイング法の系統的誤差とは何か？どのような条件下でその誤差は小さいのか？
• RQ4 非侵襲的シャドウイングアルゴリズムがシャドウイング寄与に収束することを厳密に証明できるか？
• RQ5 シャドウイングと効率的な不安定寄与の計算を組み合わせることで、全線形応答をどのように回復できるか？

主な発見

• 非侵襲的シャドウイングアルゴリズムは、線形応答の 'シャドウイング寄与' のみを計算しており、全微分ではない。
• シャドウイング法の系統的誤差は '不安定寄与' であり、これは線形応答の欠落部分である。
• 不安定次元が低く、相関が指数関数的に減少する条件下では、不安定寄与は小さい。これはシャドウイングの経験的成功を説明する。
• シャドウイング寄与のサンプリング誤差は、$ \|\tilde{e}^S\| \leq \sqrt{\frac{2C_3}{K(1 - \kappa_3)}} \|v^A\| $ で抑えられ、仮定3のもとで $ K^{-0.5} $ の速度で減少する。
• 非侵襲的シャドウイングアルゴリズムは、従来のシャドウイング法と同程度の速度でシャドウイング寄与に収束し、$ v $ の精度を確認する。
• 線形応答を近似する際の総誤差は、不安定寄与、サンプリング誤差、数値誤差の和で抑えられ、後二者は $ K \to \infty $ のとき消える。