应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 为什么我利用ShengBTE计算出来的热导率比实验的热导率低了一个数量级？
51/1返回列表
查看: 2451  |  回复: 12
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者Mr.linzh将赠送您 10 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Mr.linzh

铁虫 (初入文坛)

[求助] 为什么我利用ShengBTE计算出来的热导率比实验的热导率低了一个数量级？

我在利用 vasp + phonopy + thirdorder + ShengBTE 计算MoO3体块的热导率时，发现算出来的热导率比实验报道的热导率低了一个数量级

计算热导率所需的3个文件，1是CONTROL文件，2是2阶力常数，3是3阶力常数

2阶力常数是通过phonopy对优化好的原胞进行2*2*2扩胞，并通过vasp利用DFPT方法进行计算，最后用phonopy收集得到的。
（原胞的优化是通过vasp利用晶格弛豫（LR）的方法优化的）

3阶力常数是通过thirdorder_vasp.py对优化好的单胞扩胞2*2*2，选取第一近邻，并提交给vasp利用静态计算（ST）计算，对结果提取得到的3阶力常数。

CONTROL文件如下所示，设置了convergence=.TRUE.，并考虑上了BORN有效电荷和介电张量

&allocations
        nelements=2,
        natoms=16,
        ngrid(=4 4 1
&end
&crystal
        lfactor=1.0,
        lattvec(:,1)=3.7131244448375376    0.0000000000000000    0.0000000000000000,
        lattvec(:,2)=0.0000000000000000    3.9306023619198029    0.0000000000000000,
        lattvec(:,3)=0.0000000000000000    0.0000000000000000   14.3152604185205998,
        types(=1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2,
        elements(="Mo" "O"
        positions(:,1)=0.2500000000000000  0.9250711539947374  0.6018119741466632,
        positions(:,2)=0.7500000000000000  0.4250711539947379  0.8981880258533368,
        positions(:,3)=0.2500000000000000  0.5749288460052626  0.1018119441466677,
        positions(:,4)=0.7500000000000000  0.0749288460052624  0.3981880558533326,
        positions(:,5)=0.7500000000000000  0.9815308803785459  0.9129148960037973,
        positions(:,6)=0.2500000000000000  0.4815308803785455  0.5870851039962027,
        positions(:,7)=0.2500000000000000  0.0184691276214550  0.0870850889962059,
        positions(:,8)=0.7500000000000000  0.5184691196214541  0.4129148960037964,
        positions(:,9)=0.7500000000000000  0.4676371119733567  0.7813714805663621,
        positions(:,10)=0.2500000000000000  0.9676370819733541  0.7186285194336379,
        positions(:,11)=0.7500000000000000  0.0323628990266407  0.2813715105663571,
        positions(:,12)=0.2500000000000000  0.5323629180266459  0.2186285044336402,
        positions(:,13)=0.2500000000000000  0.4991321354305162  0.9350443894145257,
        positions(:,14)=0.2500000000000000  0.0008678675694805  0.4350444194145213,
        positions(:,15)=0.7500000000000000  0.5008678945694864  0.0649555875854757,
        positions(:,16)=0.7500000000000000  0.9991321054305136  0.5649556105854743,
        epsilon(:,1)=6.228409     0.000000     0.000000,
        epsilon(:,2)=-0.000000     5.509955     0.000000,
        epsilon(:,3)=0.000000     0.000000     4.126402,
        born(:,1,1)=   7.83082     0.00362    -0.00247,
        born(:,2,1)=  -0.00000     6.64552     0.26868,
        born(:,3,1)=   0.00000     0.54415     4.38384,
        born(:,1,2)=   7.83082     0.00362     0.00247,
        born(:,2,2)=  -0.00000     6.64552    -0.26868,
        born(:,3,2)=  -0.00000    -0.54415     4.38384,
        born(:,1,3)=   7.83082    -0.00362    -0.00247,
        born(:,2,3)=   0.00000     6.64191    -0.27171,
        born(:,3,3)=  -0.00000    -0.54415     4.38384,
        born(:,1,4)=   7.83082    -0.00362     0.00247,
        born(:,2,4)=  -0.00000     6.64191     0.27171,
        born(:,3,4)=  -0.00000     0.54415     4.38384,
        born(:,1,5)=  -1.17320    -0.00152     0.00030,
        born(:,2,5)=   0.00000    -4.87569     0.30130,
        born(:,3,5)=  -0.00000     0.45355    -0.63001,
        born(:,1,6)=  -1.17320    -0.00152    -0.00030,
        born(:,2,6)=  -0.00000    -4.87569    -0.30130,
        born(:,3,6)=  -0.00000    -0.45355    -0.63001,
        born(:,1,7)=  -1.17320     0.00152    -0.00030,
        born(:,2,7)=   0.00000    -4.87511     0.30517,
        born(:,3,7)=   0.00000     0.45355    -0.63001,
        born(:,1,8)=  -1.17320     0.00152     0.00030,
        born(:,2,8)=   0.00000    -4.87511    -0.30517,
        born(:,3,8)=  -0.00000    -0.45355    -0.63001,
        born(:,1,9)=  -0.54685    -0.00040    -0.00017,
        born(:,2,9)=   0.00000    -0.62803     0.35297,
        born(:,3,9)=   0.00000     0.30618    -2.18269,
        born(:,1,10)=  -0.54685    -0.00040     0.00017,
        born(:,2,10)=   0.00000    -0.62803    -0.35297,
        born(:,3,10)=   0.00000    -0.30618    -2.18269,
        born(:,1,11)=  -0.54685     0.00040    -0.00017,
        born(:,2,11)=   0.00000    -0.62655    -0.35523,
        born(:,3,11)=   0.00000    -0.30618    -2.18269,
        born(:,1,12)=  -0.54685     0.00040     0.00017,
        born(:,2,12)=  -0.00000    -0.62655     0.35523,
        born(:,3,12)=  -0.00000     0.30618    -2.18269,
        born(:,1,13)=  -6.11077    -0.00257    -0.00252,
        born(:,2,13)=  -0.00000    -1.14103    -0.33734,
        born(:,3,13)=   0.00000    -0.19451    -1.57114,
        born(:,1,14)=  -6.11077     0.00257    -0.00252,
        born(:,2,14)=   0.00000    -1.14102     0.33865,
        born(:,3,14)=   0.00000     0.19451    -1.57114,
        born(:,1,15)=  -6.11077     0.00257     0.00252,
        born(:,2,15)=  -0.00000    -1.14102    -0.33865,
        born(:,3,15)=  -0.00000    -0.19451    -1.57115,
        born(:,1,16)=  -6.11077    -0.00257     0.00252,
        born(:,2,16)=   0.00000    -1.14103     0.33734,
        born(:,3,16)=   0.00000     0.19451    -1.57115,
        scell(=2 2 2
&end
&parameters
        T_min=20
        T_max=400
        T_step=10
&end
&flags
        convergence=.TRUE.
        nonanalytic=.TRUE.
        isotopes=.FALSE.
        nanowires=.FALSE.
&end

楼主还尝试过扩胞3*2*2来计算2阶力常数，也尝试过第二近邻来计算3阶力常数，都没有改变热导率低一个数量级的结果

有没有哪位大神能够指导一下，计算过程中的哪个步骤可能会导致最终热导率过低的结果？有没有合理的改进建议？
回复此楼

» 猜你喜欢
1楼2020-08-06 14:51:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr.linzh

铁虫 (初入文坛)
引用回帖:
4楼: Originally posted by 123dingyc at 2020-08-06 23:39:08
至少5，到10紧邻

如果计算三阶力常数的时候算到了5近邻的话，那扩胞的大小需要增大吗？
赞一下 回复此楼
5楼2020-08-11 09:30:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr.linzh

铁虫 (初入文坛)
问题找到了，主要是Control文件里的晶格常数单位和POSCAR文件里的晶格常数单位不一致，差了一个数量级，导致了热导率计算失误
赞一下 回复此楼
7楼2020-09-15 15:24:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr.linzh

铁虫 (初入文坛)
引用回帖:
9楼: Originally posted by Sahira at 2020-09-20 09:09:03
是不是有个单位转换应该填0.1那里？
...

是的
回复此楼
11楼2021-02-25 10:50:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Mr.linzh 的主题更新
51/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 入职高校3年发表10+SCI，尽人事听天命 +6 kaoyan250 2024-05-27 6/300 2024-05-28 06:21 by gloomy6159
[硕博家园] 每到中夜，情难自抑 +3 sioc-sunj 2024-05-28 3/150 2024-05-28 06:19 by 夏天的醚
[考博] 24或25申博 +4 Jacob- 2024-05-22 7/350 2024-05-27 21:21 by 安塔瓦拉多
[基金申请] 基金上会 +13 mrKiller 2024-05-25 18/900 2024-05-27 19:51 by 阿漠先生
[基金申请] 面上基金会评专家，有回避机制吗？ +4 huang1991js 2024-05-27 4/200 2024-05-27 19:08 by 星火12
[教师之家] 女博士高校择业三天之内签合同，求支招 +38 chengmy19 2024-05-23 55/2750 2024-05-27 18:02 by 好男人￥
[教师之家] 青椒做实验吗 +3 �徐凤年 2024-05-22 3/150 2024-05-27 16:25 by otani
[硕博家园] 博士复试，申请成绩复核，有机会翻盘吗？ +23 长海二声笑 2024-05-21 31/1550 2024-05-27 12:54 by bingyulin3
[基金申请] 感觉自然基金限制通过比例就是有点扯，学学B口，化学学部，不限制比例。 +10 wsjing 2024-05-26 14/700 2024-05-27 11:57 by kanmiaolucky
[硕博家园] 求助 +3 单增李斯特21633 2024-05-25 3/150 2024-05-27 10:33 by hahamyid
[硕博家园] 周日 +6 1加油哦棒 2024-05-26 9/450 2024-05-27 10:30 by hahamyid
[材料综合] 求助，刚入门MOF合成，实验总是失败 20+3 张月半 2024-05-25 7/350 2024-05-27 08:48 by mkpooooo
[考博] 24年还有机会申博吗，求博导收留 +3 鄙人康康 2024-05-23 6/300 2024-05-27 00:29 by 鄙人康康
[硕博家园] 答辩 +14 暮色恋伊人 2024-05-22 14/700 2024-05-26 23:27 by 意缱绻·
[硕博家园] 2024博士招生 +3 big 混子 2024-05-26 3/150 2024-05-26 20:47 by 宁多缺毋滥
[教师之家] 经常觉得挺累的 +13 zylfront 2024-05-22 19/950 2024-05-25 21:08 by hjc404
[有机交流] TsCl保护羟基为什么不反应 +9 853015158 2024-05-21 29/1450 2024-05-24 19:53 by 谢志新
[考博] 大扩招 +13 version211 2024-05-21 19/950 2024-05-23 20:46 by keyaner23
[论文投稿] 期刊投稿进度 +6 jianhuang9 2024-05-21 8/400 2024-05-22 18:29 by sakuraai
[考博] 化工三篇sci只想有书读，急寻博导！ +3 8139173 2024-05-21 5/250 2024-05-22 10:56 by 安塔瓦拉多
信息提示
关闭
请填处理意见
关闭