石墨烯(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 onsurf(xlist,ylist,-E)%能量为负的部分shading flat

狄拉克锥图形

俯视图

 

 

 

 

 

 

 

 

 

 

 

参考：

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

[2]