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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Wien2k&FLAPW&ELK » 使用BoltzTrap计算的Seebeck系数结果有偏差
161/2返回列表12下一页
查看: 4030  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上提前没消息，有中的吗 +10 archvillain 2026-08-02 12/600 2026-08-02 18:46 by hitlend
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 18:44 by ca3DVAx4OB19
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 18:32 by ca3DVAx4OB19
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 9/450 2026-08-02 18:27 by ca3DVAx4OB19
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 9/450 2026-08-02 18:23 by ca3DVAx4OB19
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 13/650 2026-08-02 18:14 by ca3DVAx4OB19
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 AS6fX1tr1on4 2026-08-02 12/600 2026-08-02 18:14 by ca3DVAx4OB19
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 18:11 by ca3DVAx4OB19
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 cBGhfJFuNMca 2026-08-02 20/1000 2026-08-02 18:04 by ca3DVAx4OB19
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 20/1000 2026-08-02 18:04 by ca3DVAx4OB19
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 20/1000 2026-08-02 18:03 by ca3DVAx4OB19
[基金申请] 准备明年的基金了 +7 Tide man 2026-07-27 7/350 2026-08-02 16:43 by gltch
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +18 低垂的野花 2026-07-31 24/1200 2026-08-02 16:03 by swallowwh
[找工作] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 12/600 2026-08-02 13:43 by jHcrHMhbKCFV
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 13:34 by jHcrHMhbKCFV
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
关闭
请填处理意见
关闭