24小时热门版块排行榜    

查看: 1349  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 8/400 2026-11-02 03:24 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-02 03:18 by xzcv0Tehd8lQ
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:05 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-02 02:34 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 8/400 2026-11-02 00:37 by xzcv0Tehd8lQ
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:24 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:05 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 21:52 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 20:03 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 16:58 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
信息提示
请填处理意见