[論文レビュー] Unearthing the foundational role of anharmonicity in heat transport in glasses

Alfredo Fiorentino, Enrico Drigo+2 more|arXiv (Cornell University)|Jul 14, 2023

Thermal properties of materialsMaterials Science被引用数 3

ひとこと要約

この論文は、ガラスにおける熱伝導のアレン＝フェルドマン理論における調和近似が、赤外特異性によって本質的に発散する熱伝導率を予測することを示している。この特異性は、非調和項の効果または有限サイズのカットオフによってのみ正則化される。準調和グリーン＝キューブ理論とヘイドックの再帰的手法を用いて、著者らは、結晶およびガラス状物質の両方において、有限の熱伝導率を得るには非調和項が不可欠であることを示した。これにより、理論的欠陥にもかかわらず、実験的妥当性を示し続けるこの理論の長年のパラドックスが解消された。

ABSTRACT

The time-honored Allen-Feldman theory of heat transport in glasses is generally assumed to predict a finite value for the thermal conductivity, even if it neglects the anharmonic broadening of vibrational normal modes. We demonstrate that the harmonic approximation predicts that the bulk lattice thermal conductivity of harmonic solids inevitably diverges at any temperature, irrespective of configurational disorder, and that its ability to represent the heat-transport properties observed experimentally in most glasses is implicitly due to finite-size effects. Our theoretical analysis is thoroughly benchmarked against careful numerical simulations. Our findings thus reveal that a proper account of anharmonic effects is indispensable to predict a finite value for the bulk thermal conductivity in any solid material, be it crystalline or glassy.

研究の動機と目的

• アレン＝フェルドマン（AF）理論の調和近似が発散する熱伝導率を予測するにもかかわらず、その理論的欠陥にもかかわらず実験的妥当性を示し続ける長年のパラドックスを解消すること。
• 不規則な固体における熱輸送の調和理論が直面する赤外特異性を正則化する非調和項の役割を調査すること。
• シミュレーションにおける有限サイズ効果および界面散乱が、物理的正則化機構をどれほど模倣しているかを明確にすること。
• 大規模なガラスモデルに対してヘイドックの再帰的手法を用いた数値シミュレーションと照らし合わせ、調和的および準調和的アプローチをベンチマークすること。

提案手法

• 著者らは、準調和グリーン＝キューブ（QHGK）理論を用いて、非調和項の効果を摂動的に調和的AFフレームワークに組み込み、量子トンネル効果や人工的なカットオフに依存せずに赤外発散を正則化できるようにした。
• 大規模なガラスモデル（最大13,824原子）に対して、振動ダイナミカル構造因子（VDSF）を計算するためにヘイドックの再帰的アルゴリズムを用い、O(kN)スケーリングで調和的幅および音波減衰を効率的に評価した。
• 非晶質Si、SiO₂、および4H-SiCに対して有限サイズスケーリング解析を実施し、調和的極限における熱伝導率の発散を分離し、バルク極限への収束を評価した。
• VDSFをデバイモデルにフィットさせ、低周波数成分を抽出し、得られた熱伝導率を調和的および準調和的領域の両方でアレン＝フェルドマン式（式1）を用いて評価した。
• 直接対角化法とヘイドック法の結果を比較することで、大規模系における再帰的アプローチの正確性を検証した。
• 数値的安定性とVDSFおよび幅の計算の信頼性を保証するため、リャンツ再帰ステップ数（最大600）を段階的に増加させた収束テストを体系的に行った。

実験結果

リサーチクエスチョン

• RQ1 ガラスにおけるアレン＝フェルドマンの調和的熱輸送理論は、配置の不規則性にかかわらず、レイリー的ω⁴減衰に起因する赤外特異性により発散する熱伝導率を予測するにもかかわらず、実際のシミュレーションでは有限の値を示すのはなぜか？
• RQ2 シミュレーションにおける有限サイズ効果および界面散乱は、調和的発散をどれほど人工的に正則化しており、物理的非調和項効果と比較してどの程度の類似性を示すか？
• RQ3 量子トンネル効果や任意のカットオフに依存せずに、準調和グリーン＝キューブアプローチがバルク極限で赤外特異性を効果的に正則化できるか？
• RQ4 非調和項を適切に取り入れた場合、周波数がω_min未満のモードが熱伝導率にどの程度の寄与を示すか？
• RQ5 シミュレーションで一般的に用いられるスメアリング手順は、薄膜実験で観察される実際の界面散乱効果をどの程度正確に模倣しているか？

主な発見

• 調和的アレン＝フェルドマン理論は、配置の不規則性にかかわらず、レイリー的ω⁴減衰に起因する赤外特異性により、バルク極限で無限大の熱伝導率を予測する。
• シミュレーションにおける有限サイズ効果は、自然な低周波数カットオフ（ω_min ~ 2πc/L）を導入し、発散を人工的に正則化しており、調和的理論が有限モデルで成功を収める背後にある要因を説明している。
• 準調和グリーン＝キューブ法により、非調和項を適切に取り入れた場合、バルク極限においても赤外特異性が正則化され、量子トンネル効果や任意のカットオフを必要としない。
• 非調和項を含めた場合、ω_min未満のモードの寄与は比較的小さく、調和的発散は実際の物質では物理的に意味を持たないことが示された。
• シミュレーションで一般的に用いられる数値的スメアリング手順は、薄膜実験で観察される実際の界面散乱効果を効果的に模倣しており、隠れた物理的対応関係を示している。
• ヘイドックの再帰的手法により、大規模なガラスモデル（最大13,824原子）におけるVDSFおよび幅の正確な計算が可能となり、直接対角化法との比較により数値的安定性および収束性が確認された。