On the independence of the slow and fast scales in multiple-scale expansions, with application to Van der Pol's equation

問題設定
多重尺度法（MMS）は、弱非線形振動子のように、大きく分離された時間スケールを持つ自律微分方程式を解くために用いられる標準的な摂動技法である。MMSの基本的な指示は、高速な時間スケール（$\theta$）と低速な時間スケール（$\tau$）を、それらが共に単一の物理変数 $t$ に由来しているにもかかわらず、導出過程においては独立な変数として扱うことである。この形式的な独立性は優れた数値近似をもたらすが、厳密な正当性に欠けている。なぜなら、真の解 $Y(t; \epsilon)$ は時間並進不変性を備えており、これは時間の変数をシフトさせると、両方のスケールに結合したシフトが生じることを意味する。本論文は、$\theta$ と $\tau$ が厳密に独立であるという仮定の妥当性を、特に、高速な位相 $\theta$ のシフトが、低速なスケールに対する非自明で指数関数的に小さいシフトや、低速スケールに依存する位相補正を誘発し得るかどうかを検証することで調査している。著者らは、ホップ分岐および弱非線形振動の典型的なモデルとしてファン・デル・ポール方程式に焦点を当てている。

手法
著者らは、標準的なMMSの手順を任意の高次まで拡張することで、独立性の仮定の限界を明らかにする厳密な漸近解析を用いている。その手法は以下の通りである：

1. 高次展開： 解を小パラメータ $\epsilon$ のべき級数に展開する。解析では、$j \gg 1$ の高次までの展開係数を追跡する。
2. 階乗・べき乗アンザッツ（Factorial-Over-Power Ansatz）： 著者らは、MMSによって生成される漸近級数が発散することを示す。彼らは（ディングルの理論に着想を得た）「階乗・べき乗アンザッツ」を利用して、$j \to \infty$ における係数 $A_{j,k}$ の振る舞いを特徴付ける。これにより、級数が複素 $\tau$ 平面における特異点、具体的には $\tau_n = (2n+1)i\pi$ によって支配されていることが明らかになる。
3. ストークス現象と最適截断： 解析では、これらの特異点から放射されるストークス線を特定する。発散級数をその最小項の付近で最適に截断することにより、著者らは指数関数的に小さい剰余 $R$ を孤立させる。この剰余は、標準的なMMSが無視する「全次数を超えた（beyond-all-orders）」寄与を表している。
4. マッチングとストークス・スイッチング： 剰余は、実軸を横切るストークス線にわたって解析される。著者らは、級数の後期項の振る舞いと特異点付近の局所的な振る舞いをマッチングさせることで、指数関数的に小さい項の立ち上がり（switching on）を決定する。これには、位相シフトと低速スケールの移動の具体的な形式を見出すために、剰余に関する線形方程式を解くことが含まれる。
5. 数値的検証： 理論的な予測は、ファン・デル・ポール方程式の数値シミュレーションを通じて検証される。著者らは、様々な $\epsilon$ の値に対してリミットサイクルの初期位相に対する位相シフトを計算し、そのシフトの振幅を解析的なスケーリング則と比較する。

主要な貢献と結果
本論文の主要な結果は、低速スケールと高速スケールの仮定された独立性を打破する、非自明で指数関数的に小さい位相シフト $\Delta\theta$ の導出である。

• 変換対称性への補正： 本論文は、標準的なMMSの変換 $y(\theta, \tau) \to y(\theta + \phi, \tau)$ がシステムの対称性を完全に記述しているという仮定を無効にする。代わりに、正しい変換は指数関数的に小さい補正を含む：
  $$y(\theta, \tau) \to y(\theta + \phi + \Delta\theta(\tau, \phi, \epsilon), \tau + \Delta\tau(\tau, \phi, \epsilon))$$
  ここで、$\Delta\theta$ と $\Delta\tau$ は $\epsilon$ の指数関数的に小さい関数である。
• 位相シフトの明示的な形式： ファン・デル・ポール方程式において、長時間極限（$\tau \to \infty$）における位相シフト $\Delta\theta$ は以下のように求められる：
  $$\Delta\theta \sim 4\pi\Lambda e^{-2\pi/\epsilon} \cos(2\phi)$$
  ここで、$\Lambda \approx 0.023080$ は漸化式から導出された数値定数であり、$\phi$ は初期位相である。因子 $e^{-2\pi/\epsilon}$ は、この効果が指数関数的に小さく、$\epsilon$ のいかなるべき乗よりも速く消失することを示している。
• 非線形性と高調波の役割： 効果の大きさは、非線形性によって生成される高調波の間隔によって決定される。ファン・デル・ポール方程式は三次非線形性を含むため、奇数高調波のみが生成され、その結果、フーリエ空間における高調波の間隔は $\Delta k = 2$ となる。これにより、指数因子は $e^{-2\pi/\epsilon}$ となる。著者らは、もし二次非線形性が存在し（偶数高調波 $\Delta k = 1$ を生成）、その場合の因子が $e^{-\pi/\epsilon}$ であれば、効果が著しく大きくなることを指摘している。
• 数値的確認： 数値シミュレーションは、位相シフトの存在と、初期位相 $\phi$ に対する $\cos(2\phi)$ 依存性を裏付けている。位相シフトの振幅の振動は、$\epsilon = 0.35$ までの値に対して、予測されるスケーリング $\Delta(\epsilon) \approx 0.58 e^{-2\pi/\epsilon}$ と一致している。

意義と主張
本論文は、自律的な微分問題において、標準的な摂動論ではしばしば見落とされる現象である、低速スケールと高速スケールがどのように相互作用するかという正確なメカニズムを解明することを目的としている。

• MMSの妥当性： 著者らは、MMSにおけるスケール独立性の仮定は、多くの文脈において根本的に「誤っている」のではなく、むしろ「無害（benign）」であると主張している。スケール間の相互作用は指数関数的に弱く、ストークス線の近傍（振幅の成長が最も速い場所）でのみ発生する。これらの線から離れると、スケールは事実上デカップル（分離）するため、MMSの実用的な成功が正当化される。
• 非線形的起源： 本研究は、スケール間の結合が純粋に非線形な効果であることを確立している。線形システムでは、解はスケールが厳密に因数分解された指数関数の単純な組み合わせとなるが、独立性の崩壊は、非線形による高次高調波の生成から特異的に生じる。
• 汎用性： 計算はファン・デル・ポール方程式に対して行われているが、著者らはその結果の本質が一般的であると示唆している。補正の形式は非線形性によって生成される特定の高調波に依存するため、この手法は他のホップ分岐や弱非線形波束にも適用可能である。
• 指数漸近論： 本論文は、多重スケール問題に対して指数漸近論（特に最適截断とストークス・スイッチング）をどのように適用するかを示す詳細な実例であり、他の「全次数を超えた（beyond-all-orders）」現象を解析するためのテンプレートを提供している。