| 查看: 2188 | 回復(fù): 18 | ||
shala123456鐵桿木蟲 (知名作家)
|
[求助]
如何在vasp中TDOS費米能不歸零 已有3人參與
|
|
用split 分DOSCAR后,分出的Tdos費米能級就在0ev處,比如費米能是7.1ev時,那么其他的能量是7.1的相對值。 現(xiàn)在我的問題是:如何在DOS圖中Oev處不代表費米能,改為它的實際能量? 謝謝大家啦! ![]() 5555555555.png |
新蟲 (初入文壇)
|
把數(shù)據(jù)導(dǎo)入origin中,能量整體加上費米能,然后以新的能量作圖 發(fā)自小木蟲IOS客戶端 |
鐵桿木蟲 (知名作家)
鐵桿木蟲 (知名作家)
金蟲 (小有名氣)
|
費米能級處在0eV,是由于split_dos腳本導(dǎo)致的,VASP本身沒有參數(shù)使費米能未0eV。 不知道你用的是哪個split_dos腳本,一般在腳本開始的注釋處有以下注釋: # Modified to also align the Fermi level from # the ground state to E=0.0 具體代碼也做了使費米能出于0eV. 你可以把代碼改了。 |
鐵桿木蟲 (知名作家)
金蟲 (小有名氣)
|
用寫字板打開split_dos腳本,查找 $efermi 你應(yīng)該能找到類似這樣的幾行代碼 sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f \n", $1+(-1)*'$efermi', $2, -1*$3, $4, -1*$5, $6, -1*$7}' >> DOS$i 該行代碼是由sed和awk兩個處理文件的命令組合而成,其中 $1+(-1)*'$efermi' $1:表示指定行第一個變量,即DOSCAR中的第一列能量; $efermi:表示從OUTCAR中讀取的費米能; 你把代碼中所有$1+(-1)*'$efermi'改為$1即可。 |
金蟲 (小有名氣)
鐵桿木蟲 (知名作家)
鐵桿木蟲 (知名作家)
|
你所說的腳本我沒找到,能麻煩你發(fā)我一份嗎 這是我常用的split dos 軟件 https://yunpan.cn/cSYEDe8jMKwZE (提取碼:014c) 我一般用這個劈裂dos 還有就是VASP中有一個自帶的split_dos命令 這兩個好像都沒有腳本 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 調(diào)劑求院校招收 +4 | 鶴鯨鴿 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料求調(diào)劑 +9 | @taotao 2026-03-21 | 9/450 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
[考研] 283求調(diào)劑(080500) +4 | A child 2026-03-27 | 4/200 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 22408 359分調(diào)劑 +3 | Qshers 2026-03-27 | 3/150 |
|
|
[考研] 085600材料與化工306 +10 | z1z2z3879 2026-03-21 | 11/550 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯條 2026-03-23 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|