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

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

(* 定义哈密顿量 *)
hamiltonian[n_] := 
  Table[
   If[i == j + 1 || i == j - 1, -1, 0], {i, 1, n}, {j, 1, n}];

(* 计算能谱 *)
eigenvalues[n_] := Eigenvalues[hamiltonian[n]];

(* 输出能谱 *)
TableForm[Transpose[{Range[100], eigenvalues[100]}], 
 TableHeadings -> {None, {'Index', 'Energy'}}]

执行以上代码将会输出 100 个格点紧束缚模型的能谱，其中第一列是能级编号，第二列是能级能量。