【2022.12.9】Lammps+Python 在计算g6(r)时遇到的问题

写在前面

全部log：
【2022.11.16】Lammps+Python+MATLAB在绘制维诺图时遇到的问题

绘制g6( r )执行步骤【updated】

画了个随便看看就不对的图，老板帮忙捋了一遍步骤，发现有点小问题，现更新一下全部的步骤：
1、 导入数据
2、 利用论文中的方法找到最近邻
3、 算出每个粒子的psi6（ r ）
4、 列出所有粒子与粒子之间的距离【将单位从m改为了mm】。并对距离，即g6 （ r ）中的r，进行四舍五入
4、 去除最外圈的粒子，对所有的粒子进行两两之间的距离groupby【也可能是generate】
5、 利用公式算出g6（ r ）

如何检查图像的正确性：

1、画psi6的图
或者直接根据物理性质，观察psi6的模，如如果这个晶体更像完美的晶体，那大部分粒子的psi6应该是接近1的
2、将任意两个psi6相乘
原理同上

不是编程问题，而是数学问题的一个小bug

enmmm，虽然平时做e指数加上角度的时候默认是弧度制，但是写的时候忘记了

废稿2则：

暂定模拟盒子大小：54mmX54mm
T = 30000K的第一帧的