24小时热门版块排行榜    

查看: 1015  |  回复: 6

大瑶Alice

新虫 (小有名气)

[求助] wien2k不收敛怎么办?

我用wien2k计算了53个W原子中有一个H原子的总能量,可是能量总是在两个数之间跳动,始终也不收敛。我运行的是run_lapw命令,是不是后面加上收敛条件会成功呢?可是我算其他情况的时候也都用了这个命令,能量也都收敛了呀。。。还是初始化的参数和算法上也要适当修改呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖


感谢参与,应助指数 +1
sunyang1988: 金币+1, 谢谢交流 2012-04-06 08:31:57
如果直接run_lapw,后面没有收敛参数,实际上等价于

run_lapw -i 40 -ec 0.0001

默认最大迭代步数 40,能量收敛标准0.0001 Ry
2楼2012-04-05 21:02:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大瑶Alice

新虫 (小有名气)

引用回帖:
2楼: Originally posted by fzx2008 at 2012-04-05 21:02:07:
如果直接run_lapw,后面没有收敛参数,实际上等价于

run_lapw -i 40 -ec 0.0001

默认最大迭代步数 40,能量收敛标准0.0001 Ry

哦。。长见识了!那能量在两个数之间来回变换不收敛是怎么回事呀?把能量收敛调得大一点会不会好呢?
3楼2012-04-05 21:24:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大瑶Alice

新虫 (小有名气)

还想弱问一下:run_lapw后,出现的circle1:ETEST和CTEST是什么意思啊?他们后面的数代表什么呢?
4楼2012-04-05 21:44:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大瑶Alice

新虫 (小有名气)

没人理我。。。还有一个问题就是grep:ENE case.scf后出现能量值得行数是不是就是叠代的次数?如果是这样的话,为什么行数有时会大于默认的40次叠代次数呢?
5楼2012-04-06 19:45:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大瑶Alice

新虫 (小有名气)

请问这个wien2k的能量收敛与k个数的选取有关系吗?因为我取的是k=1(shifted),可是其他的情况我也是这样取的,有的也是可以收敛的。不知道其中的原因是什么?
6楼2012-04-06 21:13:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liyanliyan

银虫 (小有名气)


感谢参与,应助指数 +1
贺仪: 金币-1, 应助指数-1, 请勿应助! 2012-04-18 16:40:25
问一下,若设置ec是0.0001Ry,则迭代次数是默认四十次吗?
勇敢面对!
7楼2012-04-18 15:43:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 大瑶Alice 的主题更新
信息提示
请填处理意见