24小时热门版块排行榜    

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

pingguo802929

木虫 (初入文坛)


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

使用VASP对体系加压,可以使用关键字PSTRESS,经过我的实践,发现这种方法
在计算中容易出错,而且计算的时间较长。
   另一种方法是调节POSCAR文件中的放缩因子,从而改变体系的体积,大家有没有
发现这样作在OUTCAR文件中有一个‘external pressure’ 输出,可以把它记录下来,
然后将结果和使用PSTRESS的结果对比,我发现两者是一样的,而且这样做会大大缩短
计算的时间,出错的几率也很少,也不用费时费力的去拟合状态方程。
   大家有没有人知道这个 external pressure 是怎样产生出来的?这样做是不是可取?
如果这样做得当,这也许是对体系加压的好办法!
回复此楼
已阅   回复此楼   关注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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见