24小时热门版块排行榜    

查看: 3603  |  回复: 15

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专辑

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

尽自己最大的力!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

ted1277

木虫 (小有名气)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-07-29 14:18:21
gather energy.pl <case>  是把你并行计算的能量整合到一个的命令,
intrans的第三行,分别是费米能级,能量变化尺度,能量在费米能级附近变化范围,电子数目
费米能级,电子数目你要从scf文件里找着看。
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 的主题更新
信息提示
请填处理意见