24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2254  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 290调剂生物0860 +35 哇哈哈,。 2026-04-11 41/2050 2026-04-14 23:47 by Xurambo2014
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[考研] 药学求调剂 +11 喽哈加油 2026-04-14 13/650 2026-04-14 21:14 by qingfeng258
[考研] 302分求调剂 +10 凡语祈愿 2026-04-08 11/550 2026-04-14 16:50 by jiangguiquan11
[考研] 一志愿085502,267分求调剂 +19 再忙也要吃饭啊 2026-04-08 20/1000 2026-04-14 16:03 by zs92450
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 297求调剂 +23 ORCHID1 2026-04-10 26/1300 2026-04-14 13:52 by 陈皮皮
[考研] 300分求调剂 (085501机械专硕,本科扬大) +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu~
[考研] 314求调剂 +24 wakeluofu 2026-04-09 25/1250 2026-04-13 08:58 by lhj2009
[考研] 一志愿华中农微生物,288分,三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 求调剂,262机械专硕 +8 嗯yyl 2026-04-08 8/400 2026-04-12 02:31 by 秋豆菜芽
[考研] 调剂 +6 青灯不负 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 08:57 by zhq0425
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 296求调剂 +6 汪!?! 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 085404,334分,求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
信息提示
请填处理意见