matlab狄拉克锥的三维图

石墨烯(graphene)中的狄拉克锥(Dirac cone)图形。直接按照能级公式绘图，公式参考[1]中 energy of the electrons 公式。

%matlab代码
clear;
gamma0=2.8;
a=2.46;%晶格常数的值
xlist=linspace(-2,2,100);
ylist=linspace(-2,2,100);
[kx,ky]=meshgrid(xlist,ylist);
%能量与波矢kx,ky的函数关系
E=sqrt((gamma0.^2).*(1+4.*(cos(ky.*a./2).^2)+4.*(cos(ky.*a./2)).*cos(sqrt(3).*kx.*a./2)));

surf(xlist,ylist,E)%能量为正的部分
hold on
surf(xlist,ylist,-E)%能量为负的部分
shading flat

狄拉克锥图形

俯视图

参考：

[1] https://en.m.wikipedia.org/wiki/Graphene

[2] 石墨烯晶格