24小时热门版块排行榜    

查看: 1617  |  回复: 7

dxcharlary

专家顾问 (职业作家)

[交流] 【求助】isif=2是否允许改变元胞体积 已有4人参与

按照http://cms.mpi.univie.ac.at/vasp/vasp/ISIF_tag.html中的说法,isif=2的话不改变元胞体积,那么附加外界压强的话,不改变元胞形状也不改变元胞体积,只改变离子位置的话,结构优化还有意义吗?如果想要保持元胞形状,改变元胞体积和离子位置该怎么办?

[ Last edited by dxcharlary on 2010-7-29 at 21:59 ]
回复此楼

» 猜你喜欢

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

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

condensed

木虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2):谢谢指点 2010-07-30 08:28:40
isif = 3可以做到一切. isif = 2是有意义的.当isif = 3的时候,自由度很大,需要很长时间才能弛豫到平衡位置.先做isif = 2 有助于进一步的优化.
2楼2010-07-29 22:53:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
cenwanglai(金币+2):呵呵,除非 lz想找软件的bug.谢谢侯老师. 2010-07-30 09:47:36
引用回帖:
Originally posted by dxcharlary at 2010-07-29 21:58:22:
按照http://cms.mpi.univie.ac.at/vasp/vasp/ISIF_tag.html中的说法,isif=2的话不改变元胞体积,那么附加外界压强的话,不改变元胞形状也不改变元胞体积,只改变离子位置的话,结构优化还有意义吗?如 ...

施加了压强了,体积会变的,采用ISIF=3。

“那么附加外界压强的话,不改变元胞形状也不改变元胞体积,只改变离子位置的话",没有意义,原则上也是不能发生的,或者这样的过程是不物理的。不用去做,就想得到。
3楼2010-07-30 09:34:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

引用回帖:
Originally posted by condensed at 2010-07-29 22:53:27:
isif = 3可以做到一切. isif = 2是有意义的.当isif = 3的时候,自由度很大,需要很长时间才能弛豫到平衡位置.先做isif = 2 有助于进一步的优化.

一个面心立方的晶体,采用2个原子的晶胞来描述的话,需要固定晶胞的长宽比例为1:根号(2),如果采用isif=3的结果是最后优化的变形了,完全不是原来的东西了
4楼2010-07-30 14:45:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

引用回帖:
Originally posted by valenhou001 at 2010-07-30 09:34:10:


施加了压强了,体积会变的,采用ISIF=3。

“那么附加外界压强的话,不改变元胞形状也不改变元胞体积,只改变离子位置的话",没有意义,原则上也是不能发生的,或者这样的过程是不物理的。不用去做, ...

注意一下,我给的vasp的官方说明。isif=2时确实就是那么说的,不知道是不是写错了。可能是我问题没问清楚,我想说的是对于一个立方晶系而言,需要保持其高和底面的比例为1:1.414,既保持面心立方结构,而进行优化,isif应该取多少?

[ Last edited by dxcharlary on 2010-7-30 at 14:57 ]
5楼2010-07-30 14:54:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):感谢交流 2010-07-30 17:37:01
引用回帖:
Originally posted by dxcharlary at 2010-07-30 14:54:25:


注意一下,我给的vasp的官方说明。isif=2时确实就是那么说的,不知道是不是写错了。可能是我问题没问清楚,我想说的是对于一个立方晶系而言,需要保持其高和底面的比例为1:1.414,既保持面心立方结构,而进行 ...

因为VASP的优化时保持对称性的.当你优化立方晶系的时候,就算isif=3,晶体的形状(比例)是不会发生变化的.
6楼2010-07-30 17:14:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+2):谢谢解释 2010-07-30 18:35:41
引用回帖:
Originally posted by dxcharlary at 2010-07-30 14:54:25:


注意一下,我给的vasp的官方说明。isif=2时确实就是那么说的,不知道是不是写错了。可能是我问题没问清楚,我想说的是对于一个立方晶系而言,需要保持其高和底面的比例为1:1.414,既保持面心立方结构,而进行 ...

ISIF=2里面提到的是计算stress tensor,又不是施加压力。要保持高与地面的比例为1:1.414,手动自己优化就可以了,给一系列的a(a=b),并设置c=1/1.414 *a,只有优化原子位置,然后计算总能,这样得到了 E~ a的曲线,进行拟合找到最小值。
7楼2010-07-30 18:15:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+3):谢谢指导 2010-07-30 20:06:04
aylayl08(金币+1):谢谢,追加奖励 2010-07-31 19:25:23
引用回帖:
Originally posted by dxcharlary at 2010-07-29 21:58:22:
按照http://cms.mpi.univie.ac.at/vasp/vasp/ISIF_tag.html中的说法,isif=2的话不改变元胞体积,那么附加外界压强的话,不改变元胞形状也不改变元胞体积,只改变离子位置的话,结构优化还有意义吗?如 ...

你说的加外压是加PSTRESS吧,
是这样的,
你加上一个PSTRESS,
实际VASP是改通过改变晶胞体积实现的。
同时设定ISIF=2,
VASP会
1)先PSTRESS->改变体积,
2)ISIF=2,固定住这个形状,优化原子位置。

其实和PSTRESS+ISIF=3效果是一致的,

PSTRESS加外压,
如果体系发生结构相变,
HP,VP...的对应关系就不对了.

应该如valenhou所说的办法加压.
由天命,成人事!
8楼2010-07-30 18:38:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dxcharlary 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见