论文卡壳并终于续上了

• 包佳同学是郭教授的学生，研一。我们在合作写论文“Multi-partite nonlocality in the Lipkin-Meshkov-Glick model ”。

我们先前的论文都是近程模型。这次给包佳安排的题目，编程量不大，同时由于处理的是无限远程耦合模型，具有一定的创新性在里面。很适合给研究生练手，并体验完成一篇SCI的整个流程。

下面简要描述下我们发现问题并解决问题的过程。

---

包佳同学在处理LMG模型时，看见了一个没法解释的数值结果：在有限温度下，密度矩阵随外场是连续变化的，然而算出来发现nonlocality是突变的。接着就过年了，接着我各种改稿子，各种备课的事，这个问题就耽误了很久（表达下歉意）。

然后周四晚上和周五早上，终于找到了原因。

• 当时写代码的时候，为了应付服务器断电的情况，中间把优化的向量存储起来了。由于是在服务器上按照最大颗粒并行，这样操作从来没有出错。也没有意识到可能存在的风险。

• 然后放在本地计算。程序把前次计算的结果，当做了新的优化初值。这样导致一条曲线（彼此独立的数据点）都是在高维空间的一个特定的初值附近进行的优化，大错。

错误的结果：突变

• 把存储中间结果的代码注释掉，果然就正常了。

修订后的正确结果

---

• 这次问题解决了之后，按照包佳的认真劲来预计，两个月内就能得到所需的完整数据。
• 这个事也提醒了我自己：
  代码千万行，安全第一条；代码不规范，坐等泪两行~