First passage time distribution in underdamped harmonic oscillators

技術要約：弱減衰調和振動子における第一通過時間分布

問題の定義
本論文は、位置 $x$ が閾値 $x_B$ を超える際の、弱減衰調和振動子の第一通過時間分布（FPTd）$P(t_{fp})$ の計算について扱っている。FPTd は過減衰系においては確立されているが、慣性による長期的記憶効果（long-memory effects）が存在する弱減衰系（高い品質係数 $Q$ を特徴とする）においては依然として複雑である。分布の形状は $Q$ に強く依存するため、異なる減衰領域に対して異なる理論的アプローチが必要となる。著者らは、$Q \sim 10$ の中間的な値から $Q \gg 100$ の非常に大きな値に至る全範囲をカバーする包括的な $P(t_{fp})$ の導出を行い、脱出イベントを駆動する特定のノイズパターンを分析することを目的としている。

手法
著者らは、質量 $m$、剛性 $k$、粘性減衰 $\gamma$ を持ち、温度 $T$ の熱浴と結合した調和振動子のランジュバン方程式を用いてシステムをモデル化している。ダイナミクスは、無次元単位を用いた位相空間 $(x, v)$ を用いて解析される。$P(t_{fp})$ の導出は、時間スケールと初期条件に基づき、3つの異なる寄与に分割されている：

1. 即時的寄与 ($P_I$): 閾値より上で開始する軌跡 ($x(0) \ge x_B$) を考慮する。
2. 短時間寄与 ($P_{II}$): 高い初期エネルギーを持つ軌跡 ($E(0) \ge B$) が、最初の振動周期内に閾値を横切る場合を考慮する。
3. 長時間寄与 ($P_{III}$): 低い初期エネルギーを持つ軌跡 ($E(0) < B$) が、熱活性化によって障壁を越える場合を考慮する。

長時間の振る舞いを導出するために、著者らは品質係数に応じて2つの異なる形式を用いている：

• 中間的な $Q$ ($Q \sim 10$) の場合: 著者らは、吸収境界条件を持つクラマース線形微分作用素を利用する。この非自己共役作用素の固有値と固有関数を解く。長時間の FPTd は、脱出率に対応する最も遅い減衰固有値 ($\lambda_1$) によって支配される。
• 非常に大きな $Q$ ($Q \gg 100$) の場合: エネルギー拡散モデルを導出する。極座標（エネルギー $E$ と位相 $\theta$）に変換し、速い位相変数について平均化することで、システムをエネルギーのみに関するフォッカー・プランク方程式へと簡約する。平均第一通過時間は、エネルギー拡散過程に対する随伴方程式を解き、それに続く初期平衡分布による平均化によって算出される。

さらに、著者らはインスタントン理論と数値シミュレーションを用いて、脱出イベントに至るまでの位置、速度、および熱ノイズの平均軌跡を分析し、特定のノイズパターンを特定している。

主な結果

1. FPTd の分解: 全分布は、和 $P(t_{fp}) = P_I \delta(t_{fp}) + P_{II}(t_{fp}) + P_{III}(t_{fp})$ として表される。
   • $P_I$ と $P_{II}$ はハミルトン力学から導かれ、$Q$ に依存しない。
   • $P_{III}$ は $Q$ に依存し、ポテンシャルの井戸からの脱出を記述する。
2. 領域依存の脱出率:
   • 中間的な $Q$ において、脱出率 $\Gamma_B$ はクラマース作用素の支配的な固有値 $\lambda_1$ によって決定される。この作用素の数値的な対角化は、直接的なランジュバンシミュレーションと極めて良好に一致している。
   • 非常に大きな $Q$ において、エネルギー拡散アプローチは平均第一通過時間 $\tau_Z(B)$ と、それに対応する率 $\Gamma_Z$ を導出する。著者らは、大きな $Q$ において、エネルギー障壁を横切った後に閾値に到達する時間は近似的に $\pi$ （半周期）であり、全率が $\Gamma_Z = 1/(\pi + \tau_Z)$ となることを見出した。
   • 本論文は、中間的な $Q$ では、振動子が空間的閾値 $x_B$ に到達する前にエネルギー障壁を何度も横切る可能性があるため、エネルギー拡散近似が平均 FPT を過小評価することを指摘している。
3. 過渡的挙動: 大きな $Q$ において、FPTd は漸近的な指数関数的減衰に至る前に、「階段状」の減衰を示す。これらのプラトー（平坦部）は、振動周波数に関連する虚部を持つ高次固有値 ($\lambda_{n>1}$) の寄与から生じるもので、その周期性は $2\pi$ である。
4. ノイズパターンと平均軌跡: $t_{fp}$ で終了する平均軌跡の分析により、特定の形状が明らかになった：正弦波信号の指数関数的な成長と、それに続く指数関数的な減衰である。平均ノイズ $\langle \eta(t) \rangle$ は、振動子の周波数と共鳴する明確なパターンを示す。横切りの瞬間における速度分布は、$P(v) = v e^{-v^2/2}$ であり、$Q$ に依存しないことが判明した。

意義と主張
著者らは、彼らの理論的枠組みが、すべての品質係数において弱減衰調和振動子の FPTd に対する完全な記述を提供するものであると主張している。主な意義は以下の通りである：

• 手法の統一: 減衰領域に応じて、システムを正確に記述するために異なる数学的ツール（クラマース作用素の固有値 vs エネルギー拡散）が必要であることを示した点。
• 定量的一致: 理論的な予測は、ランジュバン方程式の直接的な数値シミュレーション、および $Q=7$ の振動子で以前に測定された実験データ（付随するレター [17] に参照）の両方と極めて良好に一致している。
• 物理的洞察: 脱出イベントを駆動する特定のノイズパターンの特定は、弱減衰系における平均 FPT が単なる統計的な平均ではなく、特定のノイズシグネチャを特徴とする再現可能な動的な経路（インスタントン）によって駆動されていることを示唆している。

本論文は、FPTd の指数関数的減衰はよく理解されている一方で、過渡的な挙候（プラトー）や脱出を駆動する特定のノイズパターンは、高次の効果を解析的に決定することが依然として困難な領域であるが、その存在はシミュレーションによって明確に支持されていると結論付けている。