24小时热门版块排行榜    

查看: 10456  |  回复: 57
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主

[交流] 【活动】结合WIEN2K探讨Boltztrap软件的安装编译及正确使用方法

  结合WIEN2K探讨Boltztrap软件的安装编译及正确使用方法


关于半经典输运系数计算,Boltztrap无疑使用非常广泛,下面我就简单探讨式说明下。

1.安装说明
参照帖子:http://muchong.com/html/201210/5100314.html给出的MAKEFILE文件内容,作如下说明
首先我们在linux系统中必须安装编译器及MKL库文件,现在以ifort为例说明:

Makefile文件中编译器及库路径设置必须要正确,
FC = ifort    #设置编译器为ifort                                                                                       
FOPT = -FR -mp1 -w -prec_div -pc80 -pad -ip -DINTEL_VML  #这里选择默认
LDFLAGS = -L/opt/intel/Compiler/11.0/074/mkl/lib/em64t -Vaxlib  -pthread  #这里/opt/intel/Compiler/11.0/074/mkl/lib/em64t是mkl库的路径要在你安装的机器上找到正确的 路径                                 
LIBS = -lmkl_lapack -lmkl_ia32 -lguide    #这里可以选择默认                                                                                                                                                                     
之后保存makefile文件即可,在终端make没有报错即成功。

注:关于编译器的选择,根据自己需要安装,一般我们使用intel公司出的免费软件即可,不过自13年后ifort和mkl似乎合并在一起了,导致很多人安装会出现各种问题,建议还是用11或12版分开的比较好。
2.使用说明
在使用Boltztrap计算时,需要依赖诸如WIEN2k、VASP提供输入文件,以WIEN2K为例说明之:
1)在使用时需要利用WIEN2K通过自洽计算生成分能量文件合成总的能量文件,这里是要讨论的重点,这里如果合成错误,将导致结果失败。大家可以研究下gather_energy.pl这个文件。
2)intrans也是输入的文件之一,这里面有费米能,温度,计算积分切断大小,浓度掺杂等设置,具体可参照http://arxiv.org/pdf/cond-mat/0602203.pdf所述。
3)结构文件使用计算的struct文件即可
4)准备好文件之后,有两种命令可供使用,一是/Boltztrap Boltztrap.def 这个命令需要 Boltztrap.def 文件,可从自带test copy,另外是/ xtrans Boltztrap  后面可加-T延时参数等
3.生成文件
一般我们需要的数据都会在case.trace中生成,如果有掺杂,可以在case.trace_fix_dopping文件中找到数据。还有霍尔系数,热导等均可找到。
4.其他问题等待大家探讨,欢迎大家讨论!

看到大家比较热情参与,现在追加悬赏啦,如果虫友提出或者分享一条主题相关问题或建设性意见,会得到5-20金币,金币不够还会追加,谢谢参与!


[ Last edited by ljw4010 on 2014-11-1 at 20:26 ]
回复此楼
不要说话
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ted1277

木虫 (小有名气)

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
ljw4010: 金币+5, 谢谢交流! 2015-05-06 21:20:25
ljw4010: 回帖置顶 2015-05-06 21:20:28
关于vasp+boltztrap

现在的最新的BoltzTraP 1.2.5已经很支持vasp接口了,
在压缩包“util”文件夹里有vasp2boltztrap.py和vasp2boltztrap.txt两个文件
在系统上已经安装好ase(https://wiki.fysik.dtu.dk/ase/)的基础上

根据需要修改vasp2boltztrap.txt中的
# Example 1:
#sg = Spacegroup(1)
#ao.info = {'spacegroup': sg}

# Example 2. Here the unit cell in the POSCAR file is cubic (conventional):
sg = Spacegroup(122)
ao.info = {'spacegroup': sg, 'unit_cell': 'conventional'}

并将vasp2boltztrap.txt 改名为a.py
执行命令python a.py
显示
********
None

输出的 energies.boltztrap,struct,intrans,可以用于boltztrap计算
hello,boy~
51楼2015-02-02 22:30:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 58 个回答
★
ljw4010(金币+3): 谢谢参与
7楼2014-10-30 23:54:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[职场人生] 40+中年人失业后,还能找到好工作吗? +11 xqy2009 2026-10-13 13/650 2026-10-15 21:32 by bear2007
[有机交流] 用NBS溴化时溶剂除了CCl4外还有什么 +3 18834002860 2026-10-12 6/300 2026-10-15 14:30 by 长沙摩尔信息科�
[基金申请] 申请科研项目纯靠实力能成功吗? 200+7 David1982 2026-10-13 11/550 2026-10-15 10:30 by zyqchem
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
信息提示
请填处理意见