| 查看: 2243 | 回复: 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 ] |
回复此楼» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有231人回复
-
汕头大学化学招调剂，相关内容如下。
已经有4人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
-
分子动力学模拟合作需求
已经有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