| 查看: 2253 | 回复: 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 ] |
回复此楼» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有196人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有1人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招材料相关专业调剂生
已经有10人回复
-
招收0703开头的学硕调剂生 线上复试
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
yaowenzhi
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 65.8
- 帖子: 206
- 在线: 25.4小时
- 虫号: 458345
- 注册: 2007-11-13
- 性别: MM
- 专业: 理论和计算化学
2楼2008-05-21 15:18:51
wxjbuilder
木虫 (知名作家)
- 应助: 1 (幼儿园)
- 贵宾: 0.05
- 金币: 2986.4
- 散金: 1761
- 红花: 2
- 沙发: 1
- 帖子: 5392
- 在线: 273.2小时
- 虫号: 442356
- 注册: 2007-10-27
- 性别: GG
- 专业: 理论和计算化学
3楼2008-05-21 19:12:50
4楼2008-05-23 07:47:25
红枫叶688
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4790.5
- 散金: 3
- 红花: 1
- 帖子: 227
- 在线: 59小时
- 虫号: 469001
- 注册: 2007-11-28
- 专业: 理论和计算化学
5楼2008-06-03 21:47:04
xwj007
铁杆木虫 (职业作家)
- 应助: 0 (幼儿园)
- 金币: 7454.1
- 帖子: 3957
- 在线: 1373.3小时
- 虫号: 399979
- 注册: 2007-06-12
- 性别: GG
- 专业: 高分子组装与超分子结构
6楼2008-06-03 23:05:29
python9527
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1201
- 红花: 1
- 帖子: 141
- 在线: 23.7小时
- 虫号: 391033
- 注册: 2007-06-03
- 专业: 有机分子功能材料化学
7楼2008-06-04 09:17:18
8楼2008-06-04 20:19:54
|
用略微改变后的分子几何结构做单点能计算 注意保存chk文件 假定这个单点能计算能够正常收敛 那么我们可以把这个计算中最后得到的波函数做为新的计算的初始猜测波函数 这个波函数就保存在chk文件里面 因此新的计算可以采用原来的几何构型 但是需要在输入文件里面指定guess=read （注意，不要习惯性的加上geom=check） 同时也要把前面用略微改变后的分子几何结构做单点能计算产生的chk文件包含在你的输入文件里面 比如说%chk=*.chk 如果你不加guess=read这个选项 那么程序会用默认的Harris泛函产生的波函数做为SCF的初始猜测波函数 这样做导致不收敛是因为在有些情况下是因为用Harris泛函产生的初始猜测波函数不好 所以你需要用自己产生的波函数做为新的计算的初始猜测波函数 |
9楼2008-06-05 00:48:18
10楼2008-06-06 08:45:02
5