vasp boltztrapp计算出来的数据怎么处理

小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
关于BoltzTrap-1.2.5中vasp接口程序,问题很严重啊!
我用BoltzTrap中的util目录下的vasp2boltz.py, vasp2boltz.txt的两个脚本程序转化vasp计算结果,目的是为了计算石墨烯的电导。
目前看来问题很多啊,首先vasp2boltz.txt提示如果没有装spglib要手动输入对称性,我装了sgblib,并且没有安装报错信息,结果运行脚本vasp2boltz时总是提示找不到对称性,无奈之下,我输入了对称性,
如sg = Spacegroup(6)
ao.info = {'spacegroup': sg, 'unit_cell': conventional'}
结果出来三个文件,hte.struct,hte.intrans,energies_boltrap,可是hte.struct里面的原子数目和坐标明显不对,
我的POSCAR文件是
&&Graphene
& & 14.& &0.0& & 0.0000
& & -7.& &12.& & 0.0000
& &&&0.0000& & 0.
。。。。。。。。。。。。。。
输出hte.struct
HTE output
& & 28.00574& &&&0.00000& &&&0.00000
& &-14.00287& & 24.25368& &&&0.00000
& &&&0.00000& &&&0.00000& & 32.85013
然后运行boltztrap后报错
&&1&&0&&0& && && &1
&&1&&0&&0& && && &2
error in factorization
0.003u 0.002s 0:00.00 0.0%& & & & 0+0k 0+16io 0pf+0w
2.输出文件hte.struct 中首行 GENE& && && & # use generic interface
可是这个怎么能转化成vasp界面呢?或者不许要改变这行?
3.BoltzTrap中关于电导给出的计算量是sigma(电导)/tau(电子散射时间),这就比较纠结了,总部能拿这个数据发表文章吧?至少要电导的数据,难道还要用phonopy计算声子频率,然后计算电子散射时间?再带入计算电导么?有没有更好的方法?
4.感觉BoltzTrap-1.2.5有问题啊,我的程序编译完之后,运行之后就呆在哪里,一边占用内存和cpu,一边多少个小时都不输出数据,感觉肯定程序有问题啊。
我用过BoltzTrap-1.2.0的,很好用,很快就计算完了。
摆脱各位大侠能够帮助一下,哪怕一条也行,为了用boltztrap计算电导,我装了numpy, ase, spglib,python等不下六七个程序,可是总是伺候不好它,很纠结。
单机,程序有问题,我现在还是用的老版的。
这个Boltrap没法算,要想计算估计需要计算声子散射频率,用phonopy计算吧。
程序中给出的Al的例子是用wien2k算的,那个我测试过,没问题,但是用vasp算的没有测试过。我认为你把vasp incar中设置ismear=-5, encut=400,ediff=0.00001,应该差不多,否则的话就是vasp计算的态密度与wien2k的差别较大引起来的,客观的说,我不太看好boltztrap计算结果,因为我用它算过石墨烯的电导,好像起伏较大,有点类似dos图像
用 vasp和boltztrap相接的时候,用的outcar是自洽计算得到的还是能带计算得到的?
请问是怎么解决的?我也出现了这样的问题。你的意思是用vasp2boltz.txt生成文件而不是vasp2boltz.py生成的吗
你选择用的是什么的呢?静态自洽还是能带计算? 两者的区别大吗?
你好,请问你是怎么解决的,方便指导一下嘛?谢谢,如果能留个qq,请教你一下,就更好了
研究生必备与500万研究生在线互动!
扫描下载送金币小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
求助boltztrap软件的学习资料和安装教程
这个软件找不到学习资料,求助boltztrap软件的学习资料和安装教程,小弟在此感激不尽!
谢谢,我也要一份,
恩恩太谢谢了,我的邮箱
我看了下,百度就有它的说明书,我的资料就是说明书,大家百度下就行了。
百度它的说明书,
你这有点坑啊,我还以为你有什么好的学习资料。。。
说明书谁没有…
无语。。。。
研究生必备与500万研究生在线互动!
扫描下载送金币小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
BoltzTraP如何处理二维或者低维材料载流子浓度?
求助大侠,用过BoltzTraP的虫子们知道,最后结果给出的载流子浓度单位是每个unit cell的载流子浓度,对于三维材料来说,求一下体积转化一下即可得到cm-3的浓度。
但是对于二维材料或者更低维度的材料,这个体积如何定义?对于平面的二维材料或者有褶皱的类二维材料如何解决?
有没有知道的大侠们啊,求解答!!
那平面结构不是范德瓦尔斯厚度而是0啊?
一般2D结构可以看成是从相应的块体结构剥离的,那么这个厚度就是对应的块体结构中的范德瓦尔斯厚度,例如计算石墨烯的厚度的时候,可以认为石墨烯的厚度是3.35埃。
哦哦,谢谢,明白了。
你好,如果像硅烯这类没有对应层状结构的时候,z方向的厚度应该怎么取呢?求指教
这个也是一样,你用考虑vdW的泛函优化两层(硅烯最优的堆垛方式应该是AB)硅烯结构,所求的厚度就是一层原子的厚度加上两层之间的vdW厚度
好的,谢谢啊!
我表示质疑
单位应该换一下 变成 cm-2
而不是硬给它加一个维度
它就是2维的
研究生必备与500万研究生在线互动!
扫描下载送金币小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
WIEN2k+BolTzTraP做自旋极化时,自洽总能合成问题请教!
以前没做过自旋极化的计算,现在在合成能量上出了问题:
自旋极化计算自洽分能量有case.energydn_??和case.energydn_??的分能量各20个,
1.Boltztrap 用的总能怎么合成?
2.我改进了下gather_energy.pl程序,原始代码是做耦合的代码,我在这里改动了下:
3.改动后下:
4.但是不知这样做行得通不?
prince1000
自旋极化情况下,合并成两个文件energyup和energydn,那然后呢?在Boltztrap计算时,分别计算up和dn的?还是up和dn的还要再合并一下呢?
我看到相关文献,对于磁性的化合物输运结果并没有up和dn两个结果啊,只有一条曲线,里面也没有说怎么考虑的自旋极化的情况。
研究生必备与500万研究生在线互动!
扫描下载送金币

我要回帖

更多关于 boltztrap软件包介绍 的文章

 

随机推荐