fly_elephan

銀蟲 (小有名氣)

1ST強(qiáng)帖: 1
應(yīng)助: 3 (幼兒園)
金幣: 268.4
散金: 35
紅花: 7
帖子: 219
在線: 124.2小時
蟲號: 1177420
注冊: 2010-12-25
性別: GG
專業(yè): 電路與系統(tǒng)

[交流] vasp計算結(jié)果處理態(tài)密度要用到的程序已有7人參與

在材料電學(xué)特性的研究中，能帶和態(tài)密度是我們最關(guān)心的，下面的這個程序可以用于vasp計算結(jié)果的處理，可以分離出總態(tài)密度和分態(tài)密度，這個程序是我們研究所用了很久的，絕對沒問題，有資源大家共享，歡迎下載。
#! /bin/ksh
# Script to break the VASP DOSCAR file into segments for
# each atomic site.
# JMS NRL 5/0/02

# Modified to also align the Fermi level from the ground state
# to E=0.0
# JMS NRL 6/19/02

# Executables reside here
BIN=~/tzyy/AlP/swnt10

if [ $# -ne 0 ]; then
echo "Usage: split_dos"
exit 2
fi

# Script to split the DOSCAR file into the atomic
# projections labeled by atom number
dosfile=DOSCAR
outfile=OUTCAR
infile=INCAR

# Token for splitting the files
token=$(sed -n '6 p' $dosfile | awk '{print $1}')

# Number of points
nl=$(sed -n '6 p' $dosfile | awk '{print $3}')

# Number of atoms
natom=$(sed -n '1 p' $dosfile | awk '{print $1}')

# Get the Fermi level if the OUTCAR file is present,
# else we set it to zero.
if [ -a $outfile ]; then
echo "The" $outfile "exists, we use it to get the Fermi level,"
echo "the RWIGS tag and the number of spins."
efermi=$(grep "E-fermi" $outfile | tail -1 | awk '{print $3}')
echo "Fermi level:" $efermi
nspin=$(grep "ISPIN" $outfile | tail -1 | awk '{print $3}')
if [ $nspin -eq 2 ]; then
echo "Spin polarized calculation"
else
echo "Unpolarized calculation"
fi

# 2.a
# JMS 2/3/03 Modified to accept specification by LORBIT token.
#
lorbit=$(grep "LORBIT" $outfile | tail -1 | awk '{print $3}')
if [ $lorbit -ge 10 ]; then
      echo "LORBIT > 10"
      echo "WARNING: not completely test for vasp.4.*"
echo "Use at your own risk. Please check results for consistency."
form=1
else
# 2.a
rwigs=$(grep "RWIGS" $outfile | tail -1 | awk '{print $3}' | sed 's/\.//g')
if [ $rwigs -eq -100 ]; then
      echo "RWIGS token not set"
      form=0
else
      echo "RWIGS token set"
      form=1
fi
# 2.a
fi
# 2.a

else
echo "The" $outfile "does not exist, we set the Fermi level to 0"
echo "assume an unpolarized calculation, and RWIGS not set"
form=0
nspin=1
efermi=0.0
fi

# If the outcar file is not present and you wish to set something by hand
# you should do it here. Uncomment the tokens below and set them by hand.

#efermi=2.255
#form=
#nspin=

# Get the atomic coordinates
$BIN/vp >| tmp.dat
tail -$natom tmp.dat | awk '{printf "%s %12.8f %12.8f %12.8f \n", "#", $2, $3, $4}' >| tmp.dat2

# Total is first DOS
if [ $form -eq 0 ]; then
i=0
else
i=0
fi

start=7
end=$((start+nl-1))
#echo $start $end
#exit 0
rm -f DOS0
if [ $form -eq 1 ]; then

while [ $i -le $natom ]; do

#echo $i $start $end

if [ $i -gt 0 ]; then
      sed -n ''$i' p' tmp.dat2 >| DOS$i
fi

if [ $i -gt 0 ]; then

# Atomic projected DOS
      if [ $nspin -eq 2 ]; then
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
      else
sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f %12.8f \n", $1+(-1)*'$efermi', $2, $3, $4 }' >> DOS$i
      fi

else

# Total DOS
      if [ $nspin -eq 2 ]; then
sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f \n", $1+(-1)*'$efermi', $2, -1*$3, $4, -1*$5 }' >> DOS$i
      else
sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f \n", $1+(-1)*'$efermi', $2, $3 }' >> DOS$i
      fi

      fi

start=$((end+2))
end=$((start+nl-1))
i=$((i+1))
done

else

while [ $i -le 0 ]; do

#echo $i $start $end

if [ $i -gt 0 ]; then
      sed -n ''$i' p' tmp.dat2 >| DOS$i
fi

if [ $i -gt 0 ]; then

# Atomic projected DOS
      if [ $nspin -eq 2 ]; then
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
      else
sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f %12.8f \n", $1+(-1)*'$efermi', $2, $3, $4 }' >> DOS$i
      fi

else

# Total DOS
      if [ $nspin -eq 2 ]; then
sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f %12.8f %12.8f \n", $1+(-1)*'$efermi', $2, -1*$3, $4, -1*$5 }' >> DOS$i
      else
sed -n ''$start','$end' p' $dosfile | awk '{printf "%12.8f %12.8f %12.8f \n", $1+(-1)*'$efermi', $2, $3 }' >> DOS$i
      fi

      fi

start=$((end+2))
end=$((start+nl-1))
i=$((i+1))
done

fi

exit 0

回復(fù)此樓

» 猜你喜歡

一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有8人回復(fù)
349求調(diào)劑已經(jīng)有6人回復(fù)
329求調(diào)劑已經(jīng)有9人回復(fù)
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）已經(jīng)有9人回復(fù)
327求調(diào)劑已經(jīng)有5人回復(fù)
求收留已經(jīng)有5人回復(fù)
375求調(diào)劑已經(jīng)有4人回復(fù)
340求調(diào)劑已經(jīng)有6人回復(fù)
298求調(diào)劑已經(jīng)有3人回復(fù)
085601材料工程找調(diào)劑已經(jīng)有9人回復(fù)

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

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

用vasp如何計算能帶圖和態(tài)密度圖呢？已經(jīng)有11人回復(fù)
用vasp計算態(tài)密度已經(jīng)有5人回復(fù)
VASP計算態(tài)密度出現(xiàn)問題已經(jīng)有12人回復(fù)
【求助】VASP計算電子態(tài)密度時提示錯誤WAVECAR 已經(jīng)有4人回復(fù)
【求助】求VASP計算ZnO能帶和態(tài)密度的詳細(xì)輸入設(shè)置已經(jīng)有15人回復(fù)
【求助】VASP怎樣精確計算出費米能級處的態(tài)密度值？已經(jīng)有22人回復(fù)
【求助】關(guān)于vasp態(tài)密度計算已經(jīng)有10人回復(fù)
【求助】VASP計算電子態(tài)密度已經(jīng)有34人回復(fù)
【求助】關(guān)于鋰原子vasp計算后的1s2s2p分波態(tài)密度處理已經(jīng)有12人回復(fù)

1樓 2011-05-30 10:09:48

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

gougou531

禁蟲 (正式寫手)

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

本帖內(nèi)容被屏蔽

5樓2011-11-14 14:05:43

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

查看全部 19 個回答

cenwanglai

榮譽(yù)版主 (知名作家)

老和山猥瑣派九段

1ST強(qiáng)帖: 5
應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計算化學(xué)
管轄: 仿真模擬

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
zzy870720z(金幣+2): 謝謝指教 2011-05-30 15:23:30

這個程序vtstool code中有的。

你們組看起來是做單壁納米管的，所以這個腳本可以一直用著。

如果含有f電子，就沒有辦法用了。LORBIT = 11和LORBIT = 10是有很大差異的，這個腳本需要修改才能用。

贊一下(3人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

哥是傳說

2樓2011-05-30 11:01:16

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

wujing110

禁蟲 (正式寫手)

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

本帖內(nèi)容被屏蔽

3樓2011-11-14 11:21:03

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

cenwanglai

榮譽(yù)版主 (知名作家)

老和山猥瑣派九段

1ST強(qiáng)帖: 5
應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計算化學(xué)
管轄: 仿真模擬

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

引用回帖:

3樓: Originally posted by wujing110 at 2011-11-14 11:21:03:
請問區(qū)長童鞋，LORBIT = 11時候應(yīng)該怎么改呀？我設(shè)置LORBIT = 11，一共有19列數(shù)據(jù)，但是我現(xiàn)在只想要 s+,s-,p+,p-,d+,d-，那些角動量分量要一列一列加成總的么？

是啊。第四列到九列是p的三個分量及上下旋。三個奇數(shù)列相加的p+，三個偶數(shù)列相加的p-。后面10列是d分量及自旋。上旋和下旋分別相加。

可以用excel或者origin來做列的相加。

贊一下(1人)

回復(fù)此樓

4樓2011-11-14 13:21:54

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

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

[考研] 298求調(diào)劑 +3	種圣賜 2026-03-29	3/150	2026-03-29 12:06 by longlotian
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥求調(diào)劑 +4	篤行致遠(yuǎn). 2026-03-28	4/200	2026-03-29 11:55 by longlotian
[考研] 329求調(diào)劑 +7	星野? 2026-03-26	7/350	2026-03-29 06:43 by 544594351
[考研] 289求調(diào)劑 +13	新時代材料 2026-03-27	13/650	2026-03-29 01:16 by 544594351
[考研] 316求調(diào)劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 315求調(diào)劑 +4	akie... 2026-03-28	5/250	2026-03-28 21:05 by zhq0425
[考研] 083000學(xué)碩274求調(diào)劑 +8	Li李魚 2026-03-26	8/400	2026-03-28 20:33 by 加油向未來啊
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4	小張zxy張 2026-03-26	8/400	2026-03-28 19:25 by lbsjt
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8	AZMK 2026-03-25	12/600	2026-03-28 18:15 by AZMK
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5	niko- 2026-03-27	5/250	2026-03-28 14:25 by 唐沐兒
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7	hanamiko 2026-03-27	7/350	2026-03-28 09:56 by 李上岸0921
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4	Ripcord上岸 2026-03-27	4/200	2026-03-27 17:00 by zhyzzh
[考研] 08開頭275求調(diào)劑 +4	拉誰不重要 2026-03-26	4/200	2026-03-27 14:12 by Delta2012
[考研] 314求調(diào)劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[考研] 281求調(diào)劑 +3	亞克西good 2026-03-26	5/250	2026-03-26 19:48 by 不吃魚的貓
[考研] 318求調(diào)劑 +3	plum李子 2026-03-23	3/150	2026-03-25 09:42 by 霧散后相遇lc
[論文投稿] 急發(fā)核心期刊論文 +3	賢達(dá)問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹