Numerical analysis of heat transport in classical one-dimensional systems

技術要約：古典的一次元系における熱輸送の数値解析

問題提起
本論文は、古典的な一次元（1D）系における熱輸送の数値研究において、長年存在している不一致に対処するものである。変動流体力学を中心とする理論的枠組みは、エネルギー、線運動量、および全伸張（total stretch）を保存する系は異常熱伝導性を示す（通常 $\alpha=1/3$ または $1/2$ で $L^\alpha$ に比例して発散する）と予測しているが、様々な数値シミュレーションでは、フーリエの法則を満たす有限の熱伝導率が報告されてきた。著者は、これらの「正常な」結果が真の物理的特性なのか、あるいは漸近的な発散を覆い隠してしまう有限サイズ効果によるアーティファクト（偽像）なのかを調査している。

手法
本研究では、固定境界条件を持つハミルトン力学に従う振動子鎖の非平衡定常状態（NESS）の数値シミュレーションを用いている。主要な手法の特徴は以下の通りである：

• 熱リザーバー： 鎖の両端は、速度の確率的なリセットを介して温度 $T_L$ および $T_R$ の熱浴に結合されており、これによりエネルギーフラックス $J$ を直接計算することが可能となる。
• 初期条件： 大きな系における長い過渡時間を軽減するため、特定の初期化の工夫を用いている。具体的には、長さ $L$ の鎖から得られた構成をインターリーブ（交互に配置）することで、長さ $2L$ の鎖に対して統計的に適切な初期条件を生成している。
• モデルのバリエーション： 解析は以下のモデルを対象としている：
  1. ソフト・ポイント・ガス (SPG)： 質量不均一性（$\epsilon$）を導入して可積分性を破り、平均自由行程を調整した、斥力ポテンシャル $V(u) \propto 1/u^{-2}$ を持つモデル。
  2. 非結合ポテンシャル： 有限の熱伝導率を示すと以前に主張されていた、特定の斥力ポテンシャル（片側が放物線状で、もう片側が平坦なもの）。
  3. ディン・ア・リング（Ding-a-ling）の変種： 内部運動量が保存されるモデル（調和相互作用と衝突相互作用が交互に現れる）。これは異常輸送の反例として以前に引用されたものである。
• 解析的枠組み： 有効熱伝導率 $\kappa(L)$ は、全抵抗が動力学的（拡散的）成分と流体力学的（異常的）成分の和であるとする分解モデルを用いて解析される。

主な貢献と結果

1. 熱伝導率の分解： 本論文は、準可積分系に対して元々提案された推測を、完全カオス的な非結合系へと検証・拡張している。有効熱伝導率は、以下の2つの並列な寄与の和として示される：

   • 動力学的項 ($R_{kin}$)：有限の平均自由行程 ($l$) を持つフォノン的なモードに関連し、大きな長さにおいて飽和する。
   • 流体力学的項 ($R_{an}$)：漸近的な発散 ($L^{1-\alpha}$) を引き起こす原因となる。
     有効熱伝導率は以下の形式に従う：
     $$\kappa(L) \propto \left( \frac{1}{L/l + 1/R_0} + \sigma L^\alpha \right)$$
     ここで、$R_0$ は接触抵抗、$\sigma$ は異常項の振幅である。

2. 非結合ポテンシャル (SPG およびその他)：

   • 質量が均質（$\epsilon=0$）なソフト・ポイント・ガス (SPG) のシミュレーションでは、飽和が遅いため、当初は正常な熱伝導率が示唆されていた。しかし、質量不均一性（$\epsilon > 0$）を導入することで、潜在的な発散が明らかになった。
   • データを分解モデルに適合させた結果、漸近的な領域は異常な流体力学的成分（$\alpha = 1/3$）によって支配されていることが確認された。
   • 異常な挙動が顕在化するクロスオーバー長 $L_c$ は極めて大きく（例：SPGでは $L_c \approx 2 \times 10^4$）、これは異常項の振幅 $\sigma$ が系の均質化に伴って消失するためである。これが、従来のシミュレーションが発散を検出できなかった理由である。

3. ディン・ア・リング・モデルの変種：

   • 線運動量を保存しているにもかかわらず、フーリエの法則を満たすと主張されていた、特定のディン・ア・リング・モデルの変種（文献 [19]）を再検討した。
   • シンプレクティック積分器と直接的な NESS フラックス測定を用いた新しいシミュレーションにより、熱フラックス $J$ が $L^{-2/3}$ でスケーリングすること、すなわち $\kappa \propto L^{1/3}$ であることが示された。
   • この結果は、正常な熱伝導率を示すという以前の主張と矛盾しており、以前の不一致は、局所熱フラックスの不適切な定義、あるいは不十分なシステムサイズに起因することを示唆している。

意義と主張
本論文は、一次元系における異常輸送理論の明白な違反は、その大部分が強い有限サイズ効果によって引き起こされた錯覚であると論じている。

• 異常輸送の普遍性： 著者は、エネルギー、運動量、および伸張を保存する1D系においては、漸近的な領域は常に異常な流体力学的発散によって支配されると結論付けている。
• 遮蔽メカニズム： 非結合ポテンシャルにおいて観察される「擬似的な正常」挙動は、独立した相ではなく、異常項の振幅（$\sigma$）が非常に小さく、かつ平均自由行程（$l$）が大きいことによって、動力学的な拡散成分が支配的となっているクロスオーバー領域である。
• 文献の修正： 本研究は、特定のモデル（SPGおよびディン・ア・リングの変種）の解釈を修正し、十分なシステムサイズを考慮するか、あるいは分解モデルを適用すれば、それらの漸近的挙動は理論的予測（$L^{1/3}$）と一致することを断言している。

本論文は新しい実験的セットアップを提案するものではなく、熱輸送の数値研究において、過渡的なクロスオーバー挙動と真の漸近的スケーリングを区別することの重要性を強調している。