24小时热门版块排行榜    

查看: 2262  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[基金申请] 什么时候开奖? +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 9/450 2026-08-23 10:59 by Tide man
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见