24小时热门版块排行榜    

查看: 1691  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 0703化学调剂,求各位老师收留 +7 秋有木北 2026-03-14 7/350 2026-03-15 17:30 by 小物理化学
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 求调剂 +5 鹤遨予卿 2026-03-09 5/250 2026-03-14 02:44 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-09 6/300 2026-03-14 02:14 by JourneyLucky
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 085600求调剂 +3 a邵星池 2026-03-09 3/150 2026-03-14 01:32 by JourneyLucky
[考研] 一志愿华中农业大学071010,总分三百二,求调剂 +3 困困困困坤坤 2026-03-10 3/150 2026-03-14 00:35 by JourneyLucky
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 285求调剂 +4 ytter 2026-03-12 4/200 2026-03-13 14:48 by jxchenghu
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
[考研] 一志愿山东大学,总分327,英语二79,有论文,有竞赛,已过四六级 +3 木木目目1 2026-03-09 3/150 2026-03-09 19:52 by yuningshan
信息提示
请填处理意见