24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1726  |  回复: 9
【奖励】 本帖被评价2次,作者pingguo802929增加金币 1.25
当前主题已经存档。

pingguo802929

木虫 (初入文坛)


[资源] 【讨论】使用VASP在压力下对体系的计算

使用VASP对体系加压,可以使用关键字PSTRESS,经过我的实践,发现这种方法
在计算中容易出错,而且计算的时间较长。
   另一种方法是调节POSCAR文件中的放缩因子,从而改变体系的体积,大家有没有
发现这样作在OUTCAR文件中有一个‘external pressure’ 输出,可以把它记录下来,
然后将结果和使用PSTRESS的结果对比,我发现两者是一样的,而且这样做会大大缩短
计算的时间,出错的几率也很少,也不用费时费力的去拟合状态方程。
   大家有没有人知道这个 external pressure 是怎样产生出来的?这样做是不是可取?
如果这样做得当,这也许是对体系加压的好办法!
回复此楼

» 收录本帖的淘帖专辑推荐

VASP

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf(金币+1):谢谢
没有什么用,状态方程才是正解
2楼2009-04-01 22:47:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shxincui518

铜虫 (正式写手)



wuchenwf(金币+1):谢谢
POSCAR加放缩因子一般来说对立方体系较适用,但对于非立方,例如单斜等就不行了吧!
3楼2009-04-02 08:05:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pingguo802929

木虫 (初入文坛)


楼上的专家,你好

楼上的专家: 您的 意思是说,只有用状态方程拟合出来的结果才可用?
用关键字PSTRESS加压的结果,可以相信吗?
还有我想弄明白的是OUTCAR中输出的压力是怎样得来的,是不是在VASP中也
有一套状态方程?
4楼2009-04-02 08:06:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
wuchenwf(金币+2):谢谢
这种方法不适合处理立方体系以外的体系。
因为在做结构优化时,不仅体积会变,而且形状也很容易发生变化。
如果形状发生变化了,你的方法就不对了。你得到的状态将不是最稳定的状态。
其实在给定压力后计算,虽然最后得到的压力不是绝对准确,但是是可以接受的,因为它不会变化太大。
而且vasp做的负压很不准,有些体系得到的状态方程跟本就不可信。
这一点pw等软件要不vasp好。
5楼2009-04-02 15:17:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shxincui518

铜虫 (正式写手)



spur(金币+1,VIP+0): 3Q!!谢谢参与! 4-2 23:03
vasp做负压的确糟糕的很,我的计算结果不好,就是考虑了负压。我要回去再重新计算下看看。
6楼2009-04-02 21:28:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cbFeng200881

金虫 (著名写手)


★ 一星级,一般


zxzj05(金币+1,VIP+0):3Q! 期待更多讨论! 4-10 18:29
VASP加压,搞什么呢?
7楼2009-04-03 01:29:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

化学小工

木虫 (正式写手)


★★★★★ 五星级,优秀推荐


zxzj05(金币+1,VIP+0):3Q! 期待更多讨论! 4-10 18:29
记得去年清华的一片文章就用的lz的方法。当然个人不做这方面的 研究关注一下
8楼2009-04-03 09:13:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

perrin

铜虫 (初入文坛)



zxzj05(金币+1,VIP+0):3Q! 期待更多讨论! 4-10 18:29
算一堆体积,拟合状态方程。不用多,正压5-6个,负压3个左右,状态方程拟合的很好了。
我一般用pwscf的ev.x拟合
9楼2009-04-10 02:01:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xujc1983

木虫 (著名写手)


★ ★ ★
zdhlover(金币+3,VIP+0):谢谢回帖,欢迎常来 4-17 15:59
引用回帖:
Originally posted by pingguo802929 at 2009-4-1 20:48:
使用VASP对体系加压,可以使用关键字PSTRESS,经过我的实践,发现这种方法
在计算中容易出错,而且计算的时间较长。
   另一种方法是调节POSCAR文件中的放缩因子,从而改变体系的体积,大家有没有
发现这样作在 ...

不清楚楼主要研究的问题,但是加压的话,第一种方法算完一组后,可以直接算体模量,第二种方法虽然不能直接算体模量,但是可以算出各种弹性系数,模量也就没有问题了,而且感觉第二种方法物理景象更清楚一些
推荐一篇论文,赵纪军老师关于二阶、三阶弹性系数计算方面的

http://scitation.aip.org/getabs/ ... =cvips&gifs=yes

[ Last edited by xujc1983 on 2009-4-17 at 14:25 ]
10楼2009-04-17 14:21:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pingguo802929 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初试分332,一志愿报考西北工业大学, +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[考研] 0703调剂 +18 拾玖壹 2026-04-04 20/1000 2026-04-09 20:53 by zhouxiaoyu
[考研] 调剂 化学 307 +11 73372112 2026-04-09 12/600 2026-04-09 18:52 by vgtyfty
[考研] 085501机械专硕 302分 不挑专业求调剂 +5 汪某. 2026-04-09 5/250 2026-04-09 15:38 by 蒋皓禹
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +15 努力奋斗112 2026-04-06 17/850 2026-04-08 22:27 by 凯凯要变帅
[考研] 专硕0854初试考材科基,求调剂 +7 3220548044 2026-04-06 10/500 2026-04-08 21:59 by hypershenger
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +10 努力奋斗112 2026-04-07 10/500 2026-04-08 15:01 by screening
[考研] 336材料与化工085600求调剂 +19 水星记infp 2026-04-05 22/1100 2026-04-07 21:11 by yongzhesheng
[考研] 277求调剂 数一104分 +9 瓶子PZ 2026-04-05 14/700 2026-04-07 17:52 by 蓝云思雨
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
[考研] 专硕304找调剂,一线城市最好 +3 李lsl李 2026-04-05 3/150 2026-04-06 12:16 by ffffjjjj
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 315求调剂 +5 &123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 能动调剂326专硕 +4 wan112233 2026-04-04 4/200 2026-04-04 22:47 by yu221
[考研] 321求调剂 +6 认真求上学 2026-04-03 6/300 2026-04-04 19:51 by dongzh2009
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
信息提示
请填处理意见