| 查看: 1234 | 回复: 7 | ||||
| 当前主题已经存档。 | ||||
[交流]
【求助】关于提高精度就不收敛的问题
|
||||
|
算一个优化结构时,先用hf/3-21g,正常的,再用b3lyp/3-21g,就有2个不收敛 Maximum Force 0.000002 0.000450 YES RMS Force 0.000000 0.000300 YES Maximum Displacement 0.011336 0.001800 NO RMS Displacement 0.001408 0.001200 NO 这是什么原因?我应该怎么做呢?谢谢高手指点我 |
回复此楼» 收录本帖的淘帖专辑推荐
 Gaussian |
» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有266人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有2人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
-
微纳hub，微纳加工的宝库
已经有0人回复
2楼2008-03-03 09:39:13
好好生活
木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 2883.4
- 散金: 67
- 红花: 1
- 帖子: 1204
- 在线: 351.2小时
- 虫号: 316003
- 注册: 2007-03-03
- 专业: 电化学
3楼2008-03-03 10:22:45
★ ★ ★
lei0736(金币+2,VIP+0):多谢详细解答
herosmzd(金币+1,VIP+0):多谢指教,能帮我解答下对结果的影响吗?
lei0736(金币+2,VIP+0):多谢详细解答
herosmzd(金币+1,VIP+0):多谢指教,能帮我解答下对结果的影响吗?
|
你的FORCE项已经很小了，但是DISPLACEMENT项还很大，即势能梯度很小、优化位移很大，说明势能面非常平坦。这种情况很可能是由于DFT对势能面的描述定性不正确造成的。 1楼的朋友说的是不对的，SCF=TIGHT/VERYTIGHT只能提高HF-SCF的收敛精度，对DFT的优化过程特别是优化的收敛性没有直接的帮助。 由于DFT方法和传统的从头计算法出发点不同，当DFT方法描述的势能面出现问题的情况下，建议在计算条件许可时考虑尝试用Post-HF方法尝试进行优化，如原子数不超过50个，可以尝试采用MP2方法，不超过20个，可以尝试采用CCSD方法。由于这些方法对计算机性能要求较高，望LZ慎重考虑！ |
4楼2008-03-03 13:00:22
5楼2008-03-03 17:16:18
6楼2008-03-03 19:12:08
★ ★ ★
zzgyb(金币+3,VIP+0):谢谢你的帮助源，欢迎继续关注计算模拟版！
zzgyb(金币+3,VIP+0):谢谢你的帮助源，欢迎继续关注计算模拟版！
|
个人认为，如果是由于真实的势能面非常平坦造成的，在搜索点的各方向势能梯度都已经很小的情况下，即使搜索位移仍然很大，也可以认为是优化已经完成。比较典型的例子比如乙烷分子的各种构像异构，由于旋转势能壁垒非常小，在绕碳键旋转的时候势能梯度非常小，可以认为其中任意一个构型都是比较合理的构型，可以以一定分布概率存在。 如果是由于DFT的势能面定性不正确造成的优化问题，就不能用B3LYP的方法优化了。B3LYP是一种杂化密度泛函方法，其泛函含有3个来源于HF-SCF的参数。由于这类DFT方法并非来源于量子力学第一性原理，不能简单地将此类密度泛函方法与HF-SCF和Post-HF方法等同起来。 一般来说，对于几何构型优化而言，HF-SCF的方法通常是定性正确的，衡量结果是否可以接受的标准是：键长偏差在0.02A以内，键角偏差在2度以内，LZ的结果是勉强可以接受的。在含有弱相互作用的体系里（如氢键、范德华力等），HF的精度会下降，但此类体系DFT方法通常也很难给出很好的优化结果，适宜采用CI、MP、CC等Post-HF方法 |
7楼2008-03-04 09:11:33
hairan
木虫 (著名写手)
瞌睡木虫
- QC强帖: 3
- 应助: 61 (初中生)
- 贵宾: 1.06
- 金币: 4974.3
- 散金: 8072
- 红花: 44
- 帖子: 2901
- 在线: 733.4小时
- 虫号: 432069
- 注册: 2007-08-12
- 性别: GG
- 专业: 理论和计算化学
★ ★
lei0736(金币+2,VIP+0):多谢详细解答
lei0736(金币+2,VIP+0):多谢详细解答
|
如果按照Gaussian默认的收敛标准，你的结果应该收敛了，因为当force小于0.000004， RMS force小于0.000003的时候，程序会认为force已经可以忽略不计，已经找到了stationary point。你应该是使用了scf=tight一类的关键词 我的计算经验是，遇到此种情况，先把计算停了，然后将结构调出来，重新创建输入文件，开新的计算，输入文件的关键词中不要加入guess=read或者geom=allcheck等从chk文件读内容的选项，有时后两项会很快收敛的，但并不是任何时候都有效。 |
8楼2008-03-05 15:46:45