Characteristic oscillations in frequency-resolved heat dissipation of linear time-delayed Langevin systems: Approach from the violation of the fluctuation-response relation

技術要約：線形時間遅延ランジュバン系における周波数分解された熱散逸の特性振動

問題提起
時間遅延は自然界および工学的なシステムにおいて遍在しており、非マルコフ的な記憶効果をもたらし、しばしば特有の非平衡な特徴（単一の熱浴からシステムへの正味の熱流、すなわち見かけ上の負の熱散逸など）を引き起こします。確率熱力学は、時間遅延を持つシステムに対して揺動定理や不確定性関係を拡張することに成功してきましたが、既存研究の多くは形式的な理論的枠組みや、平均熱散逸率のようなグローバルな量に焦点を当てています。時間遅延システムの熱力学的特徴について、単一の量を超えた詳細な解析は不足しています。具体的には、時間遅延の効果を実験的に特定するのに役立つ可能性にもかかわらず、異なる周波数スケールにおける熱散逸の分布を支配するメカニズムについては、未だ十分に探求されていません。

手法
著者らは、確率遅延微分方程式（SDDE）によって記述される、線形な時間遅延を持つ過減衰ランジュバン系を調査しています。核となる手法は、ハラダ・ササの等式を用いて熱散逸率を周波数成分に分解することです。この等式は、平均熱散リット率と、周波数領域における揺動・応答関係（FRR）の破れとを関連付けます。

研究の手順は以下の通りです：

1. モデルの定義： 著者らは、白色ガウスノイズによって駆動される、線形な時間局所的および時間遅延的な力を伴う一般的なクラスのシステムを定義します。まず単一の遅延時間（$\tau$）に焦点を当て、後に複数の遅延時間へと解析を拡張します。
2. 解析的導出： 線形時間遅延ランジュバン方程式の既知の解析解を用いて、定常状態の線形応答関数（$R_v$）と速度自己相関関数（$C_v$）を導出します。
3. スペクトル分解： ハラダ・ササの等式を適用することで、熱散逸のスペクトル分解 $\tilde{C}_v(\omega) - 2T\tilde{R}'_v(\omega)$ を計算します。ここで、$\tilde{C}_v$ は速度のパワースペクトル密度であり、$\tilde{R}'_v$ は応答関数の実部です。
4. 漸近解析： 低周波数および高周波数の極限におけるスペクトルの挙動を分析し、特性的なスケーリング則と振動パターンを特定します。
5. 数値的検証： 単一および複数の遅延シナリオの両方について、オイラー・丸山法を用いた数値シミュレーションを通じて解析結果を確認します。

主な貢献と結果
本論文は、線形時間遅延システムにおける熱散逸のいくつかの具体的な特性を確立しています：

• 特性振動挙動： 周波数分解された熱散逸のスペクトルは、全周波数領域にわたって持続的な振動を示します。この振動の周期は、遅延時間 $\tau$（具体的には周期 $2\pi/\tau$）によって直接決定されます。
• 高周波漸近挙動： 高周波極限（$\omega \gg |a|/\gamma, |b|/\gamma$）において、スペクトルは周波数に対して反比例する（$1/\omega$）減衰エンベロープを持つ正弦波状の振動を示します。このエンベロープの振幅は、時間遅延力の強さに比例します。この $1/\omega$ スケーリングは、マルコフ的なシステムで通常観察される $1/\omega^2$ スケーリングとは対照的です。
• 低周波の符号と大きさ： スペクトルの低周波数領域における挙動は、全熱散逸率の符号と大きさを決定します。例えば、正の遅延力（$b > 0$）は負の熱散逸（エネルギーの熱浴からシステムへの流入）をもたらす可能性があり、これは低周波数スペクトルの特定の谷の卓越性に反映されます。逆に、負の遅延力（$b < 0$）は正の熱散逸をもたらします。
• 複数遅延への一般性： 著者らはこれらの知見を複数の遅延時間を持つシステムへと拡張しています。これらの場合、スペクトルは異なる遅延時間に対応する正弦波の重ね合わせとなり、$1/\omega$ の減衰エンベロープを維持します。この解析は、振動的な性質と遅延パラメータへの依存性が、線形時間遅延システムの一般的な特徴であることを裏付けています。

意義と主張
著者らは、自らの結果が、平均的な量を見ること以上に、時間遅延システムの熱力学的特徴に関する具体的かつ詳細な解析を提供するものであると主張しています。主な意義は、その知見の実験的なアクセシビリティにあります。揺動・応答関係の破れは、軌跡解析と摂動応答を通じて実験的に測定可能であるため、特性的な振動スペクトルは、時間遅延効果を検出・分析するための特定のシグネチャーを提供します。

本論文は、このアプローチが、アクティブマター、フィードバック制御下のコロイド粒子、電気回路、機械振動子など、多様な実験システムにおいて時間遅延効果を特定するために使用できることを示唆しています。さらに、著者らは、特定の $1/\omega$ スケーリングと漸近的な振動が、マルコフ的なシステムと比較した際の、時間遅延システムの識別的な特徴として機能する可能性があると述べていますが、より広範な非線形システムにおけるこれらの挙動の検証は、今後の理論的研究の方向性であると控えめに示唆しています。本研究は新しい実験セットアップを提案するものではなく、むしろ既存のデータを用いてエネルギーコストを定量化し、記憶効果を特定するための解析の方向性を概説するものです。