 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 5039 | 回復(fù): 6 | ||
| 【獎勵】 本帖被評價3次，作者byin增加金幣 2 個 | ||
| 本帖產(chǎn)生 1 個 1ST強(qiáng)帖 ，點(diǎn)擊這里進(jìn)行查看 | ||
[資源]
用matlab腳本繪制vasp中的能帶圖
|
||
|
把EIGENVAL和.m腳本放一起，修改腳本里的一些參數(shù)“手動修改區(qū)”，然后運(yùn)行腳本即可。見笑了。 %%%%%%%%%%%%%%%%%%%% 此腳本用于VASP的能帶繪制 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc clear %%%%%%%%%%%%%%%%%%%%%%%%%%% 手動修改區(qū) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% EFermi=6.1748; %設(shè)置費(fèi)米能級 NBands=8; %一共計算的能帶數(shù)目 inter=100; %每段PATH的插值點(diǎn)數(shù)目 NPath=4; %PATH數(shù)目 Ktotal=inter*NPath; %特殊K點(diǎn)的總數(shù)目 xmin=1; %x，y坐標(biāo)軸的范圍 xmax=Ktotal; ymin=-6; ymax=7; kpath={'L','\Gamma' ,'X','K','\Gamma'}; %設(shè)置K空間的路徑 %%%%%%%%%%%%%%%%%%%%%%%%%% 圖像邊框 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% A4 = figure('PaperType','a4letter','PaperSize',[20.98 29.68]); axes1 = axes('Parent',A4,'Position',[0.2 0.1 0.6 0.8],'LineWidth',1.5,'FontSize',15.0,'FontWeight','bold'); %Position [left bottom width height] box(axes1,'on'); hold(axes1,'all'); %%%%%%%%%%%%%%%%%%%%%%%%%%% 數(shù)據(jù)讀取 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% data=textread('EIGENVAL','','headerlines',7); for i=1:1:NBands, for j=1:1:Ktotal, P(j,: )=[j,data((NBands+1)*(j-1)+1+i,2)-EFermi]; %讀出每條能帶的數(shù)據(jù)（每個i中的各個j），做出曲線 end plot(P(:,1),P(:,2),'-r.'); hold on end %%%%%%%%%%%%%%%%%%%%%%%%%%%% 坐標(biāo)軸設(shè)置 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% axis([xmin xmax ymin ymax]); set(gca,'XTick',[1,inter:inter:Ktotal]) set(gca,'XTickLabel','') j=1; %高對稱點(diǎn)符號 for i=[1,inter:inter:Ktotal] text(i,ymin-0.7,kpath(j),'FontSize',18, 'horizontalAlignment', 'center') j=j+1; end ylabel('E-E_F (eV)'); X=[xmin xmax]; %輔助線 - 水平 Y=[0 0]; line(X,Y,'LineWidth',1.5,'Color',[0 0 0]) for i=1:1: (NPath-1) %輔助線 - 垂直 X=[i*inter i*inter]; Y=[ymin ymax]; line(X,Y,'LineWidth',1.5,'Color',[0 0 0]); end print -dpng Band_Structure_of_ %輸出圖像文件 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% END %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [ Last edited by byin on 2011-5-30 at 21:27 ] |
回復(fù)此樓» 猜你喜歡
085601材料工程找調(diào)劑
已經(jīng)有11人回復(fù)
-
總分293求調(diào)劑
已經(jīng)有9人回復(fù)
-
334分 一志愿武理 材料求調(diào)劑
已經(jīng)有5人回復(fù)
-
327求調(diào)劑
已經(jīng)有6人回復(fù)
-
329求調(diào)劑
已經(jīng)有10人回復(fù)
-
275求調(diào)劑
已經(jīng)有15人回復(fù)
-
291求調(diào)劑
已經(jīng)有6人回復(fù)
-
309求調(diào)劑
已經(jīng)有7人回復(fù)
-
297求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 如何用matlab畫出石墨烯高對稱點(diǎn)的能帶結(jié)構(gòu)
已經(jīng)有3人回復(fù)
- 在VASP中的腳本
已經(jīng)有4人回復(fù)
- vasp 算能帶k點(diǎn)數(shù)目問題
已經(jīng)有10人回復(fù)
- 推薦幾個vasp處理計算結(jié)果的python腳本
已經(jīng)有66人回復(fù)
- 計算模擬需要的服務(wù)器 硬件 MS,VASP,matlab等等吧
已經(jīng)有7人回復(fù)
- NBANDS 的設(shè)置問題
已經(jīng)有6人回復(fù)
- 電機(jī)原理與設(shè)計的MATLAB分析
已經(jīng)有140人回復(fù)
- 【求助】請問如何用P4VASP在LINUX系統(tǒng)下畫能帶圖，或者導(dǎo)出數(shù)據(jù)
已經(jīng)有5人回復(fù)
- vasp_bnd的Matlab改進(jìn)版，更方便~
已經(jīng)有3人回復(fù)
- 請教VASP CI-NBE分析腳本的使用問題
已經(jīng)有6人回復(fù)
- matlab多維插值原理是什么？
已經(jīng)有5人回復(fù)
- VASP有幾個腳本文件？
已經(jīng)有8人回復(fù)
- DOSCAR畫態(tài)密度
已經(jīng)有16人回復(fù)
- VASP能帶橫坐標(biāo)怎么畫
已經(jīng)有6人回復(fù)
- 【活動or資源】vasp晶體結(jié)構(gòu)優(yōu)化
已經(jīng)有142人回復(fù)
- 【求助】vasp作業(yè)提交腳本
已經(jīng)有5人回復(fù)
- 【求助】求將vasp的輸出文件轉(zhuǎn)化成xyz坐標(biāo)的腳本
已經(jīng)有11人回復(fù)
- 【求助】Matlab 四維圖形的繪制。
已經(jīng)有12人回復(fù)
2樓2011-06-01 19:52:25
3樓2011-06-01 21:15:17
4樓2011-06-03 11:57:02
5樓2011-11-05 21:46:17
7樓2014-07-01 20:13:38
簡單回復(fù)
2014-05-01 10:57
回復(fù)
三星好評 
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 byin 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +8 | 舉個栗子oi 2026-03-24 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +5 | 景隨風(fēng) 2026-03-29 | 8/400 |
|
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +7 | WW.' 2026-03-23 | 9/450 |
|
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +8 | leaflight 2026-03-24 | 8/400 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
|
[考研] 359求調(diào)劑 +4 | 王了個楠 2026-03-25 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑！�。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|