24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2284  |  回复: 15
当前主题已经存档。

604gq

木虫 (职业作家)

学历太低不辨东西南北

[交流] 自洽场(scf)计算收敛问题及解决

自洽场(scf)计算收敛问题及解决



SCF(自洽场)不收敛则一般是L502错误

两种原因造成的不收敛:数学上不收敛
                                  化学上不合理


1,改变初猜,使用GUESS=CORE,HUCKEL,MIX关键词;
2,对于开壳层体系,先对于其相应的闭壳层离子计算到收敛,然后使用GUESS=READ关键词从chk文件中读取分子轨道初猜;
3,先用小基组得到分子轨道,然后同样使用GUESS=READ关键词从chk文件中读取;
4,使用能级提升法(Vshift),将HOMO和LUMO的能量间隔加大;
5,将SCF=SLEAZY收敛限放松,或者用SCFCON关键词改变收敛限;
6,对于使用了弥散函数的情形,需要提高积分精度,加关键词SCF=NoVarAcc;
7,略微改变分子几何结构,得到收敛的波函数后,再使用原来的几何构型读分子轨道计算;
8,考虑使用不同级别的计算方法和基组(小一些的);
9,使用强制收敛方法SCF=QC(经常用到)或者SCF=DM。
10,用iop(5/13=1)忽略不收敛,继续往下做。


[ Last edited by 604gq on 2008-5-21 at 10:45 ]
回复此楼
坚决贯彻“灌水有理,一灌到底,乱灌不上税,灌死不负责”的大无畏精神。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaowenzhi

铜虫 (小有名气)

顶,收藏了,谢谢
2楼2008-05-21 15:18:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxjbuilder

木虫 (知名作家)

不错啊,好好研究一下
悠悠天地,奈何独立苍茫?
3楼2008-05-21 19:12:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxylf

金虫 (正式写手)

顶!高,收藏了
4楼2008-05-23 07:47:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

红枫叶688

木虫 (小有名气)

好极了,Thanks
5楼2008-06-03 21:47:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwj007

铁杆木虫 (职业作家)

十分的好!谢谢分享!
6楼2008-06-03 23:05:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

python9527

金虫 (小有名气)

7,略微改变分子几何结构,得到收敛的波函数后,再使用原来的几何构型读分子轨道计算;



这一条怎么做呀?麻烦lz指点下!!!
7楼2008-06-04 09:17:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huahua1981

铜虫 (小有名气)

就是 7有问题
8楼2008-06-04 20:19:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsl6585

银虫 (初入文坛)

引用回帖:
Originally posted by python9527 at 2008-6-4 09:17:
7,略微改变分子几何结构,得到收敛的波函数后,再使用原来的几何构型读分子轨道计算;



这一条怎么做呀?麻烦lz指点下!!!

用略微改变后的分子几何结构做单点能计算 注意保存chk文件 假定这个单点能计算能够正常收敛 那么我们可以把这个计算中最后得到的波函数做为新的计算的初始猜测波函数 这个波函数就保存在chk文件里面 因此新的计算可以采用原来的几何构型 但是需要在输入文件里面指定guess=read (注意,不要习惯性的加上geom=check) 同时也要把前面用略微改变后的分子几何结构做单点能计算产生的chk文件包含在你的输入文件里面 比如说%chk=*.chk 如果你不加guess=read这个选项 那么程序会用默认的Harris泛函产生的波函数做为SCF的初始猜测波函数 这样做导致不收敛是因为在有些情况下是因为用Harris泛函产生的初始猜测波函数不好 所以你需要用自己产生的波函数做为新的计算的初始猜测波函数
9楼2008-06-05 00:48:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huahua1981

铜虫 (小有名气)

引用回帖:
Originally posted by lsl6585 at 2008-6-5 00:48:


用略微改变后的分子几何结构做单点能计算 注意保存chk文件 假定这个单点能计算能够正常收敛 那么我们可以把这个计算中最后得到的波函数做为新的计算的初始猜测波函数 这个波函数就保存在chk文件里面 因此新的 ...

哥们你怎么越说越复杂啊了啊
10楼2008-06-06 08:45:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 604gq 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 3/150 2026-04-19 21:44 by 淡雅人生27
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 297,工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 291求调剂 +12 关忆北. 2026-04-14 13/650 2026-04-19 16:50 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 6/300 2026-04-19 13:00 by Aaron_zyn
[考研] 327求调剂 +27 Xxjc1107. 2026-04-13 30/1500 2026-04-19 08:22 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 一志愿沪9,326求生物学调剂 +12 刘墨墨 2026-04-13 12/600 2026-04-18 23:31 by 路病情
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考博] 申博/考博 +3 啃面包的小书虫 2026-04-17 4/200 2026-04-17 23:54 by 阳阳阳^_^
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
信息提示
请填处理意见