24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4621  |  回复: 23

yayanoba

新虫 (小有名气)


[交流] External pressure 很疑惑!!// vasp

对于在VASP的OUTCAR里面会有External pressure 和Pullay Stress,我倒是还能找到关于pullay stress的相关的描述,在手册里面也说了是一种未达到绝对收敛而产生的一种error,但是我想知道对于这个external pressure 到底说的是一种什么外界压力,我看到过说是可以在优化的过程中通过不同的算法使得external pressure 趋于0,但是还会发现就是我在静态计算的时候并没有设置PSTRESS这样的参数,按常理普通的用实验块材下的晶格经过静态计算得到的external pressure应该是为0的,但是为什么我看到在块才下也是有值的呢,而且也不是很小,有的有十几kb之大,甚至更大,想请教一下这种情况是正常的吗,或者说在块才下静态计算可以不需要去考虑这个参数只有在优化的时候才有用或者说是我设置的PSTRESS这样的参数才有用呢,希望有知道点的任何都能为我解惑一下!! 谢谢!!
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )
★ ★ ★
yayanoba(金币+1): 谢谢参与
franch: 金币+2, 谢谢回帖交流 2013-09-23 20:49:30
你relax的时候用的ISIF是多少?如果ISIF=3的话,PSTRESS=0的话,收敛之后external pressure不会太大。
15楼2013-09-23 11:21:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wawayk

铜虫 (初入文坛)


★ ★ ★ ★
yayanoba(金币+1): 谢谢参与
franch: 金币+3, 谢谢回帖交流 2013-09-23 20:49:39
我是这样理解的:不论是直接加压还是改变体积,都应在保持给定压力或体积的情况下优化晶胞参数,从而得到总能。直接加压一般就是加等静压,各个方向压力均等;改变体积相当于一种应变,会产生类似于加压的效果;个人认为只要优化能够完成,效果应该是差不多的。块材的时候,external pressure应该是接近于0的,我查看了我的块材计算,有的是零点几KB的,有的是负一点几KB的,基本上不会太大。如果加的是compressive strain,external pressure是正值,如果加的是tensile strain,则是负值。你可以查看你算过的体系看一下是不是这样。另外,如果楼主采纳,可以打发点喽!

[ 发自手机版 http://muchong.com/3g ]
17楼2013-09-23 14:24:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youdao1

铁虫 (知名作家)



yayanoba(金币+1): 谢谢参与
Good luck
7楼2013-09-22 23:00:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qingxiang827

铁虫 (知名作家)



yayanoba(金币+1): 谢谢参与
深深祝福
8楼2013-09-22 23:02:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gitme3388

铁虫 (知名作家)



yayanoba(金币+1): 谢谢参与
Best wishes
9楼2013-09-22 23:02:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oneday1234

铁虫 (知名作家)



yayanoba(金币+1): 谢谢参与
祝福楼主
10楼2013-09-22 23:04:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

god0729

铁虫 (知名作家)



yayanoba(金币+1): 谢谢参与
愿楼主心想事成
11楼2013-09-22 23:06:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jondan

木虫 (著名写手)



yayanoba(金币+1): 谢谢参与
顶起来!!!待高手解答啊!!!

[ 发自小木虫客户端 ]
14楼2013-09-23 09:10:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yayanoba

新虫 (小有名气)


引用回帖:
15楼: Originally posted by dxcharlary at 2013-09-23 11:21:14
你relax的时候用的ISIF是多少?如果ISIF=3的话,PSTRESS=0的话,收敛之后external pressure不会太大。

如果我没有优化,仅仅静态计算的话,也会有external pressure出现吗?
16楼2013-09-23 13:06:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
16楼: Originally posted by yayanoba at 2013-09-23 13:06:41
如果我没有优化,仅仅静态计算的话,也会有external pressure出现吗?...

取决于你的结构,不过因为DFT计算的晶格常数往往和实验有些差别。比如LDA会低估晶格常数,GGA会高估。所以你用实验数据的话,会出现误差也算正常。
19楼2013-09-24 02:13:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yayanoba

新虫 (小有名气)


引用回帖:
19楼: Originally posted by dxcharlary at 2013-09-24 02:13:37
取决于你的结构,不过因为DFT计算的晶格常数往往和实验有些差别。比如LDA会低估晶格常数,GGA会高估。所以你用实验数据的话,会出现误差也算正常。...

谢谢!
20楼2013-09-24 08:49:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pmm12

铁虫 (小有名气)



yayanoba(金币+1): 谢谢参与
引用回帖:
16楼: Originally posted by yayanoba at 2013-09-23 13:06:41
如果我没有优化,仅仅静态计算的话,也会有external pressure出现吗?...

LZ,我也遇到在静态计算后的external pressre比较大的现象,relax时是用的ISIF=3,OUTCAR中的external pressre比较小,为0.几kB,但是,在之后的静态计算OUTCAR中的external pressre比较大,都十几kb了,甚至更大些,请问LZ找到问题所在了吗?亟待回复,预先谢过了~
21楼2014-02-28 17:07:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shmonray

新虫 (小有名气)



yayanoba(金币+1): 谢谢参与
external pressure就是压力xx,yy,zz分量的总和,就在external pressure上面一行。一般来说如果PSTRESS=0为零的话,fully relaxed的force的xx, yy, zz都应该接近零,然后external pressure 也接近于零
22楼2014-03-04 12:09:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卓文尔

木虫 (职业作家)



yayanoba(金币+1): 谢谢参与
我也遇到相似的问题,我切Slab表面模型加真空层后用ISIF=2优化,External pressure值4~6 kB.但是我改变缩放系数与真空层(它们的乘积为定值,相当于手动优化晶格)体系所收到的应力External pressure值接近于0,此时能量也是最低的。但是我不知道有些什么原因造成了体系应力较大。
23楼2018-05-24 22:08:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzkkkkk

新虫 (小有名气)



yayanoba(金币+1): 谢谢参与
引用回帖:
23楼: Originally posted by 卓文尔 at 2018-05-24 22:08:04
我也遇到相似的问题,我切Slab表面模型加真空层后用ISIF=2优化,External pressure值4~6 kB.但是我改变缩放系数与真空层(它们的乘积为定值,相当于手动优化晶格)体系所收到的应力External pressure值接近于0,此时 ...

你好,你的问题解决了吗?我有类似的问题,用ISIF=3优化bulk的时候,external pressure只有零点几。但是用ISIF=2优化slab的时候,external pressure到了-8.不知道这样会不会有问题。
24楼2020-03-14 15:23:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2013-09-22 22:29   回复  
yayanoba(金币+1): 谢谢参与
祝福!!!
flyxu3楼
2013-09-22 22:40   回复  
yayanoba(金币+1): 谢谢参与
haixiawu4楼
2013-09-22 22:50   回复  
yayanoba(金币+1): 谢谢参与
readytogo5楼
2013-09-22 22:50   回复  
yayanoba(金币+1): 谢谢参与
2013-09-22 22:53   回复  
yayanoba(金币+1): 谢谢参与
2013-09-22 23:35   回复  
yayanoba(金币+1): 谢谢参与
2013-09-23 08:16   回复  
yayanoba(金币+1): 谢谢参与
2013-09-23 15:40   回复  
yayanoba(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 yayanoba 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见