24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2510  |  回复: 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专辑

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

拍拍熊出隐刀

木虫 (小有名气)

自己给自己盖楼。
例子中BiTe采用自旋轨道耦合计算,输入文件energyso替换了energy文件,同时输入命令多了一项 ‘-so’,告诉程序每个能态只包含一个电子。
那么考虑自旋极化时可不可以也采用类似做法,直接每个K点的自旋向下的能态能量值并入其自旋向上的文件中去,然后输入命令是加上 '-so'项,告诉程序读入的能态只包含一个电子?
请大家多多指教。
2楼2015-11-12 16:02:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

拍拍熊出隐刀

木虫 (小有名气)

自己给自己砌砖。
我在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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

museld

金虫 (著名写手)

楼主用你的这个重复一下文献的数据看看,对比一下就知道有没有问题了。
123
4楼2015-11-23 18:23:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ouyangyulou

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 拍拍熊出隐刀 at 2015-11-12 16:02:34
自己给自己盖楼。
例子中BiTe采用自旋轨道耦合计算,输入文件energyso替换了energy文件,同时输入命令多了一项 ‘-so’,告诉程序每个能态只包含一个电子。
那么考虑自旋极化时可不可以也采用类似做法,直接每个K ...

高人,您好请问最后你是怎么处理加SOC时,计算热电的呢
创新、思路决定出路
5楼2016-07-31 21:38:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ouyangyulou

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 拍拍熊出隐刀 at 2015-11-12 16:02:34
自己给自己盖楼。
例子中BiTe采用自旋轨道耦合计算,输入文件energyso替换了energy文件,同时输入命令多了一项 ‘-so’,告诉程序每个能态只包含一个电子。
那么考虑自旋极化时可不可以也采用类似做法,直接每个K ...

高人,您好请问最后你是怎么处理加SOC时,计算热电的呢
创新、思路决定出路
6楼2016-07-31 21:40:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyzqmc

银虫 (小有名气)

引用回帖:
3楼: Originally posted by 拍拍熊出隐刀 at 2015-11-13 08:40:01
自己给自己砌砖。
我在vasp2boltz.py文件中write_bandstructure_boltztrap函数中添加了关于自旋极化情况下energy文件的转换,大家看看有没有问题啊!
        if bandstructure!=1:
             n_up=len(enk)
...

楼主,我遇到了同样的问题,这个接口转换ispin=2的能量文件转换不了,您的代码也试了不行,您最后是怎么解决啊啊
7楼2016-11-22 14:08:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

江霞

新虫 (正式写手)

你好你这个问题解决了吗?怎么做呀?可不可以告诉我呀?
8楼2017-09-22 16:32:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 拍拍熊出隐刀 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +11 熊二想上岸 2026-04-04 11/550 2026-04-05 22:21 by 醉翁wl
[考研] 0854电子信息319求调剂(接受跨专业调剂) +3 星星不眨眼喽 2026-04-05 3/150 2026-04-05 20:20 by 啵啵啵0119
[考研] 求调剂 一志愿西南交通大学085701环境工程 282分 +7 多多爱吃汉堡 2026-04-04 7/350 2026-04-05 19:47 by 福农资环_环境基
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 能动调剂326专硕 +4 wan112233 2026-04-04 4/200 2026-04-04 22:47 by yu221
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 材料调剂 +15 一样YWY 2026-04-01 15/750 2026-04-04 22:23 by hemengdong
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +11 努力奋斗112 2026-04-04 11/550 2026-04-04 20:51 by 蓝云思雨
[考研] 277求调剂 +4 12A3 2026-04-02 5/250 2026-04-04 20:28 by 蓝云思雨
[考研] 283分材料与化工求调剂 +29 罗KAKA 2026-04-02 29/1450 2026-04-03 23:56 by userper
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 材料调剂 +4 一样YWY 2026-04-03 4/200 2026-04-03 09:48 by 蓝云思雨
[考研] 295求调剂 +7 愿旅途永远坦然 2026-04-02 7/350 2026-04-03 08:22 by fangshan711
[考研] 调剂 +3 osbbx 2026-04-02 3/150 2026-04-03 07:47 by cc8418
[考研] 286分调剂 +20 Faune 2026-03-30 22/1100 2026-04-02 13:24 by clyblh
[考研] 272求调剂,接受跨专业调剂! +4 闲鱼卢 2026-03-31 4/200 2026-04-02 11:18 by guyan1000
[考研] 302求调剂一志愿北航070300,本科郑大化学 +8 圣日耳曼条 2026-04-01 11/550 2026-04-02 07:40 by chemdavid
[考研] 化学工程专硕324分,一志愿中国矿业大学求调剂 +7 耿耿1314 2026-04-01 7/350 2026-04-02 07:40 by 尚水阁主
[考研] 一志愿食品科学与工程083200求调剂 +4 XQTJZ 2026-03-30 4/200 2026-03-31 04:10 by fmesaito
信息提示
请填处理意见