用mathematic计算n+1个格点的链的哈密顿量

假设这n+1个格点分别为1, 2, ..., n+1，且相邻两个格点之间的距离为1，那么这条链的哈密顿量可以表示为：

H = -J * ∑i=1^n SiSi+1 - h * ∑i=1^n Si - h * Sn+1

其中，Si表示第i个格点上的自旋，J表示相邻两个自旋之间的交换耦合常数，h表示外加磁场的大小。

用Mathematica可以直接定义这个哈密顿量：

H = -J * Sum[Subscript[S, i] Subscript[S, i + 1], {i, 1, n}] - h * Sum[Subscript[S, i], {i, 1, n+1}]

其中，Subscript[S, i]表示第i个格点上的自旋。

然后，可以对这个哈密顿量进行求解，得到其能量本征值和能量本征态。具体方法包括对角化矩阵、求解特征方程等。