Mathematica 计算 100 个格点紧束缚模型能谱

以下是使用 Mathematica 计算 100 个格点紧束缚模型的能谱的代码：

n = 100;(*格点数*)
h = Table[If[i == j + 1 || i == j - 1, -1, 0], {i, n}, {j, n}];(*构造哈密顿矩阵*)
{eval, evec} = Eigensystem[h];(*计算哈密顿矩阵的本征值和本征向量*)
ListPlot[eval, PlotRange -> All, AxesLabel -> {'n', 'Energy'}](*绘制能谱图*)

运行代码后，将会得到一个包含 100 个本征能量的列表，并绘制出对应的能谱图。其中，x 轴为能级编号，y 轴为能量。