24小时热门版块排行榜    

查看: 2257  |  回复: 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.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 03:31 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[教师之家] 跳槽后在研项目怎么办? +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 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
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见