24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2237  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703调剂 +16 拾玖壹 2026-04-04 18/900 2026-04-07 12:49 by flydream1314
[考研] 325 调剂 +5 QQ小虾 2026-04-07 5/250 2026-04-07 11:24 by 诗与自由
[考研] 293求调剂 +3 勇远库爱314 2026-04-06 3/150 2026-04-07 11:15 by hugr
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 22408 调剂材料 +6 我叫ez 2026-04-06 7/350 2026-04-07 09:04 by MW_WTU
[考研] 302分求调剂 一志愿安徽大学085601 +12 zyx上岸! 2026-04-04 12/600 2026-04-07 02:09 by BruceLiu320
[考研] 305求调剂 +3 77Qi 2026-04-06 3/150 2026-04-06 15:36 by going home
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 315求调剂 +14 欣喜777 2026-04-04 15/750 2026-04-06 06:25 by houyaoxu
[考研] 085600,320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 计算机11408,286分求调剂 +7 木子念晞 2026-04-05 7/350 2026-04-05 19:02 by chy09050039
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 294求调剂 +6 Grey_Ey 2026-04-02 9/450 2026-04-04 22:07 by hemengdong
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-04-02 6/300 2026-04-03 10:19 by 蓝云思雨
[考研] 重庆大学材料与化工085600,初试370+,求求调剂建议 +8 shzhou_ 2026-04-01 9/450 2026-04-03 09:31 by 蓝云思雨
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8 双马尾痞老板2 2026-04-02 9/450 2026-04-02 14:45 by 5896
[考研] 生物学327,求调剂 +5 书上的梅子 2026-04-01 6/300 2026-04-02 06:47 by ilovexiaobin
[考研] 求调剂:一志愿:南京大学 专业:0705 总分320 ,本科985,四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
信息提示
请填处理意见