| 查看: 2730 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
拍拍熊出隐刀木虫 (小有名气)
|
[求助]
Vasp+Boltztrap计算自旋极化材料求助
|
||
|
有没有人了解针对自旋极化材料，怎么从vasp的输出文件转成boltztrap的输入文件？ 我查看了下vasp2boltz.py的代码，发现他在从EIGENVAL文件中已经分别读取了自旋上下各K点的能量值（energies_up和energies_dn），转换不了的原因就在于写入energies.boltztrap文件的相关代码没有： if bandstructure['nspin']!=1: print 'write_bandstructure_boltztrap: No idea what to do for nspin=%d' % bandstructure['nspin'] 所以，想问下以往做过自旋极化材料的boltztrap计算的朋友，其energies.boltztrap文件的格式是怎样的？还有网上说的需要做能量合成又是怎么一回事？ 新手，还请大家多多指教。 |
回复此楼» 收录本帖的淘帖专辑推荐
 Boltztrap专辑 |
» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP+BoltzTrap计算Al的电输运性质问题
已经有3人回复
- VASP与BolTzTrap的技术问题
已经有2人回复
- vasp2boltz.py运行出错
已经有4人回复
- vasp结果脚本转换后进行boltztrap计算，结果分析
已经有6人回复
- vasp+新版BoltzTraP计算出错
已经有6人回复
8楼2017-09-22 16:32:12
拍拍熊出隐刀
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4502.1
- 散金: 1
- 红花: 3
- 帖子: 176
- 在线: 256.5小时
- 虫号: 1238512
- 注册: 2011-03-19
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2015-11-12 16:02:34
拍拍熊出隐刀
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4502.1
- 散金: 1
- 红花: 3
- 帖子: 176
- 在线: 256.5小时
- 虫号: 1238512
- 注册: 2011-03-19
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
|
自己给自己砌砖。 我在vasp2boltz.py文件中write_bandstructure_boltztrap函数中添加了关于自旋极化情况下energy文件的转换，大家看看有没有问题啊！ if bandstructure['nspin']!=1: n_up=len(enk['energies_up']) n_dn=len(enk['energies_dn']) n_tot=n_up+n_dn e_tot=enk['energies_up']+enk['energies_dn'] e_tot.sort() fout.write("%12.8f %12.8f %12.8f %d\n" %(kp[0],kp[1],kp[2],n_tot)) line='' for i in range(n_tot): ener=e_tot if (E_Fermi_zero==True): ener=ener-bandstructure['E_Fermi'] ener=ener*yscale fout.write("%18.8f\n" %ener) 实测是可以生成energy文件，换名为XX.energyso后，输入命令加上-so项也是可以计算。但就不知道这样正不正确，还请大家多多指教。 |
3楼2015-11-13 08:40:01
4楼2015-11-23 18:23:24
40