使用 numba 加速蒙特卡洛流体模拟时，若函数依赖全局数组（如 `positions`），jit 编译会捕获其初始快照而非运行时值，导致能量计算错误、接受率异常升高——根本原因在于 numba 不支持动态全局变量引用。
在基于 Lennard-Jones 势的分子动力学或蒙特卡洛（MC）模拟中，正确使用 Numba 的核心原则是：所有被 @njit 装饰的函数必须显式接收所需数据作为参数，严禁访问模块级全局变量。原始代码中，calc_dist、calc_energy_particle 等函数直接读取全局 positions 和 L，看似简洁，实则触发了 Numba 的“编译期绑定”行为——即在首次调用时将 positions 的当前值（通常是未初始化或过时状态）固化为常量，后续位置更新完全不被 JIT 函数感知。
例如，calc_dist(i, j) 在未传入 positions 时，Numba 无法识别其应随模拟步长动态变化，导致距离计算始终基于初始晶格构型（甚至可能因内存未初始化而产生随机噪声），进而使势能 calc_LJ_potential 返回极大正值，Metropolis 判据 np.exp(-beta * delta_energy) 恒接近 1，造成虚假高接受率（如报告中的 1.999% 实为 ~100% 的整数溢出显示误差），最终能量与压强严重偏离物理预期（+114 vs −4.78）。
✅ 正确做法是重构所有 @njit 函数，将 positions 作为首参显式传递：
@njit
def calc_dist(positions, i, j):
p1, p2 = positions[i], positions[j]
dx, dy, dz = p1[0]-p2[0], p1[1]-p2[1], p1[2]-p2[2]
# 周期性边界修正（略）
return np.sqrt(dx**2 + dy**2 + dz**2)
@njit
def calc_energy_particle(positions, p):
energy = 0.0
for j in range(N):
if j != p:
dist = calc_dist(positions, p, j)
if dist <= cutoff:
energy += 4.0 * ((1.0/dist)**12 - (1.0/dist)**6)
return energy
@njit
def move_particle(positions, p, displ=0.5):
new_pos = positions[p].copy()
for dim in range(3):
new_pos[dim] += (np.random.rand() - 0.5) * displ
if new_pos[dim] >= L: new_pos[dim] -= L
elif new_pos[dim] < 0.0: new_pos[dim] += L
return new_pos
同时，主循环中需同步更新调用签名：
# ✅ 正确：每次调用均传入当前 positions prev_energy = calc_energy_particle(positions, particle_index) positions[particle_index] = move_particle(positions, particle_index) new_energy = calc_energy_particle(positions, particle_index)
⚠️ 其他关键注意事项：
- 常量也建议显式传入：虽然 L、cutoff 等标量在 Numba 中通常可安全作为闭包变量，但为明确性和可维护性，推荐统一通过参数或 @njit(fastmath=True) 配合局部常量声明；
- 避免 np.random.rand() 在 @njit 外混用：Numba 的 np.random 支持有限，应确保所有随机数生成均在 @njit 函数内完成（本例已合规）；
- 验证 JIT 编译状态：添加 print(calc_energy_total.inspect_types()) 可确认函数是否成功编译为机器码，避免隐式回退到 Python 解释执行；
- 调试技巧：临时移除 @njit 并插入 print() 日志，对比有无加速时 delta_energy 的分布，快速定位逻辑偏差点。
遵循“参数化一切”的设计范式后，Numba 不仅恢复物理结果的准确性（能量 ≈ −4.77，压力 ≈ 5.19），更带来 5–10 倍的性能提升。这印证了一个重要准则：高性能科学计算的可维护性，始于清晰的数据流契约，而非对全局状态的隐式依赖。