24小时热门版块排行榜    

查看: 1939  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

菲顶顶

铜虫 (正式写手)

[求助] 体相和表面做应变 已有1人参与

想请问大家,用vasp计算分别计算体相和表面情况的单轴和双轴应变的时候的方法是一样的吗?
我的想法是直接改变晶格常数不就行了,无论体相还是表面,如果是单轴的就改变一个方向的晶格常数,双轴的同时改变两个方向的晶格常数,然后优化。ISIF是怎么取的?体相取3?表面取2?
可是我在查一些资料的时候看到大家说的什么拉伸哪个方向就固定该方向不优化之类之类的。表示不懂,有做过应力应变的大神们可否来解答一下。
回复此楼

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

应变

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

redsnowolf

银虫 (小有名气)

vasp算应变的时候,变形是通过修改POSCAR自己加的,并且应变量要逐步增大,算大变形的时候不能把应变量一步加到位。结构优化前要修改constr_cell_relax.F文件,将施加应变的方向固定不优化,然后ISIF该怎么取就怎么取。例如只在x方向上有变形,可以在constr_cell_relax.F里加上FCELL(1,1)=0,把x方向上的力置零,重新编译,再优化时x方向的晶格常数就不变了。其它方向的、单双轴变形的情况类似,具体的在constr_cell_relax.F文件里有说明。
8楼2015-02-04 14:15:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

stractor

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
菲顶顶: 金币+10, ★★★很有帮助 2015-01-31 22:24:33
liliangfang: 金币+2, 谢谢交流 2015-02-01 09:49:00
“ISIF是怎么取的?体相取3?表面取2?”
不能直接取3和2,那样弛豫后又回到平衡状态,应变消失。
听说有人改写了VASP弛豫的代码,可以实现你的要求。
还有一种方式就是EOS拟合,比如,我在xy平面施加双轴应变,让z方向逐渐改变一点晶格常数(一些列dz),计算在该应变下,不同z的改变对于的体系总能,然后通过E~dz找到能量最低的.
2楼2015-01-31 21:57:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)

引用回帖:
2楼: Originally posted by stractor at 2015-01-31 21:57:39
“ISIF是怎么取的?体相取3?表面取2?”
不能直接取3和2,那样弛豫后又回到平衡状态,应变消失。
听说有人改写了VASP弛豫的代码,可以实现你的要求。
还有一种方式就是EOS拟合,比如,我在xy平面施加双轴应变, ...

这样确实可以的,我先试试,那ISIF到底该取什么呢,你的意思就是不优化直接计算自洽能量是吧。
3楼2015-01-31 22:24:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
菲顶顶: 金币+10, ★★★很有帮助 2015-02-01 14:08:19
引用回帖:
3楼: Originally posted by 菲顶顶 at 2015-01-31 22:24:26
这样确实可以的,我先试试,那ISIF到底该取什么呢,你的意思就是不优化直接计算自洽能量是吧。...

不用优化了,直接电子自洽计算。然后通过总能来确定。

[ 发自手机版 http://muchong.com/3g ]
4楼2015-02-01 11:21:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见