24小时热门版块排行榜    

查看: 4883  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 5+3 想要一壶桃花水 2026-03-25 7/350 2026-03-26 03:29 by sjeidudkso
[考研] 329求调剂 +3 钮恩雪 2026-03-25 3/150 2026-03-25 14:43 by 糖加冰
[考研] 085602 289分求调剂 +7 WWW西西弗斯 2026-03-24 7/350 2026-03-25 14:28 by 3Strings
[考研] 材料与化工考研调剂 +9 孅華 2026-03-22 9/450 2026-03-25 13:09 by cmz0325
[考研] 284求调剂 +15 Zhao anqi 2026-03-22 15/750 2026-03-25 12:51 by wht0531
[考研] 285求调剂 +3 AZMK 2026-03-24 3/150 2026-03-25 12:23 by userper
[考研] 299求调剂 +7 shxchem 2026-03-20 9/450 2026-03-25 10:41 by lbsjt
[考研] 0856材料专硕353求调剂 +6 NIFFFfff 2026-03-20 6/300 2026-03-24 21:02 by hello七七
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +4 曼111 2026-03-24 5/250 2026-03-24 15:44 by 星空星月
[考研] 291求调剂 +8 hhhhxn.. 2026-03-23 8/400 2026-03-23 23:15 by peike
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 一志愿070300浙大化学358分,求调剂! +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
信息提示
请填处理意见