LargeChernNumber/mplot.m at master · franklin5/LargeChernNumber · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
clear
clc
clf
EK = load('spectrum_2109.OUT');
NKX = sqrt(length(EK(:,1)));
kmax = 2;
NK = (NKX-1)/2;
akx = -kmax:2*kmax/(NKX-1):kmax;aky = akx;
[AKX,AKY]=meshgrid(akx,aky);
cutoffP = length(EK(1,:));
energy = zeros(NKX,NKX,cutoffP);
for nk = 0:NKX*NKX-1
    nkx = mod(nk,NKX);
    nky = floor(nk/NKX);
    energy(nkx+1,nky+1,:) = EK(nk+1,:);
end
Tperiod = 26.24;
figure(1)
for np = 1:cutoffP
    temp(:,:) = energy(:,:,np)/(pi/Tperiod);
    mesh(AKX,AKY,temp)
    hold on
end
hold off
xlabel('k_x/k_F')
ylabel('k_y/k_F')
zlabel('E(k_x,k_y)/(pi/T)')
title('hi=2.1,hf=0.9')
view(3)
axis([-kmax kmax -kmax kmax -3 3])
view([0 0])
%axis auto
figure(2)
nky = NK+1;
for np = 1:cutoffP
    temp = energy(:,nky,np)/(pi/Tperiod);
    plot(akx,temp)
    hold on
end
hold off
axis([-kmax kmax -3 3])