Mathematica计算格林函数矩阵元 - 10个格点的示例

以下是使用Mathematica计算10个格点的格林函数矩阵元的示例代码：

1. 定义晶格和哈密顿量

lat = Tuples[Range[-5, 5], 2];
t = 1;
h = SparseArray[{
   Band[{1, 2}, {Length[lat], Length[lat]}] -> t,
   Band[{2, 1}, {Length[lat], Length[lat]}] -> t
  }];

2. 计算格林函数矩阵元

gf = Inverse[ω + I η - h];
gfm = Table[gf[[i, j]], {i, 1, Length[lat]}, {j, 1, Length[lat]}];

其中，'ω'和'η'是自由参数，可以根据具体问题设定。'gfm'就是10个格点的格林函数矩阵元。可以通过'MatrixForm'来显示矩阵的形式：

MatrixForm[gfm]