[論文レビュー] First-principle quantum Monte-Carlo study of charge carrier mobility in organic molecular semiconductors

Johann Ostmeyer, Tahereh Nematiaram+2 more|arXiv (Cornell University)|Dec 22, 2023

Advancements in Semiconductor Devices and Circuit Design被引用数 5

ひとこと要約

本稿では、調整可能なパラメータなしで格子振動子と電荷キャリアの量子力学的ダイナミクスをシミュレートするためのハイブリッド・モンテカルロ（HMC）アルゴリズムを用いて、2次元有機分子半導体における電荷キャリア移動度について第一原理的量子モンテカルロ研究を提示する。この手法は1CPU時間未満で千分の一未満の精度を達成し、実験的移動度の傾向と一時的局在化モデルの結果を再現し、有機半導体における動的不規則性を研究する上で正確かつ効率的であることを検証する。

ABSTRACT

We present a first-principle numerical study of charge transport in a realistic two-dimensional tight-binding model of organic molecular semiconductors. We use the Hybrid Monte Carlo (HMC) algorithm to simulate the full quantum dynamics of phonons and either a single or multiple charge carriers without any tunable parameters. We introduce a number of algorithmic improvements, including efficient Metropolis updates for phonon fields based on analytic insights, which lead to negligible autocorrelation times and allow to reach sub-permille precisions at small computational cost of $O(1)$ CPU-hour. Our simulations produce charge mobility estimates that are in good agreement with experiment and that also justify the phenomenological Transient Localisation approach.

研究の動機と目的

• 量子格子振動子のダイナミクスを完全に考慮し、古典的核近似を避ける、有機分子半導体における電荷輸送を研究する第一原理的数値手法の開発。
• 従来のモデル（例：一時的局在化モデル）で用いられる緩和時間などの現象論的パrameterに依存しないことの実現。
• 現実的な2次元タイトバインディングモデルにおける電荷キャリア移動度とその温度依存性を、正確かつパラメータフリーに計算可能にする。
• 格子振動子場のHMCサンプリングにおけるアルゴリズム的改善を通じて、低い計算コストで高い統計的精度を達成すること。

提案手法

• 本研究では、三角格子上に2次元タイトバインディングハミルトニアンを用い、電子-格子振動子結合を導入することで、分子的有機半導体における電荷輸送をモデル化する。
• ハイブリッド・モンテカルロ（HMC）アルゴリズムを用いて、格子振動子と電荷キャリアの完全な量子力学的ダイナミクスをシミュレートし、調整可能なパラメータは一切使用しない。
• 解析的知見から導かれた効率的なメトロポリス更新を格子振動子場に適用することで、自己相関時間を短縮し、千分の一未満の精度を実現する。
• 現在-現在相関関数などの物理量を計算するために、熱的分配関数の虚時間経路積分表現を用いる。
• 低キャリア密度条件下では自由格子振動子近似が導かれるため、HMC結果の解析的妥当性を検証可能となる。
• 物理量と重みの間のキャンセル効果により統計誤差を最小化し、単純なサンプリング手法と比較してフラクチュエーションを低減する。

実験結果

リサーチクエスチョン

• RQ1 調整可能なパラメータを一切使用せず、第一原理的量子モンテカルロ手法が有機半導体における電荷キャリア移動度を正確に予測できるか？
• RQ2 HMCアルゴリズムは、動的不規則性を示す2次元有機半導体における量子的格子振動子ダイナミクスと電荷輸送をどの程度正確にシミュレートできるか？
• RQ3 HMC結果は、実験的移動度データおよび一時的局在化モデルの予測とどの程度一致するか？
• RQ4 提示されたHMC手法の計算コストと統計的精度は、現実的なハミルトニアンパラメータと温度範囲においてどの程度か？
• RQ5 格子振動子場の更新におけるアルゴリズム的改善は、自己相関時間と収束性にどのように影響を与えるか？

主な発見

• HMCシミュレーションは約1CPU時間で千分の一未満の統計的精度を達成し、高い計算効率を示した。
• 計算された電荷キャリア移動度の値は、さまざまな温度およびハミルトニアンパラメータにおいて、実験的測定値と良好な定量的整合性を示した。
• 一時的局在化モデルの予測と密接に一致しており、現象論的緩和時間パラメータを必要としない手法の正確性が裏付けられた。
• この手法は、完全に量子力学的な枠組みの中で、格子振動子に起因する遷移積分のフラクチュエーションに起因する動的不規則性の効果を的確に捉えた。
• 解析的メトロポリス更新を含むアルゴリズム的改善により、自己相関時間がほとんどゼロに近くなり、収束が安定した。
• 経路積分における物理量と重みのキャンセル効果により、単純なサンプリング手法と比較して統計誤差が著しく低減した。