24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

604gq

木虫 (职业作家)

学历太低不辨东西南北

应助: 0 (幼儿园)
贵宾: 0.9
金币: 2717.9
散金: 78
红花: 4
帖子: 3535
在线: 159.8小时
虫号: 190935
注册: 2006-02-20
专业: 胶体与界面化学

[交流] 自洽场（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 ]

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

坚决贯彻“灌水有理，一灌到底，乱灌不上税，灌死不负责”的大无畏精神。

1楼 2008-05-21 10:43:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lsl6585

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 443
帖子: 24
在线: 10.7小时
虫号: 565453
注册: 2008-05-29
专业: 凝聚态物性 II ：电子结构

引用回帖:

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泛函产生的初始猜测波函数不好所以你需要用自己产生的波函数做为新的计算的初始猜测波函数