版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時(shí)熱門版塊排行榜

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

返回列表

【獎(jiǎng)勵(lì)】本帖被評價(jià)3次，作者byin增加金幣 2 個(gè)

本帖產(chǎn)生 1 個(gè) 1ST強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

byin

銅蟲 (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 0 (幼兒園)
金幣: 245
帖子: 395
在線: 174.4小時(shí)
蟲號(hào): 948127

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

把EIGENVAL和.m腳本放一起，修改腳本里的一些參數(shù)“手動(dòng)修改區(qū)”，然后運(yùn)行腳本即可。見笑了。

%%%%%%%%%%%%%%%%%%%% 此腳本用于VASP的能帶繪制 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc
clear
%%%%%%%%%%%%%%%%%%%%%%%%%%% 手動(dòng)修改區(qū) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

EFermi=6.1748;                %設(shè)置費(fèi)米能級
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;                               %高對稱點(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ù)此樓

» 猜你喜歡

調(diào)劑310 已經(jīng)有12人回復(fù)
化學(xué)0703 調(diào)劑 306分一志愿211 已經(jīng)有6人回復(fù)
26申博自薦已經(jīng)有4人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有20人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有5人回復(fù)
289求調(diào)劑已經(jīng)有3人回復(fù)
085601材料工程找調(diào)劑已經(jīng)有14人回復(fù)
一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師已經(jīng)有6人回復(fù)
334分一志愿武理材料求調(diào)劑已經(jīng)有7人回復(fù)
275求調(diào)劑已經(jīng)有18人回復(fù)

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

如何用matlab畫出石墨烯高對稱點(diǎn)的能帶結(jié)構(gòu) 已經(jīng)有3人回復(fù)
在VASP中的腳本已經(jīng)有4人回復(fù)
vasp 算能帶k點(diǎn)數(shù)目問題已經(jīng)有10人回復(fù)
推薦幾個(gè)vasp處理計(jì)算結(jié)果的python腳本已經(jīng)有66人回復(fù)
計(jì)算模擬需要的服務(wù)器硬件 MS,VASP,matlab等等吧已經(jīng)有7人回復(fù)
NBANDS 的設(shè)置問題已經(jīng)有6人回復(fù)
電機(jī)原理與設(shè)計(jì)的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有幾個(gè)腳本文件？已經(jīng)有8人回復(fù)
DOSCAR畫態(tài)密度已經(jīng)有16人回復(fù)
VASP能帶橫坐標(biāo)怎么畫已經(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ù)

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

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

lrjhxh

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 3654.9
帖子: 267
在線: 371.8小時(shí)
蟲號(hào): 102835

★★★ 三星級,支持鼓勵(lì)

請問怎么會(huì)出現(xiàn)以下錯(cuò)誤：
?? 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ù)（每個(gè)i中的各個(gè)j），做出曲線

贊一下

回復(fù)此樓

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

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

byin

銅蟲 (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 0 (幼兒園)
金幣: 245
帖子: 395
在線: 174.4小時(shí)
蟲號(hào): 948127

引用回帖:

Originally posted by lrjhxh at 2011-06-01 19:52:25:
請問怎么會(huì)出現(xiàn)以下錯(cuò)誤：
?? 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

贊一下

回復(fù)此樓

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

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

lrjhxh

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 3654.9
帖子: 267
在線: 371.8小時(shí)
蟲號(hào): 102835

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

贊一下

回復(fù)此樓

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

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

lrjhxh

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 3654.9
帖子: 267
在線: 371.8小時(shí)
蟲號(hào): 102835

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

贊一下

回復(fù)此樓

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

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

白小紙bql

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 52.5
帖子: 49
在線: 18.1小時(shí)
蟲號(hào): 3009957

★★★★★ 五星級,優(yōu)秀推薦

這個(gè)星級是給lz評價(jià)么？
額，，正題，，lz您好~我是一只新蟲子，謝謝您寫的程序~好頂贊。
但是我有兩個(gè)問題，一是不同高對稱點(diǎn)之間的間距應(yīng)該會(huì)不同，這個(gè)程序處理出來時(shí)均勻的；另一個(gè)可能是精度的問題，我在用這個(gè)程序處理石墨烯能帶的時(shí)候，在Dirac錐處出現(xiàn)了平直部分；
lz能幫忙解決一下上面兩個(gè)問題么？感激不盡！

g.png

贊一下

回復(fù)此樓

7樓2014-07-01 20:13:38

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

簡單回復(fù)

zhrq7812026樓

2014-05-01 10:57 回復(fù)

三星好評

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 byin 的主題更新

返回列表

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

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

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

[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5	hanamiko 2026-03-29	5/250	2026-03-29 16:46 by 學(xué)員8dgXkO
[考研] 275求調(diào)劑 +13	Micky11223 2026-03-25	18/900	2026-03-29 16:23 by 學(xué)員8dgXkO
[考研] 298求調(diào)劑 +3	種圣賜 2026-03-29	3/150	2026-03-29 12:06 by longlotian
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥求調(diào)劑 +4	篤行致遠(yuǎn). 2026-03-28	4/200	2026-03-29 11:55 by longlotian
[考研] 317分一志愿南理工材料工程本科湖工大求調(diào)劑 +9	芋泥小鈴鐺 2026-03-28	9/450	2026-03-29 11:46 by longlotian
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +3	真名有冰 2026-03-29	4/200	2026-03-29 08:47 by qingfeng258
[考研] 材料277求調(diào)劑 +7	min3 2026-03-24	7/350	2026-03-28 11:39 by xuxiang
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7	hanamiko 2026-03-27	7/350	2026-03-28 09:56 by 李上岸0921
[考研] 考研調(diào)劑 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 085600，材料與化工321分調(diào)劑 +4	大饞小子 2026-03-27	6/300	2026-03-27 14:11 by 松花缸1201
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 292求調(diào)劑 +4	求求了收下我吧?/a> 2026-03-26	4/200	2026-03-27 10:37 by zhshch
[考研] 324求調(diào)劑 +8	hanamiko 2026-03-26	10/500	2026-03-27 08:06 by hypershenger
[考研] 一志愿南京郵電大學(xué) 288分材料考研求調(diào)劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4	gdlf9999 2026-03-24	4/200	2026-03-25 23:01 by boxking200
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 生物技術(shù)與工程 +3	1294608413 2026-03-25	4/200	2026-03-25 18:02 by 1294608413
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3	an10101 2026-03-24	7/350	2026-03-25 00:00 by 山鬼0-
[考研] 調(diào)劑 +4	13853210211 2026-03-24	4/200	2026-03-24 19:44 by ms629

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

24小時(shí)熱門版塊排行榜

byin

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

» 猜你喜歡

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

lrjhxh

★★★ 三星級,支持鼓勵(lì)

byin

lrjhxh

lrjhxh

白小紙bql

★★★★★ 五星級,優(yōu)秀推薦

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

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