如何拟合晶体的弹性常数？

珊珊来迟

请问在用band model计算迁移率的时候公式中弹性常数这个参数是怎么拟合的？我的体系是一个三斜晶系，下面是我的方法，不知是否正确：
1.先优化晶胞得到最稳定结构；
  该步我才用的方法是先用ISIF=3优化晶格参数（得到a0, b0, c0），再用ISIF=2弛豫离子位置得到最终结构（结构0），并提取该步的能量。
2.比如要得到沿a轴方向的弹性常数，固定b0, c0，将a0改变1%，2%，3%，-1%，-2%，-3%，用ISIF=2固定晶格参数进行优化，得到各晶格常数下的能量；
3.将前两步得到的共七个结构的七组数据：纵坐标能量E以及横坐标-0.03，-0.02，-0.01, 0.00, 0.01, 0.02, 0.03 作图，用二次函数拟合法计算弹性常数。

现在问题时我做出来的E-strain图是一个单调递减的变化趋势，而不是一个抛物线，也就是说应力为0的点不是能量最低的结构。我不知道是我第一步优化的结构有问题还是我拟合的方法不正确？


期待老师的回复！

卡开发发

应该是构型优化的问题，isif=3优化晶格时最好把force和stress的收敛性测试做好。

chenhong

你好，你的问题解决了吗？
我在进行计算弹性常数出现了如下错误，VERY BAD NEWS! internal error in subroutine PRICEL (probably precision problem, try to change SYMPREC in INCAR ?):
Sorry, number of cells and number of vectors did not agree.       8
我根据别人发的帖子 添加了SYMPREC=0.001，但是还是出现了如下错误！谢谢你的回复。