开边界与周期边界下的哈密顿量：物理模拟中的能量描述

在物理学中，哈密顿量是描述物理系统总能量的量。在计算机模拟中，哈密顿量通常用来描述某个系统的动力学行为。在开边界和周期边界下，哈密顿量的形式略有不同。

开边界

在开边界的情况下，物理系统的边界是固定的，因此哈密顿量不需要考虑周期性。例如，在一个立方体中，边界是固定的，因此哈密顿量可以写成：

$H = \sum_{i=1}^N \frac{p_i^2}{2m} + \sum_{i<j} V_{ij}(r_{ij})$

其中，$p_i$ 是第 $i$ 个粒子的动量，$m$ 是粒子的质量，$r_{ij}$ 是第 $i$ 和第 $j$ 个粒子之间的距离，$V_{ij}(r_{ij})$ 是粒子之间的相互作用势能。

周期边界

在周期边界的情况下，物理系统的边界是周期性的，因此哈密顿量需要考虑周期性。例如，在一个三维晶格中，哈密顿量可以写成：

$H = \sum_{i=1}^N \frac{p_i^2}{2m} + \sum_{i<j} V_{ij}(r_{ij}) + \sum_{i=1}^N \sum_{\alpha=x,y,z} W_{\alpha}(r_{i,\alpha})$

其中，$p_i$ 是第 $i$ 个粒子的动量，$m$ 是粒子的质量，$r_{ij}$ 是第 $i$ 和第 $j$ 个粒子之间的距离，$V_{ij}(r_{ij})$ 是粒子之间的相互作用势能，$W_{\alpha}(r_{i,\alpha})$ 是第 $i$ 个粒子在 $\alpha$ 方向上与边界相互作用的势能。这个势能可以写成：

$W_{\alpha}(r_{i,\alpha}) = \frac{1}{2} k_{\alpha} (r_{i,\alpha}-L_{\alpha})^2$

其中，$k_{\alpha}$ 是势能常数，$L_{\alpha}$ 是晶格在 $\alpha$ 方向上的长度。这个势能保证了边界与晶格的周期性。