关于能量-应变方法计算二维六角格子力学性质的疑惑

physics_xw

      近期在使用能量-应变方法计算一系列二维材料的力学性质，以其中的CrI3为例，其为二维六角格子，a=b，γ=120.
      按照vaspkit官网教程设置VPKIT.in，更改3D为2D，其余不变，然后使用vaspkit 201功能生成了两个计算子文件夹C11和C12，C11和C12中分别又有7个子文件夹。C12中的子文件夹分别对应施加双轴等拉伸或压缩应变后的结构，这个容易理解。而C11文件夹中的结构似乎并非是单轴应变后的结构，因为a和b同时发生了变化，并且两者变化并不相等，此外，γ角也不再严格等于120度。请问C11文件夹中施加应变的标准是什么？这个文件夹的含义是什么？
      按照VASPKIT的介绍文章(Computer Physics Communications 267 (2021) 108033)中3.1节倒数2-3段及附录B Table B.4的描述，为计算C11和C12，需要给结构施加两种应变，分别为单轴应变和双轴应变，在文中这两种应变方式分别记为1和2.对能量-应变的散点做二次拟合，根据拟合系数可以求解出C11和C12。以上是我对这种方法的理解，不知道是否正确，如有错误还请老师和同学们指正。

总结一下疑惑：
1.按照VASPKIT文章中的描述，需要施加两种模式的应变，而这两种模式与VASPKIT生成的C11和C12中的结构并无完全的对应的关系？如果是这样，那C11和C12中的应变是以何种标准施加的呢？C11和C12文件夹的命名含义是什么？通过某一个文件夹中的结构可分别单独计算出C11或C12？
2.在vaspkit为CrI3生成的计算子文件夹C11中，结构发生了并不属于文章中所描述的2种模式的变化，尤其是120度的夹角发生了变化。我原先猜想是否C11中结构是两种应变模式的结合，但显然无论怎样结合120度的夹角都不应该变化。想请教各位，这里的结构变化正常吗？是否我对这种结构变化的理解有问题？

附件中提供了CrI3原始的POSCAR及C11中0.015拉伸及压缩应变的POSCAR，以供各位参考。
所用vaspkit版本号为1.3.3，感谢王伟老师及VASPKIT开发团队创造出这么便利的工具。

vaspkit

对于二维六角格子，其应变模式可参考(Computer Physics Communications 267 (2021) 108033)中的Table B.4。不妨使用文献中的公式8和9手动计算不同应变下晶格基矢变化，并与VASPKIT生成的做比较，若还有疑惑再讨论。

yrui

你好 我也遇到了同样的问题，请问你的问题解决了吗？

vaspkit

对于二维体系，单轴拉伸是指在笛卡尔坐标下沿着x或y轴方向，由于六角基矢b在x轴有分量，因此也会变化。