vaspkit计算形变势

小菜鸡

想请问一下vaspkit目前有计算形变势的功能么，英文手册中的简单tutorials里面没有看到相应的指导，若是未提供的话，请问可以在哪里学习和复现三维材料形变势的计算呢，感谢回答

vaspkit

不能一步得到形变势，可以在计算弹性常数的基础之上再执行一次能带计算并设置LVHAR=T，对于块体材料再用VASPKIT-911提取带边并手动减去真空能级，对于二维材料直接调用VASPKIT-927可以提取带边相对于真空能级数据。

小菜鸡

vaspkit 发表于 2022-3-3 19:04
不能一步得到形变势，可以在计算弹性常数的基础之上再执行一次能带计算并设置LVHAR=T，对于块体材料再用VAS ...

首先感谢您的回答，但是我还是有一些小问题想向您咨询，您说的在计算弹性常数的基础上，是利用计算弹性常数时所生成的POSCAR作为形变后的结构进行计算，进而和未形变前计算的结果进行对比对吗？我在之前的学习中了解到块体材料的形变势计算需要在静态计算过程中添加一个icorelevel=1的参数，得到指定原子的1s能级，然后将这个能级和您讲解中提到的由带边减去真空能级所得到的真正的CBM或者VBM一起带入形变势计算公式中。请问我这样的理解正确吗？

vaspkit

小菜鸡 发表于 2022-3-3 19:46
首先感谢您的回答，但是我还是有一些小问题想向您咨询，您说的在计算弹性常数的基础上，是利用计算弹性常 ...

计算弹性常数后，执行cp CONTCAR POSCAR，能带计算，VASPKIT-911提取带边，手动减去1s能级，主要形变势计算应变变化范围小于弹性常数计算的。

一点一点会

您好，在计算三维材料的形变势的时候，通过拟合发现是二次线性的关系，请问您知道应该怎么处理三维材料的形变势才是正确的吗

vaspkit

一点一点会 发表于 2022-8-19 19:18
您好，在计算三维材料的形变势的时候，通过拟合发现是二次线性的关系，请问您知道应该怎么处理三维材料的形 ...

施加应变不能太大（典型值小于2%）；再者应变构型的带边位置可能会发生改变，因此须保证K-points包含带边带边位置。

LCN

vaspkit 发表于 2022-3-3 19:04
不能一步得到形变势，可以在计算弹性常数的基础之上再执行一次能带计算并设置LVHAR=T，对于块体材料再用VAS ...

我想问一下，生成的C22和C11_C12_C66文件，是分别对应的在y方向施加应力和在x方向施加应力的情况吗？