2D应力应变求助

mylifecg

我通过重新编译vasp code，方法如下：打开vasp文件的constr_cell_relax.F 在“REAL(q) FCELL(3,3)”下面的任意行加上FCELL(2,2)=0.0，FCELL(3,3)=0.0然后重新编译。这里实现的是只优化 x 轴，只优化 y 和 z 轴。
我的问题是，
1. 这种编译是否有问题？
2.当我对2D结构沿着Y方向（图中箭头指向）拉伸到应力突变的时候，突变后的结构，y轴不能彻底固定，优化前后，Y轴晶格矢量有微小变化，但是突变前可以实现完全固定y轴，也就是优化前后y轴完全不变。
附件是我应力突变曲线以及突变前后对应的结构图


为何突变前可以完全 固定y轴，应力突变后，Y轴不能完全固定，导致了原本矩形原胞有一定的倾斜。计算了好几次都是这个结果，是否这种现象是合理的呢？跪求高手解答！

mylifecg

问题有一处写错了，通过这种编译方法是只优化x轴，不优化y轴和z轴

卡开发发

1、FCELL(2,2)=0.0，FCELL(3,3)=0.0，那么其余非对角分量怎么处理呢？
2、可以检查一下突变前一帧结构的stress tensor是否出现了非对角分量。

刘荣

想请问不同应变下的应力值是怎么读取的呀，谢谢