[論文レビュー] Efficient simulation of multidimensional phonon transport using energy-based variance-reduced Monte Carlo formulations

Jean-Philippe M. Péraud, Nicolas G. Hadjiconstantinou|DSpace@MIT (Massachusetts Institute of Technology)|Sep 18, 2011

Thermal properties of materials参考文献 28被引用数 18

ひとこと要約

本稿では、マイクロおよびナノスケールデバイスにおける多次元フォノン輸送の効率的かつ正確なシミュレーションを実現する、エネルギーに基づく分散低減モンテカルロ法を提案する。平衡状態と非平衡状態の寄与を分離する制御変数形式を用い、個々のフォノンではなくエネルギーパケットをシミュレートすることで、正確なエネルギー保存が可能となり、統計的不確実性が著しく低減される。これにより、計算コストを一定に保ちながら、任意に小さな温度信号のシミュレーションが可能となる。その有効性は多孔質シリコンおよび熱伝導度スペクトロスコピーの研究を通じて実証された。

ABSTRACT

We present a new Monte Carlo method for obtaining solutions of the Boltzmann equation for describing phonon transport in micro and nanoscale devices. The proposed method can resolve arbitrarily small signals (e.g. temperature differences) at small constant cost and thus represents a considerable improvement compared to traditional Monte Carlo methods whose cost increases quadratically with decreasing signal. This is achieved via a control-variate variance reduction formulation in which the stochastic particle description only solves for the deviation from a nearby equilibrium, while the latter is described analytically. We also show that simulating an energy-based Boltzmann equation results in an algorithm that lends itself naturally to exact energy conservation thereby considerably improving the simulation fidelity. Simulations using the proposed method are used to investigate the effect of porosity on the effective thermal conductivity of silicon. We also present simulations of a recently developed thermal conductivity spectroscopy process. The latter simulations demonstrate how the computational gains introduced by the proposed method enable the simulation of otherwise intractable multiscale phenomena.

研究の動機と目的

• フォノン輸送における微小な温度信号をシミュレートする際の従来のモンテカルロ法の高い計算コストを低減すること。
• エネルギーに基づくボルツマン方程式を用いることで、正確なエネルギー保存を実現し、シミュレーションの忠実性を向上させること。
• 制御変数を用いた分散低減技術を導入することで、低信号領域における統計的不確実性を低減すること。
• 従来の手法では取り扱いが困難であった複雑なマルチスケールフォノン輸送現象をシミュレート可能にする。
• 多孔質シリコンおよび熱伝導度スペクトロスコピーへの応用を通じて、本手法の妥当性を検証すること。

提案手法

• 平衡状態は解析的に既知であり、非平衡状態の偏差のみを確率的にシミュレートする制御変数形式を採用する。
• 個々のフォノンではなくエネルギーパケットを記述するエネルギーに基づくボルツマン方程式を用い、正確なエネルギー保存を実現する。
• 非平衡成分のみを扱うデバイショナルモンテカルロ法により、従来手法で信号サイズの二乗に比例して増大する統計的ノイズを低減する。
• 微小な信号は平衡からの微小なずれに対応することに着目し、そのずれの低分散推定が可能となる。
• デバイュ・モデルおよびアインシュタインモデルに基づく周波数および偏光依存の群速度および緩和時間に基づくアルゴリズムを実装する。
• 一時的境界条件を組み込み、ヘヴィサイド関数を用いて空間的・時間的フォノン放出履歴を追跡する。

実験結果

リサーチクエスチョン

• RQ1 フォノン輸送のモンテカルロシミュレーションを、任意に小さな温度信号に対しても効率的に実行するにはどうすればよいか？
• RQ2 エネルギーパケットに基づくボルツマン方程式の定式化により、エネルギー保存の正確性と統計的精度はどの程度向上するか？
• RQ3 制御変数を用いた分散低減は、低信号フォノン輸送問題における計算コストをどの程度低減できるか？
• RQ4 多孔質性はナノスケールのシリコンにおける有効熱伝導度にどのように影響を与えるか？
• RQ5 提案手法は、従来では困難であった複雑で一時的かつマルチスケールの熱的現象、例えば熱伝導度スペクトロスコピーをシミュレートできるか？

主な発見

• 提案手法は、任意に小さな温度信号に対しても計算コストが一定であるのに対し、従来のモンテカルロ法では信号サイズの二乗に比例してコストが増加する。
• エネルギーに基づく定式化により正確なエネルギー保存が実現され、従来の手法で手動でエネルギー保存を補正していた方法と比較して、シミュレーションの忠実性が顕著に向上した。
• 多孔質シリコンのシミュレーションでは、配列された空孔がランダムに分布した場合よりも有効熱伝導度をより顕著に低減することが判明した。これは、バルビンフォノン散乱が強化されるためである。
• 本手法は、計算効率の高さのおかげで、従来では取り扱いが困難だった熱伝導度スペクトロスコピーの過渡応答を成功裏にシミュレートした。
• デバイショナル定式化により、信号サイズが小さくなるに従い統計的不確実性が低減され、微小な熱勾配に対しても高精度な結果が得られた。
• 一時的1次元バルビン輸送の解析解は、微小な温度変化の仮定のもとで、ΔT(x,t) = ½(1−x/Vgt)H(1−x/Vgt)ΔTl + ½(1−(L−x)/Vgt)H(1−(L−x)/Vgt)ΔTr と導出された式と一致した。