化学势计算，结果诡异。。。

周一貌似小老板被我忽悠动了，打算用 adress 算算化学势。。。但是现在结果非常诡异。。。。。。大概比正确数值小了3倍。。不知道这是不是出于某种巧合。。。。

对于结果各种解释不能，这两天各种讨论都没有什么积极结果。。。算些新例子看看。。。

写文章，非平衡态模拟缓慢进行中，寄希望于大老板

AdResS 理论文章终于接近完成。小老板确实能写啊，哈哈～～ 加上了 WCA 的例子，而决定不再另发文章。

非平衡态模拟进行中。发现能量守恒对于非平衡态模拟非常重要，直接决定 dynamics 是否正确。局部热浴的方式看上去是对的。需要验证更大的体系，计算开销是一个问题，初步估计最大的体系算一个扰动参数需要大约38天。。。啊。。。1/10年啊。。。参照体系就是全原子的，NVE的，体积很大的模拟结果。所谓“体积很大”是要检验收敛性的～～

大老板说是帮忙写相应理论的文章。。。不知道他有没有时间。。。我们真的需要更快了。。。唉。。揪心啊。。。

WCA 和非平衡态新想法

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

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