24小时热门版块排行榜    

查看: 1471  |  回复: 5

奋力拼搏

至尊木虫 (文坛精英)

[求助] 一定压力下如何计算弹性常数的问题,请各位高手同志们帮忙解决!谢谢你们啦!

由于初学VASP,不知道如何在有一定压强下计算晶体的弹性常数。更确切的说,我不明白需要在INCAR文件中设置什么参数才能做到体积(一定压强所对应的体积VO)不变的应变。也就是说,我在该压强对应的体积VO上施加形变,让形变后的体积仍然等于VO。请各位高手帮忙解决,谢谢!
回复此楼
天道酬勤
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
liliangfang: 金币+2, 应助指数+1, 谢谢指教 2012-10-13 15:16:29
优化中保持晶胞体积不变,试试ISIF=4
2楼2012-10-12 21:17:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奋力拼搏

至尊木虫 (文坛精英)

引用回帖:
2楼: Originally posted by fzx2008 at 2012-10-12 21:17:06
优化中保持晶胞体积不变,试试ISIF=4

感谢您的回复。我试过了,用ISIF=4的话好像只能保证每一个形变后的体积不变,而不能保证形变后的体积仍然等于没发生形变时的体Vo,冒昧的问一下,你们在算压力下的弹性常数是怎么设置参数的啊?
天道酬勤
3楼2012-10-12 21:27:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liliangfang: 金币+2, 应助指数+1, 谢谢提示 2012-10-13 15:16:52
奋力拼搏: 金币+10, ★★★很有帮助 2012-10-14 10:21:22
那这样,注意一下POSCAR的第二行,如果是负值,则表示体积!(正值表示universal scaling factor)

» 本帖已获得的红花(最新10朵)

4楼2012-10-12 23:14:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奋力拼搏

至尊木虫 (文坛精英)

送鲜花一朵
引用回帖:
4楼: Originally posted by fzx2008 at 2012-10-12 23:14:36
那这样,注意一下POSCAR的第二行,如果是负值,则表示体积!(正值表示universal scaling factor)

嗯,好的。谢谢您的帮助,我试试看
天道酬勤
5楼2012-10-14 10:21:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
franch: 金币+5, 鼓励交流,,呵呵 2012-10-15 12:29:34
我最近在做LZ说的这件事。

可以去看看我最近发的帖子 http://muchong.com/bbs/viewthread.php?tid=5041009&fpage=1

我的方法是这样的:
固定a/c值,保持体积不变。原来的体积是V0,现在的体积仍然是V0,只不过a/c值变化了。

由于我的体系是体心立方,就只有一个常数X,假如a/c=0.9,则通过计算
具体设置体现在POSCAR里面晶格常数那部分:
1.0
X*sqrt(0.9) 0.0 0.0
0.0              X    0.0
0.0              0.0  X/sqrt(0.9)

不知道LZ算的什么体系,根据实际情况进行计算吧。
6楼2012-10-14 20:36:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 奋力拼搏 的主题更新
信息提示
请填处理意见