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ù)此樓

» 猜你喜歡

309求調(diào)劑已經(jīng)有3人回復(fù)
070300化學(xué)354求調(diào)劑已經(jīng)有4人回復(fù)
330分求調(diào)劑已經(jīng)有3人回復(fù)
299求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑已經(jīng)有6人回復(fù)
300求調(diào)劑，材料科學(xué)英一數(shù)二已經(jīng)有8人回復(fù)
招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑已經(jīng)有5人回復(fù)
070305高分子化學(xué)與物理 304分求調(diào)劑已經(jīng)有7人回復(fù)
289求調(diào)劑已經(jīng)有13人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(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)容被屏蔽

7樓2011-11-14 14:09:06

已閱回復(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ā)表

[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6	崔wj 2026-03-26	6/300	2026-03-29 01:11 by hanserlol
[考研] 086000生物與醫(yī)藥調(diào)劑 +3	Feisty。 2026-03-28	7/350	2026-03-29 00:41 by 544594351
[考研] 348求調(diào)劑 +3	小懶蟲不懶了 2026-03-28	3/150	2026-03-29 00:39 by 544594351
[考研] 295材料工程專碩求調(diào)劑 +7	1428151015 2026-03-27	7/350	2026-03-28 19:58 by S240
[考研] 材料與化工272求調(diào)劑 +9	阿斯蒂芬2004 2026-03-28	9/450	2026-03-28 15:21 by 1018329917
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 266分，求材料冶金能源化工等調(diào)劑 +7	哇呼哼呼哼 2026-03-27	9/450	2026-03-28 12:22 by zllcz
[考研] 0856，材料與化工321分求調(diào)劑 +12	大饞小子 2026-03-27	13/650	2026-03-28 10:56 by self2008
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4	哇呼哼呼哼 2026-03-26	4/200	2026-03-27 17:02 by zhyzzh
[考博] 26申博 +3	加油沖啊！ 2026-03-26	3/150	2026-03-27 15:38 by cls512
[考研] 343求調(diào)劑 +4	贈我一本書 2026-03-23	4/200	2026-03-27 00:40 by wxiongid
[考研] 321求調(diào)劑 +6	wasdssaa 2026-03-26	6/300	2026-03-26 20:57 by sanrepian
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4	JKSOIID 2026-03-26	4/200	2026-03-26 15:58 by 不吃魚的貓
[考研] 081700 調(diào)劑 267分 +11	迷人的哈哈 2026-03-23	11/550	2026-03-26 15:41 by zzll406
[考研] 281求調(diào)劑 +6	Koxui 2026-03-24	7/350	2026-03-26 15:37 by 無際的草原
[考研] 材料專碩 335 分求調(diào)劑 +4	拒絕冷暴力 2026-03-25	4/200	2026-03-25 18:45 by haxia
[考研] 296求調(diào)劑 +4	汪�。�！ 2026-03-25	7/350	2026-03-25 16:41 by 汪��？！
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5	超贊的 2026-03-24	5/250	2026-03-24 21:07 by 星空星月
[考研] 材料專碩找調(diào)劑 +5	哈哈哈吼吼吼哈 2026-03-23	5/250	2026-03-24 19:07 by 了了了了。。
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte