求讲解vasp利用改变晶格常数进行加压的方法

碧雪霜花剑

最近在进行二维材料的研究，需要加压去研究它的性质，但在大师兄learn vasp教程中没找到具体的方法，在小*虫中搜索也没找到具体的参数设置，好像有改变pstress参数和改变晶格常数两种方法？找到加压的文献但是没有具体的过程，所以想问下通过改变晶格常数加压的具体操作是怎么实现的？有没有什么文献教程是讲这个的，希望各位解答

欧拉角

优化的时候设置ISIF=3，可以得到压强接近零的晶格常数，然后设置ISIF=2，改变晶格常数（加压的话晶格常数往小改），然后计算opt收敛之后，在OUTCAR里面搜索带有“external pressure”的行，就可以找到相应晶格常数下晶体受到的压强。可以通过插值来得到想要的压强所对应的晶格常数。

注意ISIF=3的计算需要设置PREC=high，并且截断能要比平时上调30%。
详情可以看这里
https://www.vasp.at/wiki/index.p ... ns_and_Pulay_stress

碧雪霜花剑

欧拉角 发表于 2021-1-24 17:45
优化的时候设置ISIF=3，可以得到压强接近零的晶格常数，然后设置ISIF=2，改变晶格常数（加压的话晶格常数往 ...

感谢回答，这两天我在看文献还有侯老师的相关教程。
如果我不需要计算具体加压的数值，只是计算成倍数的应力如1%~5%的改变，那么仅改变晶格常数，例如施加1%的拉力，只需要a，b晶格常数*1.01，也不需要添加pstress参数，这样的方法是可行的吗？
另外所谓的单轴应变和双轴应变之类的是只改变a和同时改变a，b的区别吗？
感谢

欧拉角

碧雪霜花剑 发表于 2021-1-25 17:43
感谢回答，这两天我在看文献还有侯老师的相关教程。
如果我不需要计算具体加压的数值，只是计算成倍数的 ...

那样的话，施加的就不是拉力而是施加的应变。如果你需要的就是看应变带来的影响，那当然就可以。这种情况下不需要设置pstress参数。
关于单轴应变和双轴应变，你说的是对的。可以参考这里：
单轴应变：https://csmbrannon.net/2012/08/0 ... nd-uniaxial-strain/
双轴应变：https://eng.libretexts.org/Books ... s/11.04:_Section_4-

碧雪霜花剑

欧拉角 发表于 2021-1-26 09:52
那样的话，施加的就不是拉力而是施加的应变。如果你需要的就是看应变带来的影响，那当然就可以。这种情况 ...

感谢

sunjk

isif=3结构优化后，改变晶格常数，应该用isif=4？还是isif=2？？

sunjk

sunjk 发表于 2022-1-18 19:31
isif=3结构优化后，改变晶格常数，应该用isif=4？还是isif=2？？

哦哦，就是isif=2合理

飞翔的猪

欧拉角 发表于 2021-1-24 17:45
优化的时候设置ISIF=3，可以得到压强接近零的晶格常数，然后设置ISIF=2，改变晶格常数（加压的话晶格常数往 ...

请问老师，您这里的“优化的时候设置ISIF=3”指的是优化bulk还是优化slab？

人心向背

欧拉角 发表于 2021-1-24 17:45
优化的时候设置ISIF=3，可以得到压强接近零的晶格常数，然后设置ISIF=2，改变晶格常数（加压的话晶格常数往 ...

请问利用改变晶格常数的方法模拟加压是否必须采用分数坐标，对笛卡尔坐标来说还算等比例压缩吗

zzffzz33

欧拉角 发表于 2021-1-24 17:45
优化的时候设置ISIF=3，可以得到压强接近零的晶格常数，然后设置ISIF=2，改变晶格常数（加压的话晶格常数往 ...

你好想请教下，这种改变压力的方法和AIMD改变压力的方法有什么不同吗？另外在DFT计算中是不是还可以调整PSTRESS参数来改变压力？他们之间有啥区别呢？