版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

拍拍熊出隱刀

木蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4502.1
散金: 1
紅花: 3
帖子: 176
在線: 256.5小時(shí)
蟲(chóng)號(hào): 1238512
注冊(cè): 2011-03-19
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[求助] Vasp+Boltztrap計(jì)算自旋極化材料求助

有沒(méi)有人了解針對(duì)自旋極化材料，怎么從vasp的輸出文件轉(zhuǎn)成boltztrap的輸入文件？
我查看了下vasp2boltz.py的代碼，發(fā)現(xiàn)他在從EIGENVAL文件中已經(jīng)分別讀取了自旋上下各K點(diǎn)的能量值（energies_up和energies_dn），轉(zhuǎn)換不了的原因就在于寫(xiě)入energies.boltztrap文件的相關(guān)代碼沒(méi)有：
if bandstructure['nspin']!=1:
print 'write_bandstructure_boltztrap: No idea what to do for nspin=%d' % bandstructure['nspin']
所以，想問(wèn)下以往做過(guò)自旋極化材料的boltztrap計(jì)算的朋友，其energies.boltztrap文件的格式是怎樣的？還有網(wǎng)上說(shuō)的需要做能量合成又是怎么一回事？
新手，還請(qǐng)大家多多指教。

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

Boltztrap專輯

» 猜你喜歡

276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷已經(jīng)有8人回復(fù)
接收2026碩士調(diào)劑(學(xué)碩+專碩) 已經(jīng)有3人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有3人回復(fù)
081700 調(diào)劑 267分已經(jīng)有3人回復(fù)
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑已經(jīng)有3人回復(fù)
336化工調(diào)劑已經(jīng)有4人回復(fù)
08工學(xué)調(diào)劑已經(jīng)有10人回復(fù)
263求調(diào)劑已經(jīng)有6人回復(fù)
298求調(diào)劑已經(jīng)有8人回復(fù)
考研化學(xué)308分求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

VASP+BoltzTrap計(jì)算Al的電輸運(yùn)性質(zhì)問(wèn)題已經(jīng)有3人回復(fù)
VASP與BolTzTrap的技術(shù)問(wèn)題已經(jīng)有2人回復(fù)
vasp2boltz.py運(yùn)行出錯(cuò) 已經(jīng)有4人回復(fù)
vasp結(jié)果腳本轉(zhuǎn)換后進(jìn)行boltztrap計(jì)算，結(jié)果分析已經(jīng)有6人回復(fù)
vasp+新版BoltzTraP計(jì)算出錯(cuò) 已經(jīng)有6人回復(fù)

1樓 2015-11-12 14:40:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

拍拍熊出隱刀

木蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4502.1
散金: 1
紅花: 3
帖子: 176
在線: 256.5小時(shí)
蟲(chóng)號(hào): 1238512
注冊(cè): 2011-03-19
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

自己給自己蓋樓。
例子中BiTe采用自旋軌道耦合計(jì)算，輸入文件energyso替換了energy文件，同時(shí)輸入命令多了一項(xiàng) ‘-so’，告訴程序每個(gè)能態(tài)只包含一個(gè)電子。
那么考慮自旋極化時(shí)可不可以也采用類似做法，直接每個(gè)K點(diǎn)的自旋向下的能態(tài)能量值并入其自旋向上的文件中去，然后輸入命令是加上 '-so'項(xiàng)，告訴程序讀入的能態(tài)只包含一個(gè)電子？
請(qǐng)大家多多指教。

贊一下(1人)

回復(fù)此樓

2樓2015-11-12 16:02:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

拍拍熊出隱刀

木蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 4502.1
散金: 1
紅花: 3
帖子: 176
在線: 256.5小時(shí)
蟲(chóng)號(hào): 1238512
注冊(cè): 2011-03-19
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

自己給自己砌磚。
我在vasp2boltz.py文件中write_bandstructure_boltztrap函數(shù)中添加了關(guān)于自旋極化情況下energy文件的轉(zhuǎn)換，大家看看有沒(méi)有問(wèn)題�。�
      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)
實(shí)測(cè)是可以生成energy文件，換名為XX.energyso后，輸入命令加上-so項(xiàng)也是可以計(jì)算。但就不知道這樣正不正確，還請(qǐng)大家多多指教。

贊一下

回復(fù)此樓

3樓2015-11-13 08:40:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

museld

金蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 7 (幼兒園)
金幣: 813.7
散金: 2415
紅花: 8
帖子: 1719
在線: 451.2小時(shí)
蟲(chóng)號(hào): 1782594
注冊(cè): 2012-04-27
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

樓主用你的這個(gè)重復(fù)一下文獻(xiàn)的數(shù)據(jù)看看，對(duì)比一下就知道有沒(méi)有問(wèn)題了。

贊一下

回復(fù)此樓

123

4樓2015-11-23 18:23:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ouyangyulou

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 254
紅花: 12
帖子: 124
在線: 95.9小時(shí)
蟲(chóng)號(hào): 2850255
注冊(cè): 2013-12-05
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

2樓: Originally posted by 拍拍熊出隱刀 at 2015-11-12 16:02:34
自己給自己蓋樓。
例子中BiTe采用自旋軌道耦合計(jì)算，輸入文件energyso替換了energy文件，同時(shí)輸入命令多了一項(xiàng) ‘-so’，告訴程序每個(gè)能態(tài)只包含一個(gè)電子。
那么考慮自旋極化時(shí)可不可以也采用類似做法，直接每個(gè)K ...

高人，您好請(qǐng)問(wèn)最后你是怎么處理加SOC時(shí)，計(jì)算熱電的呢

贊一下

回復(fù)此樓

創(chuàng)新、思路決定出路

5樓2016-07-31 21:38:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ouyangyulou

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 254
紅花: 12
帖子: 124
在線: 95.9小時(shí)
蟲(chóng)號(hào): 2850255
注冊(cè): 2013-12-05
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

高人，您好請(qǐng)問(wèn)最后你是怎么處理加SOC時(shí)，計(jì)算熱電的呢

贊一下

回復(fù)此樓

創(chuàng)新、思路決定出路

6樓2016-07-31 21:40:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wyzqmc

銀蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 963.2
散金: 100
帖子: 212
在線: 107.6小時(shí)
蟲(chóng)號(hào): 4323453
注冊(cè): 2015-12-30
性別: GG
專業(yè): 半導(dǎo)體物理

引用回帖:

3樓: Originally posted by 拍拍熊出隱刀 at 2015-11-13 08:40:01
自己給自己砌磚。
我在vasp2boltz.py文件中write_bandstructure_boltztrap函數(shù)中添加了關(guān)于自旋極化情況下energy文件的轉(zhuǎn)換，大家看看有沒(méi)有問(wèn)題啊！
if bandstructure!=1:
n_up=len(enk)
...

樓主，我遇到了同樣的問(wèn)題，這個(gè)接口轉(zhuǎn)換ispin=2的能量文件轉(zhuǎn)換不了，您的代碼也試了不行，您最后是怎么解決啊啊

贊一下

回復(fù)此樓

7樓2016-11-22 14:08:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

江霞

新蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 49.2
散金: 150
紅花: 2
帖子: 718
在線: 108.6小時(shí)
蟲(chóng)號(hào): 4465223
注冊(cè): 2016-03-04
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

你好你這個(gè)問(wèn)題解決了嗎？怎么做呀？可不可以告訴我呀？

贊一下

回復(fù)此樓

8樓2017-09-22 16:32:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主拍拍熊出隱刀的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 070300化學(xué)求調(diào)劑 +5	苑豆豆 2026-03-20	5/250	2026-03-23 09:39 by guoliang1816
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11	楓橋ZL 2026-03-18	13/650	2026-03-22 20:26 by edmund7
[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +4	zxts12 2026-03-21	7/350	2026-03-22 09:57 by zxts12
[考研] 材料學(xué)碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 306求0703調(diào)劑一志愿華中師范 +5	紙魚(yú)ly 2026-03-21	5/250	2026-03-21 17:11 by 學(xué)員8dgXkO
[考研] 268求調(diào)劑 +9	簡(jiǎn)單點(diǎn)0 2026-03-17	9/450	2026-03-21 15:37 by lature00
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過(guò) 本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 332求調(diào)劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機(jī)用戶 2026-03-17	9/450	2026-03-20 14:24 by 無(wú)懈可擊111
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 283求調(diào)劑 +3	聽(tīng)風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠