根据平衡条件，即稳定结合时
求得 则可以求得每一摩尔氢分子晶体的结合能为
计算中没有考虑零点能的量子修正，这是造成理论和实验值之间巨大差别的原因。
是1.5的图 是3.2的图
是3.3的图
3.2 讨论N 个原胞的一维双原子链（相邻原子间距为a），其2N 个格波解，当M = m时与一维单原子链的结果一一对应。
解：如图所示，质量为M 的原子位于2n ?1, 2n +1, 2n + 3?? 质量为m 的原子位于2n, 2n +
2,2n + …. 牛顿运动方程为
每个原胞有两个，共有2N 个形式相同的独立方程。形式解为：
代回运动方程有
这是一个以 A 、B 为未知量的齐次线性方程组，有解的条件是系数行列式为零
有两组不同的解：
q 的取值范围是：对应于每个 q 值，有两个格波，共计2N 个格波。
当M = m时，两组解变为
初看似乎仍为双值函数，
但是由于原来取布里渊区为为实际区域大小的一半，所以
当我们把布里渊区扩展为时，就不必用双值表示了，变为
这时当然就没有光学波了
3.3 考虑一双原子链的晶格振动，链上最近邻原子间力常数交替为c 和10c 。令两种原子质量相同，且最近邻间距为a / 2。求在k = 0和k =π / a处的ω(k)。大略地画出色散关系。此问题模拟如2 H 这样的双原子分子晶体。
解：可以这样考虑这个问题， H2分子组成一维晶体，分子内部的相互作用较强，力常数为 10c ，相邻的原子间作用较弱，力常数为c ，第s 个分子中的两个原子的位移分别用
表示：
将试探解
代入上式 有
是u ，ν 的线性齐次方程组，存在非零解的条件为
当k = 0时，
当k =π / a 时
ω2与k 的关系如下图所示，这是一个双原子(例如2 H )晶体。令
3.6 求出一维单原子链的频率分布函数。 解：设单原子链长度 L = Na
频率分布函数