版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

【獎勵】本帖被評價3次，作者byin增加金幣 2 個

本帖產(chǎn)生 1 個 1ST強帖，點擊這里進行查看

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

byin

銅蟲 (正式寫手)

1ST強帖: 1
應助: 0 (幼兒園)
金幣: 245
帖子: 395
在線: 174.4小時
蟲號: 948127

[資源] 用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 ]

回復此樓

» 猜你喜歡

343求調(diào)劑085601 已經(jīng)有3人回復
08工科求調(diào)劑286 已經(jīng)有3人回復
337求調(diào)劑已經(jīng)有3人回復
329求調(diào)劑，一志愿西北工業(yè)大學，材料工程（085601）已經(jīng)有10人回復
一志愿北京化工大學材料與化工（085600）296求調(diào)劑已經(jīng)有16人回復
291求調(diào)劑已經(jīng)有6人回復
299求調(diào)劑已經(jīng)有10人回復
面上5B能上會嗎？已經(jīng)有4人回復
322求調(diào)劑：一志愿湖南大學材料與化工（085600），已過六級。已經(jīng)有4人回復
085600 材料與化工 329分求調(diào)劑已經(jīng)有14人回復

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

如何用matlab畫出石墨烯高對稱點的能帶結(jié)構(gòu) 已經(jīng)有3人回復
在VASP中的腳本已經(jīng)有4人回復
vasp 算能帶k點數(shù)目問題已經(jīng)有10人回復
推薦幾個vasp處理計算結(jié)果的python腳本已經(jīng)有66人回復
計算模擬需要的服務器硬件 MS,VASP,matlab等等吧已經(jīng)有7人回復
NBANDS 的設(shè)置問題已經(jīng)有6人回復
電機原理與設(shè)計的MATLAB分析已經(jīng)有140人回復
【求助】請問如何用P4VASP在LINUX系統(tǒng)下畫能帶圖，或者導出數(shù)據(jù) 已經(jīng)有5人回復
vasp_bnd的Matlab改進版，更方便~ 已經(jīng)有3人回復
請教VASP CI-NBE分析腳本的使用問題已經(jīng)有6人回復
matlab多維插值原理是什么？已經(jīng)有5人回復
VASP有幾個腳本文件？已經(jīng)有8人回復
DOSCAR畫態(tài)密度已經(jīng)有16人回復
VASP能帶橫坐標怎么畫已經(jīng)有6人回復
【活動or資源】vasp晶體結(jié)構(gòu)優(yōu)化已經(jīng)有142人回復
【求助】vasp作業(yè)提交腳本已經(jīng)有5人回復
【求助】求將vasp的輸出文件轉(zhuǎn)化成xyz坐標的腳本已經(jīng)有11人回復
【求助】Matlab 四維圖形的繪制。已經(jīng)有12人回復

1樓 2011-05-30 21:22:10

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lrjhxh

木蟲 (小有名氣)

應助: 2 (幼兒園)
金幣: 3654.9
帖子: 267
在線: 371.8小時
蟲號: 102835

看來問題還是不能解決。請問VASP計算之后的輸出結(jié)果中，EIGENVAL中的第一~第三行表示什么意思？畫圖時各高對稱K點的坐標是如何定義的？怎樣轉(zhuǎn)換各K點的坐標來畫圖？你所提供的腳本中，讀出每條能帶的數(shù)據(jù)算法是怎樣的？從CHGCAR中獲得電荷密度圖的Matlab腳本有沒有？你是通常使用什么工具進行輸出結(jié)果的分析？謝謝！

贊一下

回復此樓

5樓2011-11-05 21:46:17

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 7 個回答

lrjhxh

木蟲 (小有名氣)

應助: 2 (幼兒園)
金幣: 3654.9
帖子: 267
在線: 371.8小時
蟲號: 102835

★★★ 三星級,支持鼓勵

請問怎么會出現(xiàn)以下錯誤：
?? Attempted to access data(731,2); index out of bounds because
size(data)=[728,4].

Error in ==> bandstructure at 31
P(j,: )=[j,data((NBands+1)*(j-1)+1+i,2)-EFermi];
%讀出每條能帶的數(shù)據(jù)（每個i中的各個j），做出曲線

贊一下

回復此樓

2樓2011-06-01 19:52:25

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

byin

銅蟲 (正式寫手)

1ST強帖: 1
應助: 0 (幼兒園)
金幣: 245
帖子: 395
在線: 174.4小時
蟲號: 948127

引用回帖:

Originally posted by lrjhxh at 2011-06-01 19:52:25:
請問怎么會出現(xiàn)以下錯誤：
?? Attempted to access data(731,2); index out of bounds because
size(data)=[728,4].

Error in ==> bandstructure at 31
P(j,: )=[j,data((NBands+1)*(j-1)+1+i,2 ...

對不起，可能有bug，請將EIGNVAL發(fā)送到我郵箱，我試試看。

byin87@gmail.com

贊一下

回復此樓

3樓2011-06-01 21:15:17

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lrjhxh

木蟲 (小有名氣)

應助: 2 (幼兒園)
金幣: 3654.9
帖子: 267
在線: 371.8小時
蟲號: 102835

重新試過了，可以了，但圖不好看，還需要不斷改進才行。謝謝

贊一下

回復此樓

4樓2011-06-03 11:57:02

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 7 個回答

☆ 無星級 ★ 一星級 ★★★ 三星級 ★★★★★ 五星級

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

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

[考研] 296求調(diào)劑 +5	彼岸t 2026-03-29	5/250	2026-03-29 16:54 by 熊一刀
[考研] 總分293求調(diào)劑 +7	加一一九 2026-03-25	10/500	2026-03-29 16:53 by 楊楊楊紫
[考研] 327求調(diào)劑 +6	汲亦昊 2026-03-29	6/300	2026-03-29 13:40 by peike
[考研] 352分-085602-一志愿985 +5	海納百川Ly 2026-03-29	5/250	2026-03-29 09:57 by Sjndkwm
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6	烤麥芽 2026-03-27	8/400	2026-03-28 10:40 by 烤麥芽
[考研] 295求調(diào)劑 +5	1428151015 2026-03-27	6/300	2026-03-28 04:04 by fmesaito
[考研] 266求調(diào)劑 +11	陽陽哇塞 2026-03-27	12/600	2026-03-27 17:56 by yu221
[考研] 279 分求調(diào)劑 +4	睡個好覺_16 2026-03-24	4/200	2026-03-27 15:05 by 醉在風里
[考研] 調(diào)劑 +3	李嘉圖·S·路 2026-03-27	3/150	2026-03-27 11:19 by wangjy2002
[考研] 0703化學338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 材料調(diào)劑 +8	匹克i 2026-03-23	8/400	2026-03-27 08:11 by hypershenger
[考研] 294分080500材料科學與工程求調(diào)劑 +4	柳溪邊 2026-03-26	4/200	2026-03-26 21:14 by XPU李慶
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 材料與化工304求B區(qū)調(diào)劑 +3	邱gl 2026-03-26	6/300	2026-03-26 18:03 by 邱gl
[考研] 347求調(diào)劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper
[考研] 340求調(diào)劑 +5	話梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin
[有機交流] 有機合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

byin

[資源] 用matlab腳本繪制vasp中的能帶圖

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

lrjhxh

lrjhxh

★★★ 三星級,支持鼓勵

byin

lrjhxh

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: