WCA 和非平衡态新想法

WCA 实验成功，发现 g^3 在 HY I 结果有明显区别。

非平衡态模拟关键是提出合适的边界条件。在全原子区域足够大时吸收边界条件就应该可以。写代码做验证，嗯。