當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

fly_elephan

銀蟲 (小有名氣)

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

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

在材料電學(xué)特性的研究中，能帶和態(tài)密度是我們最關(guān)心的，下面的這個(gè)程序可以用于vasp計(jì)算結(jié)果的處理，可以分離出總態(tài)密度和分態(tài)密度，這個(gè)程序是我們研究所用了很久的，絕對(duì)沒問題，有資源大家共享，歡迎下載。
#! /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ù)此樓

» 猜你喜歡

一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷已經(jīng)有8人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
求化學(xué)調(diào)劑已經(jīng)有5人回復(fù)
085600，材料與化工321分求調(diào)劑已經(jīng)有10人回復(fù)
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分已經(jīng)有5人回復(fù)
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑已經(jīng)有3人回復(fù)
085600，專業(yè)課化工原理，320分求調(diào)劑已經(jīng)有4人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有16人回復(fù)
290求調(diào)劑已經(jīng)有3人回復(fù)
327求調(diào)劑已經(jīng)有3人回復(fù)

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

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

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

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

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

wujing110

禁蟲 (正式寫手)

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

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

9樓2011-11-14 16:40:10

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

查看全部 19 個(gè)回答

cenwanglai

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

老和山猥瑣派九段

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

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

這個(gè)程序vtstool code中有的。

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

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

贊一下(3人)

回復(fù)此樓

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

哥是傳說

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

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

wujing110

禁蟲 (正式寫手)

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

本帖內(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小時(shí)
蟲號(hào): 537452
注冊(cè): 2008-04-01
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 仿真模擬

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

引用回帖:

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

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

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

贊一下(1人)

回復(fù)此樓

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

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

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

[考研] 調(diào)劑310 +12	溫柔的晚安 2026-03-25	13/650	2026-03-29 20:01 by 無際的草原
[考研] 考研調(diào)劑 +7	小蠟新筆 2026-03-29	7/350	2026-03-29 19:00 by 學(xué)員8dgXkO
[考研] 277跪求調(diào)劑 +6	1915668 2026-03-27	10/500	2026-03-29 16:03 by 王亮_大連醫(yī)科大
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6	芝士雪baoo 2026-03-28	8/400	2026-03-29 08:16 by 松花缸1201
[考研] 調(diào)劑考研 +3	王杰一 2026-03-29	3/150	2026-03-29 08:09 by fmesaito
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4	IceGuo 2026-03-26	5/250	2026-03-29 01:25 by griffith2014
[考研] 一志愿華北電力大學(xué)能動(dòng)專碩，293，求調(diào)劑 +3	15537177284 2026-03-23	5/250	2026-03-28 16:11 by xxxsssccc
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9	丹青奶蓋 2026-03-26	10/500	2026-03-28 07:45 by barnett0632
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)
[考研] 324求調(diào)劑 +5	hanamiko 2026-03-26	5/250	2026-03-27 10:33 by wangjy2002
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 0703化學(xué)求調(diào)劑 +3	丹青奶蓋 2026-03-26	5/250	2026-03-26 20:11 by macy2011
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4	JKSOIID 2026-03-26	4/200	2026-03-26 15:58 by 不吃魚的貓
[考研] 生物技術(shù)與工程 +3	1294608413 2026-03-25	4/200	2026-03-25 18:02 by 1294608413
[考研] 求b區(qū)院校調(diào)劑 +4	周56 2026-03-24	5/250	2026-03-25 17:12 by yishunmin
[考研] 0854AI CV方向招收調(diào)劑 +4	章小魚567 2026-03-23	4/200	2026-03-25 17:04 by CoderLoser