 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 5044 | 回復(fù): 6 | ||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)3次，作者byin增加金幣 2 個(gè) | ||
| 本帖產(chǎn)生 1 個(gè) 1ST強(qiáng)帖 ，點(diǎn)擊這里進(jìn)行查看 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
byin銅蟲(chóng) (正式寫(xiě)手)
|
[資源]
用matlab腳本繪制vasp中的能帶圖
|
|
|
把EIGENVAL和.m腳本放一起，修改腳本里的一些參數(shù)“手動(dòng)修改區(qū)”，然后運(yùn)行腳本即可。見(jiàn)笑了。 %%%%%%%%%%%%%%%%%%%% 此腳本用于VASP的能帶繪制 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc clear %%%%%%%%%%%%%%%%%%%%%%%%%%% 手動(dòng)修改區(qū) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% EFermi=6.1748; %設(shè)置費(fèi)米能級(jí) NBands=8; %一共計(jì)算的能帶數(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ù)（每個(gè)i中的各個(gè)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; %高對(duì)稱(chēng)點(diǎn)符號(hào) 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ù)此樓» 猜你喜歡
317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑
已經(jīng)有11人回復(fù)
-
332求調(diào)劑
已經(jīng)有12人回復(fù)
-
材料專(zhuān)碩 085600求調(diào)劑
已經(jīng)有5人回復(fù)
-
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可
已經(jīng)有14人回復(fù)
-
283求調(diào)劑（080500）
已經(jīng)有14人回復(fù)
-
材料工程專(zhuān)碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專(zhuān)碩調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑 有機(jī)化學(xué)考研356分
已經(jīng)有11人回復(fù)
-
環(huán)境科學(xué)與工程334分求調(diào)劑
已經(jīng)有8人回復(fù)
-
考研調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 如何用matlab畫(huà)出石墨烯高對(duì)稱(chēng)點(diǎn)的能帶結(jié)構(gòu)
已經(jīng)有3人回復(fù)
- 在VASP中的腳本
已經(jīng)有4人回復(fù)
- vasp 算能帶k點(diǎn)數(shù)目問(wèn)題
已經(jīng)有10人回復(fù)
- 推薦幾個(gè)vasp處理計(jì)算結(jié)果的python腳本
已經(jīng)有66人回復(fù)
- 計(jì)算模擬需要的服務(wù)器 硬件 MS,VASP,matlab等等吧
已經(jīng)有7人回復(fù)
- NBANDS 的設(shè)置問(wèn)題
已經(jīng)有6人回復(fù)
- 電機(jī)原理與設(shè)計(jì)的MATLAB分析
已經(jīng)有140人回復(fù)
- 【求助】請(qǐng)問(wèn)如何用P4VASP在LINUX系統(tǒng)下畫(huà)能帶圖，或者導(dǎo)出數(shù)據(jù)
已經(jīng)有5人回復(fù)
- vasp_bnd的Matlab改進(jìn)版，更方便~
已經(jīng)有3人回復(fù)
- 請(qǐng)教VASP CI-NBE分析腳本的使用問(wèn)題
已經(jīng)有6人回復(fù)
- matlab多維插值原理是什么？
已經(jīng)有5人回復(fù)
- VASP有幾個(gè)腳本文件？
已經(jīng)有8人回復(fù)
- DOSCAR畫(huà)態(tài)密度
已經(jīng)有16人回復(fù)
- VASP能帶橫坐標(biāo)怎么畫(huà)
已經(jīng)有6人回復(fù)
- 【活動(dòng)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ù)
byin
銅蟲(chóng) (正式寫(xiě)手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 0 (幼兒園)
- 金幣: 245
- 帖子: 395
- 在線: 174.4小時(shí)
- 蟲(chóng)號(hào): 948127
3樓2011-06-01 21:15:17
2樓2011-06-01 19:52:25
重新試過(guò)了，可以了，但圖不好看，還需要不斷改進(jìn)才行。謝謝 4樓2011-06-03 11:57:02
5樓2011-11-05 21:46:17
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +6 | j的立方 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 085601材料工程找調(diào)劑 +16 | oatmealR 2026-03-29 | 17/850 |
|
|
|
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +6 | 來(lái)財(cái)aa 2026-03-25 | 6/300 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
|
[考研] 0703化學(xué) +13 | 妮妮ninicgb 2026-03-27 | 13/650 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+11
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +7 | 26要上岸 2026-03-28 | 7/350 |
|
|
|
[考研] 289求調(diào)劑 +13 | 新時(shí)代材料 2026-03-27 | 13/650 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚(yú) 2026-03-26 | 9/450 |
|
|
|
[考研] 一志愿北化085600材料專(zhuān)碩275|有文章專(zhuān)利｜求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|