[交流] 【求助】如何在平衡状态下进行晶格参数的优化？

采用状态方程拟合得到平衡状态下的体积、体弹模量。
在该体积下，重复上面第二步和第三步，得到平衡状态下的其他的晶胞参数。这一步也就是：在得到了E(V)曲线以后，通过状态方程拟合得到平衡状态下的体积，计算出上面脚本中变量$i的值，并改变$i的循环值，再运行runcell计算一次，得到其他的结构参数c和位置u。

      ALATT0   VOL0 (ULA**3)  B0 (MBAR)     B0PRIME       E0 (EV)      E0 (RY)

alat=  11.63799 b0=  155.6274 b0p=   854.216 E0=  -5.095122
       11.6380    394.07058    155.62739    854.21573   -138.64846    -10.19037
        6.1586