24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2306  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +5 yexuqing 2026-04-19 5/250 2026-04-22 10:01 by easeheart
[论文投稿] 急需审稿人!!! +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
请填处理意见