计算化学公社

标题: 求讲解vasp利用改变晶格常数进行加压的方法 [打印本页]

作者
Author:
碧雪霜花剑    时间: 2021-1-24 11:39
标题: 求讲解vasp利用改变晶格常数进行加压的方法
最近在进行二维材料的研究,需要加压去研究它的性质,但在大师兄learn vasp教程中没找到具体的方法,在小*虫中搜索也没找到具体的参数设置,好像有改变pstress参数和改变晶格常数两种方法?找到加压的文献但是没有具体的过程,所以想问下通过改变晶格常数加压的具体操作是怎么实现的?有没有什么文献教程是讲这个的,希望各位解答

作者
Author:
欧拉角    时间: 2021-1-24 17:45
优化的时候设置ISIF=3,可以得到压强接近零的晶格常数,然后设置ISIF=2,改变晶格常数(加压的话晶格常数往小改),然后计算opt收敛之后,在OUTCAR里面搜索带有“external pressure”的行,就可以找到相应晶格常数下晶体受到的压强。可以通过插值来得到想要的压强所对应的晶格常数。

注意ISIF=3的计算需要设置PREC=high,并且截断能要比平时上调30%。
详情可以看这里
https://www.vasp.at/wiki/index.p ... ns_and_Pulay_stress
作者
Author:
碧雪霜花剑    时间: 2021-1-25 17:43
欧拉角 发表于 2021-1-24 17:45
优化的时候设置ISIF=3,可以得到压强接近零的晶格常数,然后设置ISIF=2,改变晶格常数(加压的话晶格常数往 ...

感谢回答,这两天我在看文献还有侯老师的相关教程。
如果我不需要计算具体加压的数值,只是计算成倍数的应力如1%~5%的改变,那么仅改变晶格常数,例如施加1%的拉力,只需要a,b晶格常数*1.01,也不需要添加pstress参数,这样的方法是可行的吗?
另外所谓的单轴应变和双轴应变之类的是只改变a和同时改变a,b的区别吗?
感谢
作者
Author:
欧拉角    时间: 2021-1-26 09:52
碧雪霜花剑 发表于 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-
作者
Author:
碧雪霜花剑    时间: 2021-1-26 13:16
欧拉角 发表于 2021-1-26 09:52
那样的话,施加的就不是拉力而是施加的应变。如果你需要的就是看应变带来的影响,那当然就可以。这种情况 ...

感谢
作者
Author:
sunjk    时间: 2022-1-18 19:31
isif=3结构优化后,改变晶格常数,应该用isif=4?还是isif=2??
作者
Author:
sunjk    时间: 2022-1-18 19:39
sunjk 发表于 2022-1-18 19:31
isif=3结构优化后,改变晶格常数,应该用isif=4?还是isif=2??

哦哦,就是isif=2合理
作者
Author:
飞翔的猪    时间: 2022-9-22 16:19
欧拉角 发表于 2021-1-24 17:45
优化的时候设置ISIF=3,可以得到压强接近零的晶格常数,然后设置ISIF=2,改变晶格常数(加压的话晶格常数往 ...

请问老师,您这里的“优化的时候设置ISIF=3”指的是优化bulk还是优化slab?
作者
Author:
人心向背    时间: 2023-2-4 17:36
欧拉角 发表于 2021-1-24 17:45
优化的时候设置ISIF=3,可以得到压强接近零的晶格常数,然后设置ISIF=2,改变晶格常数(加压的话晶格常数往 ...

请问利用改变晶格常数的方法模拟加压是否必须采用分数坐标,对笛卡尔坐标来说还算等比例压缩吗
作者
Author:
zzffzz33    时间: 2023-12-2 22:25
欧拉角 发表于 2021-1-24 17:45
优化的时候设置ISIF=3,可以得到压强接近零的晶格常数,然后设置ISIF=2,改变晶格常数(加压的话晶格常数往 ...

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3