[论文解读] 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）算法对二维有机分子半导体中的载流子迁移率进行了研究，该方法在无调整参数的情况下模拟了声子与载流子的量子动力学。该方法在不到一CPU小时内实现了千分之零点一以内的精度，并复现了实验中的迁移率趋势以及瞬态局域化模型的结果，验证了其在研究有机半导体中动态无序问题上的准确性与高效性。

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.

研究动机与目标

• 开发一种用于研究有机分子半导体中载流子输运的第一性原理数值方法，全面考虑量子声子动力学并避免经典核近似。
• 消除现有模型（如瞬态局域化模型）中对经验参数（如弛豫时间）的依赖。
• 实现在真实二维紧束缚模型中有机半导体载流子迁移率及其温度依赖关系的精确、无参数计算。
• 通过HMC采样中的算法改进，实现在低计算成本下获得高统计精度。

提出的方法

• 研究采用三角晶格上的二维紧束缚哈密顿量，引入电子-声子耦合，以模拟分子有机半导体中的载流子输运。
• 使用混合蒙特卡罗（HMC）算法模拟声子与载流子的完整量子动力学，且无可调参数。
• 基于解析洞察推导出高效的声子场Metropolis更新，显著降低自相关时间，实现千分之零点一以内的精度。
• 采用虚时路径积分形式化热力学配分函数，以计算电流-电流关联函数等可观测量。
• 在低载流子密度下推导出自由声子近似，为HMC结果提供解析验证。
• 通过可观测量与权重之间的抵消效应最小化统计误差，相比简单采样方法波动更小。

实验结果

研究问题

• RQ1 是否能够通过第一性原理的量子蒙特卡罗方法，在无需经验参数的情况下准确预测有机半导体中的载流子迁移率？
• RQ2 HMC算法在模拟具有动态无序的二维有机半导体中量子声子动力学与载流子输运方面的表现如何？
• RQ3 HMC结果与实验迁移率数据及瞬态局域化模型的符合程度如何？
• RQ4 所提出的HMC方法在真实哈密顿量参数与温度条件下的计算成本与统计精度如何？
• RQ5 声子场更新中的算法改进对自相关时间与模拟收敛性有何影响？

主要发现

• HMC模拟在约一CPU小时内实现了千分之零点一以内的统计精度，展现出极高的计算效率。
• 计算得到的载流子迁移率在不同温度与哈密顿量参数下与实验测量值具有良好的定量一致性。
• 结果与瞬态局域化模型的预测高度吻合，验证了该方法的准确性，且无需引入经验弛豫时间参数。
• 该方法成功在完全量子力学框架下捕捉了声子介导的跃迁积分涨落所引起的动态无序效应。
• 算法改进（包括解析Metropolis更新）使自相关时间可忽略不计，且收敛稳定。
• 由于路径积分中可观测量与权重之间的抵消效应，HMC方法在统计误差降低方面优于简单采样方法。