24小时热门版块排行榜    

查看: 1809  |  回复: 5
当前主题已经存档。

fish2282

铁虫 (初入文坛)

[交流] 【求助】请教这是怎么一回事Density matrix is not changing but DIIS error

Density matrix is not changing but DIIS error= 1.20D-06 CofLast= 4.45D-01.
The SCF is confused.
这是怎么一回事啊?怎么样解决类似的问题啊?请高手指教一下。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqx

木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
erylingjet(金币+2,VIP+0):感谢交流! 9-17 09:31
DIIS(direct inverse iteration subspace method)是一种用来改进收敛的手段,使用过程中需要求解一个线性方程组。此处的问题应该是求解线性方程组出现错误(可能是处于数值计算稳定性方面的原因)。你可尝试SCF=NoDiis关键字,如果没有其他问题出现就ok了。
PS: NoDiis意味着在迭代的过程中不采用DIIS方法,带来的后果是计算时间会变长,这一点对于较大的任务会更明显。
2楼2009-09-17 09:16:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yflchx

金虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+2,VIP+0): 9-17 10:58
用SCF=QC可以避免“The SCF is confused”错误,因为这种错误仅仅出现在DIIS计算过程中。考虑到QC方法比较耗费计算时间,可以采用SCF=(XQC,MAXCONVENTIONALCYCLES=N)。N要取一个合理的数值,在DIIS方法出现错误之前跳转为QC方法计算。用#P打印出来DIIS迭代的具体信息就可以确定合理的N值。
3楼2009-09-17 09:47:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fish2282

铁虫 (初入文坛)

引用回帖:
Originally posted by yflchx at 2009-9-17 09:47:
用SCF=QC可以避免“The SCF is confused”错误,因为这种错误仅仅出现在DIIS计算过程中。考虑到QC方法比较耗费计算时间,可以采用SCF=(XQC,MAXCONVENTIONALCYCLES=N)。N要取一个合理的数值,在DIIS方法出现错误之 ...

我用过这个关键词了,不太管用。
4楼2009-09-22 11:07:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fish2282

铁虫 (初入文坛)

谢谢各位!
5楼2009-09-22 11:08:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yflchx

金虫 (正式写手)

★ ★
yjcmwgk(金币+2,VIP+0):斑竹评分,欢迎常来 9-22 18:56
scf=(xqc,maxconventionalcycles=N) 其中的N你取的多少?
一定要先找到具体SCF迭代几个CYCLES后出错,然后才能确定N值。
比如,DIIS到第35步出错,那么N可以取为34。
祝你好运。
引用回帖:
Originally posted by fish2282 at 2009-9-22 11:07:




我用过这个关键词了,不太管用。

6楼2009-09-22 15:12:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fish2282 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 CyTWnZiM0cS6 2026-08-02 5/250 2026-08-02 20:52 by 3VGOm6paSgqv
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 CyTWnZiM0cS6 2026-08-02 3/150 2026-08-02 20:44 by 3VGOm6paSgqv
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 CyTWnZiM0cS6 2026-08-02 5/250 2026-08-02 20:43 by 3VGOm6paSgqv
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 AS6fX1tr1on4 2026-08-02 9/450 2026-08-02 20:32 by 3VGOm6paSgqv
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 5/250 2026-08-02 20:31 by 3VGOm6paSgqv
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 20:24 by 3VGOm6paSgqv
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 AS6fX1tr1on4 2026-08-02 6/300 2026-08-02 20:21 by 3VGOm6paSgqv
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 12/600 2026-08-02 20:21 by drbart
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 AS6fX1tr1on4 2026-08-02 6/300 2026-08-02 20:14 by 3VGOm6paSgqv
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 AS6fX1tr1on4 2026-08-02 14/700 2026-08-02 20:11 by 3VGOm6paSgqv
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 7/350 2026-08-02 20:11 by 3VGOm6paSgqv
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 AS6fX1tr1on4 2026-08-02 13/650 2026-08-02 20:04 by 3VGOm6paSgqv
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 13/650 2026-08-02 20:04 by 3VGOm6paSgqv
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 19:54 by 3VGOm6paSgqv
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 cBGhfJFuNMca 2026-08-02 21/1050 2026-08-02 19:52 by 3VGOm6paSgqv
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 9/450 2026-08-02 18:27 by ca3DVAx4OB19
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 18:24 by ca3DVAx4OB19
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
请填处理意见