我来帮你梳理这个问题的完整解决思路，从波矢的取值规则到离散求和、连续近似再到数值实现，一步步拆解：

1. 核心公式与系统参数 #

首先明确所有关键的物理量和约束：

• 系统：长度为$L$的一维谐振子链，含$N-1$个谐振子，采用固定边界条件
• 波矢取值：固定边界下，波矢$k$必须满足$kL = m\pi$（$m=1,2,...,N-1$），即$k_m = \frac{m\pi}{L}$，相邻波矢的间隔$\Delta k = \frac{\pi}{L}$
• 色散关系：$\omega(k) = 2\left|\sin\left(\frac{k}{2}\right)\right|$（由$\omega^2(k)=4\sin2(k/2)$推导而来，频率取正值）
• 玻色-爱因斯坦分布：$n(k) = \frac{1}{\exp\left(\frac{\hbar\omega(k)}{k_B T}\right)-1}$，其中$\hbar$为约化普朗克常数，$k_B$为玻尔兹曼常数
• 内能表达式：$U = \sum_{k} \hbar\omega(k)n(k)$，即对所有允许的波矢$k$求和

2. 离散求和的具体形式 #

把波矢的离散取值代入内能公式，得到：
$$
U = \sum_{m=1}^{N-1} \hbar \cdot 2\sin\left(\frac{m\pi}{2L}\right) \cdot \frac{1}{\exp\left(\frac{2\hbar \sin\left(\frac{m\pi}{2L}\right)}{k_B T}\right)-1}
$$
这个形式适合有限$N$（短链）的精确计算，但当$N$很大时，离散求和可以近似为积分，计算效率更高。

3. 连续极限近似（长链情况，$N\to\infty$） #

当链足够长（$N$很大，$L=(N-1)a$，$a$为晶格常数），离散波矢可以看作连续变量，求和转化为积分：

• 波矢的密度为$\frac{L}{\pi}$（因为$k$的范围是$0\sim\pi$，总共有$N-1$个点，$\Delta k=\frac{\pi}{L}$，所以单位波矢区间的点数为$\frac{1}{\Delta k}=\frac{L}{\pi}$）
• 内能的积分形式：
  $$
  U \approx \frac{2\hbar L}{\pi} \int_{0}^{\pi} \sin\left(\frac{k}{2}\right) \cdot \frac{1}{\exp\left(\frac{2\hbar \sin\left(\frac{k}{2}\right)}{k_B T}\right)-1} dk
  $$
  可以通过变量替换简化积分，比如令$x=\frac{k}{2}$，则$dk=2dx$，积分限变为$0\sim\frac{\pi}{2}$，得到：
  $$
  U = \frac{4\hbar L}{\pi} \int_{0}^{\pi/2} \sin x \cdot \frac{1}{\exp\left(\frac{2\hbar \sin x}{k_B T}\right)-1} dx
  $$

4. 数值计算实现（以Python为例） #

如果需要对有限$N$的系统做数值计算，用Python可以快速实现，代码示例如下：

import numpy as np

# 自定义系统参数
L = 10.0          # 链的长度（单位：m）
N = 1000          # 链的总格点数，谐振子数量为N-1
hbar = 1.0546e-34 # 约化普朗克常数（单位：J·s）
kb = 1.3806e-23   # 玻尔兹曼常数（单位：J/K）
T = 300.0         # 温度（单位：K）

# 生成所有允许的波矢
k_values = np.linspace(np.pi/L, (N-1)*np.pi/L, N-1)

# 计算每个波矢对应的频率
omega = 2 * np.sin(k_values / 2)

# 计算玻色-爱因斯坦分布
n_k = 1 / (np.exp(hbar * omega / (kb * T)) - 1)

# 计算总内能
total_U = np.sum(hbar * omega * n_k)

print(f"系统的总内能为：{total_U:.4e} J")

5. 极端温度下的近似结果 #

• 高温极限（$k_B T \gg \hbar\omega(k)$）：此时$\exp\left(\frac{\hbar\omega(k)}{k_B T}\right)\approx1+\frac{\hbar\omega(k)}{k_B T}$，所以$n(k)\approx\frac{k_B T}{\hbar\omega(k)}$，代入内能得$U\approx(N-1)k_B T$，符合经典能量均分定理（每个谐振子贡献$k_B T$的能量）。
• 低温极限（$k_B T \ll \hbar\omega(k)$）：只有低频模式（$k$很小，$\omega(k)\approx k$）被激发，内能近似为$U\propto T^{3/2}$，这是一维玻色系统的典型低温行为。

内容的提问来源于stack exchange，提问作者Joseph Doggo