First passage time for an underdamped harmonic oscillator and application to the power of an information engine

技術要約：不足減衰調和振動子の第一通過時間と情報エンジンへの応用

問題の定義
第一通過時間（First Passage Time, FPT）分布の推定は、化学反応速度から情報エンジンに至るまで、統計物理学における基礎的な問題である。FPT分布は、ランジュバン方程式によって記述される過減衰系については確立されているが、不足減衰（underdamped）領域においては依然として困難な課題である。不足減衰系では、慣性によって変位の二階微分が確率運動方程式に導入されるため、プロセスが非マルコフ的になる。その結果、FPT確率の計算は歴史的に「困難な作業」であり、閉形式の解（closed-form solution）は特定の極限（例：長時間の自由粒子）においてのみ一般的に利用可能であった。本論文は、有限の品質係数 $Q$ を持つブラウン型調和振動子が、位置閾値 $x_B$ を横切る際のFPT分布 $P(t_{fp})$ の計算に取り組むものである。

手法
著者らは、全時間スケールおよび品質係数をカバーするために、3つの異なるアプローチを組み合わせた統一的な理論フレームワークを提案している。

1. PDFの分解： 全確率密度関数 $P(t_{fp})$ は、以下の3つの寄与に分解される：

   • $P_I \delta(t_{fp})$：初期位置 $x(0) \geq x_B$ である場合の、瞬時的なトリガーを表す。
   • $P_{II}(t_{fp})$：最初の振動周期内（$0 < t_{fp} \leq 2\pi$）のダイナミクスを記述する。
   • $P_{III}(t_{fp})$：長時間挙動（$t_{fp} > 2\pi$）を捉える。

2. 短時間領域（$t_{fp} \leq 2\pi$）： 低減衰（$Q \gg 1$）の場合、著者らはハミルトニアン近似を採用する。ここでは、全エネルギー $E$ が1周期にわたってほぼ一定であると仮定する。

   • 瞬時トリガー ($P_I$)：初期位置のボルツマン分布を用いて計算され、$P_I = \frac{1}{2}\text{erfc}(\sqrt{B})$ となる。ここで $B = \frac{1}{2}x_B^2$ はエネルギー障壁である。
   • 第1周期 ($P_{II}$)：ダイナミクスは決定論的な正弦波運動 $x(t) = x_E \cos(\theta_0 + t)$ として扱われる。確率は、初期エネルギー $E$ と位相 $\theta_0$ の結合分布にわたって積分することによって導出される。これは、時刻 $t$ においてアクセス可能な最大エネルギー $E^\dagger(t)$ によって決定されるプラトー（平坦部）と、それに続く急速な減衰をもたらす。

3. 長時間領域（$t_{fp} > 2\pi$）： 振動周期を超える時間において、系の挙動は熱雑音によるエネルギー拡散によって支配される。

   • 高 $Q$ 極限：この問題は、エネルギー障壁を越えるクラマース脱出問題へとマッピングされる。脱出率 $\Gamma(t)$ は、時間依存のエネルギー閾値 $E^\dagger(t)$ に基づいて計算され、最終的に一定の率 $\Gamma_B$ に収束する。
   • 有限の $Q$：著者らは、吸収境界条件 $x_B$ を持つ確率方程式に関連するクラマース微分演算子の最も遅い固有値 $-\lambda_1(x_B)$ を計算することによって、脱出率 $\Gamma$ を決定する。
   • 結果として得られるテイル（裾）$P_{III}$ は、指数関数的減衰 $e^{-\Gamma_B t_{fp}}$ に従う。

主な結果
論文は、エネルギー障壁 $B$ と品質係数 $Q$ のみによってパラメータ化された、FPT分布の閉形式の解析解（式2）を提示している。

• 実験的検証： 理論的予測は、室温において $Q \approx 7$ の不足減衰マイクロカンチレバーを用いた実験によってテストされた。位置は2 MHzのサンプリングレートでピコメートル分解能で測定された。
• 一致性： 実験データは、以下の3つの異なる領域において、理論モデルと「極めて良好に一致」した：
  1. 瞬時トリガーの確率 ($P_I$)。
  2. $0 < t < \pi$ における分布のプラトー値 ($P_{II}$)。
  3. 長時間のテイルにおける指数減衰率 ($P_{III}$)。
• 情報エンジンへの応用： 導出されたFPT分布は、情報エンジンの出力（パワー）を推定するために適用された。このセットアップでは、振動子が閾値 $x_B$ を横切る際にポテンシャルの井戸がシフトすることで、仕事が抽出される。平均出力 $\langle P \rangle$ は、平均抽出仕事と、平均FPTおよび滞留時間の和の比として計算された。理論的な出力曲線は実験データと一致し、閾値 $x_B \approx 2$ およびシフト $L \approx 1$ という最適動作点を特定した。

意義と主張
著者らは、本研究が、短時間のハミルトニアン力学と長時間の拡散的脱出の間のギャップを埋め、全範囲の品質係数をカバーする解を提供することにより、不足減答系のFPT問題に対して「新たな知見」を与えるものであると主張している。

• 理論的進展： 論文は、導出された閉形式の解（式2）が、不足減衰系のFPT統計における進歩であることを主張している。これまで、不足減衰系における既知の閉形式の解は、自由粒子（ランダム加速過程）の長時間極限のみであったことが指摘されている。
• クラマース理論との区別： 著者らは、彼らの記述が、通常は大きな障壁の拡散的横断に適用される標準的なクラマース率理論とは異なることを強調している。彼らの分布では、指数関数的なテイルは、1回の振動周期よりも有意に長い時間においてのみ出現し、短時間の挙動は慣性振動によって支配される。
• 実用的有用性： この手法により、慣性を伴う情報エンジンの出力を、分布の全容を知ることなく精密に推定することが可能になる。著者らは、これが、FPT統計を不足減衰ダイナミクスを含む他の分野（確率熱力学など）へ適用する道を開くものであると述べている。

論文は、品質係数の影響に関する計算と数値シミュレーションの詳細については、付随する論文 [22] に記載されていると述べて締めくくられている。