计算化学公社

标题: 2D应力应变求助 [打印本页]

作者
Author: mylifecg 时间: 2017-4-15 22:27
标题: 2D应力应变求助
我通过重新编译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轴不能完全固定，导致了原本矩形原胞有一定的倾斜。计算了好几次都是这个结果，是否这种现象是合理的呢？跪求高手解答！

作者
Author: mylifecg 时间: 2017-4-15 22:28
问题有一处写错了，通过这种编译方法是只优化x轴，不优化y轴和z轴

作者
Author: 卡开发发 时间: 2017-4-16 11:06
1、FCELL(2,2)=0.0，FCELL(3,3)=0.0，那么其余非对角分量怎么处理呢？

2、可以检查一下突变前一帧结构的stress tensor是否出现了非对角分量。

作者
Author: 刘荣 时间: 2024-7-2 16:32
想请问不同应变下的应力值是怎么读取的呀，谢谢

欢迎光临计算化学公社 (http://bbs.keinsci.com/)