Kinetic temperatures and inertial effects in a nonequilibrium bead-spring model

技術要約：非平衡ビーズ・スプリング・モデルにおける運動学的温度と慣性効果

問題提起
本論文は、異なる温度の複数の熱浴に結合された系における、非平衡定常状態（NESS）の特性評価という課題に取り組んでいる。この領域における中心的な困難は、温度の定義が一意ではないことである。なぜなら、NESSは一般にボルツマン分布に従わず、異なる自由度が単一の値に平衡しない可能性があるからである。様々な温度のような量（例：揺動散逸関係による有効温度、構成温度など）が存在するものの、慣性系におけるエネルギー貯蔵および不可逆性に対するそれらの役割は依然として不明確である。具体的には、エネルギー貯蔵と不可逆性における慣性の明確な役割、および特異な極限（質量ゼロの過減衰限界または無限質量限界）における観測量の振る舞いについて、さらなる調査が必要とされている。

手法
著者らは、任意の質量（$m_1, m_2$）を持つ2つの結合されたビーズからなる最小モデルを調査している。これらのビーズはバネ（$\kappa$）によって接続されており、それぞれが摩擦係数（$\gamma_1, \gamma_2$）とバネ（$k_1, k_2$）を介して、異なる温度（$T_1, T_2$）の熱浴に結合されている。系はアンダーダンプド（低減衰）ランジュバン方程式によって記述される。

分析には、以下の二重のアプローチが用いられている：

1. 解析的手法： 定常状態の確率分布が、共分散行列を用いたアプローチによって導出される。系はオルンシュタイン＝ウーレンベック過程として扱われ、共分散行列はリアプノフ方程式を解くことで決定される。これにより、定常状態の揺動、相関、および熱力学的観測量を厳密に計算することが可能となる。
2. 数値的手法： 著者らは、オイラー・マルヤマ・メソッドを用いて、アンダーダンプド・ランジュバン動力学の数値シミュレーションを行う。エルゴード性と収束特性を検証するために、複数の軌跡にわたるアンサンブル平均（EA）と、単一の軌跡に沿った時間平均（TA）を比較している。

主要な貢献と結果

• 不可逆性の不変性とエネルギーの経路依存性：
  本研究は、輸送・不可逆性の観測量とエネルギー貯蔵の観測量との間の根本的な区別を明らかにしている。対称的なセットアップ（バネ定数 $k$ と摩擦係数 $\gamma$ が同一の場合）において、ビーズの質量を交換（$m_1 \leftrightarrow m_2$）しても、エントロピー生成と熱流は不変である。しかし、エネルギー貯蔵の観測量（平均エネルギー、有効温度）は、質量の配置に明示的に依存する。これは、慣性が時間非対称な不可逆性の尺度を変えることなく、時間対称な観測量に対して選択的に影響を与えることを示している。

• 組織変数としての運動学的温度：
  著者らは、各ビーズの平均運動エネルギーによって局所的に定義される運動学的温度（$T_{v\alpha} = m_\alpha \langle v_\alpha^2 \rangle$）が、系の熱力学における自然な変数であることを示している。

  • 輸送： 運動学的温度の差（$T_{v1} - T_{v2}$）が、熱輸送とエントロピー生成を制御する。
  • エネルギー： 運動学的温度の和が、一般化された分配則 $\langle E \rangle = T_{v1} + T_{v2}$ を通じて平均エネルギーを決定する。
  • 構成関係： 運動学的温度の差に基づいて熱伝導率（$\lambda$）を定義すると、その構成関係は無限質量限界においても意味を成し続ける。一方で、浴の温度（$T_1 - T_2$）に基づく定義は、この極限において一貫した駆動力を提供できない。

• 特異極限における挙動：

  • 過減衰限界（$m \to 0$）： 系は過減衰限界へと滑らかに収束する。エントロピー生成はこの領域で最大となる。著者らは、この極限が速度を積分して消去する粗視化手順とは異なり、後者では「隠れたエントロピー」の寄与が生じ得ることを指摘している。
  • 無限質量限界（$m \to \infty$）： 一般に、動力学が凍結するため、熱流とエントロピー生成は消失する。しかし、最適な質量比（式19）によって定義される特定の軌跡に沿っては、有限の熱伝導率が残存する。この極限においては、運動学的温度を用いて表現された構成関係のみが、明確に定義されたままとなる。
  • 経路依存性： 有効温度の式（$T_e = C_1 T_1 + (1-C_1)T_2$）における係数 $C_1$ は、無限質量限界および過減速限界において経路依存的であり、これは有効温度が相対的な質量分布の記憶を保持していることを意味する。対照的に、輸送係数は、極限への接近方法に依存しない一意の極限値を有する。

• エルゴード性と収束：
  系はエルゴード的であることが示されている。数値比較によれば、時間平均された観測量は、同一の計算コストにおいてアンサンブル平均よりも大幅に速く収束するが、アンサンブル平均は軌跡の数が増えるにつれて変動が小さくなる。

意義と主張
本論文は、非平衡系におけるエネルギー貯蔵と不可逆性を区別する統一的な枠組みを提供することを主張している。運動学的温度を支配変数として特定することで、著者らは、対称的なセットアップにおける質量交換の下で、輸送と不可逆性は頑健である一方で、エネルギー観測量は経路依存性を示すことを示した。

著者らは、これらの結果が線形ランジュバン系およびガウス型定常状態に特有のものであることを強調している。エネルギー貯蔵と輸送の分離、および特定の不変性は、確率的浴を持たない非線形またはハミルトニアン系では成立しない可能性があると控えめに述べている。本研究は、運動学的温度が、従来の浴温度の定義が破綻する極端な極限において、一貫した熱力学的関係を定義するために不可欠であることを示唆している。これらの知見は、より複雑な不均一ネットワークの熱力学を探求するための基礎を提供し、光学トラップされたコロイドのようなメゾスコピックなプラットフォームで検証可能である。