如何推导量子修正方法对模拟温度进行修正的公式
我采用分子动力学的软件lammps对氮化硼求导热导率,后发现计算结果偏大,想通过量子化修正对结果进行后处理。参考石墨烯类似计算,发现通常假设模拟温度下声子总能量等于量子修正温度下的声子总能量。该假设由公式表达为图1所示,最后推导为图2所示公式。请问有做过量子化修正的朋友吗?是如何推导的呢?没有学过量子化学感觉好吃力,感激不尽。
图1
图2 返回小木虫查看更多
今日热帖
我采用分子动力学的软件lammps对氮化硼求导热导率,后发现计算结果偏大,想通过量子化修正对结果进行后处理。参考石墨烯类似计算,发现通常假设模拟温度下声子总能量等于量子修正温度下的声子总能量。该假设由公式表达为图1所示,最后推导为图2所示公式。请问有做过量子化修正的朋友吗?是如何推导的呢?没有学过量子化学感觉好吃力,感激不尽。
图1
图2 返回小木虫查看更多
这要用到级联方程修正了
嗨,这不是量子化学,是固体物理中的晶格动力学。这种修正或许简单化,而且不是必须的。经典MD算的结果还是蛮准的,我算过室温下无限大系统的热导率,应该是600W/mK左右,与晶格动力学方法的结果差不多。
好的,谢谢樊老师,但是经典MD和晶格MD具体是什么我不明白,想请教您如何判断经典MD和晶格MD呢
,
推导过程已经解决了,需要用到你计算材料的声子态密度和声子态波色方程,这两个方程是关于声子频率、修正温度、声子群速度的函数,运用数学中积分换元替换,巧妙的将包含声子频率、声子群速度在内的修正项化简成德拜温度,最终得到最后的公式。
请问能不能用ms做氮化硼热导率?