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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Wien2k&FLAPW&ELK » 使用BoltzTrap计算的Seebeck系数结果有偏差
161/2返回列表12下一页
查看: 4100  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiaofeiyu0

木虫 (正式写手)

[求助] 使用BoltzTrap计算的Seebeck系数结果有偏差 已有2人参与

我用BoltzTrap软件计算Mg2Si的seebeck系数。k-point从10*10*10到38*38*38 计算得出的seebeck系数，数量级都与其他文献相差1个数量级。不知道什么原因？

此外，我运算了BoltzTrap软件给的CoSb3测试例子，完全用它的输入文件，运算出来的*.trace文件
#       Ef[Ry] T [K]            N         DOS(Ef)           S             s/t               R_H        kappa0         c                 chi
   0.07494   50.0000     72.00000000  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
   0.07494  100.0000     72.00000000  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
   0.07494  150.0000     72.00000000  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
   0.07494  200.0000     72.00000000  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00  0.00000000E+00
。。。。。。
费米能级处显示0.07而 *.intrans文件中费米能级输入的是0.55475。 不明白是什么决定*.trace文件的费米能级。

自己也按照文献计算了Bi2Te3的seebeck系数发现数量级也是不对。

请大家帮我解答下，非常感谢！
回复此楼

» 收录本帖的淘帖专辑推荐

Boltztrap专辑

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
尽自己最大的力！
1楼 2014-07-07 20:33:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

ted1277

木虫 (小有名气)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-07-29 14:18:21
gather energy.pl <case>  是把你并行计算的能量整合到一个的命令，
intrans的第三行，分别是费米能级，能量变化尺度，能量在费米能级附近变化范围，电子数目
费米能级，电子数目你要从scf文件里找着看。
赞一下(1人) 回复此楼
hello,boy~
15楼2014-07-26 10:12:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

xiaofeiyu0

木虫 (正式写手)
希望用BoltzTrap软件的人，帮我解答一下。这个问题我想了很多办法，都不知道怎么解决。
赞一下 回复此楼
尽自己最大的力！
2楼2014-07-10 08:41:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
xiaofeiyu0: 金币+5 2014-07-11 16:06:36
把你的intrans文件设置贴出来看看，或者找个其他版本的试试。
赞一下 回复此楼
不要说话
3楼2014-07-10 10:29:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaofeiyu0

木虫 (正式写手)
引用回帖:
3楼: Originally posted by ljw4010 at 2014-07-10 10:29:34
把你的intrans文件设置贴出来看看，或者找个其他版本的试试。

intrans:
WIEN                      # Format of DOS                                                                             
0 0 0 0.0                 # iskip (not presently used) idebug setgap shiftgap                                             
0.23921 0.0005 0.15 96.   # Fermilevel (Ry), energygrid, energy span around Fermilevel, number of electrons                     
CALC                    # CALC (calculate expansion coeff), NOCALC read from file                                     
5                         # lpfac, number of latt-points per k-point                                                                
BOLTZ                     # run mode (only BOLTZ is supported)                                                             
.15                       # (efcut) energy range of chemical potential                                                     
1000. 50.                  # Tmax, temperature grid                                                                     
-1.            # energyrange of bands given individual DOS output sig_xxx and dos_xxx (xxx is band number)
HISTO

我试了1.2.5版本  这里面我算了下LiZnSb的例子 这个结果和参考的结果类似，至少数量级是对的。
赞一下 回复此楼
尽自己最大的力！
4楼2014-07-11 16:10:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主
引用回帖:
4楼: Originally posted by xiaofeiyu0 at 2014-07-11 16:10:33
intrans:
WIEN                      # Format of DOS                                                                             
0 0 0 0.0                 # iskip (not presently used) idebug setgap shiftgap                                             
0.23921 0.0005 0.15 96.   # Fermilevel (Ry ...

1.2.5版本比以前版本功能提升了
赞一下 回复此楼
不要说话
5楼2014-07-11 16:14:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaofeiyu0

木虫 (正式写手)
引用回帖:
5楼: Originally posted by ljw4010 at 2014-07-11 16:14:32
1.2.5版本比以前版本功能提升了...

但我计算的Mg2Si数量级还是不对的
赞一下 回复此楼
尽自己最大的力！
6楼2014-07-11 16:26:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主
引用回帖:
6楼: Originally posted by xiaofeiyu0 at 2014-07-11 16:26:26
但我计算的Mg2Si数量级还是不对的...

没道理啊，1.2.5版本，我现在也在用，很正常。你算软件包自带的例子正常，算你的就有问题了，是不是你能量文件有问题，1.2.5版本设置不合理会报错的
赞一下 回复此楼
不要说话
7楼2014-07-11 16:38:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaofeiyu0

木虫 (正式写手)
引用回帖:
7楼: Originally posted by ljw4010 at 2014-07-11 16:38:28
没道理啊，1.2.5版本，我现在也在用，很正常。你算软件包自带的例子正常，算你的就有问题了，是不是你能量文件有问题，1.2.5版本设置不合理会报错的...

我也觉得没道理。能量文件是我用wien2k直接计算出来的。
赞一下 回复此楼
尽自己最大的力！
8楼2014-07-11 17:38:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主
引用回帖:
8楼: Originally posted by xiaofeiyu0 at 2014-07-11 17:38:36
我也觉得没道理。能量文件是我用wien2k直接计算出来的。...

那你用其他的wien2k生成的energy试试
赞一下 回复此楼
不要说话
9楼2014-07-11 18:07:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaofeiyu0

木虫 (正式写手)
引用回帖:
9楼: Originally posted by ljw4010 at 2014-07-11 18:07:35
那你用其他的wien2k生成的energy试试...

我用并行算出的能量文件。你知道怎么把并行文件 生成一个Boltztrap能用的能量文件吗？我怀疑我是这一步弄错了
赞一下 回复此楼
尽自己最大的力！
10楼2014-07-11 22:07:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaofeiyu0 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 00:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 00:40 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 00:28 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 8/400 2026-08-22 23:51 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 10/500 2026-08-22 23:42 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 3/150 2026-08-22 22:03 by OEbVnUOu01ol
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 让我中一个面上吧！ +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭