24小时热门版块排行榜    

查看: 2031  |  回复: 7

kule150

木虫 (小有名气)

[求助] 使用ISIF=7进行结构优化,总是出问题,是否需要另设一些特定参数呢 已有1人参与

现在需要对现在的体系进行优化,想固定一下体积做一次结构弛豫(ISIF=7),看结构的体积变化:
但试了多次总是有这样的问题:ZBRENT: fatal error in bracketing  please rerun with smaller EDIFF, or copy CONTCAR to POSCAR and continue
或者是: ZBRENT: fatal internal in brackting system-shutdown; contact gK immediat

我看网上帖子说:你所遇到的错误是因为计算误差或最低点附近能量曲线过于平坦之类的原因无法找到此路径上的最低点。如果确认输入文件没有问题的话,没有关系,从当前CONTCAR继续就行了(即停止寻找最低点,直接跳出循环外找下一个路径)如果反复出现这样的问题,不妨换IBRION=1试试。

这些解决办法我都试过,还是不行,我在想是不是在ISIF=7的情况下,需要另设一些参数呢,vasp新手,经验不足,望虫友们多指教
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哈士狼

新虫 (著名写手)

算什么?居然用7,你咋不用∞

发自小木虫Android客户端
2楼2017-10-17 12:18:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小熊尼尼

木虫 (正式写手)

ISIF的含义找一下说明书呀,固定体积的驰豫是ISIF=4

发自小木虫Android客户端
3楼2017-10-17 13:22:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kule150

木虫 (小有名气)

引用回帖:
2楼: Originally posted by 哈士狼 at 2017-10-17 12:18:34
算什么?居然用7,你咋不用∞

能设无穷大的话,我还真想试一下
4楼2017-10-17 13:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kule150

木虫 (小有名气)

引用回帖:
3楼: Originally posted by 小熊尼尼 at 2017-10-17 13:22:58
ISIF的含义找一下说明书呀,固定体积的驰豫是ISIF=4

嗯,说错了,是固定晶格结构和原子位置,只改变体积,我就想纯粹的看一下体积的改变,如果ISIF=3的话,不只体积变,晶格形状也变了。话说,一般ISIF=7一般什么场合下用,我算的是合金体系
5楼2017-10-17 13:48:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hs0218

金虫 (正式写手)

DDr.

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
kule150: 金币+10, ★★★很有帮助 2017-10-18 20:48:04
ISIF7找平衡体积没错
但是一般确实会出现各种问题吗,试着先用低kpoint粗略的跑一次然后再用CONTCAR接着跑
我用isif7跑过200多个的无序合金都没问题。
路漫漫其修远兮
6楼2017-10-17 14:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kule150

木虫 (小有名气)

引用回帖:
6楼: Originally posted by hs0218 at 2017-10-17 14:19:24
ISIF7找平衡体积没错
但是一般确实会出现各种问题吗,试着先用低kpoint粗略的跑一次然后再用CONTCAR接着跑
我用isif7跑过200多个的无序合金都没问题。

嗯,好嘞,这么一说,顿时心里有底了
7楼2017-10-17 16:22:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huihui521

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by kule150 at 2017-10-17 13:48:38
嗯,说错了,是固定晶格结构和原子位置,只改变体积,我就想纯粹的看一下体积的改变,如果ISIF=3的话,不只体积变,晶格形状也变了。话说,一般ISIF=7一般什么场合下用,我算的是合金体系...

所以请问楼主不想改变晶体形状进行优化,可以用ISIF=7吗,因为看多数人都不会用7,所以很疑惑
8楼2018-04-15 13:10:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kule150 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了,不知道这次命运如何 +3 雨打竹帘 2026-08-30 3/150 2026-08-30 09:43 by KNIGHT2011
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 08:33 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:17 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:35 by ZPa0EcMwuECS
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 06:00 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:51 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 面上意见出来了 +8 黄鸟于飞Chao 2026-08-29 16/800 2026-08-29 20:45 by 黄鸟于飞Chao
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
信息提示
请填处理意见