partial charge可以很方便地輸出每一條能帶上每乙個波矢所對應的電荷密度，有利於我們更深層次得理解體系的電子結構性質。下面將介紹如何利用第一性原理軟體vasp去計算體系的partial charge。
計算partial charge要用到以下幾個開關：
lpard，該開關有兩個取值，t和f。t代表計算partial charge。如果僅僅只設定lpard，那麼總的電荷密度就會被寫入chgcar
iband，用於指定需要計算partial charge的能帶。例如設定iband= 20 21 22 23，那麼vasp就會計算從第20到第23四條帶的電荷密度
eint，用於指定需要計算partial charge的能帶的能量範圍，一般設定為兩個實數。如果只設定乙個實數，那麼第二個實數預設為費公尺能（也就是說計算從指定的能量值到費公尺能之間能帶的partial charge）；如果eint設定了，但nbmod沒設定，那麼vasp會自動將nbmod設定為-2
nbmod，設定為整數，如果nbmod>0，那麼該值就是iband對應帶的數目，如果iband已經設定了，nbmod也會被調整成語iband相應的值；如果nbmod=0，那麼vasp會計算所有能帶的電荷密度（包括被填充的和未被填充的）；如果nbmod=-1，那麼vasp會計算總的電荷密度，-1是nbmod的預設值；如果nbmod=-2，那麼vasp會計算eint指定的能量區間的partial charge；nbmod=-3和-2基本相同，差別僅在於能量區間取到了eint設定的值和費公尺能之間
kpuse，用於指定計算partial charge的k點，例如設定kpuse= 12 3 4，vasp將計算前四個k點的partial charge
lsepb，設定為t表明分開計算每條帶的電荷密度，並將其輸出到檔案parchg.nb.⋆；如果設定為f，vasp會將選中的帶的電荷密度合併到乙個檔案parchg.allb.⋆或者parchg進行輸出
lsepk，設定為t表明分開計算每個k點的電荷密度，並將其輸出到檔案
計算時，首先進行一次一般的靜態自洽，incar的設定如下：
system=***
istart=0
icharg=2
nsw=0
prec=n
encut=400
ismear=0
sigma=0.01
ediff=0.1e-04
靜態自洽結束後，然後加上上面給出的開關再跑一次，這裡給出乙個incar設定的例項：
system=***
istart=1
icharg=1
nsw=0
prec=n
encut=400
ismear=0
sigma=0.01
ediff=0.1e-04
lpard=t
iband=568 569 570
lsepb=t
執行結束後會出現三個檔案，因為在這裡我設定iband=568 569 570，所以後面會出現parchg.0568.allk、parchg.0569.allk、parchg.0570.allk三個檔案，將其改名成vesta可以識別的檔案，例如chgcar568、chgcar569、chgcar570，然後將這三個檔案拖進vesta中去。
最後，需要說的是，在vesta中可以選擇objects→properties→isosu***ces在其中選中電荷密度等值面（見下圖），然後在isosu***ce level中通過調整等值面的值來調整圖中partial charge的形狀。
歡迎關注計算材料學
如何計算 軟體如何計算疊合板？
中國建築裝配化，已經進入國家發展規劃，所以，現在越來越多的建築設計都是裝配化設計，我們自己接到的工程，裝配的圖紙也越來越多。今天的案例，就講乙個我們最近遇到的乙個裝配式工程中疊合板的案例。圖1 所謂的疊合板，就是下面是預製板，上面是現澆板，大家看一下圖1，就是施工人員正在放置疊合板的下部預製部分，圖...
如何計算idf
我們這道乙個單詞的idf的計算公式可以如下 上面這個公式需要先要有乙個語料庫，如果此時你恰好沒有這樣的乙個語料庫，那你要怎麼去估計這個idf呢？現在假定在一篇 中國的蜜蜂養殖 文章當中，假定該文長度為1000個詞，中國 蜜蜂 養殖 各出現了20次，則這三個詞的詞頻 tf 都為0.02.為了計算這些詞...
如何計算素數
題目 如何計算素數？思路 素數只可以被1和它本身相除，所以我打算把素數從2除到它減去1 解答 for i in range 100,201 for j in range 2,i 1 if i j 0 break else print d i,end 但是這個是行不通的，像101會被除99次，也就會出...