經(jīng)驗分享

計算材料熱輸運性質(zhì)時需要計算模型原子間力常數(shù)，包括可處理獲得聲子色散的二階力常數(shù)，和三階力常數(shù)甚至四階。通常使用phono3py或者thirdorder生成三階力常數(shù)的計算文件。

在phono3py中可以通過調(diào)整原子位移控制產(chǎn)生的任務(wù)數(shù)量，而在thirdorder中則主要控制臨近數(shù)值。

在有關(guān)資料中有說在thirdorder擴包生成計算任務(wù)時臨近數(shù)盡量取10以上，才能計算得到收斂的晶格熱導(dǎo)率。

詳細(xì)請參考鏈接：

https://www.bilibili.com/video/BV1rL41147Rg/?spm_id_from=333.999.0.0

但對于不同的體系和不同的擴包倍數(shù)來講，在使用thirdorder擴包時會有臨近數(shù)值的上限。

比如筆者使用一個2原子的晶胞（空間群Fm-3m）進行擴胞，擴胞倍數(shù)為2×2×2（此擴包倍數(shù)較小，無法計算得到有效和精確的結(jié)構(gòu)，僅為案例需要），當(dāng)臨近數(shù)值大于5時所產(chǎn)生的位移POSCAR就會達到上限88個，并提示warning: supercell too small to find n-th neighbours，無法找到所設(shè)置的近鄰。后續(xù)即使增大臨近數(shù)值，也不會增加三階力常數(shù)計算的任務(wù)數(shù)量，即該擴胞倍數(shù)下最大臨近數(shù)值為5。

當(dāng)增大擴胞倍數(shù)至4×4×4，當(dāng)臨近數(shù)值達到12才得到達到任務(wù)上限；

當(dāng)增大擴胞倍數(shù)至5×5×5，當(dāng)臨近數(shù)值達到17才得到達到任務(wù)上限；

所以即便是同一體系，擴胞倍數(shù)不同的情況下，使用臨近數(shù)值也需要測試和研究。

在可使用的計算資源有限時，使用合適的臨近數(shù)值是非常有必要的。

愿有所成

IEchoQ

引喻失義? ?妄自菲薄