24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者sky423将赠送您 50 个金币

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

sky423

新虫 (著名写手)

应助: 0 (幼儿园)
金币: 1595.1
散金: 352
红花: 8
帖子: 1048
在线: 103.6小时
虫号: 5800207
注册: 2017-03-01
专业: 有机合成

[求助] 高斯不收敛，如何优化已有1人参与

#p b3lyp/lanl2dz opt pop=(full,npa)

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB+HF-LYP) =  -1714.45115901    A.U. after  129 cycles
         Convg  = 0.8068D-05          -V/T =  2.0529
         S**2 = 0.9172
KE= 1.628346096785D+03 PE=-1.275059909267D+04 EE= 5.005782636124D+03
Annihilation of the first spin contaminant:
S**2 before annihilation    0.9172, after    0.7568
Convergence failure -- run terminated.
Error termination via Lnk1e in d:program filesg03l502.exe at Fri Jun 18 18:59:12 2021.
Job cpu time:  1 days  5 hours 55 minutes 14.0 seconds.
File lengths (MBytes):  RWF= 206 Int=    0 D2E=    0 Chk=    18 Scr=    1

请问高手，我这个从报错看应该是没有收敛，但是应该如何优化或者修改？救救孩子吧。。。。@jiewei

回复此楼

» 猜你喜欢

香港理工大学计算材料与物理课题组招聘博后已经有15人回复
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有131人回复
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？已经有2人回复
英国利物浦大学屠昕和雍雪课题组招收CSC博士生已经有27人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）已经有0人回复
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人已经有45人回复
英国贝尔法斯特女王大学招收 2026 年 CSC/BBSRC phd 机器学习催化剂设计方向已经有0人回复
法国里昂大学招生CSC 资助博士生已经有0人回复
法国里昂大学接收CSC 资助访问学者已经有0人回复
北京工业大学材料科学与工程学院《院士团队》招收工程专业博士研究生（计算优先）已经有1人回复

1楼 2021-06-21 16:30:55

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

paramecium86

版主 (著名写手)

QC强帖: 5
应助: 586 (博士)
贵宾: 0.01
金币: 26316.1
红花: 348
帖子: 2021
在线: 1440.5小时
虫号: 2466443
注册: 2013-05-16
性别: GG
专业: 理论和计算化学
管辖: 量子化学

引用回帖:

15楼: originally posted by sky423 at 2021-06-28 00:44:00
你好，我能再麻烦您问您一个问题吗？就是我输入文件结构是正常的，每个键都连接很好，并且中间也没有死点，完成了，但是结构完成之后部分键就断开了，这属于正常现象吗？我就单纯的是普通优化，这样完全有背我原始 ...

如果基组设置没有问题。电荷和多重度也设置的没有什么问题。结构本身原子间的距离并没有太大异常而仅仅是gaussview不显示某些键那一般不会有什么问题。
见我这个教程3.3节 https://muchong.com/t-14242735-1

3.3 化学键的显示问题。
诸如为什么优化完成后“缺了什么什么键“，“哪个键变成了单/双/三键”，“那个键不见了”等问题。
这类问题的回答是：没有关系。gaussview怎么显示的，都没什么所谓。键显示的符合预期也好或者不符合预期也罢都不能说明计算结果本身是对还是错。
用高斯做的量化计算需要注意的是量化计算中，或者说薛定谔方程中根本没有“键”这种东西。一切都是波函数。也就是一个一个的分子轨道。优化出来的就是原子间的相对位置，因为薛定谔方程需要的只是原子的位置信息。
而gaussview在显示分子的时候会根据一个内在的键长标准（gaussview内置了一套原子的半径信息，会拿原子间的距离相应的两个原子的半径之和作比较，比如键长超过了原子半径之和的百分之多少就属于单键，之类之类的。），完全是根据原子间的距离来显示键的种类的。跟电子真正在分子中的排布一毛钱关系都没有。所以要分析原子间的成键需要进行波函数分析。算算键级之类的去判断。gaussview显示的完全就是为了好看，但是没有任何化学意义的。并且也是可以自己手动更改键的种类的。