版塊導(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í)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

下一頁(yè)

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

» 猜你喜歡

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ù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有6人回復(fù)
調(diào)劑求院校招收已經(jīng)有5人回復(fù)
調(diào)劑310 已經(jīng)有8人回復(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的回帖

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中有的。

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

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

贊一下(3人)

回復(fù)此樓

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

哥是傳說(shuō)

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來(lái)做列的相加。

贊一下(1人)

回復(fù)此樓

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

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

gougou531

禁蟲 (正式寫手)

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

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

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

已閱回復(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è)紅包，謝謝回帖

引用回帖:

5樓: Originally posted by gougou531 at 2011-11-14 14:05:43:
區(qū)長(zhǎng)，我i怎么覺得這些數(shù)據(jù)是在procar中，不是在doscar中，設(shè)置lorbit=11只是影響doscar，doscar影響小吧，只是當(dāng)lorbit=11時(shí)doscar中像O原子會(huì)有d軌道的分量。謝謝

POSCAR和PROCAR中都有。至于二者有什么區(qū)別，我沒有研究過(guò)~

贊一下(1人)

回復(fù)此樓

6樓2011-11-14 14:06:49

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

gougou531

禁蟲 (正式寫手)

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

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

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

已閱回復(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è)紅包，謝謝回帖

引用回帖:

7樓: Originally posted by gougou531 at 2011-11-14 14:09:06:
我的DOSCAR中，只有s，p，d的上下自旋，沒有你說(shuō)的將P軌道投影到px，py，pz啊。這是怎么回事，是我的設(shè)置有問題

要保證所用的勢(shì)是pw91。

千萬(wàn)千萬(wàn)要仔細(xì)看手冊(cè)。這是千百年來(lái)，學(xué)習(xí)vasp和其他軟件的人的忠告~

贊一下(1人)

回復(fù)此樓

8樓2011-11-14 14:11:09

已閱回復(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的回帖

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è)紅包，謝謝回帖

引用回帖:

9樓: Originally posted by wujing110 at 2011-11-14 16:40:10:
謝謝區(qū)長(zhǎng)童鞋
我還要一個(gè)問題就是，我看到文獻(xiàn)中有做出t2g和eg軌道的圖，t2g由dxy,dxz,dyz組成么？那我把對(duì)應(yīng)的這幾列相加就是t2g軌道么？那 eg軌道呢？求區(qū)長(zhǎng)給講解講解

這個(gè)是d五重兼并軌道（能級(jí)）的分裂。分裂成上兩組，兩上三下或者三上兩下。三個(gè)一組的是t2g，兩個(gè)一組的是eg。這兩組兼并能級(jí)為什么要加和呢？

加不加隨你吧，我也不清楚你要干什么用。要注意的是t2g和eg所對(duì)應(yīng)的是什么。dxy之類的d后面的下標(biāo)與你的晶胞坐標(biāo)軸的選取有關(guān)。

贊一下(1人)

回復(fù)此樓

10樓2011-11-14 16:50:58

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 fly_elephan 的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 332求調(diào)劑 +6	蕉蕉123 2026-03-28	6/300	2026-03-29 00:37 by 544594351
[考研] 學(xué)碩274求調(diào)劑 +9	Li李魚 2026-03-26	9/450	2026-03-28 21:42 by bymhappy
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7	小小甜面團(tuán) 2026-03-25	7/350	2026-03-28 20:25 by 唐沐兒
[考研] 085600，材料與化工321分求調(diào)劑 +9	大饞小子 2026-03-28	9/450	2026-03-28 14:56 by 神馬都不懂
[考研] 329求調(diào)劑 +6	星野? 2026-03-26	6/300	2026-03-28 14:14 by 唐沐兒
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-27	6/300	2026-03-28 14:10 by 唐沐兒
[考研] 0856求調(diào)劑 +11	zhn03 2026-03-25	12/600	2026-03-28 13:32 by 唐沐兒
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 070300求調(diào)劑306分 +4	26要上岸 2026-03-27	4/200	2026-03-28 13:06 by 唐沐兒
[考研] 316求調(diào)劑 +5	Pigcasso 2026-03-24	5/250	2026-03-27 12:10 by zhshch
[考研] 359求調(diào)劑 +4	王了個(gè)楠 2026-03-25	4/200	2026-03-27 08:43 by 不吃魚的貓
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 0703化學(xué)求調(diào)劑 +3	丹青奶蓋 2026-03-26	5/250	2026-03-26 20:11 by macy2011
[考研] 求調(diào)劑 +8	Auroracx 2026-03-22	8/400	2026-03-26 19:55 by 不吃魚的貓
[考研] 334分一志愿武理-080500 材料求調(diào)劑 +4	李李不服輸 2026-03-25	4/200	2026-03-25 21:26 by 星空星月
[考研] 生物技術(shù)與工程 +3	1294608413 2026-03-25	4/200	2026-03-25 18:02 by 1294608413
[考研] 0854人工智能方向招收調(diào)劑 +4	章小魚567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325
[考研] 材料調(diào)劑 +3	iwinso 2026-03-23	3/150	2026-03-25 11:29 by greychen00
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060

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

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

» 猜你喜歡

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

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

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

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

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