2025-11-10T02:50:53.145942

Two approaches to average stochastic perturbations of integrable systems

Kuksin

We discuss two approaches to study the long-time behaviour and infinite-time behaviour of solutions for integrable hamiltonian systems under small stochastic perturbations. Then we compare these results with those for deterministic perturbations of integrable systems.

academic

可積系統の平均確率摂動に対する2つのアプローチ

基本情報

論文ID: 2509.17159
タイトル: Two approaches to average stochastic perturbations of integrable systems
著者: Sergei Kuksin
分類: math.DS（力学系）
発表日時: 2025年10月11日（arXiv v2）
論文リンク: https://arxiv.org/abs/2509.17159

要旨

本論文は、小確率摂動を受ける可積ハミルトン系の解の長時間挙動および無限時間挙動を研究する2つの方法について論じ、これらの結果を可積系の確定的摂動の結果と比較している。

研究背景と動機

問題の背景

古典的問題: 可積ハミルトン系の小摂動の大時間スケールでの挙動を研究することは力学系の古典的問題であり、18世紀末に遡る
従来の方法: 既存研究は主に2つの場合に焦点を当てている：
- A) 摂動が小さい滑らかなベクトル場である場合
- B) 摂動が小さいハミルトン場である場合
新しい課題: 本論文は第3の場合——確率摂動を考察する。これは理論的にも応用的にも重要な意義を持つ

研究の動機

確率摂動は確定的摂動よりも実際の物理系におけるノイズの影響をより適切に反映する
既存の平均化理論は主に確定的系を対象としており、確率的な場合の理論は相対的に未発達である
確率摂動下の長時間動力学挙動を扱うための新しい数学的ツールの開発が必要である

核心的貢献

可積系の確率摂動の平均化問題を研究する2つの異なる方法を提案
作用-角度座標を用いた高速-低速系の方法を確立し、問題を低速変数の研究に変換
有効方程式の方法を開発し、元の座標系で直接平均化方程式を構成
収束定理を証明し、確率摂動下の作用変数の極限挙動を与える
確定的摂動の結果との詳細な比較を提供し、確率性の安定化作用を明らかにする

方法の詳細

基本的な設定

シンプレクティック空間 $(R^{2n}_{x,y}, dx \wedge dy)$ におけるBirkhoff可積ハミルトン系を考える。ハミルトン関数は作用変数 $I_j = (x_j^2 + y_j^2)/2$ のみに依存する。複素座標 $v_j = x_j + iy_j$ を導入すると、系は以下のように書ける：

$\frac{\partial}{\partial t}v_k(t) = i\nabla_k H(I)v_k, \quad k = 1, \ldots, n$

ε-小確率摂動

以下の確率摂動系を考える（低速時間 $\tau = εt$ ）：

$\dot{v}_k(\tau) = iε^{-1}\nabla_k H(I)v_k + P_k(v) + \sum_{j=1}^{n_1} B_{kj}(v)\dot{\beta}^c_j(\tau)$

ここで $\{\beta^c_k(t)\}$ は標準独立複Wiener過程である。

基本的な仮定

(A1) Anosov非退化条件: ほぼすべての $I \in R^n$ に対して、周波数ベクトル $\nabla H(I)$ の成分は整数上で線形独立である (A2) 局所Lipschitz条件と多項式増大条件 (A3) 方程式の適切性とモーメントの一様有界性 (A4) ノイズの非退化性：行列 $B(v)$ のランクは $n$ に等しい

第1の方法：高速-低速系

作用-角度座標変換

作用-角度座標 $(I, \phi)$ を導入する：

$I_j = \frac{1}{2}|v_j|^2$
$\phi_j = \arg v_j$

低速変数方程式

作用変数は以下を満たす： $\dot{I}^ε(\tau) = F(v^ε) + G(v^ε)\dot{\beta}(\tau)$

これは高速-低速系であり、作用 $I^ε_k$ は低速変数、角度 $\phi^ε_k$ は高速変数である。

平均化方程式

高速変数上で平均化することにより以下を得る： $\dot{I}(\tau) = \langle\langle F \rangle\rangle(I) + \langle\langle G \rangle\rangle(I)\dot{\beta}(\tau)$

ここで平均化は以下のように定義される： $\langle\langle F \rangle\rangle(I) = \frac{1}{(2\pi)^n} \int_{T^n} F(I, \phi) d\phi$

主要な結果

定理2.1: 解 $v^ε(\tau)$ の作用曲線 $I^ε(\tau)$ は分布の意味で相対コンパクトであり、各極限過程は平均化方程式の弱解である。

制限事項: 行列 $B(v)$ が定数でない場合、平均化係数はHölder-1/2滑らかさのみを持つ可能性があり、解の一意性の問題につながる。

第2の方法：有効方程式

直接平均化

Krylov-Bogolyubov平均化に着想を得て、複素座標で直接有効方程式を構成する：

$\dot{a}(\tau) = \langle\langle P \rangle\rangle(a) + \langle\langle B \rangle\rangle(a)\dot{\beta}^c(\tau)$

ここで平均化は回転演算子により定義される： $\langle\langle P \rangle\rangle(a) = \frac{1}{(2\pi)^n} \int_{T^n} \Phi_\omega \circ P(\Phi_{-\omega}a) d\omega$

主要な利点

作用-角度座標の特異性を回避する
元の方程式の滑らかさを保持する
常に一意解を持つ

核心定理

定理3.1: 仮定A1)-A4)の下で、有効方程式は一意解を持ち、 $ε \to 0$ のとき、 $I^ε(\tau)$ は分布の意味で $a(\tau)$ の作用に収束する。

定理3.3: さらに仮定A5)（長時間推定）を仮定すれば、収束性は時間に関して一様に成立する。

定理3.4: $P(v) = P_1(v) + P_2(v)$ であり、 $P_2$ がハミルトン場である場合、有効方程式は非ハミルトン部分のみを平均化するよう修正できる。

実験的設定と応用

応用シーン

5.1 振子チェーン: 非線形振子チェーンを考える： $\ddot{q}_k = -Q(q_k), \quad k = 1, \ldots, n$ ここで位置エネルギー $Q$ は滑らか、凸、奇関数であり、零点近傍で $Q(q) = αq + βq^3 + o(|q|^4)$ である。

5.2 減衰/駆動ハミルトン系: 方程式を考える： $\dot{v}_k(\tau) = iε^{-1}\nabla_k H(I(v))v_k + 2i\frac{\partial}{\partial \bar{v}_k}h(v) - ν_k v_k + b_k \dot{\beta}^c_k(\tau)$