python - 费米曲面图

Question

本质上，我正在尝试绘制二维费米曲面图。即 K 空间中某个平面的 f(n,vec_k)=e_f 的 2D 切割，并进行插值。具体来说，我有一个 numpy 数组:Eigen，具有形状，

Eigen.shape = (100,100,100,10)，其中前三个索引位于向量 vec_k 上，第三个索引是频带索引“n”。如何获得任意曲面 Eigen == e_f 的切割，其中 e_f 是标量？

Answer 1

这通常使用行进立方体算法来解决。您应该在此处查看 MLab 的 contour3d() 函数:http://docs.enthought.com/mayavi/mayavi/auto/mlab_helper_functions.html#mayavi.mlab.contour3d 。这正是你想要的。

contour3d(Eigen[:,:,:, bandNo], contours=e_f)

应该为你做。您还可以为所有带指定多个表面，如下所示:

for bandNo in range(totalBands):
    contour3d(Eigen[:,:,:, bandNo], contours=[eF1, eF2, eF3], opacity=0.4)

您应该得到如下图:http://docs.enthought.com/mayavi/mayavi/_images/enthought_mayavi_mlab_contour3d1.jpg