关于VASP计算二维材料应力-应变曲线的问题
本人最近刚开始接触二维材料的力学性能计算，对二维材料的应力应变曲线怎么求不太清楚，看到一篇文献说在a方向上施加一个应变（比如拉伸1%），然后弛豫b方向上的晶格直到应力分量（conjugate stress components）小于0.01GPa。请问这个是怎么实现的呢？按理说vasp里是没办法固定某一个晶格参数不变而弛豫另一个晶格参数的吧，因为ISIF=2应该是同时弛豫晶格参数和晶胞体积的。我在猜是不是这样的，在a方向上拉伸1%之后，对b方向做状态方程计算，找到能量最低的点，此时认为这是在该拉伸下的稳定晶格常数，然后取此时的OUTCAR里的应力（in kB 那行后面的）作为最终应力，然后再拉伸到2%，以此类推，最后画出应力应变曲线。个人理解啊，希望大神能够解答。。。。还有对于二维材料力学性能的计算一般都需要做哪些工作呢，小白求指导。。。。。
给出文献关于这一段的描述截图：
QQ截图20170809162914.jpg
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
他们自己修改代码，可以只优化一个方向
请问是要怎么修改呢？是需要向脚本一样在计算的时候去赋予权限然后计算吗？
修改constr_cell_relax.F，重新编译，以实现
，
我感觉这个可以依靠poscar的selective dynamic来处理，就是对原子只优化a方向，固定b方向或者反过来。
瞎扯
如果材料是二维 unit cell是长方形 应变加在某一个坐标轴方向 可以通过selective dynamics 将T -> F fix住lattice。如果不是这种情况 就要改vasp的代码 来实现了。