应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp自洽计算过程中不收敛的解决办法交流
293/3返回列表上一页123
查看: 8503  |  回复: 28
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

obaica
引用回帖:
18楼: Originally posted by Feynmanux at 2016-09-04 13:21:31
帮顶，前几天恰好出现了这么个情况。

做优化到 只要一个离子步就达到精度结束。isif =2 7 都是如此。

然而到了自洽，却怎么也不收敛，就在一个能量值附近起伏。

IALGO 48 38 都试过了，不收敛；改动 sigma ...

棒

发自小木虫Android客户端
21楼2016-09-04 13:28:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk
引用回帖:
20楼: Originally posted by Feynmanux at 2016-09-04 13:23:51
把没收敛（比如默认60步结束后）写出来的波函数，读一遍吗？...

是
22楼2016-09-04 13:46:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
obaica: 金币+5 2016-09-04 19:21:25
如果难收敛首要得考虑模型是否合理，其次看看计算参数选择是否合适、是否经过严格测试，最后实在不行再从算法下手，自洽中的电子步算法换几种之后都不行的话，就得考虑是不是结构本身的问题了。

检查一下之前离子弛豫过程中每一步的电子步是否都收敛了，如果离子弛豫过程中有电子步是没收敛的，之后自洽可能没法正常收敛。

不排除结构不合理的情况，如果结构不合理，离子弛豫前几步可能电子步不是正常收敛的，但是对于vasp，电子步不收敛计算并不会停止（类似于Gaussian的IOp(5/13=1)），这样的结果可能会有些问题。

我遇到过极少情况确实会遇到特定的泛函优化得到的晶格然后切面计算确实收敛不了的，不过这可能是很特殊的情况了。
赞一下(1人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
23楼2016-09-04 14:45:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happinesslyl

新虫 (小有名气)
引用回帖:
20楼: Originally posted by Feynmanux at 2016-09-04 13:23:51
把没收敛（比如默认60步结束后）写出来的波函数，读一遍吗？...

是啊，你可以试试

发自小木虫Android客户端
回复此楼

» 本帖已获得的红花（最新10朵）

Feynmanux
24楼2016-09-05 08:41:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Feynmanux
送红花一朵
引用回帖:
24楼: Originally posted by happinesslyl at 2016-09-05 08:41:25
是啊，你可以试试
...

3x
25楼2016-09-05 22:39:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

金天赐

木虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
引用回帖:
10楼: Originally posted by happinesslyl at 2016-08-25 08:48:02
设置istart=1, icharg=11,读取波函数和电荷密度再算一遍，一般几步就能收敛

前辈，我是结构优化时不收敛，设置istart=1, icharg=1,读取电荷密度(波函数是空白，我就没有复制到新的计算文件里)，把contcar改成poscar再算一遍，可不可以？计算结果可信吗？有时候要反复这样算2,3次才能收敛，结果可信吗？另外问一下，icharg一定要等于11吗，我设的等于1,也能读取电荷密度呀。

发自小木虫Android客户端
赞一下 回复此楼
26楼2016-12-07 11:00:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

不久又一年

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
23楼: Originally posted by 卡开发发 at 2016-09-04 14:45:36
如果难收敛首要得考虑模型是否合理，其次看看计算参数选择是否合适、是否经过严格测试，最后实在不行再从算法下手，自洽中的电子步算法换几种之后都不行的话，就得考虑是不是结构本身的问题了。

检查一下之前离子 ...

您好，我想请教一下，在结构优化过程中大部分电子步没有收敛应该如何解决？

发自小木虫Android客户端
赞一下 回复此楼
27楼2017-08-10 15:53:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

追枫少年

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
23楼: Originally posted by 卡开发发 at 2016-09-04 14:45:36
如果难收敛首要得考虑模型是否合理，其次看看计算参数选择是否合适、是否经过严格测试，最后实在不行再从算法下手，自洽中的电子步算法换几种之后都不行的话，就得考虑是不是结构本身的问题了。

检查一下之前离子 ...

感谢分享见解！
我在在做切面计算时，静态计算很难收敛，有的人说界面的计算就是比较困难，可否分享随着这类计算的见解
赞一下 回复此楼
28楼2017-12-22 16:26:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

莹莹的狗

铁虫 (小有名气)
引用回帖:
18楼: Originally posted by Feynmanux at 2016-09-04 13:21:31
帮顶，前几天恰好出现了这么个情况。

做优化到 只要一个离子步就达到精度结束。isif =2 7 都是如此。

然而到了自洽，却怎么也不收敛，就在一个能量值附近起伏。

IALGO 48 38 都试过了，不收敛；改动 sigma ...

缩小k点可以吗。。
回复此楼
29楼2019-06-27 22:34:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 obaica 的主题更新
293/3返回列表上一页123
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭