24小时热门版块排行榜    

查看: 1688  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 313分生物学求调剂 +6 Yyt杨1 2026-03-09 8/400 2026-03-14 03:00 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 7/350 2026-03-13 17:28 by xujiaoszu
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 293求调剂,一志愿陕师大生物学 +3 ??????.?.??? 2026-03-09 3/150 2026-03-11 10:02 by 学员8dgXkO
[考研] 一志愿山东大学,总分327,英语二79,有论文,有竞赛,已过四六级 +3 木木目目1 2026-03-09 3/150 2026-03-09 19:52 by yuningshan
信息提示
请填处理意见