Pulay Stress
Pulay Stress 产生的原因是当晶胞体积发生变化时，由于计算时的平面波基组不全（ENCUT 不够大）使得计算的收敛性存在问题，进而导致应力张量计算的错误，最终导致体系的形状、体积以及能量计算结果不正确。
防止Pulay Stress问题出现的两个常见方法：
(1) 选取最佳晶格参数，即BM方程拟合的方法；
(2) 增加ENCUT，允许更多的平面波基组。
ENCUT取值
当体系体积变化时，VASP官方文档给出的ENCUT值设置建议为ENCUT=1.3*max(ENMAX)，或者也可以更为保险地设置为600或700eV。
在INCAR参数的汇总章节中强调同一个计算中的ENCUT参数需要保持一致，那此处增大ENCUT值的操作是否和一致性的要求冲突呢？
其实，本步计算的目的就是获取稳定的晶格参数，再利用获得的晶格参数来保证下一步计算的准确性。可以将晶格参数理解成一个计算中的超参数，这一个步骤就是一个对超参数进行寻优的过程。要实现寻优就必须要有变化的自变量和对应的待优化应变量，这样的类比之下，计算参数的调整也就可以理解了。
