 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 5042 | 回復(fù): 6 | ||
| 【獎勵】 本帖被評價3次，作者byin增加金幣 2 個 | ||
| 本帖產(chǎn)生 1 個 1ST強帖 ，點擊這里進行查看 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[資源]
用matlab腳本繪制vasp中的能帶圖
|
||
|
把EIGENVAL和.m腳本放一起，修改腳本里的一些參數(shù)“手動修改區(qū)”，然后運行腳本即可。見笑了。 %%%%%%%%%%%%%%%%%%%% 此腳本用于VASP的能帶繪制 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc clear %%%%%%%%%%%%%%%%%%%%%%%%%%% 手動修改區(qū) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% EFermi=6.1748; %設(shè)置費米能級 NBands=8; %一共計算的能帶數(shù)目 inter=100; %每段PATH的插值點數(shù)目 NPath=4; %PATH數(shù)目 Ktotal=inter*NPath; %特殊K點的總數(shù)目 xmin=1; %x，y坐標軸的范圍 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 %%%%%%%%%%%%%%%%%%%%%%%%%%%% 坐標軸設(shè)置 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% axis([xmin xmax ymin ymax]); set(gca,'XTick',[1,inter:inter:Ktotal]) set(gca,'XTickLabel','') j=1; %高對稱點符號 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ù)此樓» 猜你喜歡
面上5B能上會嗎？
已經(jīng)有4人回復(fù)
-
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）
已經(jīng)有9人回復(fù)
-
322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級。
已經(jīng)有4人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有14人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有5人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有21人回復(fù)
-
26申博自薦
已經(jīng)有5人回復(fù)
-
296求調(diào)劑
已經(jīng)有5人回復(fù)
-
總分293求調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 如何用matlab畫出石墨烯高對稱點的能帶結(jié)構(gòu)
已經(jīng)有3人回復(fù)
- 在VASP中的腳本
已經(jīng)有4人回復(fù)
- vasp 算能帶k點數(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ù)
- 電機原理與設(shè)計的MATLAB分析
已經(jīng)有140人回復(fù)
- 【求助】請問如何用P4VASP在LINUX系統(tǒng)下畫能帶圖，或者導(dǎo)出數(shù)據(jù)
已經(jīng)有5人回復(fù)
- vasp_bnd的Matlab改進版，更方便~
已經(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能帶橫坐標怎么畫
已經(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坐標的腳本
已經(jīng)有11人回復(fù)
- 【求助】Matlab 四維圖形的繪制。
已經(jīng)有12人回復(fù)
7樓2014-07-01 20:13:38
2樓2011-06-01 19:52:25
3樓2011-06-01 21:15:17
重新試過了，可以了，但圖不好看，還需要不斷改進才行。謝謝 4樓2011-06-03 11:57:02
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩調(diào)劑 +3 | 椰椰。 2026-03-29 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
|
[考研] 309求調(diào)劑 +7 | 誰不是少年 2026-03-29 | 7/350 |
|
|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。，。，。，。i 2026-03-23 | 6/300 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
|
[考研] 上海電力大學(xué)材料防護與新材料重點實驗室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|