版塊導(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)劑公告

返回列表

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

» 猜你喜歡

一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑已經(jīng)有15人回復(fù)
085601材料工程找調(diào)劑已經(jīng)有15人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有15人回復(fù)
085600 286分材料求調(diào)劑已經(jīng)有10人回復(fù)
南京大學(xué)化學(xué)調(diào)劑已經(jīng)有14人回復(fù)
一志愿北京工業(yè)大學(xué)，324分求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑已經(jīng)有6人回復(fù)
26申博自薦已經(jīng)有6人回復(fù)
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分已經(jīng)有10人回復(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中有的。

你們組看起來是做單壁納米管的，所以這個(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的回帖

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ū)別，我沒有研究過~

贊一下(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的上下自旋，沒有你說的將P軌道投影到px，py，pz啊。這是怎么回事，是我的設(shè)置有問題

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

千萬千萬要仔細(xì)看手冊(cè)。這是千百年來，學(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 的主題更新

返回列表

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

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

[考研] 總分293求調(diào)劑 +8	加一一九 2026-03-25	11/550	2026-03-29 19:53 by 無際的草原
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4	幸運(yùn)的三木 2026-03-29	5/250	2026-03-29 14:49 by Miko19
[考研] 309求調(diào)劑 +7	誰不是少年 2026-03-29	7/350	2026-03-29 13:17 by mumin1990
[考研] 本科雙非材料，跨考一志愿華電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
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5	0857求調(diào)劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7	小小甜面團(tuán) 2026-03-25	7/350	2026-03-28 20:25 by 唐沐兒
[考研] 材料求調(diào)劑一志愿哈工大324 +7	閆旭東 2026-03-28	9/450	2026-03-28 08:51 by Xu de nuo
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8	曼111 2026-03-24	9/450	2026-03-28 07:58 by YYYYX1234
[考研] 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
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚的貓
[考研] 286求調(diào)劑 +4	lim0922 2026-03-26	4/200	2026-03-27 10:28 by guoweigw
[考研] 327求調(diào)劑 +7	prayer13 2026-03-23	7/350	2026-03-26 20:48 by 不吃魚的貓
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-26	3/150	2026-03-26 19:57 by nihaoar
[考研] 340求調(diào)劑 +5	話梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin
[考研] 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 了了了了。。
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta

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ì)您同樣有幫助: